信号与系统实验报告（川大版）

title('幅值频谱'); xlabel('频率(HZ)'); ylabel('幅值'); subplot(2,1,2); plot(f,angle(Y12)); title('相位频谱'); xlabel('频率(HZ)'); ylabel('相位'); figure; subplot(2,1,1); plot(f,abs(Y13)); title('幅值频谱'); xlabel('频率(HZ)'); ylabel('幅值'); subplot(2,1,2); plot(f,angle(Y13)); title('相位频谱'); xlabel('频率(HZ)'); ylabel('相位'); 9

10

fs=128; t=1/fs; L1=256; n=0:1:L1-1; h2=[-0.0238,0.0562,-0.0575,-0.1302,0.5252, -0.6842,-0.3129,5.6197,5.6197,-0.3129,-0.6842,0.5252,-0.1302,-0.0575,0.0562,-0.0238]; x1=cos(0.25*pi*n); x2=cos(1.25*pi*n); x3=x2+x1; N = 2^nextpow2(L1); H2=fft(h2,N)/N*2; Y21=fft(x1,N)/N.*2.*fft(h2,N)/N.*2; Y22=fft(x2,N)/N.*2.*fft(h2,N)/N.*2; Y23=fft(x3,N)/N.*2.*fft(h2,N)/N.*2; f=fs/N*(0:1:N-1); figure; subplot(2,1,1); plot(f,abs(H2)); title('幅值频谱'); xlabel('频率(HZ)'); ylabel('幅值'); subplot(2,1,2); plot(f,angle(H2)); title('相位频谱'); xlabel('频率(HZ)'); ylabel('相位'); figure; subplot(2,1,1); plot(f,abs(Y21)); title('幅值频谱'); xlabel('频率(HZ)'); ylabel('幅值'); subplot(2,1,2); plot(f,angle(Y21)); title('相位频谱'); xlabel('频率(HZ)'); ylabel('相位'); figure; subplot(2,1,1); plot(f,abs(Y22)); title('幅值频谱'); xlabel('频率(HZ)'); ylabel('幅值'); 11

subplot(2,1,2); plot(f,angle(Y22)); title('相位频谱'); xlabel('频率(HZ)'); ylabel('相位'); figure; subplot(2,1,1); plot(f,abs(Y23)); title('幅值频谱'); xlabel('频率(HZ)'); ylabel('幅值'); subplot(2,1,2); plot(f,angle(Y23)); title('相位频谱'); xlabel('频率(HZ)'); ylabel('相位'); 12