f_bode.m

控制系统计算机辅助设计——MATLAB语言与应用(源代码）

function H=f_bode(a,eta,b,gamma,w,Gc)
for i=1:length(w)
    P=b*((sqrt(-1)*w(i)).^gamma.'); Q=a*((sqrt(-1)*w(i)).^eta.');
    H(i)=P/Q;  % 注意上面的 .'，表示直接转置，不能是共轭转置
end
H=frd(H,w); if nargin==6, H=H*Gc; end
if nargout==0, bode(H); end