
2074
.pdf
3.3. Обработка модулем «Множественная регрессия»
Multiple Regression
Вернемся к окну на рис 3.27 и выберем модуль Multiple Regression (строка, расположенная над выделенной). Определим переменные в модели (зависимые и предикты-независимые). Используем опцию Пошаговая с включением. Имеется возможность использовать пошаговый метод как с включением предиктов, так и с их исключением. На каждом шаге добавления предиктов рассчитывают коэффициент множественной корреляции. Его квадрат – коэффициент детерминации – показывает качество построенной модели ОК. В диалоговом окне результатов указываются стандартизованные коэффициенты регрессии.
Рис.3.27. Окно результатов оценивания параметров модели
Таблица 3 . 9 Результаты оценивания параметров переменных
Regression Summary for Dependent Variable: Y (new.sta)
Множественный коэффициент корреляции – R= 0,99318636 Коэффициент детерминации – RІ= 0,98641915 – очень хорошее значение
Скорректированный [1-(1-R2)(n/n-КолКоэфРегр)] – Adjusted RІ= 0,98224043
F(4,13)=236,06 – значимость регрессии, но не адекватность |
Уровень значимости ошибки |
||||||
Станд.ошибка оценки – Std.Error of estimate: 0,00011 |
|
p<0,00000 – крайне незначителен |
|||||
|
Стандартиз. |
Ст.ошибка |
Не стандарт. |
Ст.ошибка |
Критерий |
Ур.знач. |
|
|
коэф регрес |
БЕТА |
коэф.регресс |
|
В |
Стьюдента |
ошибки |
Своб.чл. |
BETA |
St.Err.ofBETA |
B |
St. Err.of B |
t(13) |
p-level |
|
Intercpt |
|
|
0,004915 |
2,66E-05 |
184,7333 |
1,29E-23 |
|
X1 |
-0,05037 |
0,033912 |
-4,70E-05 |
3,14E-05 |
-1,48535 |
0,161293 |
|
X2 |
0,7083 |
0,033565 |
0,000655 |
3,10E-05 |
21,1022 |
1,93E-11 |
|
X3 |
-0,5467 |
0,033481 |
-0,00048 |
2,92E-05 |
-16,33 |
4,82E-10 |
|
X4 |
0,026688 |
0,032667 |
2,39E-05 |
2,92E-05 |
0,81698 |
0,428659 |
171
Например, вводим значения независимых переменных Х1, Х2, Х3, Х4 и зависимой Y. В результате обработки получаем результат (см. рис.3.24). Значимые коэффициенты выделяются красным цветом. Сравните значения выделенных показателей с другими и сделайте вывод о степени их достоверности и влияния: коэффициентов регрессии, критериев Стьюдента, уровней значимости ошибки. Адекватность модели в данных модулях проверяется величиной и распределением остатков. Величина БЕТА показывает влияние факторов в сопоставимых единицах (как бы в «кодированном» виде).
Левые кнопки позволяют определить статистики (сверху вниз): Итоговая таблица регрессии (табл.3.9); Дисперсионный анализ (табл.3.10); Ковариации коэффициентов (табл.3.11); Текущая матрица выметания (табл.3.12); Частные корреляции (табл.3.13) более 0,8 свидетельствуют о «дублировании» параметров и должны быть исключены. Низкие значения ковариации между коэффициентами В говорят об отсутствии дублирования показателей (корреляции между ними).
Таблица 3 . 1 0 Анализ дисперсии зависимой переменной
|
Analysis of Variance; DV: Y (new.sta) |
|
|
|
|
|
Суммаквадратов |
Числостепеней |
Среднее |
|
Ур. знач. |
|
остатков |
свободы |
квадр. |
|
ошибки |
|
Sums of Squares |
df |
Mean Squares |
F |
p-level |
Regress. |
1,18E-05 |
4 |
2,95E-06 |
236,0576 |
5,42E-12 |
Residual |
1,63E-07 |
13 |
1,25E-08 |
|
|
Total |
1,2E-05 |
|
|
|
|
|
|
|
Таблица 3 . 1 1 |
||
|
Ковариация коэффициентов регрессии В |
|
|||
|
|
|
|
|
|
|
Covariances of Regression Coefficients B; DV: (new.sta) |
|
|||
|
X1 |
X2 |
X3 |
|
X4 |
X1 |
9,83E-10 |
1,78E-10 |
-1,5E-10 |
|
-1,3E-10 |
X2 |
1,78E-10 |
9,63E-10 |
1,49E-10 |
|
-2,3E-11 |
X3 |
-1,5E-10 |
1,49E-10 |
8,5E-10 |
|
2,69E-11 |
X4 |
-1,3E-10 |
-2,3E-11 |
2,69E-11 |
|
8,55E-10 |
|
|
|
|
Таблица 3 . 1 2 |
||
|
Текущий статус размаха переменных матрицы |
|
||||
|
|
|
|
|
|
|
|
Current Status of Sweep Matrix; DV: Y (new.sta) |
|
||||
|
X1 |
X2 |
X3 |
X4 |
|
Y |
X1 |
-1,10083 |
-0,19892 |
0,183163 |
0,153536 |
|
-0,05037 |
X2 |
-0,19892 |
-1,07842 |
-0,17734 |
0,026221 |
|
0,708296 |
X3 |
0,183163 |
-0,17734 |
-1,07301 |
-0,03309 |
|
-0,54672 |
X4 |
0,153536 |
0,026221 |
-0,03309 |
-1,02147 |
|
0,026688 |
Y |
-0,05037 |
0,708296 |
-0,54672 |
0,026688 |
|
0,013581 |
172

