数字信号处理课程设计(上传).docVIP

二○一*～二○一*学年第 一 学期 信息科学与工程学院 课程设计报告书 课程名称： 数字信号处理课程设计 班 级： 学 号： 姓 名： 指导教师： 二○一* 年 一月 一、需求分析：的正弦型干扰信号。语音信号的采样频率为。要求通过计算机仿真对滤波器的性能进行分析 (2) 针对含有5Hz、15Hz和30Hz的混合正弦波信号，试设计滤波系统，5Hz和30Hz的正弦分量，阻带的最小衰减不小于50dB。要求：确定采样频率、采样长度，对进行采样；数字滤波器； 通过计算机仿真对滤波器的性能进行分析设计设计滤波系统，混合正弦波信号5Hz和30Hz的正弦分量 解：①确定采样频率、采样长度，对进行采样 详细设计x=sin(2*pi*t*5)+sin(2*pi*t*15)+sin(2*pi*t*30); figure(1) subplot(211) plot(t,x);title(原始信号的时域波形图); axis([0 1 -3 3]);grid on; xlabel(时间/s); subplot(212) Y=fft(x,N); F=[0:N-1]*fs/N; plot(F,abs(Y));%原始信号的频谱 title(原始信号的频谱); grid on; xlabel(频率/Hz) Wp = [14 16]/(fs/2);%带通滤波器的通带截止频率 Ws = [12 18]/(fs/2);%带通滤波器的阻带截止频率 Ap = 3;%带通滤波器的通带最大衰减 As = 50;%带通滤波器的阻带最小衰减 [N1 Wn] = buttord(Wp,Ws,Ap,As);%产生滤波器阶数和截止频率 [B,A] = butter(N1,Wn);%产生系统函数系数 [h,w]=freqz(B,A);%滤波器的频率响应 figure(2); subplot(211); plot(w*fs/(2*pi),20*log(h)); grid on; title(IIR滤波器幅频特性); xlabel(频率/Hz); subplot(212); plot(w*fs/(2*pi),angle(h)); grid on; title(IIR滤波器相频特性); xlabel(频率/Hz); y = filter(B,A,x);%原始信号经滤波器后的信号 figure(3); subplot(211); plot(t,y);%经滤波后的信号的时域波形图 title(IIR滤波后信号的频谱); grid on; axis([0 5 -1.5 1.5]); xlabel(时间/s); Y1=fft(y,N); subplot(212); plot(F,abs(Y1));%经滤波后的信号的频谱 title(IIR滤波后信号的频谱); xlabel(频率/Hz);grid on; Wc=[14 16]/(fs/2);%FIR截止频率 b=fir1(127,Wc);%采用默认哈明窗，经计算127阶，产生滤波器系数 figure(4) [h1,w1]=freqz(B,A)%滤波器的频率响应 y2=filter(b,1,x); Y2=fft(y2,N);%滤波信号的FFT figure(4) subplot(211) plot(w1*fs/(2*pi),20*log(abs(h1))) title(FIR滤波器幅频特性) xlabel(频率/Hz); axis([-5 115 -1200 250]) grid on; subplot(212) plot(w1*fs/(2*pi),angle(h1)) title(FIR滤波器相频特性) xlabel(频率/Hz) grid on; figure(5) subplot(211) plot(t,y2);grid on; title(FIR滤波后信号的时域波形图) axis([0 5 -1.5 1.5]); subplot(212) plot(F,abs(Y2)) title(FIR滤波后信号的频谱图) grid on; 四、调试分析调试存在问题总结包括课程设计过程的收获、课程的认识等内容需求分析：设计详细设计调试分析测试数据测试输出结果调试存在问题