[問題] convolution的結果產生了shift已回收

看板MATLAB作者 (查無暱稱)時間16年前 (2010/02/27 21:55), 編輯推噓1(102)
留言3則, 2人參與, 最新討論串1/2 (看更多)
※ [本文轉錄自 Math 看板] 作者: olys (查無暱稱) 看板: Math 標題: [微積] convolution的結果產生了shift 時間: Sat Feb 27 14:57:44 2010 大家好,第一次在這裡發問。 我最近在作訊號平滑的研究, 手上有一維訊號f(x),還有一高斯常態分布函數g(x)(mu值設為0,sigma由演算法算出來) 平滑的內容就是將f和g作convolution,然後取中間1/2段的部分當作結果h(x) 但是不知道為什麼作出來的h(x)的y值和原來訊號f(x)的y值差很遠, 找不到是什麼數學原因讓h(x)產生了類似shift的現象。 但平滑化倒是成功了 底下是圖,紅色是f(x),藍色是平滑化後的h(x) http://olys.myweb.hinet.net/gsmooth.jpg
-- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 114.47.154.205

02/27 15:52,
你用的Gaussian filter的方程式是怎麼樣?
02/27 15:52
matlab的normpdf http://www.mathworks.com/access/helpdesk/help/toolbox/stats/normpdf.html ※ 編輯: olys 來自: 114.47.154.205 (02/27 16:01)

02/27 16:02,
先感謝關注,如果需要補充清楚的部分請盡量提出
02/27 16:02

02/27 16:13,
你應該是用MATLAB吧 先把g normalize=> g = g/sum(g)
02/27 16:13

02/27 19:27,
剛剛去試, 還是一樣shift到下面
02/27 19:27

02/27 21:55,
我轉到matlab版也問問看好了, 有結果會回來報
02/27 21:55
-- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 140.127.208.68

02/28 02:57, , 1F
shift? 你先說說看戈為什麼你取中間1/2段吧?
02/28 02:57, 1F

02/28 13:49, , 2F
我是參考別人寫的程式,而且convn後矩陣長度是m+n-1也無法和
02/28 13:49, 2F

02/28 13:50, , 3F
直接作比較。另外從圖中也能看到取中間的確波型是被平滑成功
02/28 13:50, 3F
※ 編輯: olys 來自: 140.127.208.68 (02/28 13:51)
文章代碼(AID): #1BYIHMqp (MATLAB)
文章代碼(AID): #1BYIHMqp (MATLAB)