Реферат Рішення зворотного завдання динаміки

trix (Aw8, 'Оператор Aw8' );

disp ( 'Оператор Aw8' );

disp (Aw8 (1:8,1:8));

% оператор правою лінійної частини

Aw67 = inv (Ae + Aw7 * Aw6) * Aw6;

Aw671 = inv (Ae + Ksh * Aw67) * Aw67;

Aw_r = Kz * inv (Ae + Aw8 * Aw671 * Aw5) * (Aw671 * Aw5);

% pr_matrix (Aw_r, 'Оператор правій частині Aw_r' );

disp ( 'Оператор правій частині Aw_r' );

disp (Aw_r (1:8,1:8));

5. Спектральна характеристика вхідного сигналу

В 

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

fprintf ('5. Спектральна характеристика вхідного сигналу n ');

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

u = zeros (1, Nt) +1;

Cu = fwht (u ');

% pr_matrix (Cu, 'Cu');

disp ( 'Оператор Cu' );

disp (Cu (1:8));

6. Розрахунок вихідного сигналу методом послідовних наближень

В 

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

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;