发布时间 : 星期一 文章海南大学信号与系统实验报告更新完毕开始阅读a6189518c5da50e2524d7f50
7-1 clear syms t;
f=exp(-2*t)*sym('Heaviside(t)')+exp(2*t)*sym('Heaviside(-t)'); subplot(3,1,1);
ezplot(f);ylabel('f(t)'); axis([-3,3,0,1.2]); F=fourier(f); subplot(3,1,2);
ezplot(abs(F));title('幅度频谱'); im=imag(F); re=real(F);
phase=atan(im/re); subplot(3,1,3);
ezplot(phase);title('相位频谱');
7-2 clear syms t w;
F=2*sin(w)/w;
subplot(2,2,1);ezplot(F); f=ifourier(F,t);
subplot(2,2,2);ezplot(f);title('f(t)');axis([-1 1.0 0 1.5]) subplot(2,2,3);ezplot(abs(F));title('幅度频谱'); im=imag(F); re=real(F);
phase=atan(im/re);
subplot(2,2,4);ezplot(phase);title('相位频谱');
7-3
w=linspace(0,10,200); b=[1];
a=[0.08 0.4 1]; [h,w]=freqs(b,a,w); subplot(2,1,1); plot(w,abs(h));
set(gca,'xtick',[0 1 2 3 4 5 6 7 8 9 10]); set(gca,'ytick',[0 0.4 0.707 1]); grid on;
subplot(2,1,2); plot(w,angle(h));
set(gca,'xtick',[0 1 2 3 4 5 6 7 8 9 10]); grid on;
7-4
b=[1 0 22500]; a=[1 200 20000]; [h,w]=freqs(b,a,100); subplot(211); plot(w,abs(h)); hold on; grid on;
subplot(212); plot(w,angle(h)); grid on; hold on;