>> clear >> a=[1]; >> b=[1 ];
>> w=linspace(-2*pi,2*pi,512); >> h=freqs(a,b,w);
>> subplot(2,1,2),plot(w/pi,angle(h)); >>subplot(2,1,1),plot(w/pi,abs(h));
二
>> clear >> a=[2 0];
>> b=[1 2,10001];
>> w=linspace(-2*pi,2*pi,512); >> h=freqs(a,b,w);
>> subplot(211),plot(w/pi,abs(h)); >> subplot(212),plot(w/pi,angle(h));
>> clear >> syms t;
>> v1=(1+cos(t))*cos(100*t); >> V=laplace(v1); >> syms s;
>> H=2*s/(s^2+2*s+10001); >> V2=V*H;
>> v2=ilaplace(V2); >> ezplot(v2,[*pi,*pi])
因篇幅问题不能全部显示，请点此查看更多更全内容