淘豆網
1/11
下載文檔
文檔分類:高等教育 > 實驗設計

數字信號處理實驗報告.doc


下載后只包含 1 個 DOC 格式的文檔,沒有任何的圖紙或源代碼,查看文件列表 我要舉報
0/100
您的瀏覽器不支持進度條
更多>>該用戶其他文檔
下載所得到的文件列表
數字信號處理實驗報告.doc
文檔介紹:
數字信號處理實驗報告.docEvaluationWarning:ThedocumentwascreatedwithSpire..重慶交通大學學生實驗報告實驗課程名稱 開課實驗室 數學實驗室 學 院理學院年級口專業班信息與計算科學2學生姓名王家南學號632222020203開課時間2013至2014學年第1學期優良中差排版的規范性報吿的完抽方法的正確性內容的創新性總成績指導教師韓君§1實驗一用卷積運算求濾波器的響應實驗內容:給定一個FIR濾波器,它的單位沖激響應為/j(n)=0.95n, 0SW24.現有一輸入T(n)=6(n)+24(n一40)+26(n-70)+6(n-80),宛=0,1.2,???,120求其輸出t/(n),并在同一張圖匕畫出x(n)和y(n).程序代碼:Convexam.m:nx=l:120x=(impseq(1,0,119)+impseq(2,40,120)+impseq(2,70,120)+impseq(1,80,120));nh=0:24;h=0.95?Anh;[y/ny]=my_conv(x,h);subplot(2,1,1),stem(nx,x);titie(,輸入信號")subplot(2,1,2)rstem(nyzy);title{'卷積濾波')impseq.m:function[x]=impseq(nO,nl,n2)if(nl>n2)errorendx=n0?*[zeros(1,(nl-1)),1zzeros(lz(n2-nl))];my_conv.m:function[y,ny]=my_conv(x,h)ny=[0:(length(x)+length(h)-2)];y=conv(x,h);實驗結果:輸入信號21.51©0.520°;60卷積濾泯4080100120心得體會:通過本次實驗,了解到了conv函數的用法;同時了解到函數名不得和白己取的m文件不能重名否則會出現錯誤,同時學會了stem函數。§2實驗二級聯型到直接型的形式轉換實驗內容:給是一個級聯形式表示的濾波器、其傳遞函數為2-zi 1+2Z-1+Z-21-0.75z一11+0.81Z-2將其轉換為直接型,并比較級聯型和直接型的單位脈沖響應.程序代碼Casfiltr.m:functiony=casfiltr(bO,B,A,x)[KzL]=size(B);N=length(x);w=zeros(K+lfN);w(lz:)=x;fori=l:1:Kw(i+l,:)=filter(B(iz:),A(i,:),w(i,:));endy=bO*w(K+lr:);cas2direin:function[b,a]=cas2dir(bO,B,A)[K,L]=size(B);b=[l];a=[l];fori=l:1:Kb=conv(b,B(i,:));a=conv(azA(iz:));endb=b*bO;casfiltrexam.m:b0=2;B=[l/-0.5/0;l/2/l];A=[l/-0.75/0;l/0/0.81];N=64;n=l:N;x=[(n-1)==0];h_cas=casfiltr(bO,B,A,x);subplot(2,1,1)stem(nzh_cas,1filled1,1k');axis([0,70,-4,61);titie( &彳^iio|1)[bza]=cas2dir(bO,B,A);hdir=filter(b#a,x);subplot(2,1,2)stemCrizhdir,'filled','k1);axis([0,70,-4,61);title(,6±^ODip¥i»Ad3aii6|1)實驗結果:級聯型單位脈沖響應6 1 1 1 1 4-f0 10 20 30 40 50 60 70直接型單位脈沖響應6 1 1 1 1 4-f0 10 20 30 40 50 60 70心得體會:通過本次實驗了解到了級聯型轉直接性的轉變方法,在函數層血上了解了級聯轉肓連的過程,對這兩種脈沖響應加深了印象。§3實驗三用窗函數法設計FIR濾波器實驗內容:#1漢明窗i殳計一個FIR低通濾波器,該濾波器截II:頻率他=0?45碼窗的丨、隻為M取121.諸在圖屮価出:加漢明窗前后低通濾波器的沖激響應以及濾波器的幅頻特性.程序代碼:Ideal」p-m:functionhd=ideal_lp(wc,M)alpha=(M-l)/2;n=[0:1:(M-l)];m=n-alpha+eps;%加一個小數以避免零做除數沒有意義hd=sin(wc*m)./(pi*m);[db,mag,pha,grd,w]=freqz_m(hd,[1])Freqz_m.m:function[dbzmag,pha,grd,w]=freqz_m(b,a)[H,w]=freqz(b 內容來自淘豆網www.eawpme.live轉載請標明出處.
北京单场开奖