※ 引述《Maddulin (what else do u focus?)》之銘言： : ※ 引述《flybow (yyyyyyyy)》之銘言： : : 對不起，我對訊號處理真的不是很懂，有些問題想請問一下 : : 看了書，還是有很多不懂的地方 : : 我的系統是一個在水中的散射體 : : 時間t 0：0.1：20 : : 我的input: sin(0.5pi x t) : : 進入一個系統 output為 0.3sin(0.5pi x t) : : 能量和振幅有衰減，我想算這個系統的脈衝響應h(t) : : 我是先求轉移函數H（w)= fft(output,2048)/fft(input,2048) : : 在逆傅立業轉換可得 h(t)=ifft(H(w)) : : 我驗證 input ＊ h(t) ＝output 成功 : : 那麼我再模擬另一種輸入訊號 兩個不同頻率弦波的組合 : : 我的input: sin(0.3pi x t) ＋ sin(0.6pi x t) : : 進入系統後 output: 0.5*sin(0.3pi x t) + 0.24*sin(0.6pi*t) : : 越高頻的衰減越多 : : （看書上說，如果abs(H(w))不等於1，會有振幅失真） : : 同樣利用上述的方法 : : H（w)= fft(output,2048)/fft(input,2048) : : h(t)=ifft(H(w)) : : 我的h(t)在大的脈衝值後面接一堆對時間的漣漪震盪，真的很奇怪 : : 文獻說這就是Gibbs 現象？可是並沒有說原因 : : 請問這是為什麼，那些漣漪的產生的原因為何？ : : 有沒有辦法克服這個後面的漣漪呢？？ : : 請問各位一下，不好意思 : : （我看到文獻是用窗口函數截下這個大的振幅值，後面漣漪去掉，這樣截函數， : : 不是會對頻域造成影響嗎？看書是說會有頻譜洩漏的問題？） : 首先一個簡單的算術問題 0*a=0 請問如何得a : 你懂我問的問題嗎? 我只想先和你說 你這樣算的H 並非真正系統的 H : 除了特定的頻率外,你了解嗎? : 另外以DFT計算FT有一件條件,我想我就不多說,即然你用就代表你懂 : 所以你的輸入經過fft轉換後應當是spark,但數值計算僅能近似 : 然後再回到我第一個問題,你能不能明白我的point在那? : 所以你最好先確定你這樣的運算沒有問題,剩下的問題才有意義 第一個問題我真的回答不出來， 不過我確實是要算特定頻率的散射體的h(t) 我求出來確實有個SPARK，我翻了很多書找了許久始終沒找到解決我這個問題的 關鍵。 想了想，我確實還是想不出第一個問題的核心跟關於這個問題的相關性， 能不能相煩說的詳細些。 不過真的謝謝你的回答 --

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