《數字信號處理》上機全部源代碼調試通過,完整版

1、《數字信號處理數字信號處理》上機全部源代碼調試通過，完整版上機全部源代碼調試通過，完整版（高西全，第四版（高西全，第四版）實驗一實驗一%實驗1：系統(tǒng)響應及系統(tǒng)穩(wěn)定性closeallclearall%調用fliter解差分方程，由系統(tǒng)對un的響應判斷穩(wěn)定性%內容1：調用filter解差分方程，由系統(tǒng)對u(n)的響應判斷穩(wěn)定性A=[10.9]B=[0.050.05]x1n=[11111111zeros(150)]x2n=ones(1128)

2、hn=impz(BA58)subplot(221)y=h(n)tstem(hny)title((a)系統(tǒng)單位脈沖響應h(n))y1n=filter(BAx1n)subplot(222)y=y1(n)tstem(y1ny)title((b)系統(tǒng)對R8(n)的響應y1(n))y2n=filter(BAx2n)subplot(224)y=y2(n)tstem(y2ny)title((c)系統(tǒng)對u(n)的響應y2(n))y1n=filter(B

3、Ax1n)subplot(222)y=y1(n)tstem(y1ny)title((b)系統(tǒng)對R8(n)的響應y1(n))y2n=filter(BAx2n)subplot(224)y=y2(n)tstem(y2ny)title((c)系統(tǒng)對u(n)的響應y2(n))%內容2：調用conv函數計算卷積x1n=[11111111]%產生信號x1n=R8nh1n=[ones(110)zeros(110)]h2n=[12.52.51zeros(

4、110)]y21n=conv(h1nx1n)y22n=conv(h2nx1n)figure(2)subplot(221)y=h1(n)tstem(h1ny)%調用函數tstem繪圖title((d)系統(tǒng)單位脈沖響應h1(n))subplot(222)y=y21(n)tstem(y21ny)title((e)h1(n)與R8(n)的卷積y21(n))subplot(223)y=h2(n)tstem(h2ny)%調用函數tstem繪圖tit

5、le((f)系統(tǒng)單位脈沖響應h2(n))實驗二實驗二%時域采樣理論驗證程序exp2a.mTp=641000%觀察時間Tp=64微秒%產生M長采樣序列x(n)%Fs=1000T=1FsFs=1000T=1FsM=TpFsn=0:M1A=444.128alph=pi502^0.5omega=pi502^0.5xnt=Aexp(alphnT).sin(omeganT)Xk=Tfft(xntM)%M點FFT[xnt)]yn=xa(nT)subp

6、lot(321)tstem(xntyn)%調用自編繪圖函數tstem繪制序列圖boxontitle((a)Fs=1000Hz)k=0:M1fk=kTpsubplot(322)plot(fkabs(Xk))title((a)TFT[xa(nT)]Fs=1000Hz)xlabel(f(Hz))ylabel(幅度)axis([0Fs01.2max(abs(Xk))])%=================================%Fs=

7、300Hz和Fs=200Hz的程序與上面Fs=1000Hz的程序完全相同。%%%%%%%%%%%%fs=300HzTp=641000%觀察時間Tp=64微秒%產生M長采樣序列x(n)%Fs=1000T=1FsFs=300T=1FsM=TpFsn=0:M1A=444.128alph=pi502^0.5omega=pi502^0.5xnt=Aexp(alphnT).sin(omeganT)Xk=Tfft(xntM)%M點FFT[xnt)]y

8、n=xa(nT)figure(2)subplot(321)tstem(xntyn)%調用自編繪圖函數tstem繪制序列圖boxontitle((b)Fs=300Hz)k=0:M1fk=kTpsubplot(322)plot(fkabs(Xk))title((b)TFT[xa(nT)]Fs=1000Hz)xlabel(f(Hz))ylabel(幅度)axis([0Fs01.2max(abs(Xk))])%%%%%%%%%%%%%%%%%%