實驗1 離散時間信號的MATLAB實現(xiàn)-_第1頁
實驗1 離散時間信號的MATLAB實現(xiàn)-_第2頁
實驗1 離散時間信號的MATLAB實現(xiàn)-_第3頁
實驗1 離散時間信號的MATLAB實現(xiàn)-_第4頁
實驗1 離散時間信號的MATLAB實現(xiàn)-_第5頁
全文預覽已結(jié)束

下載本文檔

版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進行舉報或認領(lǐng)

文檔簡介

1、 圖 1.10 線性時不變系統(tǒng)的輸出 第四部分 序列的傅里葉變換與信號重構(gòu) MATLAB 實現(xiàn) 例 1-11 xa(t = e 1000|t| ,求其傅立葉變換 X a( j 。 解:嚴格說,在 MATLAB 中不使用 symbolic 工具箱是不能分析模擬信號的,但是當采樣時 間間隔充分小的時候,可產(chǎn)生平滑的圖形,當時間足夠長,可顯示出所有的模型,也就是可 以近似的分析。此例中, xa(t 為 fh=2000Hz 的帶限信號,因此取 1 T = 5 × 10 5 << = 25 × 10 5 2 × 2000 程序清單如下。 Dt=0.00005;t

2、=-0.005:Dt:0.005; xa=exp(-1000*abs(t; Wmax=2*pi*2000;K=500;k=0:1:K;W=k*Wmax/K; %連續(xù)時間傅立葉變換 Xa=xa*exp(-j*t'*W*Dt;Xa=real(Xa; W=-fliplr(W,W(2:501; Xa=fliplr(Xa,Xa(2:501; figure(1 subplot(2,1,1;plot(t*1000,xa,'.' xlabel('t in msec'ylabel('xa(t'gtext('模擬信號' subplot(2,1

3、,2;plot(W/(2*pi*1000,Xa*1000,'.' xlabel('Frequence in KHz'ylabel('Xa(jw*1000'gtext('連續(xù)時間傅立葉變換' 程序運行結(jié)果如圖 1.11 所示。 %模擬信號 圖 1.11 模擬信號及傅立葉變換曲線 例 1-12 以 例 1-11 中的 xa(t 說明采 樣頻率對頻域 特 性的 影 響 ,分 別 取采樣頻率 為 fs=5000Hz 和 fs=1000Hz,繪出 X a( j 曲線。 解:程序清單如下。 Dt=0.00005;t=-0.005:Dt:0.0

4、05; xa=exp(-1000*abs(t; Ts=0.0002;n=-25:1:25; x=exp(-1000*abs(n*Ts; K=500;k=0:1:K;w=pi*k/K; X=x*exp(-j*n'*w;X=real(X; w=-fliplr(w,w(2:501; X=fliplr(X,X(2:501; figure subplot(2,1,1;plot(t*1000,xa,'.' ylabel('xa1(t' gtext('離散信號'hold on stem(n*Ts*1000,x;hold off subplot(2,1,

5、2;plot(w/pi,X,'.'ylabel('X1(jw' gtext('離散時間傅立葉變換' gtext('Ts=0.2msec' 程序運行結(jié)果如圖 1.12 所示。 %離散時間傅立葉變換 %離散時間信號 %模擬信號 圖 1.12 fs=5000Hz 時的離散信號及傅立葉變換曲線 對上面的 MATLAB 程序稍加修改,就可以得到 fs=1000Hz 時的 X 2( j ,程序運行結(jié)果如 圖 1.13 所示??梢钥闯?,當 fs=5000Hz 時, X 1( j 的曲線和例 1.14 中的 X a( j 曲線完全一 致, 沒有發(fā)

6、生混疊現(xiàn)象, 而當 fs=1000Hz 時,X 2( j 的曲線與例 1.14 中的 X a( j 曲線不同, 產(chǎn)生了混疊的現(xiàn)象。 圖 1.13 fs=1000Hz 時的離散信號及傅立葉變換曲線 例 1-13 對例 1-11 中產(chǎn)生的離散序列 x1(n和 x2(n,采用 sinc 函數(shù)進行內(nèi)插重構(gòu)。 解:程序清單如下 Ts1=0.0002;Fs1=1/Ts1;n1=-25:1:25;nTs1=n1*Ts1; x1=exp(-1000*abs(nTs1; Ts2=0.001;Fs2=1/Ts2;n2=-5:1:5;nTs2=n2*Ts2; x2=exp(-1000*abs(nTs2; Dt=0.

7、00005;t=-0.005:Dt:0.005; %模擬信號重構(gòu) %離散時間信號 xa1=x1*sinc(Fs1*(ones(length(nTs1,1*t-nTs1'*ones(1,length(t; xa2=x2*sinc(Fs2*(ones(length(nTs2,1*t-nTs2'*ones(1,length(t; subplot(2,1,1;plot(t*1000,xa1,'.'ylabel('xa1(t' title('從 x1(n重構(gòu)模擬信號 xa(t'hold on stem(n1*Ts1*1000,x1;hold off subplot(2,1,2;plot(t*1000,xa2,'.'ylabel('xa2(t' title('從 x2(n重構(gòu)模擬信號 xa(t'hold on stem(n2*Ts2*1000,x2;hold off 程序運行結(jié)果如圖 1.37 所示。 將圖 1.37 中的重構(gòu)后的模擬信號曲線與圖 1.34 中的原始模擬信號曲線進行比較, 可以 看出用離散信號 x1(n重構(gòu)出的模擬信號與原始信號誤差很小,而用離散信號 x2

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責。
  • 6. 下載文件中如有侵權(quán)或不適當內(nèi)容,請與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

最新文檔

評論

0/150

提交評論