资源描述：
《合工大自动化信号分析与处理实验报告》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、合肥工业大学电气与自动化工程学院实验报告自动化专业班级学号姓名口期—指导教师黄云志共—页第—页成绩实验一含噪声语音信号频谱分析及滤波实验报告要求：1、实验报告包括四部分：实验原理、实验内容、实验程序、结果分析；分别占实验报告总成绩的10%,10%,40%,40%；2、实验程序及结果分析如内容雷同，均不给分；一、实验原理含噪声语音信号通过低通滤波器，高频的噪声信号会被过滤掉，得到清晰的无噪声语音信号。二、实验内容录制一段个人口己的语音信号，并对录制的信号进行采样；画出采样后语音信号的时域波形和频谱图；在语咅信号中增加正弦噪声信号(口己设置儿个频率
2、的正弦信号)，对加入噪声信号后的语咅信号进行频谱分析；给定滤波器的性能指标，采用窗函数法和双线性变换设计数字滤波器，并画出滤波器的频率响应；然示用白己设计的滤波器对采集的信号进行滤波，画出滤波后信号的时域波形和频谱，并对滤波前后的信号进行对比试听，分析信号的变化。给出数字低通滤波器性能指标:如，通带截止频率fp=10000Hz,阻带截止频率fs=12000Hz(可根据自己所加入噪声信号的频率进行阻带截止频率设置)，阻带最小衰减Rs=50dB,通带最大衰减Rp=3dB(也可白己设置)，采样频率根据自己语音信号采样频率设定。三、实验程序1.原始信号
3、fs二22050;%设逍采样频率xl=wavread(?C:UsersmDocumentsFFOiitpiit无标题.wav');%读取声音sound(xl,fs);%播放原始声音yl=fft(xl,length(xl));%FFTf=fs*(0:511)/1024;figure(1);%建立图形plot(xl);%原始信号时域图形title('原始语音信号时域波形')xlabelC时间’)；ylabel(，幅值');axis([0,5*10^4,-2,21);figure(2);plot(f,abs(yl(l:512)))%做原始语
4、音信号的FFT频谱图xlabel('频率');ylabel('幅度');titleC原始语音信号FFT频谱幅频')谡音苗号时■渎博1.加噪声clearallfs二22050;%设置采样频率xl=wavread(，C:UsersmDocumentsFFOutput无标题.wav');%读取声音f=fs*(0:511)/1024;t=0:1/22050:(length(xl)-l)/22050;d二[0.5*cos(2*pi*10000*t)]';%加噪声x2二xl+d;sound(x2,22050);%播放加噪音后的信号y2=fft(x
5、2);figure(3)Plot(t,x2)%加噪信号时域图形titleC加噪声后的时域信号波形')；xlabelC时间');ylabel('幅值');figure(4)subplot(211);plot(abs(xl));titleC原始语音信号频谱')；xlabel('频率’);ylabelC幅值');subplot(212);plot(abs(y2));titleC加噪声后的信号频谱')；xlabel('频率’);ylabelC幅值');6000wIfII—厂1L-厂「1111■—-——ii.1——,亠」QCLZaJICL5Q-S1IN”
6、41-ttL8er司1.设计MR滤波器clearallfs=22050;wp=O.4*pi;ws二0.7*pi;rp=3;rs=50;xl=wavread('C:UscrsmDocumcntsFFOutput无标题.wav’);%读取声iff=fs*(0:22499)/22500:%F=([l:N]-1)*Fs/N;%换算成实际的频率值t=0:1/22050:(length(xl)-1)/22050;%>ir义噪声信号d二[0.7*sin(3*pi*10000*t)]';x2二xl+d;%加噪声[n,wn]=buttord(wp/pi
7、,ws/pi,rp,rs);%设il'llR滤波器[bz,az]=butter(n,wn/pi);[h,w]=freqz(bz,az,512,fs);%滤波器幅频响应x3二fi1ter(bz,az,x2);%滤波y0=fft(x2,22500);yl=fft(x3,22500);sound(x3,22500)期播放滤波后的声曾figure(1)plot(abs(h));%滤波器幅频响应图title(s滤波器幅频响应图’)figure(2)subplot(2,1,1)plot(t,x2)%画出滤波而的时域图titleC滤波前的时域波形')；sub
8、plot(2,1,2)plot(t,x3);titleC滤波后的时域波形');figure(3)subplot(211)plot(f,abs(yO(1
