函数图（complex notation）

kidarllim

如果有复数的话怎么出图？？
比如说
有个这样的函数

G=(1-r^2-i2s*r)/((1-r^2)^2+(2s*r^2))

use complex notation

friendchj

可以把G的实部与虚部分开画，或者画G的幅度和相位

kidarllim

figure(1);clf
zeta=[0.05;0.1;0.15;0.25;0.5;1.00;1.25]; r=[0:0.01:3];
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 编辑 ]

friendchj

可以这样画在一个图：
plot3(r,ReG,ImG)
or
plot(r,abs(G),r,angle(G))

ChaChing

LZ这个看来好像是要画FRF!? 直觉LZ程序满多问题的
FRF好像很少人画real/imag!?
看看函数subplot, real, imag, freqresp, bode, nyquist ...等, 有无用

[ 本帖最后由 ChaChing 于 2010-3-25 09:12 编辑 ]

kidarllim

恩，谢谢~  
貌似是虚数画图的话也可以求虚数的大小图。
意思就是说√（a^2+b^2）=lGl