fig6_41.m

来自「包含了控制理论的各种仿真程序」· M 代码 · 共 32 行

% Fig. 6.41   Feedback Control of Dynamic Systems, 4e 
%             Franklin, Powell, Emami
%

clear all;
close all;

num=86*conv([1 1],[1 2 43.25]);
p=conv([1 2 82],[1 2 101]);
den=conv([1 0 0],p);
w=logspace(-1,2,1000);
[mag,phas]=bode(num,den,w);
[re,im]=nyquist(num,den,w);

ii=200:1000;
plot(re(ii),im(ii));
xlabel('Real');
ylabel('Imaginary');
grid;
axis([-1.2 1.2 -1.2 1.2])
axis('square')
hold on;

% add circle at mag = 1
w2=0:.05:2*pi;
x=cos(w2);
y=sin(w2);
plot(x,y,'r--')
title('Figure 6.41 Nyquist plot of a complex system.');
hold off;