資源描述:
《dsp實驗matlab程序?qū)嵗酚蓵T上傳分享,免費在線閱讀,更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。
1、2.離散時間信號(序列)的產(chǎn)生利用MATLAB產(chǎn)生和繪制下列有限長序列:單位脈沖序列,單位階躍序列,矩形序列x=-10:10;y=[];fori=1:21ifx(i)==0y(i)=1;elsey(i)=0;endendsubplot(3,1,1);%在子圖畫單位沖激序列stem(x,y,'*');title('沖激');xlabel('n');holdon;fori=1:21ifx(i)<0y(i)=0;elsey(i)=1;endendsubplot(3,1,2);%在子圖畫單位階躍序列stem(x,y,'*');title('階躍');xlabel('n');hold
2、on;fori=1:21ifx(i)<-4y(i)=0;elseifx(i)>4y(i)=0;elsey(i)=1;endendendsubplot(3,1,3);%在子圖畫矩形脈沖序列stem(x,y,'*');title('矩形');xlabel('n');holdon在一幅圖上繪出曲線,,A=2,α=0.5,f=2Hz,A1=1,A2=0.5,A3=0.2,f=2Hz。t=[0:0.5:720]*pi/180;x1=5*sin(pi*t/5+pi/3);x2=5*cos(t-pi/3);x3=2*sin(3*t).*cos(2*t-pi/3);subplot(3,1,
3、1);%在子圖中畫出一系列正弦余弦曲線plot(t,x1,'-red',t,x2,'-gr',t,x3,'-bl');title('正弦余弦曲線');xlabel('t');holdon;A=2;a=0.5;f=2;x4=A*exp(-a*t).*sin(2*pi*f*t);subplot(3,1,2);%在子圖中畫出正弦衰減信號plot(t,x4)title('正弦衰減信號');xlabel('t');holdon;A1=1;A2=0.5;A3=0.2;f=2;x5=A1*sin(2*pi*f*t)+A2*sin(2*pi*2*f*t)+A3*sin(2*pi*3*f*t
4、);subplot(3,1,3);%在子圖中畫出諧波信號plot(t,x5)title('諧波信號');xlabel('t');holdon;3.序列的運算生成下列序列:(1)利用MATLAB編程完成上述兩序列的卷積,并繪制運算后序列的波形。x=[12345];h=[1212];n1=[01234];n2=[2345];y=conv(x,h);%計算序列x與h的卷積和序列yn0=n1(1)+n2(1);%計算卷積和序列y的起點位置n3=length(x)+length(h)-2+n0;%計算卷積和序列y的終止位置n=n0:1:n3;%確定卷積和y非零樣值的時間向量subpl
5、ot(3,1,1);stem(n1,x,'.');%在子圖繪x(n)title('x(n)');xlabel('n');ylabel('x(n)');subplot(3,1,2);stem(n2,h,'.');%在子圖繪h(n)title('h(n)');xlabel('n');ylabel('h(n)');subplot(3,1,3);stem(n,y,'.');%在子圖繪卷積和y(n)title('y(n)=x(n)*h(n)');xlabel('n');ylabel('y(n)');(2)編寫一個完成卷積的函數(shù),輸入是兩個序列,輸出是結(jié)果。function[y,n]=
6、jsjuanji(x,n1,h,n2)%計算序列卷積和y(n)=x(n)*h(n)%y:卷積和y(n)對應(yīng)的非零樣值向量%n:卷積和y(n)對應(yīng)的時間向量%x:x(n)對應(yīng)的非零樣值向量%n1:x(n)對應(yīng)的時間向量%h:h(n)對應(yīng)的非零樣值向量%n2:h(n)對應(yīng)的時間向量n3=n1(1)+n2(1);%卷積和序列y的起始位置n4=length(x)+length(h)-1;%卷積和序列y的長度m=0;fori=1:n4%求卷積和yifi>length(x)x(i)=0;endforj=1:iifj>length(h)h(j)=0;endm=m+x(i+1-j)*h(j
7、);endy(i)=m;m=0;endn1=n1(1):1:n1(1)+n4-1;n1:%x(n)對應(yīng)的新的時間向量n2=n2(1):1:n2(1)+n4-1;%h(n)對應(yīng)的新的時間向量n=n3:1:n3+n4-1;%y(n)對應(yīng)的時間向量subplot(3,1,1);stem(n1,x,'.');%在子圖繪x(n)title('x(n)');xlabel('n');ylabel('x(n)');subplot(3,1,2);stem(n2,h,'.');%在子圖繪h(n)title('h(n)');xlabel