Таблица 3 . 1 3
Частные корреляции
Variables currently in the Equation; DV: Y (new.sta)
|
|
Частнаякорр. |
Общаякорр. |
ТолерантностьКоэф. детерм. |
Стьюдент |
Ур.знач. ош. |
|
|
Beta in |
Partial Cor. |
Semipart Cor. |
Tolernce |
R-square |
t(13) |
p-level |
X1 |
-0,05037 |
-0,38091 |
-0,04801 |
0,908408 |
0,091592 |
-1,48535 |
0,161293 |
X2 |
0,708296 |
0,985715 |
0,682056 |
0,927281 |
0,072719 |
21,10224 |
1,93E-11 |
X3 |
-0,54672 |
-0,97648 |
-0,5278 |
0,93196 |
0,06804 |
-16,3296 |
4,82E-10 |
X4 |
0,026688 |
0,220987 |
0,026406 |
0,978982 |
0,021018 |
0,81698 |
0,428659 |
Максимальная корреляция с откликом у второго и третьего фактора, т.е. именно они обеспечивают его значение. Высокая толерантность отдельных факторов (как и низкое значение их R2) говорит о недостаточном их личном вкладе в результат (возможно толерантность и R2 перепутаны местами!).
По центру вверху кнопка – Предсказать зависимую переменную. Ни-
же ее кнопка Redundansy – избыточность – даст табл. 3.14.
Таблица 3 . 1 4 Показатели независимых переменных
Redundancy of Independent Variables; DV: Y (new.sta)
R-square column contains R-square of respective
variable with all other independent variables
|
Toleran. |
R-square |
Partial Cor. |
Semipart Cor. |
X1 |
0,908408 |
0,091592 |
-0,38091 |
-0,04801 |
X2 |
0,927281 |
0,072719 |
0,985715 |
0,682056 |
X3 |
0,93196 |
0,06804 |
-0,97648 |
-0,5278 |
X4 |
0,978982 |
0,021018 |
0,220987 |
0,026406 |
Левая нижняя кнопка Корреляции и описательные статистики – Ковариации даст табл. 3.15. Кнопка над ней – Анализ остатков – откроет окно на рис.3.28.
Таблица 3 . 1 5
Ковариация переменных
Covariances (new.sta)
|
X1 |
X2 |
X3 |
X4 |
Y |
X1 |
0,823529 |
-0,17647 |
0,176471 |
0,117647 |
-0,00024 |
X2 |
-0,17647 |
0,823529 |
-0,17647 |
1,54E-18 |
0,000632 |
X3 |
0,176471 |
-0,17647 |
0,928105 |
-0,00654 |
-0,00057 |
X4 |
0,117647 |
1,54E-18 |
-0,00654 |
0,879085 |
1,86E-05 |
Y |
-0,00024 |
0,000632 |
-0,00057 |
1,86E-05 |
7,04E-07 |
173

