資源描述:
《數(shù)字信號處理_共軛對稱、共軛反對稱.doc》由會員上傳分享,免費在線閱讀,更多相關內(nèi)容在行業(yè)資料-天天文庫。
1、......xxxx大學實驗報告學生姓名_xxx_學號_xxxxxxx_年級班級_xxxxxxx_實驗項目_xxxxxxxx_實驗時間_xxxxxxxxx_實驗二一、實驗目的:1.充分熟悉復指數(shù)函數(shù)find、sigshift、sigfold函數(shù)的使用;2.熟悉序列的加、減、乘、除、移位、折疊的計算;3.能夠畫出結(jié)果的圖形。二、實驗步驟:1.用help查找find、sigshift、sigfold函數(shù)的使用情況;2.編輯并生成函數(shù)sigadd.m(序列相加)function[y,n]=sigadd(
2、x1,n1,x2,n2)%實現(xiàn)y(n)=x1(n)+x2(n)%[y,n]=sigadd(x1,n1,x2,n2)%y=在包含n1和n2的n點上求序列和%x1=在n1上的第一序列%x2=在n2上的第二序列(n2可與n1不等)n=min(min(n1),min(n2)):max(max(n1),max(n2));%y(n)的長度y1=zeros(1,length(n));y2=y1;%初始化y1(find((n>=min(n1))&(n<=max(n1))==1))=x1;%具有y的長度的x1y2(
3、find((n>=min(n2))&(n<=max(n2))==1))=x2;%具有y的長度的x2.專業(yè)word可編輯.......y=y1+y2;%序列相加.1.編輯并生成函數(shù)sigmult.m(序列相乘)function[y,n]=sigmult(x1,n1,x2,n2)%實現(xiàn)y(n)=x1(n)*x2(n)%[y,n]=sigmult(x1,n1,x2,n2)%y=在n區(qū)間上的乘積序列,n包含n1和n2%x1=在n1上的第一序列%x2=在n2上的第二序列(n2可與n1不等)n=min(min
4、(n1),min(n2)):max(max(n1),max(n2));%y(n)的長度y1=zeros(1,length(n));y2=y1;%初始化y1(find((n>=min(n1))&(n<=max(n1))==1))=x1;%具有y的長度的x1y2(find((n>=min(n2))&(n<=max(n2))==1))=x2;%具有y的長度的x2y=y1.*y2;%序列相乘2.編輯并生成函數(shù)sigshift.m(序列移位)function[y,n]=sigshift(x,m,n0)%實現(xiàn)
5、y(n)=x(n-n0)%[y,n]=sigshift(x,m,n0)n=m+n0;y=x;3.編輯并生成函數(shù)sigfold.m(序列折疊).專業(yè)word可編輯.......function[y,n]=sigfold(x,n)%實現(xiàn)y(n)=x(-n)%[y,n]=sigfold(x,n)y=fliplr(x);n=-fliplr(n);1.編輯并生成實現(xiàn)兩序列相乘和相加程序clc;clear;x1=[0,1,2,3,4,3,2,1,0];n1=-2:6;x2=[2,2,0,0,0,-2,-2];
6、n2=2:8;[y1,n]=sigmult(x1,n1,x2,n2);[y2,n]=sigadd(x1,n1,x2,n2);subplot(2,2,1);stem(n1,x1);title('序列x1')xlabel('n');ylabel('x1(n)');subplot(2,2,2);stem(n2,x2);title('序列x2')xlabel('n');ylabel('x2(n)');subplot(2,2,3);stem(n,y1);title('兩序列相乘')xlabel('n');y
7、label('y1(n)');subplot(2,2,4);stem(n,y2);title('兩序列相加')xlabel('n');ylabel('y2(n)');運行以上程序得到的圖形:.專業(yè)word可編輯.......體會:相加或相乘時,兩序列尺度要保持一致。1.改變坐標尺度重新畫圖,使圖形更加直觀;.專業(yè)word可編輯.......1.在命令窗口打出x1,x2,y1,y2的值,并貼近實驗報告;2.編輯并生成實現(xiàn)序列移位和折疊程序x1=[0,1,2,3,4,3,2,1,0];n1=-2:6;
8、[y1,n2]=sigshift(x1,n1,2);[y2,n3]=sigfold(x1,n1);subplot(3,1,1);stem(n1,x1);title('序列x1')xlabel('n');ylabel('x1(n)');subplot(3,1,2);stem(n2,y1);title('序列移位')xlabel('n');ylabel('y1(n)');subplot(3,1,3);stem(n3,y2);title('序列折疊').專業(yè)word可編輯.......x