用matlab實現(xiàn)功率譜分析_第1頁
用matlab實現(xiàn)功率譜分析_第2頁
用matlab實現(xiàn)功率譜分析_第3頁
全文預(yù)覽已結(jié)束

下載本文檔

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

文檔簡介

1、EXERCISE 208051302功率譜分辨率(A )分別產(chǎn)生兩個離散時間(以10赫茲的采樣率，100個時間點)正弦函數(shù)的抽樣 函數(shù)，表達(dá)式為 X(t)=2.6sin(4.2nt + 和 Y(t) = 2.1 sin(4.4nt + 0),其中 U(0, 4n)，0U(-n,n)且獨立，分別將函數(shù)保存在數(shù)據(jù)文件m10_3x. dat、 m10_y. dat,中；Z(t) = X(t) + Y(t),讀取兩個數(shù)據(jù)文件中的記錄數(shù)據(jù)，添加相應(yīng)的條件產(chǎn)生Z(t) 的記錄數(shù)據(jù)并保存在m10_3z.dat ；使用pwelch估測和繪制Z(t)的功率譜； 重復(fù)上述，用10赫茲的采樣率，1000個時間點，調(diào)

2、整序列長度；根據(jù)以上結(jié)果討論Z(t)的頻譜分辨率。程序如下：1、以10赫茲的采樣率，100個時間點時fai= unifrnd (0,6.28,1, 1);seita= unifrnd (-3.14,3.14,1, 1);t=0:0.1:9.9;a=4.2*3.14*t + fai;X=2.6*sin(a); %產(chǎn)生 X(t)savefile = .m10_3x.dat;save(savefile, X);Y= 2.1*sin(4.4*3.14*t+seita); %產(chǎn)生 Y( t)savefile = .m10_3y. dat;save(savefile, Y);Z= X+ Y; %產(chǎn)生 Z(

3、 t)savefile = .m10_3z.dat;save(savefile, Z);pwelch(Z);繪得功率譜如下:20o o o o o o 1-1-2-3-4 a-dEB3._-w2.mp) JCJMrlbs_-堂者CL2、用10赫茲的采樣率，1000個時間點時fai= unifrnd (0,6.28,1, 1);seita= unifrnd (-3.14,3.14,1, 1);t=0:0.1:99.9;a=4.2*3.14*t + fai;X=2.6*sin(a); %產(chǎn)生 X(t)savefile = .m10_3x.dat;save(savefile, X);Y= 2.1*sin(4.4*3.14*t+seita); %產(chǎn)生 Y( t)savefile = .m10_3y. dat;save(savefile, Y);Z= X+ Y; %產(chǎn)生 Z (t)savefile = .m10_3z.dat;save(savefile, Z);pwelch(Z);功率譜如下

溫馨提示

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

評論

0/150

提交評論