Рис.3.28. Окно результатов оценивания остатков модели
В левой половине (сверху вниз) кнопки: Описательные статистики;
Итоги регрессии; Остатки и предсказания; Статистики ДарбинаУотсона (табл.3.16); Сохранять остатки и предсказания; Построенные графики (табл. 3.17) – Графики остатков, Графики выбросов, Графики предсказанных значений (расчетные по модели).
|
|
|
Таблица 3 . 1 6 |
|
|
Статистики Дарбина-Уотсона |
|
|
|
|
|
|
Durbin-Watson d (new.sta) and serial correlation of residuals |
||
|
|
Durbin-Watson d |
Serial Corr. |
Estimate |
|
1,805276 |
0,087239 |
Таблица 3 . 1 7 Графики и значения опытных, расчетных показателей
Plots of residuals (A) – Raw Residual |
Опытные |
Расчетные |
Разница |
|
|
|
Observed |
Predictd |
Residual |
|
1 |
2 |
3 |
4 |
1 . . . |
* . . . . |
0,004611 |
0,004666 |
-5,5E-05 |
2 . . . |
. * . . . |
0,006139 |
0,006069 |
6,97E-05 |
3 . . . |
.* . . . |
0,003831 |
0,003807 |
2,34E-05 |
4 . . . |
. * . . . |
0,005083 |
0,005024 |
5,96E-05 |
5 . . . |
*. . . . |
0,004681 |
0,004714 |
-3,4E-05 |
6 . . . |
.* . . . |
0,006139 |
0,006117 |
2,2E-05 |
7 . . . |
. . * . . |
0,004014 |
0,003855 |
0,000159 |
8 . . . |
. * . . . |
0,005133 |
0,005071 |
6,18E-05 |
9 . . . |
* . . . . |
0,004875 |
0,004915 |
-4E-05 |
174

Окончание табл. 3 . 1 7
|
1 |
|
2 |
3 |
4 |
10 . . |
. . |
* . . . |
0,006139 |
0,00607 |
6,85E-05 |
11 . |
. * |
. . . . |
0,005289 |
0,005414 |
-0,00013 |
12 . . |
* . |
. . . . |
0,004217 |
0,004417 |
-0,0002 |
13 . . |
. . |
* . . . |
0,006139 |
0,006093 |
4,59E-05 |
14 . . |
* . |
. . . . |
0,004317 |
0,00451 |
-0,00019 |
15 . . |
. * . . . . |
0,005022 |
0,00507 |
-4,8E-05 |
|
16 . . |
. . |
* . . |
0,00385 |
0,003738 |
0,000112 |
17 . . |
. . |
* . . . |
0,004294 |
0,004238 |
5,65E-05 |
18 . . |
. . |
* . . . |
0,003778 |
0,003761 |
1,72E-05 |
Minimum . |
.* . . . . . |
0,003778 |
0,003738 |
-0,0002 |
|
Maximum . |
. . . . * . . |
0,006139 |
0,006117 |
0,000159 |
|
Mean . |
. |
. * . . . |
0,004864 |
0,004864 |
9,05E-11 |
Median |
. . . .* . . . |
0,004778 |
0,004815 |
2,27E-05 |
По центру кнопки: Предсказанные остатки; Предсказанные и квадратичные остатки; Предсказанные и наблюдаемые (рис.3.29);
Наблюдаемые и остатки; Наблюдаемые и квадрат остатков; Остатки и удаленные остатки. Гистограммы: График наблюдаемых; График предсказанных; График остатков.
Observed Values
Predicted vs. Observed Values Dependent variable: Y
0,006
0,006
0,006
0,005
0,005
0,004
0,004 |
|
|
|
|
|
|
|
0,004 |
|
|
|
|
|
Regression |
|
|
|
|
|
|
95% confid. |
||
0,003 |
0,004 |
0,005 |
0,005 |
0,006 |
0,006 |
||
|
Predicted Values
Рис.3.29. Соответствие расчетных и исходных значений 95 %-му доверительному интервалу
Высокое значение статистики Дарбина-Уотсона говорит о независимости остатков, поэтому адекватность модели принимается (не
отвергается) на уровне =р.
По результатам опытных и расчетных значений (табл.3.17), используя программу Excel, найдено значение F-тест=0,97783, что говорит об одинаковой степени плотности (разброса) показателей, а с учетом незначительной величины значений Разности – о достоверности модели. Кроме того, значение критерия Фишера, определенное программой Mathcad
(пример), говорит об адекватности модели: F=1,258 Fкр=2,272.
175

