資源描述：
《數(shù)字信號(hào)處理(matlab版)上機(jī)實(shí)驗(yàn)操作》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫(kù)。
1、實(shí)驗(yàn)一離散時(shí)間信號(hào)與系統(tǒng)一、實(shí)驗(yàn)?zāi)康模?、熟悉常見(jiàn)離散時(shí)間信號(hào)的產(chǎn)生方法；2、熟悉離散時(shí)間系統(tǒng)的單位脈沖響應(yīng)和單位階躍響應(yīng)的求解方法；3、熟悉離散時(shí)間信號(hào)經(jīng)過(guò)離散時(shí)間系統(tǒng)的響應(yīng)的求解方法。二、實(shí)驗(yàn)內(nèi)容：已知離散時(shí)間系統(tǒng)差分方程為y(n)-0.5y(n-1)+0.06y(n-2)=x(n)+x(n-1)，求1、該系統(tǒng)的單位脈沖響應(yīng)并繪圖；2、該系統(tǒng)的單位階躍響應(yīng)并繪圖；3、已知x(n)=可自己指定用filter函數(shù)經(jīng)過(guò)系統(tǒng)的響應(yīng)并繪圖；4、用conv_m函數(shù)求系統(tǒng)響應(yīng)并繪圖。三、實(shí)驗(yàn)平臺(tái)：MATLAB集成系統(tǒng)四、
2、設(shè)計(jì)流程：此處寫(xiě)個(gè)人自己的設(shè)計(jì)流程五、程序清單：此處寫(xiě)程序內(nèi)容六、調(diào)試和測(cè)試結(jié)果：此處寫(xiě)程序的執(zhí)行結(jié)果和實(shí)驗(yàn)過(guò)程中的調(diào)試經(jīng)過(guò)、出現(xiàn)的錯(cuò)誤和對(duì)應(yīng)的解決方法七、教師評(píng)語(yǔ)與成績(jī)?cè)u(píng)定此處由老師填寫(xiě)上機(jī)操作：實(shí)驗(yàn)一離散時(shí)間信號(hào)與系統(tǒng)實(shí)驗(yàn)內(nèi)容：1．脈沖響應(yīng)>>b=[1,1];a=[1,-0.5,0.06];n=[-10:25];>>impz(b,a,n);>>title('ImpulseResponse');xlabel('n');ylabel('h(n)')172.單位階躍響應(yīng)>>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)響應(yīng)');xlabel('m'
5、);ylabel('y(m)');17實(shí)驗(yàn)二離散信號(hào)與系統(tǒng)的連續(xù)頻域分析一、實(shí)驗(yàn)?zāi)康模?、掌握離散時(shí)間信號(hào)的DTFT的MATLAB實(shí)現(xiàn)；2、掌握離散時(shí)間系統(tǒng)的DTFT分析；3、掌握系統(tǒng)函數(shù)和頻率相應(yīng)之間的關(guān)系。二、實(shí)驗(yàn)內(nèi)容：1、自定義一個(gè)長(zhǎng)度為8點(diǎn)的信號(hào)，信號(hào)幅度值也由自己任意指定，對(duì)該信號(hào)作DTFT，分別畫(huà)出幅度譜和相位譜2、已知離散時(shí)間系統(tǒng)差分方程為y(n)-0.5y(n-1)+0.06y(n-2)=x(n)+x(n-1)，求出并畫(huà)出其頻率響應(yīng)3、求該系統(tǒng)系統(tǒng)函數(shù)，并畫(huà)極零點(diǎn)圖，并通過(guò)freqz函數(shù)求頻率響
6、應(yīng)。三、實(shí)驗(yàn)平臺(tái)：MATLAB集成系統(tǒng)四、設(shè)計(jì)流程：此處寫(xiě)個(gè)人自己的設(shè)計(jì)流程五、程序清單：此處寫(xiě)程序內(nèi)容六、調(diào)試和測(cè)試結(jié)果：此處寫(xiě)程序的執(zhí)行結(jié)果和實(shí)驗(yàn)過(guò)程中的調(diào)試經(jīng)過(guò)、出現(xiàn)的錯(cuò)誤和對(duì)應(yīng)的解決方法七、教師評(píng)語(yǔ)與成績(jī)?cè)u(píng)定此處由老師填寫(xiě)上機(jī)操作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
