Читать курсовая по физике: "Решение обратной задачи динамики" Страница 3

(Назад) (Cкачать работу)

fprintf('6. Расчет выходного сигнала\n');

fprintf('---------------------------\n');

Cd_old = zeros(Nt,1);

Ce = Cu-Cd_old;

Cx = Aw_1*Ce;

x = iwht(Cx)';

xf = egsp_f(x,xm);

Cxf = fwht(xf');

Cd_new = Aw_r*Cxf;

Ceps = Cd_new-Cd_old;

Eeps = sqrt(Ceps'*Ceps);

fprintf('Нулевое приближение\n');

fprintf('Eeps = %e;\n',Eeps);

d = iwht(Cd_new)';

figure; clf;

plot(t,d);

xlabel('t, c');

ylabel('delta(t)');

Niter = 0;

while Eeps >= 1e-8

Niter = Niter+1;

Cd_old = Cd_new;

Ce = Cu-Cd_old;

Cx = Aw_1*Ce;

x = iwht(Cx)';

xf = egsp_f(x,xm);

Cxf = fwht(xf');

Cd_new = Aw_r*Cxf;

Ceps = Cd_new-Cd_old;

Eeps = sqrt(Ceps'*Ceps);

end

fprintf('%i-ое приближение:\n',Niter);

fprintf('Eeps = %e;\n',Eeps);

d = iwht(Cd_new)';

%my_plot2(t,d,'t, c','delta(t)');

plot(t,d);

xlabel('t, c');

ylabel('delta(t)');

grid on; toc;

Похожие работы