关于课程设计1.docVIP

巴特沃斯模拟地通数字滤波器 %声音信号 [x,fs,bits]=wavread(D:\wode\shangbuqi);%采集语音信号 sound(x,fs,bits); N=length(x); X=fft(x,N);%求原始信号频谱 n=0:N-1; subplot(421);plot(x);grid; title(原波形图); subplot(422);plot(abs(X));grid; title(原频谱图) %加噪音 y=x+0.01rand(N,1);%加信噪比为13的白噪声 sound(y,fs,bits); subplot(423);plot(y);grid; title(加噪声后的波形图); Y=fft(y,N);%对加噪后的信号求频谱 subplot(424);plot(abs(Y));grid; title(加噪声后频谱图); %滤波器 Fs=8000; %滤波器指标 fp=1200;wp=2*pi*fp/Fs; fs=1600;ws=2*pi*fs/Fs; wap=2*Fs*tan(wp/2); was=2*Fs*tan(ws/2); Rp=1;As=100; [M,wc]=buttord(wap,was,Rp,As,s);%求阶数 [z,p,k]=buttap(M);%求零极点 [bp,ap]=zp2tf(z,p,k);%求模拟滤波器系数 [bs,as]=lp2lp(bp,ap,wap);%去归一化 [bz,az]=bilinear(bs,as,Fs/2);%用双线性变换法转为数字系数 [h,w]=freqz(bz,az);%计算数字滤波器频响 subplot(425); plot(w,abs(h));grid; title(巴特沃斯数字低通滤波器特性); %滤波 l=filter(bz,az,y);%对原声信号进行过滤 sound(l,fs,bits); subplot(426);plot(l);grid; title(滤波后波形图); L=fft(l,N);%对过滤后的声音信号进行傅里叶变换求幅频响应 subplot(427);plot(abs(L));grid; title(滤波后图谱); 二．椭圆数字高通滤波器 %声音信号 [x,fs,bits]=wavread(D:\wode\shangbuqi);%采集原始语音信号 sound(x,fs,bits); N=length(x); X=fft(x,N);%对原始信号进行傅里叶变换求频谱 n=0:N-1; subplot(421); plot(x);grid; title(原信号波形图); subplot(422); plot(abs(X));grid; title(原信号频谱图) %加噪音 y=x+0.01*rand(N,1);%对原始语音信号加噪 sound(y,fs,bits); subplot(423); plot(y);grid; title(加噪声后的波形图); Y=X+0.01*(N,1); subplot(424); plot(abs(Y));grid; title(加噪声后频谱图); %加滤波器 fp=5000;Rp=1;%指标 fs=4800;As=100; Fs=16000; wp=2*fp/Fs; ws=2*fs/Fs; [n,wn]=ellipord(wp,ws,Rp,As);%求椭圆数字滤波器阶数和通带截止频率 subplot(425); [B,A]=ellip(n,Rp,As,wn,high);%求滤波器系统函数的系数 freqz(B,A);%幅频响应 %滤波 h=filter(B,A,y);%对原始信号滤波 sound(h,fs,bits); subplot(426); plot(h);grid; title(滤波后信号波形图); H=fft(h,N);%对滤波后的原始信号求频响 subplot(427); plot(abs(H));grid; title(滤波后频谱图); 声音信号 三．布莱克曼窗函数滤波器 %声音信号 [x,fs,bits]=wavread(D:\wode\shangbuqi);%采集语音信号 sound(x,fs,bits); N=length(x); X=fft(x,N);%原信号求频谱 n=0:N-1; subplot(421); plot(x);grid; title(原信号波形图); subplot(422); plot(abs(X));grid; title(原信号频谱图) %加噪音 y=x+0.01rand(N,1);%原信号加噪 sound(y,fs,bits); subplot(423); plot(y);grid; title(加噪声后的波形图);