Пример
Определение существенности различий выборок (листинг программы Mathcad)
В результате обработки опытных данных программой Statistica получены две выборки результатов: первая – расчетные значения по математической модели, вторая – результаты опытов.
Количество значений в выборках (число строк плана): N |
18 i 1 N |
|
|||
Результаты опыта: |
|
|
|
|
|
y1 0.004611 |
y2 0.005083 |
y3 0.004014 |
y4 0.006139 |
y5 0.00385 |
y6 0.003778 |
y7 0.006139 |
y8 0.004681 |
y9 0.005133 |
y10 0.00431 |
y11 0.00429 |
y12 0.00423 |
y13 0.003831 |
y14 0.006139 |
y15 0.004875 |
y16 0.005022 |
y17 0.00377 |
y18 0.00376 |
Результаты расчетов: |
|
|
|
|
|
x1 0.004666 |
x2 0.005024 |
x3 0.003855 |
x4 0.00607 |
x5 0.006093 |
x6 0.003738 |
x7 0.006069 |
x8 0.004714 |
x9 0.005071 |
x10 0.005414 |
x11 0.00451 |
x12 0.00423 |
x13 0.003807 |
x14 0.006117 |
x15 0.004915 |
x16 0.004417 |
x17 0.00507 |
x18 0.00376 |
Среднее значение результата: |
|
|
|
|
|
yi |
xi |
|
|
|
|
i |
|
i |
|
|
|
N |
|
|
N |
|
|
Дисперсия каждой выборки значений: |
|
||||
S2y |
1 |
yi |
Ycp |
2 |
|
|
N 1 |
|
|||
|
|
|
i |
|
….. |
|
S2y |
6.881 10 7 |
|
Ycp 0.00465 |
Xcp 0.00486 |
||||
S2x |
1 |
|
xi Xcp |
2 |
|
N 1 |
|
||||
|
|
|
|
i |
|
S2x 6.942 10 7
Уровень критерия проверки гипотезы (значимости ошибки): 0.05 |
|
||||
F-критерий Фишера: |
S0 S2x |
S1 |
S2y |
F max(S) |
F 1.258 |
|
|
|
|
min(S) |
|
Число степеней свободы остаточной дисперсии: N 1 |
|
||||
Критическое значение F-критерия Фишера: |
Fkp qF 1 |
||||
|
Fkp 2.272 |
F |
Fkp 1 |
|
Поскольку условие оказалось не ложным (равным 1, а не 0), то гипотезу следует принять (не отвергнуть), то есть мы должны сделать вывод о равноточности или равнорассеянности дисперсий (различия представленных выборок незначимы), следовательно, полученная модель адекватно описывает результаты опытов.
Квантиль распределения Стьюдента (критическое значение):
|
|
|
|
|
|
|
T qt |
1 |
|
2 |
2N 2 |
|
T 2.032 |
Проверка на существенность различия между замеряемыми значениями объектов:
T1 |
T |
|
T1 2.389 10 3 |
|
|
|
|
|
S2y S2x |
|
|
|
|
|
|||
T0 |
Xcp Ycp |
T0 2.148 10 4 |
|
Xcp Ycp |
|
T |
S2y S2x 0 |
|
|
|
Разница принимается существенной, если выполняется условие: значение равно 1, а не нулю. Так как условие не выполнено, разница между выборками не существенна.
176

