Решение задач оптимального управления с использованием математической системы MATLAB и пакета имитационного моделирования SIMULINK. Сивохин А.В - 111 стр.

UptoLike

111
if (R<sqrt(2)./2 && R>0)
disp('Слабо-демпфированное звено')
end
if (R==0)
disp('Недемпфированное звено')
end
%
%-- 4.Частотная функция:
%
Kiw=1./(A.*(i.*w).^2+B.*i.*w+C);
%
%-- 5.Модуль частотной функции:
%
subplot(2,2,1)
plot(w,abs(Kiw),'r')
xlabel('w')
ylabel('abs(Kiw)')
%
%-- 6.Фаза частотной функции:
%
subplot(2,2,2)
plot(w,angle(Kiw),'r')
xlabel('w')
ylabel('angle(Kiw)')
%
%-- 7.Годограф на плоскости:
%
subplot(2,2,3)
plot(real(Kiw),imag(Kiw),'r')
xlabel('real(Kiw)')
ylabel('image(Kiw)')
%
%-- 8.Годограф в трехмерном пространстве(комета):
%
subplot(2,2,4)
comet3(real(Kiw),imag(Kiw),w)
xlabel('real(Kiw)')
ylabel('image(Kiw)')
zlabel('w')
pause
%
%-- 9.Конец функции RKiw.
   if (R0)
      disp('Слабо-демпфированное звено')
   end
   if (R==0)
      disp('Недемпфированное звено')
   end
%
%-- 4.Частотная функция:
%
   Kiw=1./(A.*(i.*w).^2+B.*i.*w+C);
%
%-- 5.Модуль частотной функции:
%
   subplot(2,2,1)
   plot(w,abs(Kiw),'r')
   xlabel('w')
   ylabel('abs(Kiw)')
%
%-- 6.Фаза частотной функции:
%
   subplot(2,2,2)
   plot(w,angle(Kiw),'r')
   xlabel('w')
   ylabel('angle(Kiw)')
%
%-- 7.Годограф на плоскости:
%
   subplot(2,2,3)
   plot(real(Kiw),imag(Kiw),'r')
   xlabel('real(Kiw)')
   ylabel('image(Kiw)')
%
%-- 8.Годограф в трехмерном пространстве(комета):
%
   subplot(2,2,4)
   comet3(real(Kiw),imag(Kiw),w)
   xlabel('real(Kiw)')
   ylabel('image(Kiw)')
   zlabel('w')
   pause
%
%-- 9.Конец функции RKiw.




                                    111