信号与线性系统课程设计报告

信号与线性系统课程设计报告课题名称：数字语音信号分析系统设计%---ExecutesonbuttonpressinFFT.functionpushbutton3_Callback(hObject,eventdata,handles)figure(2);[x,fs,nbits]=wavread('G:\design\qian1.wav');subplot(1,1,1)N=100;n=0:N-1;t=0.01*n;q=n*2*pi/N;X=fft(x,N);plot(q,abs(X));title('FFT')%---Executesonbuttonpressinjiaruzaosheng.functionpushbutton4_Callback(hObject,eventdata,handles)figure(3);[x,fs,nbits]=wavread('G:\design\qian1.wav');subplot(2,1,1)N=100;n=0:N-1;t=0.01*n;q=n*2*pi/N;X=fft(x,N);plot(q,abs(X));title('FFT')subplot(2,1,2)N=100;n=0:N-1;t=0.01*n;M=length(x);k=0:M-1;q=n*2*pi/N;y=x+0.05*cos(1.5*pi*k)';Y=fft(y,N);plot(q,abs(Y));title('FFT(withnoise)')%---ExecutesonbuttonpressinButterworthFilter.functionpushbutton7_Callback(hObject,eventdata,handles)figure(4);%ButterworthFilterwp=0.25*pi;ws=0.35*pi;Rp=1;As=17;T=1;OmegaP=wp/T;OmegaS=ws/T;[cs,ds]=afd_butt(OmegaP,OmegaS,Rp,As);[b,a]=imp_invr(cs,ds,T)[db,mag,pha,w]=freqz_m(b,a);subplot(2,1,1);plot(w/pi,mag);title('MagnitudeResponse')axis([0,1,0,1.1])subplot(2,1,2);plot(w/pi,db);title('MagnitudeinDB')axis([0,1,-40,5]);%---ExecutesonbuttonpressinChebyshevFilter.functionpushbutton8_Callback(hObject,eventdata,handles)figure(5);%ChebyshevFilterwp=0.25*pi;ws=0.35*pi;Rp=1;As=17;T=1;OmegaP=(2/T)*tan(wp/2);OmegaS=(2/T)*tan(ws/2);ep=sqrt(10^(Rp/10)-1);Ripple=sqrt(1/(1+ep*ep));Attn=1/(10^(As/20));A1=1/Attn;a1=sqrt(A1*A1-1)/ep;a2=OmegaS/OmegaP;N=ceil(logm(a1+sqrt(a1*a1-1))/logm(a2+sqrt(a2*a2-1)));fprintf('\nChebyshevFilterOrder=%2.0f\n',N)[cs,ds]=u_chb1ap(N,Rp,OmegaP);[B,A]=bilinear(cs,ds,T)[db,mag,pha,w]=freqz_m(B,A);subplot(2,1,1);plot(w/pi,mag);title('MagnitudeResponse');axis([0,1,0,1.1])subplot(2,1,2);plot(w/pi,db);title('MagnitudeinDB');axis([0,1,-40,5]);%---ExecutesonbuttonpressinButterworth.functionpushbutton9_Callback(hObject,eventdata,handles)figure(6);[x,fs,nbits]=wavread('G:\design\qian1.wav');subplot(3,1,1)N=100;n=0:N-1;t=0.01*n;q=n*2*pi/N;X=fft(x,N);plot(q,abs(X));title('FFT')subplot(3,1,2)N=100;n=0:N-1;t=0.01*n;M=length(x);l=0:M-1;q=n*2*pi/N;y=x+0.05*cos(1.5*pi*l)';Y=fft(y,N);plot(q,abs(Y));title('FFT(withnoise)')%ButterworthFilterwp=0.25*pi;ws=0.35*pi;Rp=1;As=17;T=1;OmegaP=wp/T;OmegaS=ws/T;[cs,ds]=afd_butt(OmegaP,OmegaS,Rp,As);[b,a]=imp_invr(cs,ds,T)[db,mag,pha,w]=freqz_m(b,a);subplot(3,1,3)y2=filter(b,a,y)N=100;n=0:N-1;t=0.01*n;M=length(x);q=n*2*pi/N;Y2=fft(y2,N);plot(q,abs(Y2));%---ExecutesonbuttonpressinChebyshev.functionpushbutton10_Callback(hObject,eventdata,handles)figure(7);[x,fs,nbits]=wavread('G:\design\qian1.wav');subplot(3,1,1)N=100;n=0:N-1;t=0.01*n;q=n*2*pi/N;X=fft(x,N);plot(q,abs(X));title('FFT')subplot(3,1,2)N=100;n=0:N-1;t=0.01*n;M=length(x);l=0:M-1;q=n*2*pi/N;y=x+0.05*cos(1.5*pi*l)';Y=fft(y,N);plot(q,abs(Y));title('FFT(withnoise)')%ChebyshevFilterwp=0.25*pi;ws=0.35*pi;Rp=1;As=17;T=1;OmegaP=(2/T)*tan(wp/2);OmegaS=(2/T)*tan(ws/2);ep=sqrt(10^(Rp/10)-1);Ripple=sqrt(1/(1+ep*ep));Attn=1/(10^(As/20));A1=1/Attn;a1=sqrt(A1*A1-1)/ep;a2=OmegaS/OmegaP;N=ceil(logm(a1+sqrt(a1*a1-1))/logm(a2+sqrt(a2*a2-1)));fprintf('\nChebyshevFilterOrder=%2.0f\n',N)[cs,ds]=u_chb1ap(N,Rp,OmegaP);[b,a]=bilinear(cs,ds,T)[db,mag,pha,w]=freqz_m(b,a);y2=filter(b,a,y)subplot(3,1,3)N=100;n=0:N-1;t=0.01*n;M=length(x);l=0:M-1;q=n*2*pi/N;Y2=fft(y2,N);plot(q,abs(Y2));%---Executesonbuttonpressinshiyucanshujisuan.functionpushbutton14_Callback(hObject,eventdata,handles)[x,fs,nbits]=wavread('G:\design\qian1.wav');N=100;n=0:N-1;t=0.01*n;M=length(x);%duanshinenglianghanshua=0;y2=0;fori=20:50y=x(i);a=a+y2;y2=y^2;i=i+1;endE=ab=0;y3=0;fori=20:50y=x(i);y3=abs(y);b=b+y3;i=i+1;endM=bc=0;y4=0;yf=0;yl=0;sgnf=0;sgnl=0;fori=20:50yf=x(i);ifyf>=0sgnf=1;elsesgnf=0;end;yl=x(i-1);ifyl>=0sgnl=1;elsesgnl=0;end;y4=y4+abs(sgnf-sgnl);c=c+y4/2;i=i+1;endZ=cset(handles.edit1,'string',E);set(handles.edit2,'string',M);set(handles.edit3,'string',Z);functionedit1_Callback(hObject,eventdata,handles)functionedit1_CreateFcn(hObject,eventdata,handles)ifispcset(hObject,'BackgroundColor','white');elseset(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor'));endfunctionedit2_Callback(hObject,eventdata,handles)functionedit2_CreateFcn(hObject,eventdata,handles)ifispcset(hObject,'BackgroundColor','white');elseset(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor'));endfunctionedit3_Callback(hObject,eventdata,handles)functionedit3_CreateFcn(hObject,eventdata,handles)ifispcset(hObject,'BackgroundColor','white');elseset(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor'));end公式：%afd_butt.mfunction[b,a]=afd_butt(Wp,Ws,Rp,As)N=ceil((log10((10^(Rp/10)-1)/(10^(As/10)-1)))/(2*log10(Wp/Ws)));%上条语句为求滤波器阶数：N为整数；%ceil朝正无穷大方向取整；fprintf('\nButterworthFilterOrder=%2.0f\n',N)OmegaC=Wp/((10^(Rp/10)-1)^(1/(2*N)))%求对应于N的3db截止频率；[b,a]=u_buttap(N,OmegaC);%freqz_m.mfunction[db,mag,pha,w]=freqz_m(b,a);[H,w]=freqz(b,a,1000,'whole');%在0-2*pi之间选取N个点计算频率响应H=(H(1:501))';%频率响应w=(w(1:501))';%频率mag=abs(H);%响应幅度db=20*log10((mag+eps)/max(mag));%增益pha=angle(H);%相位%imp_invr.mfunction[b,a]=imp_invr(c,d,T)[R,p,k]=residue(c,d);%部分分式展开p=exp(p*T);%从模拟到数字极点对应关系,部分分式系数相同[b,a]=residuez(R,p,k);%将部分分式的形式变换成多项式之比的形式b=real(b');%求出数字滤波器系数a=real(a');%u_buttap.mfunction[b,a]=u_buttap(N,Omegac);[z,p,k]=buttap(N);%归一化巴特沃思模拟低通滤波器原形%传输函数用极点形式表示p=p*Omegac;%将代入上式，相当于分子乘以，极点乘以k=k*Omegac^N;B=real(poly(z));%poly为构造具有指定根的