2014-06-08 31 views
0

各種距離的時候,我有一個信號x(t)與它的價值計算信號的PSD只是想:如何在MATLAB

t = 0, 1, 3, 4, 7, 8, 10, 15, ... 
x = 3, 4, 5, 4, 6, 7, 4, 8, ... 

,你可以看到,這個問題herre之間的時間距離2任何不相等的樣本。你能告訴我如何使用matlab來繪製它的psd嗎? 非常感謝。

+0

感覺就像一門功課來證明FFT的理解。與其僅僅調用一個fft函數,學生應該理解這個基本思想,即它是用正弦曲線計算數據乘積積分的快速方法。您可以在任何點評估正弦曲線,但問題變成如何計算離散積分。定義矩形或梯形來表示在每個點處集成的體積。此時您可能會認識到重新採樣或插值輸入的相似性。 interp1(有些人可能會考慮零填寫。) – dennis

+0

謝謝你回答我的問題。這不是理解fft的學生作業。我工作在我的projet上,結果是一個數據文件,這是我在該文件上看到的,我不知道它。 –

+0

如果您對信號一無所知,您可以嘗試對'x'進行插值,使其在統一網格中進行採樣。查看['interp1'](http://www.mathworks.com/help/matlab/ref/interp1.html)(或者在matlab中鍵入'help interp1') – ThP

回答

0

首先,我會重新採樣插值數據以獲得一個有規律的採樣信號。對於插值算法,我使用樣條曲線,但也可以使用interp1,這會更快;

t2 = t(1):step:t(end); %with this I have a t signal with every sample separated by step 
x2 = spline(t, x, t2); 

然後你可以用fftpwelch PSD您的X2信號...