基于matlab语音处理.zip

clear all close all %----------------------- fs1=10;%fs1=10hz fs2=50;%fs2=50hz; fs=500;%采样率1000hz wp=2*pi*fs1*2;%转换数字频率 ws=2*pi*fs2; rp=3; rs=60; [n,wn]=buttord(wp,ws,rp,rs,'s'); %截止频率wn,阶数n [b,a]=butter(n,wn,'s'); [bz,az]=impinvar(b,a,fs); [h,w]=freqz(bz,az);%连续时间频率响应 figure(1); plot(w*fs*0.5/pi,abs(h)); axis([0,50,0,1.5]); %----------------------- w1=2*pi*fs1/fs;%转换数字频率 w2=2*pi*fs2/fs; N=1024; n=0:N-1; y1=sin(w1.*n);%10hz信号 y2=sin(w2.*n); z=y1+y2; t=n./fs; figure(2); subplot(2,2,1); plot(t,y1,'--b');%10hz信号 axis([0,0.2,-1,1]); xlabel('time'); ylabel('range'); title('原始时间信号1'); grid on; subplot(2,2,2); plot(t,y2,':r');%50hz信号 axis([0,0.1,-1,1]); xlabel('time'); ylabel('range'); title('原始时间信号2'); grid on; subplot(2,2,3); plot(t,z,'--r');%叠加信号 axis([0,0.2,-2,2]); xlabel('time'); ylabel('range'); title('叠加时间信号'); grid on; fz=filter(bz,az,z);%过滤波器 subplot(2,2,4); plot(t,fz); axis([0,0.5,-1,1]); xlabel('time'); ylabel('range'); title('滤波后信号'); grid on; %------频谱分析------ fy1=fft(y1); fy2=fft(y2); fz1=fft(z); fz2=fft(fz); ffs=n./N.*fs; figure(3); subplot(2,2,1); plot(ffs,abs(fy1)); axis([0,0.2*fs,0,1024]); title('信号1频谱'); xlabel('fre/HZ'); ylabel('range'); subplot(2,2,2); plot(ffs,abs(fy2)); axis([0,0.2*fs,0,1024]); title('信号2频谱'); xlabel('fre/HZ'); ylabel('range'); subplot(2,2,3); plot(ffs,abs(fz1)); axis([0,0.2*fs,0,1024]); title('合成信号频谱'); xlabel('fre/HZ'); ylabel('range'); subplot(2,2,4); plot(ffs,abs(fz2)); axis([0,0.2*fs,0,1024]); title('滤波后信号频谱'); xlabel('fre/HZ'); ylabel('range');