资源描述：
《数字信号处理matlab课程设计》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、课题一数字信号处理系统设计一、项目要求用本课程所学的数字信号处理理论知识，设计一个具有信号的采集、处理、传输、显示和存储等功能的系统，内容如下：１、录制一段语音信号，并对录制的语音信号进行采样（采样频率可取fs=22050Hz）;２、画出采样后的语音信号的时域波形和频谱图；３、滤波器的性能指标：低通滤波器：通带边界频率ｆｐ＝１ｋＨｚ，通带最大衰减Ａｐ＝１ｄB； 阻带边界频率ｆｐ＝１.2ｋＨｚ，阻带最小衰减Ａｐ＝１00ｄB；高通滤波器：通带边界频率ｆｐ＝5ｋＨｚ，通带最大衰减Ａｐ＝１ｄB； 阻带边界频率ｆｐ＝4.
2、8ｋＨｚ，阻带最小衰减Ａｐ＝１00ｄB；带通滤波器：通带上限截止频率fp2=3kHz,通带下限截止频率fp1=1.2kHz;阻带上限截止频率fs2=3.2kHz,通带下限截止频率fp1=1kHz;通带最大衰减Ap=1dB,阻带最小衰减As=100dB;采用双线性变换法设计滤波器，并画出滤波器的频率响应；４、用自己设计的滤波器对采样的信号进行滤波，画出滤波后信号的时域波形和频谱图，并对滤波前后的信号进行对比，分析信号的变化；５、回放语音６、用ＧＵＩ设计一个信号系统的用户界面。二、实验所要用到的ＭＡＴＬＡＢ函数１、语音信号的采
3、样与播放wavread();[y，fs，bite]=wavread（）；%语音信号的采样sound（y，fs，bite）；%播放语音2.滤波器：IIR：butte（）；%巴特沃思滤波器cheby1();%切比雪夫I滤波器elliptical（）；%椭圆滤波器3.频率响应：[h，f]=freqz（b,a,n,fs）freqz(b,a,n,fs)5.快速傅里叶变换fft(x,n)6.画曲线plot(x,y)stem(x,y)7.在MATLAB中，设计辅助低通原型巴特沃思和切比雪夫滤波器的阶数和截止频率；1）利用buttord和
4、cheblord确定阶数；2）[num，den]=butter（N，Wn），[num，den]=cheby1（N，Wn）3）lp2hp,lp2bp,lp2bs可以完成低通滤波器到高通，带通，带阻的转换4）使用biliner对模拟滤波器进行双线性变换，求得数字滤波器的传输函数系数三、数字滤波器（代码）1.IIR低通滤波器fp=1000;ft=5500;fs=1200;wp=2*pi*fp/ft;ws=2*pi*fs/ft;Fp=2*ft*tan(wp/2);Fs=2*ft*tan(ws/2);[n,Omgc]=buttord
5、(Fp,Fs,1,100,'s');[z,p,k]=buttap(n);B=k*real(poly(z));A=real(poly(p));[b1,a1]=lp2lp(B,A,Omgc);[ba1,aa1]=bilinear(b1,a1,ft);%模拟转数字[Ha,w]=freqz(ba1,aa1);%求频率相应plot(w*ft/(2*pi),20*log10(abs(Ha)));title('IIR低通滤波器');xlabel('频率/HZ');ylabel('幅值');1.IIR高通滤波器fp=5000;ft=250
6、00;%取抽样频率 fs=4800;Rp=1;As=100;wp=2*pi*fp/ft;%通带频率 ws=2*pi*fs/ft;%截止频率 Fp=2*ft*tan(wp/2);Fs=2*ft*tan(ws/2);[n,Omgc]=ellipord(Fp,Fs,Rp,As,'s');%计算阶数n和截止频率 [z,p,k]=ellipap(n,Rp,As);B3=k*real(poly(z));A3=real(poly(p));[b3,a3]=lp2hp(B3,A3,Omgc);[b
7、a3,aa3]=bilinear(b3,a3,ft);[Ha,w]=freqz(ba3,aa3);plot(w*ft/(2*pi),20*log10(abs(Ha)));title('IIR高通滤波器');xlabel('频率/HZ');ylabel('幅值');1.IIR带通滤波器fp1=1200;fp2=3000;fs1=1000;fs2=3200;ft=10000;As=100;Rp=1;wp1=2*pi*fp1/ft;wp2=2*pi*fp2/ft;ws1=2*pi*fs1/ft;ws2=2*pi*fs2/ft;F
8、p1=2*ft*tan(wp1/2);Fp2=2*ft*tan(wp2/2);Fp=[Fp1,Fp2];Fs1=2*ft*tan(ws1/2);Fs2=2*ft*tan(ws2/2);Fs=[Fs1,Fs2];bw=Fp2-Fp1;w0=sqrt(Fp1*Fp2);%通带宽和中心频率
