[問題] FFT的相位已回收
關於FFT的問題
如果設定一個gauss函數y且相位是0
做Y=ifft(y)後取unwrap(angle(Y))的值應該也要是0
但運算出來的值卻有偏差
請問是程式的問題，還是內建angle函數的問題?
程式碼如下
clc;clear;close all
t=0:1e-15:2e-12;
sig=70e-15;
c0=1e-12;
y=gaussmf(t,[sig c0]);
Y=fft(y);
Y=fftshift(Y);
figure
subplot(2,2,1)
plot(abs(y));
title('temporal amplitude')
subplot(2,2,2)
plot(abs(Y));
title('spectrum amplitude')
subplot(2,2,3)
plot(unwrap(angle(y)));
title('temporal phase')
subplot(2,2,4)
plot(unwrap(angle(Y)));
title('spectrum phase')
謝謝
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.117.162.104
推
09/29 23:27, , 1F
09/29 23:27, 1F
→
09/29 23:27, , 2F
09/29 23:27, 2F
→
10/01 13:44, , 3F
10/01 13:44, 3F
討論串 (同標題文章)