基于MATLAB的数字信号滤波处理

1课程设计内容求 1 1 12课程设计理论基础 1 1 2 33课程设计的具体实现及仿真结果分析 43.1语音信号采集 4 4 4 5 5 6 7 8 8 8 9 104设计总结 1121课程设计内容及要求集，将语音信号转换成计算机能够运算的有限长序列，加一定的噪声，并选②运用wavread()函数对采集的语音信号进行时频分析③运通rand()或randn()函数给语音信号加一定的噪声信号，并且对加fp2=3000Hzfc2=3200Hz⑤运用fftfilt()函数对加噪语音信号进行滤波处理有限长单位冲激响应（FIR）滤波器有以下3dj负)ddj负)j负n4d5这种窗函数是近似于给定旁办电平，使主办具有最大能量意义下的最佳窗函3课程设计的具体实现及仿真结果分析3.1语音信号采集用录音设备录下一段音频以yr.wav格式存储在电脑的中，以便接下的过给某一向量。再对其进行采样，记住采样频率和采样点数。接下来，对语音Y=fft(y,N);%傅里叶plot(y);title('原始信号plot(abs(Y));6Noise=0.01*randn(N,2);plot(Si);title('加噪语音信S=fft(Si);7），p=1-10.^(-rp/20);fpts=[wpws];8b21=fir1(n21,wn21,Kaplot(w/pi,abs(h));fp2=3000Hzfc2=3200Hz[n,wn,bta,ftype]=kaiserb22=fir1(n,wn,ftype,kaiser(n+1,bta),'noscale');plot(w1/pi,20*log10(abs(hh1)))xlabel('归一化频率w');ylabel('幅度/db');9p=1-10.^(-rp/20);%通带阻带波纹wp=0.9;ws=0.7;fpts=[wswp];[n23,wn23,beta,ftype]=kaiserord(fpts,mag,dev);b23=fir1(n23,wn23,'highplot(w*12000*0.5/pi,abs(h));noise=0.01*randn(n,2);z21=fftfilt(b21,s);m21=fft(z21);plot(abs(S),'g');%noise=0.01*randn(n,2);z22=fftfilt(b22,s);m22=fft(z22);plot(abs(S),'g');%高通滤波器*****************************noise=0.01*randn(n,2);z23=fftfilt(b23,s);m23=fft(z23);本次课程设计完成了基于matlab的语音信号滤波处理，用wavread()函行加噪处理并且进行频谱分析，最后用窗函数法设计FIR低通、高通、带通数字滤波器，分别对加噪后的语音信号进行滤波处理，得到滤波处图。本次课程设计基本完成了设计的指标和要求，在课程设计过程 ()函数进行采集后，发现没有时频分析图，初步判断可能是音乐文件过大，后经剪切音乐文件后，再用wavread()函数采集，就可以得到时频分析图了，在设计FIR滤波器的过程中，部分程序参考了网上的设计程序，但是当用文件，以及书本的