examp4_14.m

来自「matlab的大量实例代码」· M 代码 · 共 6 行

[Pxx,f]=psd_estm(y,2028,T);
num=G.num{1}.*(-1).^[length(G.num{1})-1:-1:0];
den=G.den{1}.*(-1).^[length(G.den{1})-1:-1:0];
GG=3*G*tf(num,den); % 计算 GG(s)=P G(s)G(-s)
[mag,p]=bode(GG,f); % 计算 GG(jw)
semilogx(f,20*log10(mag(:)),f,20*log10(Pxx))