Расчетные значения (см. рис. 3.29) хорошо сходятся с исходными параметрами, что подтверждает достоверность модели. По отдельности каждый фактор (например, Х1 по рис. 3.30) не обеспечивает результата, так как велик разброс значений относительно 95 %-го доверительного интервала.
X1 vs. Y
Y = ,00486 - ,0003 * X1
Correlation: r = -,3088
0,006 |
|
|
|
|
|
|
|
|
0,006 |
|
|
|
|
|
|
|
|
0,006 |
|
|
|
|
|
|
|
|
0,005 |
|
|
|
|
|
|
|
|
Y |
|
|
|
|
|
|
|
|
0,005 |
|
|
|
|
|
|
|
|
0,004 |
|
|
|
|
|
|
|
|
0,004 |
|
|
|
|
|
|
|
|
0,004 |
|
|
|
|
|
|
Regression |
|
|
|
|
|
|
|
95% confid. |
||
-1,2 |
-0,8 |
-0,4 |
0,0 |
0,4 |
0,8 |
1,2 |
||
|
||||||||
|
|
|
X1 |
|
|
|
|
Рис.3.30. Наличие линейной связи между зависимой и частной независимой переменной
Справа кнопки (см. рис. 3.28): Графики остатков – Нормальный,
Полунормальный, Без тренда; Диаграммы рассеяния, включающие переменную – Две переменные, Остатки и переменная, Предсказанные и переменная, Частные остатки.
Изучая данные графики и полученные значения, выполняем анализ и делаем необходимый вывод.
Контрольные вопросы
1.Как осуществляется ввод данных в программе STATISTICA?
2.Как осуществляется выборочная корреляция факторов?
3.Как определяется критерий наименьшей значимой разности?
4.Как осуществляется построение гистограмм?
5.Как выполняется обработка модулем «Множественная регрессия» ?
6.Как определить соответствие расчетных и исходных значений 95 %-му доверительному интервалу?
7.Определение результатов оценивания параметров модели.
177
4. ПРИМЕР ВЫПОЛНЕНИЯ ДИССЕРТАЦИИ
ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ «ПЕНЗЕНСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ АРХИТЕКТУРЫ И СТРОИТЕЛЬСТВА»
На правах рукописи
ТАРАСОВ Александр Иванович
ОЦЕНКА ТЕХНИЧЕСКОГО СОСТОЯНИЯ ДВИГАТЕЛЕЙ В ПРОЦЕССЕ ЭКСПЛУАТАЦИИ И РЕМОНТА АВТОМОБИЛЕЙ
Специальность 05.22.10 – Эксплуатация автомобильного транспорта
Диссертация на соискание ученой степени кандидата технических наук
Научный руководитель – кандидат технических наук, доцент Лянденбурский В.В.
Пенза – 2013
178
СОДЕРЖАНИЕ |
|
ВВЕДЕНИЕ ...................................................................................................... |
181 |
Глава 1. СОВРЕМЕННОЕ СОСТОЯНИЕ ВОПРОСА................................. |
185 |
1.1. Общие сведения о неисправностях и их влияние |
|
на состояние системы питания дизельного двигателя.................... |
185 |
1.2. Роль диагностирования в повышении эффективности технической |
|
эксплуатации автомобильных дизелей............................................. |
189 |
1.3. Сравнительный анализ разработок в области диагностирования |
|
автомобильных дизелей..................................................................... |
191 |
1.5.Анализ методов поиска неисправностей .......................................... |
199 |
1.5. Классификация контрольно-измерительных приборов................. |
204 |
1.6. Выводы и задачи исследования........................................................ |
215 |
Глава 2. ТЕОРЕТИЧЕСКИЕ ИССЛЕДОВАНИЯ......................................... |
216 |
2.1. Описание объекта исследований...................................................... |
216 |
2.2. Режимы обеспечения работоспособности автомобиля и его |
|
топливной системы............................................................................. |
217 |
2.3. Целевая функция исследований........................................................ |
218 |
2.4. Вероятностно-логический подход к выявлению неисправностей |
|
автомобилей......................................................................................... |
224 |
2.4.1. Вероятностно-логическая модель как способ диагностирования |
|
автомобиля.................................................................................... |
224 |
2.4.2. Система коэффициентов вероятностно-логической модели |
|
поиска неисправностей автомобилей......................................... |
228 |
2.5. Выводы по второй главе:................................................................... |
240 |
Глава 3. ЭКСПЕРИМЕНТАЛЬНЫЕ ИССЛЕДОВАНИЯ............................ |
241 |
3.1. Общая методика исследований......................................................... |
241 |
3.2. Обоснование плана и объема исследований.................................... |
242 |
3.3. Общее назначение и цели выполнения экспериментальных |
|
исследований....................................................................................... |
244 |
3.4. Характеристика объекта исследований и методика проведения |
|
исследований....................................................................................... |
246 |
3.5. Анализ статистических данных по отказам элементов дизельной |
|
топливной системы при эксплуатации автомобилей КАМАЗ в |
|
условиях Пензы и Рязани................................................................... |
247 |
3.6. Сравнительный анализ результатов исследования поиска |
|
неисправностей................................................................................... |
253 |
3.6.1. Определение времени простоя автомобилей................................ |
253 |
179 |
|
3.6.2. Анализ удельных затрат на поиск неисправностей автомобилей |
|
КАМАЗ.......................................................................................... |
256 |
3.6.3. Анализ удельных затрат и эффективности применения |
|
встроенной системы диагностирования для топливной системы |
|
автомобилей КАМАЗ................................................................... |
260 |
3.7. Выводы по третьей главе................................................................... |
263 |
Глава 4. РЕАЛИЗАЦИЯ РЕЗУЛЬТАТОВ ИССЛЕДОВАНИЙ................... |
265 |
4.1. Структура и описание работы программы по диагностированию |
|
технического состояния автомобиля................................................ |
265 |
4.2. Макетный образец встроенной системы диагностирования |
|
автомобильных дизелей..................................................................... |
273 |
4.3. Оценка экономической эффективности внедрения системы |
|
диагностирования автомобильных дизелей..................................... |
276 |
ЗАКЛЮЧЕНИЕ................................................................................................ |
280 |
ПРИЛОЖЕНИЯ ............................................................................................... |
293 |
180