資源描述：
《數(shù)字信號處理(matlab版)上機實驗操作》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫。
1、實驗一離散時間信號與系統(tǒng)一、實驗目的：1、熟悉常見離散時間信號的產(chǎn)生方法；2、熟悉離散時間系統(tǒng)的單位脈沖響應和單位階躍響應的求解方法；3、熟悉離散時間信號經(jīng)過離散時間系統(tǒng)的響應的求解方法。二、實驗內(nèi)容：已知離散時間系統(tǒng)差分方程為y(n)-0.5y(n-1)+0.06y(n-2)=x(n)+x(n-1)，求1、該系統(tǒng)的單位脈沖響應并繪圖；2、該系統(tǒng)的單位階躍響應并繪圖；3、已知x(n)=可自己指定用filter函數(shù)經(jīng)過系統(tǒng)的響應并繪圖；4、用conv_m函數(shù)求系統(tǒng)響應并繪圖。三、實驗平臺：MATLAB集成系統(tǒng)四、
2、設(shè)計流程：此處寫個人自己的設(shè)計流程五、程序清單：此處寫程序內(nèi)容六、調(diào)試和測試結(jié)果：此處寫程序的執(zhí)行結(jié)果和實驗過程中的調(diào)試經(jīng)過、出現(xiàn)的錯誤和對應的解決方法七、教師評語與成績評定此處由老師填寫上機操作：實驗一離散時間信號與系統(tǒng)實驗內(nèi)容：1．脈沖響應>>b=[1,1];a=[1,-0.5,0.06];n=[-10:25];>>impz(b,a,n);>>title('ImpulseResponse');xlabel('n');ylabel('h(n)')172.單位階躍響應>>x=stepseq(0,-10,25);
3、s=filter(b,a,x);Warning:Functioncallstepseqinvokesinexactmatchd:MATLAB7workSTEPSEQ.M.>>stem(n,s)>>title('StepResponse');xlabel('n');ylabel('s(n)')3.>>a=[1,-0.5,0.06];b=[1,1];>>n=-20:120;>>x1=exp(-0.05*n).*sin(0.1*pi*n+pi/3);>>s1=filter(b,a,x1);>>stem(n,s1
4、);;xlabel('n');ylabel('s1(n)');174．>>a=[1,-0.5,0.06];b=[1,1];>>n=-20:120;>>h=impz(b,a,n);>>x1=exp(-0.05*n).*sin(0.1*pi*n+pi/3);>>[y,m]=conv_m(x1,n,h,n);Warning:Functioncallconv_minvokesinexactmatchd:MATLAB7workCONV_M.M.>>stem(m,y);title('系統(tǒng)響應');xlabel('m'
5、);ylabel('y(m)');17實驗二離散信號與系統(tǒng)的連續(xù)頻域分析一、實驗目的：1、掌握離散時間信號的DTFT的MATLAB實現(xiàn)；2、掌握離散時間系統(tǒng)的DTFT分析；3、掌握系統(tǒng)函數(shù)和頻率相應之間的關(guān)系。二、實驗內(nèi)容：1、自定義一個長度為8點的信號，信號幅度值也由自己任意指定，對該信號作DTFT，分別畫出幅度譜和相位譜2、已知離散時間系統(tǒng)差分方程為y(n)-0.5y(n-1)+0.06y(n-2)=x(n)+x(n-1)，求出并畫出其頻率響應3、求該系統(tǒng)系統(tǒng)函數(shù)，并畫極零點圖，并通過freqz函數(shù)求頻率響
6、應。三、實驗平臺：MATLAB集成系統(tǒng)四、設(shè)計流程：此處寫個人自己的設(shè)計流程五、程序清單：此處寫程序內(nèi)容六、調(diào)試和測試結(jié)果：此處寫程序的執(zhí)行結(jié)果和實驗過程中的調(diào)試經(jīng)過、出現(xiàn)的錯誤和對應的解決方法七、教師評語與成績評定此處由老師填寫上機操作1>>n=0:7;>>x=(0.9*exp(j*pi/3)).^n;>>w=-pi:pi/200:pi;>>X=dtft(x,n,w);>>magX=abs(X);angX=angle(X);>>subplot(2,1,1);plot(w/pi,magX);xlabel('w/
7、pi');ylabel('幅度
8、X
9、');>>subplot(2,1,2);plot(w/pi,angX);xlabel('w/pi');ylabel('相位(rad/π)');172>>a=[1,-0.5,0.06];b=[1,1,0];>>m=0:length(b)-1;l=0:length(a)-1;>>w=0:pi/500:pi;>>num=b*exp(-j*m'*w);>>den=a*exp(-j*l'*w);>>H=num./den;>>magH=abs(H);angH=angle(H);>>H1=
10、freqz(b,a,w);>>magH1=abs(H1);angH1=angle(H1);>>subplot(2,2,2);plot(w/pi,angH/pi);grid;xlabel('w(frequencyinpiunits)');ylabel('相位(rad/π)');>>subplot(2,2,1);plot(w/pi,magH);grid;xlabel('w(frequenc
