海南大学信号与系统实验报告

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;