mr.m

文件是编辑好的m文件

function [Mr,Pr,Wr]=mr(G)
% calculate Mr Pr Wr

% Copyright 2001-2002 ECUST. 
%   $Revision: 1.2 $  $Date: 2001/05/23 21:22:41 $

[mag,phase,w]=bode(G);magg(1,:)=mag(1,:);pha(1,:)=phase(1,:);
[M,i]=max(magg);Mr=20*log10(M);Pr=pha(1,i);Wr=w(i,1);