f_bode.m

来自「包含大量遗传算法程序」· M 代码 · 共 8 行

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