《现代信号分析与处理》课终作业(20090709)

《《现代信号分析与处理》课终作业(20090709)》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、《现代信号分析与处理》课终作业姓名：司爱威年级：2009级博士研究生学号：111200911005专业：车辆工程2009-7-1610一、线性卷积程序1.1子函数程序（matlab语言）functiony=saw_fun_conv(x,h)%x,h分别表示两个要卷积的信号Nx=length(x);Nh=length(h);Ny=Nx+Nh-1;fori=Nx+1:Nyx(i)=0;endfori=Nh+1:Nyh(i)=0;endfori=1:Nysum=0;forj=1:isum=sum+x(j)*h(i-j+1);endy(i)=sum;end1.2验证程序

2、x=[1234];h=[11];y=saw_fun_conv(x,h)subplot(2,2,1);stem(x,'bo');xlabel('x(k)');gridon;subplot(2,2,3);stem(h,'bo');axis([1405]);xlabel('h(k)');gridon;subplot(1,2,2);stem(y,'bo');xlabel('y(n)');gridon;1.3验证波形如图1.1所示。其中x=[1234];h=[11];y=[13574].10图1.1线性卷积程序验证波形10二、DFT程序2.1子函数程序functionX=

3、saw_fun_dft(x)%x表示原信号N=length(x);W=exp(-j*2*pi/N);fork=1:Nsum=0;forn=1:Nsum=sum+x(n)*W.^((n-1)*(k-1));endXa(k)=abs(sum);endfori=1:N/2X(i)=Xa(i);end2.2验证程序N=256;f1=0.1;f2=0.2;fs=1;a1=5;a2=3;w=2*pi/fs;x=a1*sin(w*f1*(0:N-1))+a2*sin(w*f2*(0:N-1))+randn(1,N);X=saw_fun_dft(x);subplot(2,1,1

4、);plot(x);xlabel({'n','(a)'});ylabel('x(n)');gridon;f=0:1/N:fs/2-1/N;subplot(2,1,2);plot(f,X);xlabel({'f/Hz','(b)'});ylabel('X(f)');gridon;2.3验证波形如图2.1所示。其中x=5sin(0.2n)+3sin(0.4n)+randn(1,256),n=0,2,…255图（a）为原信号波形，图（b）为DFT变换后的波形。10图2.1DFT程序验证波形10三、窗函数法设计FIR滤波器程序3.1子函数程序functionh=saw_

5、fun_firwin(N,fc)%低通滤波器，选用汉明窗，N表示滤波器阶数，fc表示截止频率。ifrem(N,2)==0n2=N/2-1;flag=1;elsen2=N/2;flag=0;endwc1=2*pi*fc;fori=0:n2s=i-N/2;h(i+1)=(sin(wc1*s)/(pi*s))*(0.54-0.46*cos(2*i*pi/(N+1)));h(N-i+1)=h(i+1);endif(flag==1)h(N/2+1)=wc1/pi;end3.2验证程序N=60;fc=0.15;f1=0.1;f2=0.2;fs=1;a1=5;a2=3;w=2

6、*pi/fs;x=a1*sin(w*f1*(0:N-1))+a2*sin(w*f2*(0:N-1))+randn(1,N);X=saw_fun_dft(x);%利用第2题DFT程序N2=20;h=saw_fun_firwin(N2,fc)x2=saw_fun_conv(x,h);%利用第1题线性卷积程序X2=saw_fun_dft(x2);%利用第2题DFT程序subplot(321);stem(h,'.');xlabel({'n','(a)'});ylabel('h(n)');gridon;subplot(322);f11=0:0.5/N:0.5-0.5/N;

7、plot(f11,abs(freqz(h,1,N)));xlabel({'f','(b)'});ylabel('H(f)');gridon;subplot(323);plot(x);xlabel({'n','(c)'});ylabel('x(n)');gridon;subplot(324);10f12=0:1/N:0.5-1/N;plot(f12,X);xlabel({'f','(d)'});ylabel('X(f)');gridon;subplot(325);plot(x2);xlabel({'n','(e)'});ylabel('x2(n)');gridon

8、;subplot(326