Re: [請益] 算是訊號與系統的問題
你是用 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
討論串 (同標題文章)
本文引述了以下文章的的內容:
完整討論串 (本文為第 2 之 3 篇):