你是用 Matlab 模擬嗎? 如果我猜的沒錯, 問題可能是出在你沒有做 coherent sampling 假設輸入信號 v(t) 是一個頻率為 fin 的 sine wave fft 取 N 個點, 每個點之間的時間差為 T (i.e.取樣頻率為1/T) coherent sampling 必須滿足下面的條件: fin * T = m / N (m 為整數, 能選個質數更好!) 為何 delta function 不會有這個問題? 因為它的 Fourier transform 是一個常數, 就不會有取樣週期的問題啦~ ※ 引述《jockerv (嗯嗯)》之銘言： : 先是簡單的v(t)，r，c串聯的一個迴路 : i(t)是流過這個迴路的電流 : 先以Laplace分析：s-domain : I(s)=(1/r) x V(s) x (s/(s+a)) a=1/(r*c),x為乘法 : 轉回time-domain： : i(t)=(1/r) x V(t) * (-a x exp(-at)) * 為convolution : 驗証I(s)是否等同於v(t)在time-domain上進行convolution : 即I(s)=(1/r) x [V(s) x H(s)] : i(t)=(1/r) x [v(t) * h(t)] : 驗証方法： : 輸入訊號v(t)為非連續函數，故V(s)=fft(v(t)).......matlab語法 : H(s)=fft(h(t)) : I(s)為V(s)與H(s)點乘結果再除以r : 將I(s)以ifft轉換回time-domain後，與i(t)進行比較 : 問題： : 1.比較後，兩者圖形有差(差不少)，目前我還找不出原因= = : 2.只有當v(t)為δ(t)形式的訊號，才會使結果相差不大 : 從其它地方extract出來的訊號都不行 : 原因： : 1.會不會是因為h(t)是以laplace得來，所以在進行H(s)=fft(h(t)) : 會有不符合的情況？(fourier跟laplace的積分範圍不同) : 2.想的方法跟本上就錯了？@@ : 希望有人能給我些建議，或是討論 : 謝謝 --

※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 67.186.57.226