ex4_20.m

来自「华东理工大学自动化系《控制系统分析、设计和应用》教材配套程序实例,有关说明： 」· M 代码 · 共 31 行

% ex4_20
% 已知多变量系统的传递函数矩阵，绘制格希哥仁圆和格希哥仁带
clear;num={[1],[0.75];[0.75],[1]};
den={[1 1],[2 1];[2 1],[0.5 1]};
G=tf(num,den);
b=10;
for i=1:2;
    for k=1:2
        w=[0,logspace(-1,0.2,b)];
        gg(i,k,:)=polyval(deconv(den{i,k},num{i,k}),j*w);
    end;
end;
figure('pos',[100,100,400,210],'color','w');
    g(1,:)=gg(1,1,:);
    r(1,:)=abs(gg(1,2,:));
    g(2,:)=gg(2,2,:);
    r(2,:)=abs(gg(2,1,:));
for l=1:2
    axes('pos',[0.05+0.5*(l-1) 0.2 0.4 0.7]);
    for n=1:b+1
        t=0:.01:2*pi;
        x(n,:)=r(l,n).*cos(t)+real(g(l,n));
        y(n,:)=r(l,n).*sin(t)+imag(g(l,n));
        plot(x(n,:),y(n,:));
        hold on;
    end;
grid;
plot(g(l,:)); 
end;