freqm1.m

华东理工大学自动化系《控制系统分析、设计和应用》教材配套程序实例,有关说明： 1. 安装MATLAB软件。安装时

function [Gmm,Pmm,Wcgm,Wcpm,S]=freqm1(G,m,varargin);
% 输入信号：G：系统；
%           m：衰减度；m=0.221表示4：1（约定值）；m=0.433表示10：1；
%           VARARGIN：频率响应图类型；不输入时绘制伯德图
%           类型有：nyquist,bode,nichols,margin
% 输出信号：无：绘制有关频率响应图；nyquist,bode,nichols图。
%           四个：得到稳定裕度和交角频率；VARARGIN必须是'margin'。
% 示例：G=zpk([],[-1 -2 -3],6);freqm(G,0.221,'nyquist');

% Copyright 2001-2002 ECUST. 
%   $Revision: 1.2 $  $Date: 2002/02/12 15:38:45 $

syms s nums dens denm numm den num w;
if nargin==1,disp('Bode Diagram,with m=0.221');plottype='bode';m=0.221;end;
if nargin==2,plottype='bode';disp('Bode Diagram');end; 
if nargin> 2,plottype= varargin {1};end;
G=tf(G);num=poly2sym(G.num{1},s);den=poly2sym(G.den{1},s);
[MAG,PHASE,ww] = bode(G);DT=get(G,'inputdelay');
wM=max(ww);wm=min(ww);w=logspace(log10(wm),log10(wM),100);
numm=subs(num,s,(j-m)*w);denm=subs(den,s,(j-m)*w);
Gg=numm./denm.*exp(m*w*DT);Gm=Gg./(1+Gg);sysm=frd(Gm,w);set(sysm,'inputdelay',DT);
[Gmm,Pmm,Wcgm,Wcpm]=margin(sysm);