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

p>
В 

Графік вихідного сигналу при навантаженні = 573

В 

На малюнках 1 і 2 представлені результати аналізу системи з використанням методу матричних операторів і з використанням функцій Уолша для вхідного сигналу і для порівняння наведені графіки необхідного вихідного сигналу, а також сигналу, який може забезпечити дана система при значенні навантаження = 573.

Лістинг програм:

В 

Програма аналізу електрогідравлічного слідкуючого приводу (Рульової машинки як окремого елемента системи самонаведення) з використанням спектрального методу (базис функцій Уолша)

close all;

clear all;

clc;

warning off;

tic;

1. Параметри системи і інтервал дослідження

В 

egsp_data;

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

fprintf ('1. Інтервал дослідження n ');

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

fprintf ('tmin =% e, c; n', tmin);

fprintf ('tmax =% e, c; n', tmax);

fprintf ('Nt =% i; n', Nt);

fprintf (' n');

2. Формування системи базисних функцій

В 

settime (T);

setsize (Nt);

Ai = mkint;

Ad = inv (Ai);

Ae = eye (Nt);

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

fprintf ('2. Формування системи функцій Уолша n ');

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

% pr_matrix (Ai, 'Оператор інтеграції Ai');

disp ( 'Оператор інтеграції Ai');

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

% pr_matrix (Ad, 'Оператор диференціювання Ad');

disp ( 'Оператор диференціювання Ad');

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

3. Розрахунок операторів лівої лінійної частини

В 

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

fprintf ('3. Оператори лівої лінійної частини n ');

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

% оператор ПФ W1 (s) - електричної частини

Aw1 = inv (RS * (Ty * Ae + Ai)) * (Ky1 * KU * Ai);

% pr_matrix (Aw1, 'Оператор Aw1 ');

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

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

% оператор ПФ W2 (s) - електромагнітного перетворювача і частина витрат

Aw2 = inv (CS * (Tem ^ 2 * Ae +2 * Tem * dzem * Ai + Ai ^ 2)) * (KFi * Kqh * Ai ^ 2);

% pr_matrix (Aw2, 'Оператор Aw2' );

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

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

% оператор ПФ W3 (s) - руху золотника і частина витрат

Aw3 = inv (Kqp1 * (Cp + Cg) * (Tz ^ 2 * Ae +2 * Tz * dzz * Ai + Ai ^ 2)) * (Az * Ai ^ 2);

% pr_matrix (Aw3, 'Оператор Aw3' );

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

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

% оператор ПФ W4 (s) - місцевої зворотного зв'язку

Aw4 = Az * Ad;

% pr_matrix (Aw4, 'Оператор Aw4' );

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

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

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

Aw34 = inv (Ae + Aw4 * Aw3) * Aw3;

Aw_1 = Aw34 * Aw2 * Aw1;

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

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

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

4. Розрахунок операторів правою лінійної частини

В 

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

fprintf ('4 . Оператори правою лінійної частини n ');

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

% оператор ПФ W5 (s) - рівняння витрат

Aw5 = inv (Kqp * (Tg * Ae + Ai)) * (Ap * l * Ai);

% pr_matrix (Aw5, 'Оператор Aw5' );

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

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

% оператор ПФ W6 (s) - навантаження

Aw6 = J * Ai ^ 2;

% pr_matrix (Aw6, 'Оператор Aw6' );

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

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

% оператор ПФ W7 (s) - тертя

Aw7 = Kf * Ad;

% pr_matrix (Aw7, 'Оператор Aw7' );

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

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

% оператор ПФ W8 (s) - місцева зворотній зв'язок

Aw8 = Ap * l * Ad;

% pr_ma...