数字信号处理实验报告96186new

《数字信号处理实验报告96186new》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、数字信号处理——王家南重庆交通大学学生实验报告实验课程名称数字信号处理开课实验室数学实验室学院理学院年级11专业班信息与计算科学2学生姓名王家南学号631122020203开课时间2013至2014学年第1学期9数字信号处理——王家南程序代码：Convexam.m：nx=1:120x=(impseq(1,0,119)+impseq(2,40,120)+impseq(2,70,120)+impseq(1,80,120));nh=0:24;h=0.95.^nh;[y,ny]=my_conv(x,h);subplot(2,1,1),stem(nx,x);title('输入信号’)subplo

2、t(2,1,2),stem(ny,y);title('卷积滤波')impseq.m：function[x]=impseq(n0,n1,n2)if(n1>n2)errorendx=n0.*[zeros(1,(n1-1)),1,zeros(1,(n2-n1))];my_conv.m:function[y,ny]=my_conv(x,h)ny=[0:(length(x)+length(h)-2)];y=conv(x,h);9数字信号处理——王家南实验结果：心得体会：通过本次实验，了解到了conv函数的用法；同时了解到函数名不得和自己取的m文件不能重名否则会出现错误，同时学会了stem函数。9

3、数字信号处理——王家南程序代码：Casfiltr.m:functiony=casfiltr(b0,B,A,x)[K,L]=size(B);N=length(x);w=zeros(K+1,N);w(1,:)=x;fori=1:1:Kw(i+1,:)=filter(B(i,:),A(i,:),w(i,:));endy=b0*w(K+1,:);cas2dir.m:function[b,a]=cas2dir(b0,B,A)[K,L]=size(B);b=[1];a=[1];fori=1:1:Kb=conv(b,B(i,:));a=conv(a,A(i,:));endb=b*b0;casfilt

4、rexam.m:b0=2;B=[1,-0.5,0;1,2,1];A=[1,-0.75,0;1,0,0.81];N=64;n=1:N;x=[(n-1)==0];9数字信号处理——王家南h_cas=casfiltr(b0,B,A,x);subplot(2,1,1)stem(n,h_cas,'filled','k');axis([0,70,-4,6]);title('¼¶ÁªÐ͵¥Î»Âö³åÏìÓ¦')[b,a]=cas2dir(b0,B,A);hdir=filter(b,a,x);subplot(2,1,2)stem(n,hdir,'filled','k');axis([0,70,-4

5、,6]);title('Ö±½ÓÐ͵¥Î»Âö³åÏìÓ¦')实验结果：心得体会：通过本次实验了解到了级联型转直接性的转变方法，在函数层面上了解了级联转直连的过程，对这两种脉冲响应加深了印象。9数字信号处理——王家南程序代码：Ideal_lp.m:functionhd=ideal_lp(wc,M)alpha=(M-1)/2;n=[0:1:(M-1)];m=n-alpha+eps;%加一个小数以避免零做除数没有意义hd=sin(wc*m)./(pi*m);[db,mag,pha,grd,w]=freqz_m(hd,[1])Freqz_m.m:function[db,mag,pha,gr

6、d,w]=freqz_m(b,a)[H,w]=freqz(b,a,1000,'whole');H=(H(1:1:501))';w=(w(1:1:501))';mag=abs(H);db=20*log10((mag+eps)/max(mag));pha=angle(H);grd=grpdelay(b,a,w);chuangexam.mM=121;n=[0:1:M-1];wc=0.45*pi;hd=ideal_lp(wc,M);w_rec=(boxcar(M))';w_ham=(hamming(M))';h_rec=hd.*w_rec;h_ham=hd.*w_ham;[db1,mag1,p

7、ha1,grd1,w1]=freqz_m(h_rec,[1]);[db2,mag2,pha2,grd2,w2]=freqz_m(h_ham,[1]);figure(1)subplot(2,2,1);stem(n,hd,'k');gridtitle('加窗前脉冲响应')axis([0,M-1,-0.2,0.5]);ylabel('hd(n)')subplot(2,2,2);stem(n,hd,'k');gridtitle('加窗后脉冲响应')9数