DSP课程大作业英文版图文

LabReportofDSPLab1Discrete-TimeFourierTransform(1Generateandplotthediscrete-timesignal(2DetermineandplottheDTFTofthesignal,observing(iiiDetermineandplottheoutputofthesystem,iftheinputis.thewillbe，wecangetthreepolepoints:，CharacteristicequationC“”ereweusethefunctionfilter,wecangettheimpulseresponse:H“”Chart1-3UnitStepResponseWeget22.1234,sothesystemisstable.ans=(3iftheinputis,theinputsequenceandoutputsequenceareshowed:n=[0:30]x=(0.93.^n.*cos(0.14*pi*n+pi/3;stem(n,x;title('x(n';xlabel('n';ylabel('x(n'1-2.n=[0:30];x=(0.89.^n.*(stepseq(0,30,0w=[0:500]*2*pi/500;X=exp(j*w./(exp(j*w-0.89*ones(1,501;mag=abs(X;ang=angle(Xfigure(1;plot(w/pi,mag;gridylabel('|X(w|';xlabel('w(inpiunit';title('DTFTofx(n-magnitude'figure(2;plot(w/pi,ang;gridylabel('argX(w/pi';xlabel('w(inpiunit';title('DTFTofx(n-angle'1-3.(2num=[9,-4,8,-1];den=[1,-0.65,0.35,0]x=impseq(0,0,100;n=[0:100]h=filter(num,den,xstem(n,h;xlabel('n';ylabel('y(n';title('ImpluseResponse';sum(abs(h(3num=[9,-4,8,-1];den=[1,-0.65,0.35,0]n=[0:200];x=(2+5*cos(0.27*pi*n+3*sin(0.77*pi*n.*stepseq(0,200,0figure(1;stem(n,x;xlabel('n';ylabel('x(n';title('InputSequence';h=filter(num,den,x;figure(2;stem(n,h;xlabel('n';ylabel('y(n';title('OutputResponse'Functions:(1function[x,n]=impseq(n0,n1,n2n=[n1:n2];x=[(n-n0==0];(2function[x,n]=stepseq(n1,n2,n0ifnargin~=3n=n1:n2;:Withoutlossofgenerality,supposethattheamplitudesofthetwofrequencycomponentsaredifferent,theyare2and1.5.Sothesignalwithoutnoiseis:.wecanusethefunction“awgn”toaddaGaussianwhiteNoise,andhereweassumetheradiois-2dB.Sowecangetthecharts:.Supposethatthedatarecordlengthis.figure(1;plot(t,x;title('原始信号';xlabel('t';y=awgn(x,-2;figure(2;plot(t,y;title('加噪信号';xlabel('t';X1=fft(x1,200;X2=fft(x2,200;figure(1;figure(2;Lab3DigitalFilterDesignSolution:‘’WiththefunctionFilter,wecangettheoutputresponsesofunitimpulseandstepsignal.Theoutputresponsesare：xMyclassnumberisF0903018,mystudentIDnumberis5090309481.SoX=（9+3+1+8+5+9+3+9+4+8+1）mod(20=0.Bn’DesignanFIRbandstopfilterwiththefollowingspecifications:and.ofSdesignrequirementsWindow,actualunitimpulseresponse,ofBandelimination:amplitudefrequencycharacteristicscurveAddendum：MatlabProgramList481,0x2=stepseq(0,50,0;y2=filter(b,a,x2;0wp1=0.4*pi;ws1=0.45*pi;wp2=0.7*pi;ws2=0.65*pi;tr_width=min((ws1-wp1,(wp2-ws2;M=ceil(6.2*pi/tr_width+1n=0:1:M-1;wc1=(ws1+wp1/2;wc2=(ws2+wp2/2;hd=ideal_bs(wc1,wc2,M;w_hanning=(hanning(M';h=hd.*w_hanning;[db,mag,pha,w]=freqz_m(h,[1];delta_w=2*pi/1000;Ap=-(min(db(1:1:wp1/delta_w+1As=-round(max(db(ws1/delta_w+1:1:ws2/delta_w+1subplot(2,2,1;stem(n,hd,'.';title('IdealImpulseResponse';subplot(2,2,2;stem(n,w_hanning,'.';title('HanningWindow';subplot(2,2,3;stem(n,h,'.';title('ActualImpulseresponse';subplot(2,2,4;plot(w/pi,db;title('ManitudeResponseindB';axis([01-10010];Functions:function[db,mag,pha,w]=freqz_m(b,a;[H,w]=freqz(b,a,1000,'whole';H=(H(1:1:501';w=(w(1:1:501';mag=abs