數(shù)字信號處理實驗報告一

姓名：彭浩 班別：電氣2班 學號：3113001318實驗一時域中的離散信號一、實驗目的1．熟悉matlab程序。2．掌握時域中的離散信號基本運算。二、實驗內(nèi)容和要求了解和掌握時域中的離散信號加、乘、移位、卷積等基本運算；三、實驗主要儀器設(shè)備和材料1．計算機2．Matlab軟件四、實驗方法、步驟及結(jié)果測試（1）連續(xù)信號和離散信號問題：請計算連續(xù)信號離散化后前五個點的絕對誤差和相對誤差。程序如下：dt=0.001;t=0:dt:6; %建立連續(xù)自變量向量xa=sqrt(t)+cos(t); %原始的連續(xù)時間信號xa(t)T=0.5;n=0:6/T; %T為采樣周期，建立離散自變量向量x=sqrt(n*T)+cos(n*T); %采樣周期為T的離散時間信號x(n)deltax=0.5; %deltax為x的量化單位xq=round(x/deltax)*deltax; %舍入量化后的數(shù)字信號xqn=1;whilen<=5 %循環(huán)五次分別前五個點的絕對誤差和相對誤差輸出ea=xq(n)-x(n)er=ea/x(n)n=n+1;end輸出：(為了美觀需要對輸出格式做了調(diào)整)ea=0 er=0ea=-0.0847 er=-0.0534ea=-0.0403 er=-0.0262ea=0.2045 er=0.1579ea=0.0019 er=0.0019（2)離散信號基本序列及運算問題：編程實現(xiàn)：1）程序如下：>> a=impseq(-1,-3,3); b=impseq(3,-3,3); x1=1.5*a-b輸出：x1=001.5000000-1.00002)令x=[1-246-5810],產(chǎn)生并畫出下列序列的樣本：X2(n)=3x(n-4)+x(n-4)+2x(n)程序如下：>> x=[1-246-5810]; nx=1:7; [a,na]=seqshift(x,nx,4); b=4*a; nb=na; c=2*x; nc=nx; [x2,n2]=seqadd(b,nb,c,nc) stem(x2)輸出：x2=2-4812-683624-203240n2=12345678910113)計算x1+x2,x1*x2,,(x1、x2同1),2))，繪制輸出中間結(jié)果及最終結(jié)果。程序如下：>> x1=[001.5000-1]; nx1=-3:3; x2=[2-4812-683624-203240]; nx2=1:11; [X1,N1]=seqadd(x1,nx1,x2,nx2)%X1為兩個序列的和 [X2,N2]=seqmult(x1,nx1,x2,nx2)%X2為兩個序列的積 [X3,N3]=convwthn(x1,nx1,x2,nx2)%X3為兩個序列的卷積輸出：（(為了美觀需要對輸出格式做了調(diào)整)X1=Columns1through13001.500002.0000-4.00007.000012.0000-6.00008.000036.000024.0000-20.0000Columns14through1532.000040.0000X2=000000-800000000X3=003-61218-11164624-244024-2420-32-40五、實驗感想通過這次實驗，我熟悉了matlab的相關(guān)操作，以及掌握了在時域中的離散信號加、乘、移位、卷積等基本運算。在實驗中我也遇到了不少麻煩：首先之前對matlab接觸較少，對其還是比較生疏，這直接影響了我在在課堂上完成實驗的速度和準確率，其次，對于序列的運算，由于開始沒注意且沒調(diào)用函數(shù)，每次序列運算我都忽略了序列的位置向量問題，在運算前沒有使兩個序列的位置向量對齊就進行運算，還有就是序列乘法運算的時沒有使用點乘，這都導致了我每次都得不出正確結(jié)果，后來調(diào)用了老師文件中的函數(shù)，使我解決