Расчетное значение кривой Xp и значение искаженной помехами кривой Xizm=Xp+shum
Нужно найти такие параметры TT1 и TT2, которые обеспечивают минимум критерия близости рассчитанного как (((Xizm-Xp)^2)/N)^0.5
Эта задача у меня выполнена.
А задание состоит в следующем: Уменьшать шаг (shag) до тех пор пока не будет достигнута заданная точность ez (среднеквадратичная ошибка), равная e=((Xe-Xp)*(Xe-Xp)')/N^0.5, где Xe рассчитанное с данными параметрами T1 и T2, а Xp -расчетное значение, но уже с найденными TT1 и TT2, которые обеспечивают минимальный критерий близости.
Помогите а то уже замучился с ней. Код рабочей программы.
Код:
Код: Выбрать все
K=100; T1=4.58; T2=1.72; N=50; i=1:1:N; t=(20/N)*i;
shum=40*rand(1,50);
Xp=K*(1+(T1/(T2-T1))*exp(-t/T1)+(-T2/(T2-T1))*exp(-t/T2));
Xizm=Xp+shum-20;
shag=0.1; ez=5;
for m=1:50;
T10=shag*m;
for l=0:(m-1);
T20=shag*l;
Xp=K*(1+(T10/(T20-T10))*exp(-t/T10)+(-T20/(T20-T10))*exp(-t/T20));
F=(Xizm-Xp)*((Xizm-Xp)');
if m==1 & l==0 Fmin=F; end;
if F<Fmin Fmin=F; T1min=T10; T2min=T20; end;
end;
end;
T0=(Fmin/(N-2))^0.5;
TT1=T1min;
TT2=T2min;
Xe=K*(1+(T1/(T2-T1))*exp(-t/T1)+(-T2/(T2-T1))*exp(-t/T2));
Xp=K*(1+(TT1/(TT2-TT1))*exp(-t/TT1)+(-TT2/(TT2-TT1))*exp(-t/TT2));
e=((Xe-Xp)*(Xe-Xp)')/(N-1)^0.5;
plot(t,Xp,t,Xe,t,Xizm);