<P>我画出来了,大家可以看看,呵呵<BR>w=logspace(-2,4,500);<BR>Y = 0.1:0.005:0.6;<BR>magnitude=zeros(101,500);<BR>for k=1:101,<BR>r=0.005*(k-1)+0.1;<BR>b=[(9e-004*sin(5*r)+0.007*sin(10*r)),(0.0056*sin(10*r)+26.1e-004*sin(5*r)),(216*sin(5*r)+126*sin(10*r))];<BR>a=1.0e+009 * [0.0000, 0.0000, 0.0003, 0.0002, 4.3200];<BR>SYS=tf(b,a);<BR>[mag, phase]=bode(SYS,w);<BR>magdB = 20*log10(mag);<BR>magnitude(k,:) = magdB(1,1,:);<BR>end<BR>w=log10(w);<BR>[w1,Y1] = meshgrid(w,Y);<BR>mesh(w1,Y1,magnitude)</P> |