函数图(complex notation)
如果有复数的话怎么出图??比如说
有个这样的函数
G=(1-r^2-i2s*r)/((1-r^2)^2+(2s*r^2))
use complex notation 可以把G的实部与虚部分开画,或者画G的幅度和相位 figure(1);clf
zeta=; r=;
for i=1:length(zeta),
ReG=(1-r.^2)./((1-r.^2).^2+(2*zeta(i)*r).^2);
ImG=-(2*zeta(i).*r)./((1-r.^2).^2+(2*zeta(i)*r).^2);
figure(1); plot(r,ReG); plot(r,ImG);hold on
end
figure(1)
title('Magnitude of the Frequency Response')
xlabel('\omega/\omega_n'); ylabel('|G(i\omega)|'); grid
是这样么?但是分开画的话不是共14个图么? 需要的是7个图线,分开画之后画不出来要求的曲线,在线求解~~~
[ 本帖最后由 ChaChing 于 2010-3-29 00:32 编辑 ] 可以这样画在一个图:
plot3(r,ReG,ImG)
or
plot(r,abs(G),r,angle(G)) LZ这个看来好像是要画FRF!? 直觉LZ程序满多问题的
FRF好像很少人画real/imag!?
看看函数subplot, real, imag, freqresp, bode, nyquist ...等, 有无用
[ 本帖最后由 ChaChing 于 2010-3-25 09:12 编辑 ] 恩,谢谢~
貌似是虚数画图的话也可以求虚数的大小图。
意思就是说√(a^2+b^2)=lGl
页:
[1]