從歐拉公式看三角波的單邊譜與雙邊譜

1、在有限區間上,凡滿足狄利赫裏條件的周期函數都可以展開爲傅里葉級數。其三角函數形式的頻譜爲單邊譜復指數函數形式的頻譜爲雙邊譜

在這裏插入圖片描述
參考:https://wenku.baidu.com/view/af2ea95b7dd184254b35eefdc8d376eeaeaa17de.html

實數系列的FFT後的X(k)是共軛對稱的,即實部偶對稱,虛部奇對稱,matlab fft()得到的是雙邊譜,前部分是正頻率部分,後部分是負頻率部分。

2、看看fft()後的頻率

x=[1,1.8,1.8,1];
n=length(x);
fs=8;
dt=1/fs;
ti=0:dt:(n-1)*dt;
nfft=8;
y=fft(x,nfft);
df=fs/nfft;
fi=0:df:(nfft-1)*df;

subplot(312)
plot(fi,real(y))
title('real')
xlabel('f(Hz)')
subplot(313)
plot(fi,imag(y))
title('imag')
xlabel('f(Hz)')
x_ifft=ifft(y);
x_ifft_output=(x_ifft(1:n));
subplot(311)
plot(ti,x_ifft_output)
title('x(t)')
xlabel('t(s)')

figure
%Determin Nyquist frequency:
Nyq=fs/2;
%Save frequency array:
fi_array=(-Nyq:df:Nyq-df);%fftshift()後的頻率 (0點右移Nyq個點!!)
y_shift=fftshift(y);
subplot(312)
plot(fi_array,real(y_shift))
title('real')
xlabel('f(Hz)')
subplot(313)
plot(fi_array,imag(y_shift))
title('imag')
xlabel('f(Hz)')
subplot(311)
plot(ti,x)
title('x(t)')
xlabel('t(s)')

在這裏插入圖片描述
可見,matlab fft()後的頻率爲:[0:df:Nqy-df,-Nqy:df:-df]是其正負頻率。

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章