Анализ переходных процессов при исследовании динамических моделей технических систем. Навыки выполнения анализа переходных процессов в динамических моделях с графической интерпретацией полученных результатов, страница 2

F=F0*sin(w.*t);   %..гармоническое_воздействие

for i=1:length(t) %..цикл

if (t(i)<=0.5)    %..оператор_выбора

H(i)=0;

else H(i)=1;

end;

end;

figure(1);                          %..создаем_графическое_окно

clf();                              %..очистка_графического_окна

subplot(3,1,1);

plot(t,F,'g',t,H);                  %..график_гармонической_возмущающей_силы_F(t)_и_ступенчатого_воздействия

xlabel('t');                        %..служит_для_вывода_название_координаты_X

ylabel('F..H');                     %..служит_для_вывода_название_координаты_Y

title('График F(t) и H(t)');        %..служит_для_вывода_название_графика

grid;                               %..сетка

subplot(3,1,2);

plot(t,F,'r');                      %..график_гармонической_возмущающей_силы_F(t)

xlabel('t');                        %..служит_для_вывода_название_координаты_X

ylabel('F');                        %..служит_для_вывода_название_координаты_Y

title('График F(t)');               %..служит_для_вывода_название_графика

grid;                               %..сетка

subplot(3,1,3);

plot(t,H);                          %..график_ступенчатого_воздействия

xlabel('t');                        %..служит_для_вывода_название_координаты_X

ylabel('H');                        %..служит_для_вывода_название_координаты_Y

title('График H(t)');               %..служит_для_вывода_название_графика

grid;                               %..сетка

%..2.Рассчитать значение функции перемещения динамической системы без воздействия начальных значений перемещения и скорости с учетом гармонической возмущающей силы F(t). Построить график этой функций.

c=G*d^4/(8*D^3*i);               %..жесткость_пружины

p=((c*l+m*10)/(2*m*l))^1/2;      %..частота_собственных_колебаний_системы

n=a/(2*m*l);                     %..коэффициент_вязкого_сопротивления

y0=[0 0];                        %..начальные_значения

[T,Y]=ode45(@vid1,[0 t_k],y0,[],n,p,F0,w);   %..вызов_функции

figure(2);                       %..создаем_графическое_окно

clf()                            %..очистка_графического_окна

subplot(3,1,1);

plot(T,Y(:,1),'r');              %..график_перемещения_с_учетом_гармонической_возмущающей_силы_F(t)

xlabel('t');                %..служит_для_вывода_название_координаты_X

ylabel('S');                %..служит_для_вывода_название_координаты_Y

title('График S(t) с учетом гармонической возмущающей силы F(t)')        %..служит_для_вывода_название_графика

grid;                       %..сетка

subplot(3,1,2);

plot(T,Y(:,2));             %..график_скорости_с_учетом_гармонической_возмущающей_силы_F(t)

xlabel('t');                %..служит_для_вывода_название_координаты_X

ylabel('v');                %..служит_для_вывода_название_координаты_Y

title('График v(t) с учетом гармонической возмущающей силы F(t)')        %..служит_для_вывода_название_графика

grid;                       %..сетка

clear F           %..очистка_переменных

F=F0*sin(w.*T);   %..гармоническое_воздействие

subplot(3,1,3);

plot(T,F,'g');                      %..график_гармонической_возмущающей_силы_F(t)

xlabel('t');                        %..служит_для_вывода_название_координаты_X

ylabel('F');                        %..служит_для_вывода_название_координаты_Y

title('График возмущающей силы гармонического воздействия F(t)')         %..служит_для_вывода_название_графика

grid;                               %..сетка

%..3.Рассчитать значение функции перемещения динамической системы без воздействия начальных значений перемещения и скорости с учетом ступенчатого воздействия. Построить график этой функций.

y0=[0 0];                               %..начальные_значения

[T,Y]=ode45(@vid2,[0 t_k],y0,[],n,p);   %..вызов_функции

figure(3);                              %..создаем_графическое_окно

clf()                                   %..очистка_графического_окна

subplot(3,1,1);

plot(T,Y(:,1),'r');         %..график_перемещения_с_учетом_ступенчатого_воздействия

xlabel('t');                %..служит_для_вывода_название_координаты_X

ylabel('S');                %..служит_для_вывода_название_координаты_Y

title('График S(t) с учетом ступенчатого воздействия')        %..служит_для_вывода_название_графика

grid;                       %..сетка

subplot(3,1,2);

plot(T,Y(:,2));             %..график_скорости_с_учетом_ступенчатого_воздействия

xlabel('t');                %..служит_для_вывода_название_координаты_X

ylabel('v');                %..служит_для_вывода_название_координаты_Y

title('График v(t) с учетом ступенчатого воздействия')        %..служит_для_вывода_название_графика

grid;                       %..сетка

clear H           %..очистка_переменных