东南大学Matlab作业3.doc

《东南大学Matlab作业3.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、.word可编辑.MatlabWorksheet3PartA1.Usingfunctionconv_m.mtomakeconvolutionbetweenthefollowingtofunctions(xandh):x=[3,11,7,0,-1,7,-5,0,2];h=[11,9,0,-7,-3,2,0-1];nx=[-2:6];nh=[0:7];Plotthefunctionsandconvolutionresults.Answer:home_work_3_A_1.mx=[3,11,7,0,-1,7,-5,0,2];h=[11,9,0,-7,-3,2,0-1];nx=[-2:

2、6];nh=[0:7];[y,ny]=conv_m(x,nx,h,nh);subplot(3,1,1);stem(nx,x);ylabel('x[n]');axis([-610-2020]);subplot(3,1,2);stem(nh,h);ylabel('h[n]');axis([-610-2020]);subplot(3,1,3);stem(ny,y);xlabel('n');ylabel('y[n]');axis([-610-5050]);conv_m.mfunction[y,ny]=conv_m(x,nx,h,nh)%convolutionfunction.专业.专注

3、..word可编辑.%conv_m.mnyb=nx(1)+nh(1);nye=nx(length(x))+nh(length(h));ny=[nyb:nye];y=conv(x,h);end1.Plotthefrequencyresponseoverforthefollowingtransferfunctionbyletting,whereisthefrequency(rad/sample).,withappropriatelabelsandtitle..Answer:home_work_3_A_2.mdelta=0.01;Omega=0:delta:pi;z=exp(j.*O

4、mega);H=z./(z.^2+1.6*z+0.9);subplot(2,1,1);plot(Omega,abs(H),'m','LineWidth',2);xlabel('0

5、H(Omega)

6、');axis([0pi0max(abs(H))]);.专业.专注..word可编辑.subplot(2,1,2);plot(Omega,atan2(imag(H),real(H)),'k','LineWidth',2);xlabel('0

7、uencyresponse');axis([0pi-pipi]);1.Useffttoanalysefollowingsignalbyplottingtheoriginalsignalanditsspectrum..Answer:home_work_3_A_3.m%Examinesignalcomponents%samplingrate1024HzN=1024;dt=1/N;t=0:dt:1-dt;%****definethesignal**x=sin(2*pi*32.*t)+sin(2*pi*137.*t)+sin(2*pi*467.*t);.专业.专注..word可编辑.%

8、***********************subplot(2,1,1);plot(t,x);axis([01-11]);xlabel('nT(seconds)');ylabel('x[n]');X=fft(x);df=1;f=0:df:N-1;subplot(2,1,2);plot(f,abs(X));axis([0N/20120]);xlabel('k(Hz)');ylabel('

9、X[k]

10、');1.UsethefastFouriertransformfunctionffttoanalysefollowingsignal.Plottheoriginalsignal,an

11、dthemagnitudeofitsspectrumlinearlyandlogarithmically.ApplyHammingwindowtoreducetheleakage..ThehammingwindowcanbecodedinMatlabas.专业.专注..word可编辑.forn=1:Nhamming(n)=0.54+0.46*cos((2*n-N+1)*pi/N);end;whereNisthedatalengthintheFFT.Answer:home_work_3_A_4