MATLAB 板


LINE

剛剛無聊推導了一下 線性回歸方程式 令理想的直線方程式為 y=ax+b 而我們手上有的資料點為 yi=ax+b+c 此時c為描述實際資料與模型的誤差 我們目標是令誤差越小越好,有很多種求法 最簡單的就是最小平方法 所以我們的目標函數建立如下 Q = Σ(yi-axi-b)^2 做偏微分運算 ∂/∂a (Q) = 2*Σ[(yi-axi-b)xi] = 0 ∂/∂b (Q) = 2*Σ(yi-axi-b) = 0 整理方程式 Σ(xi^2) * a + Σ(xi) * b = Σ(xiyi) Σ(xi) * a + n * b = Σ(yi) 解聯立,克拉瑪簡單解一下 |Σ(xi^2) Σ(xi)| | a | | Σ(xiyi)| | | | | = | | | Σ(xi) n | | b | | Σ(yi) | Δ = n*Σ(xi^2) - Σ(xi)Σ(xi) Δa = n*Σ(xiyi) - Σ(xi)Σ(yi) Δb = Σ(xi^2)*Σ(yi) - Σ(xiyi)*Σ(xi) a = [n*Σ(xiyi) - Σ(xi)Σ(yi)] / [n*Σ(xi^2) - Σ(xi)Σ(xi)] b = [Σ(xi^2)*Σ(yi) - Σ(xiyi)*Σ(xi)] / [n*Σ(xi^2) - Σ(xi)Σ(xi)] 有更簡單的形式,不過原則點到此 --------------------------------------- 好der,剛剛只是複習,接下來回到你的問題 Q = Σ(yi-axi)^2 ∂/∂a (Q) = 2*Σ[(yi-axi)xi] = 0 Σ(xiyi) - aΣ(xi^2) = 0 a = Σ(xiyi)/Σ(xi^2) 故得證 ※ 引述《fxxkjoe1231 (糟糕)》之銘言: : 大家好,我手邊有兩筆資料。 : 按照理論的推導,這兩筆資料應為y=ax的關係, : 然而實際在進行回歸的時候仍會有截距的存在, : 請問該如何強迫截距為零呢? : 謝謝大家 --



※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.113.125.75
※ 文章網址: http://webptt.com/m.aspx?n=bbs/MATLAB/M.1416681060.A.461.html
1F:推 fxxkjoe1231: 非常感謝! 11/23 10:16







like.gif 您可能會有興趣的文章
icon.png[問題/行為] 貓晚上進房間會不會有憋尿問題
icon.pngRe: [閒聊] 選了錯誤的女孩成為魔法少女 XDDDDDDDDDD
icon.png[正妹] 瑞典 一張
icon.png[心得] EMS高領長版毛衣.墨小樓MC1002
icon.png[分享] 丹龍隔熱紙GE55+33+22
icon.png[問題] 清洗洗衣機
icon.png[尋物] 窗台下的空間
icon.png[閒聊] 双極の女神1 木魔爵
icon.png[售車] 新竹 1997 march 1297cc 白色 四門
icon.png[討論] 能從照片感受到攝影者心情嗎
icon.png[狂賀] 賀賀賀賀 賀!島村卯月!總選舉NO.1
icon.png[難過] 羨慕白皮膚的女生
icon.png閱讀文章
icon.png[黑特]
icon.png[問題] SBK S1安裝於安全帽位置
icon.png[分享] 舊woo100絕版開箱!!
icon.pngRe: [無言] 關於小包衛生紙
icon.png[開箱] E5-2683V3 RX480Strix 快睿C1 簡單測試
icon.png[心得] 蒼の海賊龍 地獄 執行者16PT
icon.png[售車] 1999年Virage iO 1.8EXi
icon.png[心得] 挑戰33 LV10 獅子座pt solo
icon.png[閒聊] 手把手教你不被桶之新手主購教學
icon.png[分享] Civic Type R 量產版官方照無預警流出
icon.png[售車] Golf 4 2.0 銀色 自排
icon.png[出售] Graco提籃汽座(有底座)2000元誠可議
icon.png[問題] 請問補牙材質掉了還能再補嗎?(台中半年內
icon.png[問題] 44th 單曲 生寫竟然都給重複的啊啊!
icon.png[心得] 華南紅卡/icash 核卡
icon.png[問題] 拔牙矯正這樣正常嗎
icon.png[贈送] 老莫高業 初業 102年版
icon.png[情報] 三大行動支付 本季掀戰火
icon.png[寶寶] 博客來Amos水蠟筆5/1特價五折
icon.pngRe: [心得] 新鮮人一些面試分享
icon.png[心得] 蒼の海賊龍 地獄 麒麟25PT
icon.pngRe: [閒聊] (君の名は。雷慎入) 君名二創漫畫翻譯
icon.pngRe: [閒聊] OGN中場影片:失蹤人口局 (英文字幕)
icon.png[問題] 台灣大哥大4G訊號差
icon.png[出售] [全國]全新千尋侘草LED燈, 水草

請輸入看板名稱,例如:Tech_Job站內搜尋

TOP