信號檢測課程實驗報告

一、實驗目的熟悉MATLAB軟件的相關函數(shù)的使用方法、各種信號的時域波形。實驗內容1、用MATLAB軟件繪制各類典型信號的時域波形并且給出程序源代碼。(1)指數(shù)信號、指數(shù)序列；(2)單位階躍函數(shù)、單位階躍序列；(3)單位脈沖信號；單位抽樣序列；(4)矩形信號，矩形序列信號，矩形寬度自定；(5)正弦信號、正弦序列，其中振幅和頻率自選。2、假設存在一個單邊離散指數(shù)序列，參數(shù)自選；為隨機噪聲，方差和均值自選；繪制，和的圖形。3、繪制你所采集的電信號的時域波形（部分即可），進行必要的說明和分析。4、選作提高：編寫一個程序，運行后可以在鍵盤上依次任意輸入兩個序列，計算兩個序列的線性卷積，并輸出計算結果。（相關函數(shù)：input，conv）三、實現(xiàn)方案1、用MATLAB軟件繪制各類典型信號的時域波形并且給出程序源代碼。(1)指數(shù)信號、指數(shù)序列

(2)單位階躍函數(shù)、單位階躍序列單位脈沖信號、單位抽樣序列矩形信號、矩形序列信號正弦信號、正弦序列

程序源代碼closeallclearallclfq1=-5;q2=15;p1=0;p2=5;%繪制指數(shù)信號x1=q1:0.01:q2;y1=1.2.^x1;subplot(5,2,1);plot(x1,y1);grid;title('指數(shù)函數(shù)y1=1.2^x1');xlabel('x1');ylabel('y1');axis([q1,q2,-1,10]);%繪制指數(shù)序列n1=q1:q2;z1=1.2.^n1;subplot(5,2,2);stem(n1,z1);grid;title('指數(shù)序列z1=1.2^n1');xlabel('n1');ylabel('z1');axis([q1,q2,-1,10]);%繪制單位階躍信號x2=q1:0.01:q2;y2=(x2-p1)>=0;subplot(5,2,3);plot(x2,y2);grid;title('單位階躍信號');xlabel('x2');ylabel('y2');axis([q1,q2,-1,10]);%繪制單位階躍序列n2=q1:q2;z2=(n2-p1)>=0;subplot(5,2,4);stem(n2,z2);grid;title('單位階躍序列');xlabel('n2');ylabel('z2');axis([q1,q2,-1,10]);%繪制單位脈沖信號x3=linspace(q1,q2,100*(q2-q1+1));length_x3=length(x3);y3=zeros(1,length_x3);forj=1:100*(q2-q1+1)ifj==100*(p2-q1)+1y3(j)=1;endendsubplot(5,2,5);plot(x3,y3);grid;title('單位脈沖信號');xlabel('x3');ylabel('y3');axis([q1,q2,-1,10]);%繪制單位抽樣序列n3=linspace(q1,q2,q2-q1+1);length_n3=length(n3);z3=zeros(1,length_n3);forj=1:q2-q1+1ifj==p2-q1+1z3(j)=1;endendsubplot(5,2,6);stem(n3,z3);grid;title('單位抽樣序列');xlabel('n3');ylabel('z3');axis([q1,q2,-1,10]);%繪制矩形信號x4=q1:0.01:q2;y41=(x4-p1)>=0;y42=(x4-p2)<=0;y4=3*y41.*y42;subplot(5,2,7);plot(x4,y4);grid;title('矩形信號');xlabel('x4');ylabel('y4');axis([q1,q2,-1,10]);%繪制矩形序列n4=-5:15;z41=(n4-p1)>=0;z42=(n4-p2)<=0;z4=3*z41.*z42;subplot(5,2,8);stem(n4,z4);grid;title('矩形序列');xlabel('n4');ylabel('z4');axis([q1,q2,-1,10]);%繪制正弦信號omig=2;x5=-5:0.01:15;y5=3*sin(omig*x5);subplot(5,2,9);plot(x5,y5);grid;title('正弦信號');xlabel('x5');ylabel('y5');axis([q1,q2,-5,5]);%繪制正弦序列omig=2;n5=q1:0.5:q2;z5=3*sin(omig*n5);subplot(5,2,10);stem(n5,z5);grid;title('正弦序列');xlabel('n5');ylabel('z5');axis([q1,q2,-5,5]);假設存在一個單邊離散指數(shù)序列，參數(shù)自選；為隨機噪聲，方差和均值自選；繪制，和的圖形。假設存在一個右單邊的離散指數(shù)序列，n=[0,19];存在一個隨機噪聲，方差為0.9，均值為2。繪制，和的圖形，如下所示：程序源代碼：closeallclearallN=20;n=0:N-1;x1=exp(0.1*n);x2=randn(1,N);x2=x2-mean(x2);x2=x2/std(x2);a=2;b=0.9;x2=a+b*x2;x3=x1+x2;subplot(2,2,1);stem(n,x1);title('s[n]');gridsubplot(2,2,2);stem(n,x2);title('d[n]');gridsubplot(2,1,2);stem(n,x3);title('s[n]+d[n]');gridmean(x2)var(x2)繪制你所采集的電信號的時域波形（部分即可），進行必要的說明和分析。將一段wav格式的音頻信號tada.wav讀入到MATLAB中，利用plot函數(shù)畫出該電信號的波形，如下圖所示：程序源代碼：%繪制采集的電信號的時域波形closeallclearallfs=20000;[y,fs,nbits]=wavread('tada.wav');sound(y,fs);plot(y);grid;title('信號波形');xlabel('時間t');ylabel('幅值y');選作提高：編寫一個程序，運行后可以在鍵盤上依次任意輸入兩個序列，計算兩個序列的線性卷積，并輸出計算結果。（相關函數(shù)：input，conv）程序源代碼：%兩序列線性卷積closeallclearallx1=input('Pleaseenteraarray\nx1=');x2=input('Pleaseenteranotherarray\nx2=');y=conv(x1,x2);Y運行結果：Pleaseenteraarrayx1=[123]Pleaseenteranotherarrayx2=[1234]y=1410161712結果