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/cn.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灯, 水草

请输入看板名称,例如:Soft_Job站内搜寻

TOP