Организация, проведение и обработка результатов испытаний на основе положений пассивного эксперимента, построение автокорреляционной и взаимнокорреляционной функций входных переменных, получение уравнения регрессии, страница 3

c2=zeros(kon-1,1);

c5=zeros(kon-1,1);

c6=zeros(kon-1,1);

a1=(1/((kon-1)*(rez_disp*soc_disp)^0.5));

a2=(1/((kon-1)*(rez_disp*envir_disp)^0.5));

a5=(1/((kon-1)*(soc_disp*soc_disp)^0.5));

a6=(1/((kon-1)*(envir_disp*envir_disp)^0.5));

for i=1:kon-1

b1=(soc(i,1)-soc_sred)*(rez(i+1,1)-rez_sred);

b2=(envir(i,1)-envir_sred)*(rez(i+1,1)-rez_sred);

b5=(soc(i,1)-soc_sred)*(soc(i+1,1)-soc_sred);

b6=(envir(i,1)-envir_sred)*(envir(i+1,1)-envir_sred);

c1(i,1)=a1*b1;                                                                                              %взаимная корреляция выходного параметра и влияния социальной среды

c2(i,1)=a2*b2;                                                                                              %взаимная корреляция выходного параметра и влияния окружающей среды

c5(i,1)=a5*b5;                                                                                              %автокорреляционная функция влияния социальной среды

c6(i,1)=a6*b6;                                                                                              %автокорреляционная функция влияния окружающей среды

end

figure(1);

plot(c1(:,1));

figure(2);

plot(c2(:,1));

figure(3);

plot(rez(:,1));

figure(4);

plot(prise_new(:,1));

figure(5);

plot(c5(:,1));

%-------------------------------------------------Расчёт коэффициентов уравнения регрессии--------------------------------------------------------b1=zeros(2,1);

b2=zeros(2,1);

Ey_rez=zeros(kon,1);

Ey_prise=zeros(kon,1);

Ex_soc=zeros(kon,1);

Ex_envir=zeros(kon,1);

Ex1=zeros(kon,2);

for i=1:kon

Ey_rez(i)=(rez(i)-rez_sred)/rez_disp;                                                          %переход к нормированному представлению данных

Ex_soc(i)=(soc(i)-soc_sred)/soc_disp;                                                        %переход к нормированному представлению данных

Ex_envir(i)=(envir(i)-envir_sred)/envir_disp;                                             %переход к нормированному представлению данных

end 

for i=1:kon  

Ex1(i,1)=Ex_soc(i);

Ex1(i,2)=Ex_envir(i);

end

b1=(Ex1'*Ex1)^-1*Ex1'*Ey_rez;                                                                 %расчёт коэффициентов первого уравнения регрессии

b1

k=kon/(3600/shag);

Ex_2=zeros(k,2);

for i=1:k

Ex_2(i,1)=soc_k(i);

Ex_2(i,2)=envir_k(i);

end

b2=(Ex_2'*Ex_2)^-1*Ex_2'*s                                                                      %расчёт коэффициентов второго уравнения регрессии

%-------------------------------------------------------------------------------------Yrasch1=zeros(kon,1);    

Yrasch1=Ex1*b1;                                                                                          % расчетные значения выходной переменной по температуре помещения

Yrasch2=zeros(k,1);

Yrasch2=Ex1*b2;                                                                                          % расчетные значения выходной переменной по стоимости

Yrasch2_=zeros(k,1);

for i=2:k

Yrasch2_(1)=Yrasch2(1);

Yrasch2_(i)=Yrasch2(10*i);

end

%----------------------------------------------Расчёт коэффициентов по критерию  Стьюдента--------------------------------------------------------c1=(soc'*soc)^-1

c2=(envir'*envir)^-1

t1=zeros(2,1);

t2=zeros(2,1);

t1(1)=abs(b1(1))/(soc_disp*sqrt(c1));

t1(2)=abs(b1(2))/(envir_disp*sqrt(c2));

t2(1)=abs(b2(1))/(soc_disp*sqrt(c1));

t2(2)=abs(b2(2))/(envir_disp*sqrt(c2));

t1

t2

%---------------------------------------------------------Проверка по критерию Фишера----------------------------------------------------figure(7)

plot(soc);

hold on;

figure(7)

plot(envir);

hold off;

figure(8);

plot(rez);

hold on;

figure(8);

plot(Yrasch1);

hold off;

figure(9);

plot(Yrasch2_);

hold on;

figure(9);

plot(s);

hold off;

Y_1=zeros(kon,1);

Y_2=zeros(kon,1);

Y_1=Ex_soc*b1';

Y_2=Ex_envir*b2';

rez_=zeros(kon,1);

s_=zeros(k,1);

for i=1:kon

rez_(i)=rez(i)-Y_1(i);

end

for i=1:k

s_(i)=s(i)-Y_2(i);

end

S_ost1=(rez_'*rez_)/(kon-2);

S_ost2=(s_'*s_)/(k-2);

F1=rez_disp/S_ost1

F2=s_disp/S_ost2

if F1<1.2

disp('1-ое уравнение регрессии адекватно эксперименту')

end

if F1>=1.2

disp('Уравнение регрессии НЕ адекватно эксперименту')

end

if F2<1.2

disp('2-ое уравнение регрессии адекватно эксперименту')

end

if F2>=1.2

disp('Уравнение регрессии НЕ адекватно эксперименту')

end

Приложение 3.

Модель объекта испытаний.