matlab噪声滤波分析,应用Matlab对含噪声语音信号进行频谱分析及滤波

2、加入噪声

y=wavread('C:\\Users\\acer\\Desktop\\');

y=y(:,1); %矩阵维度置换成1 subplot(2,2,1); plot(y);

title('加噪前的时域曲线');

N=length(y)-1; n=0:1/fs:N/fs;

x=*sin(40*pi*n); %设置正弦噪声信号 x=x.'; %置换成矩阵 z=y+x; %添加噪声 subplot(2,2,2); plot(z);

title('加噪后的时域曲线'); sound(z,40000)

3、IIR滤波器设计

clc;clear;close all;

fs=22050;x1=wavread('C:\\Users\\acer\\Desktop\\'); x1=x1(:,1)

t=0:1/22050:(size(x1)-1)/22050; %设置并添加噪声信号 d=*sin(40*pi*t); d=d.'; x2=x1+d;

wp=*pi;ws=*pi;Rp=1;Rs=15; %通阻带截止和通阻带衰减 Fs=22050;Ts=1/Fs;

wp1=2/Ts*tan(wp/2); %将模拟指标转换成数字指标 ws1=2/Ts*tan(ws/2);

[N,Wn]=buttord(wp1,ws1,Rp,Rs,'s'); %选择滤波器的最小阶数 [Z,P,K]=buttap(N); %创建butterworth模拟滤波器

[Bap,Aap]=zp2tf(Z,P,K); [b,a]=lp2lp(Bap,Aap,Wn);