24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

loujing

铁杆木虫 (正式写手)

应助: 19 (小学生)
金币: 5835.5
散金: 2028
红花: 8
帖子: 322
在线: 110.2小时
虫号: 58180
注册: 2005-02-20
性别: GG
专业: 计算机应用技术

[求助] 关于用Maple或Mathematica实现符号多项式拟合已有1人参与

最近在看李航老师的统计学习方法，想用Maple或Mathematica验证一下第一章中的多项式拟合。

看起来结果都不太对。
这是我异想天开了（这问题用Maple和Mathematica是无法解决的），还是哪里写错了，万望大家指教，十分感谢。

Mathematica：

CODE:

f[Subscript[w, j]] = 1/2 \!\(\*UnderoverscriptBox[\(\[Sum]\), \(i = 1\), \(n\)]\*SuperscriptBox[\((\*UnderoverscriptBox[\(\[Sum]\), \(j = 0\), \(m\)]\*SubscriptBox[\(w\), \(j\)] \*SubsuperscriptBox[\(x\), \(i\), \(j\)] - \*SubscriptBox[\(y\), \(i\)])\), \(2\)]\)

Maple:

CODE:

L := proc (w__j) options operator, arrow; (1/2)*(sum((sum(w__j*x__i^j, j = 0 .. m)-y__i)^2, i = 1 .. n)) end proc

回复此楼

» 猜你喜欢

Cas 72-43-5需要30g，定制合成，能接单的留言已经有8人回复
求助:我三月中下旬出站，青基依托单位怎么办？已经有6人回复
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？已经有8人回复
磺酰氟产物，毕不了业了！已经有5人回复
论文终于录用啦！满足毕业条件了已经有25人回复
2026年机械制造与材料应用国际会议（ICMMMA 2026）已经有3人回复
自荐读博已经有3人回复
不自信的我已经有5人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有8人回复

1楼 2016-01-08 11:33:37

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

xzczd

木虫 (小有名气)

应助: 72 (初中生)
金币: 1941.4
红花: 11
帖子: 290
在线: 144.2小时
虫号: 3305534
注册: 2014-07-03

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
loujing: 金币+30, ★★★很有帮助 2016-01-08 23:56:20

引用回帖:

17楼: Originally posted by loujing at 2016-01-08 17:28:16
原书截图

...

……刚发现我为了凑出和你的顶楼相似的结果犯了个错误。盯着四楼的图看了半天（真伤眼），现在我基本确定，正确的推导是：

CODE:

f[i_, m_] := Sum[w[j]*x[i]^j, {j, 0, m}]

l = 1/2 Sum[(f[i, m] - y[i])^2, {i, n}];

rulej1[w_] := Sum[expr_, {i_, downup__}] /; ! FreeQ[w, i] :> sum[{i, downup}][expr]

rulej2 = sum[_]'[_] :> 1;

rulej3 = sum[i_][expr_] :> Sum[expr, i];

ruleExpand = Sum[(a_ + b_) c_, i_] :> Sum[(a + b) c // Expand, i];

ruleDistribute = Sum[expr1_ + expr2_, i_] :> Sum[expr1, i] + Sum[expr2, i];

ruleExtract = Sum[a_ b_., {i_, downup__}] /; FreeQ[a, i] :> a Sum[b, {i, downup}];

(D[l /. rulej1[w@j], w[j]] /. {rulej2, rulej3} //. {ruleExpand, ruleDistribute, 

     ruleExtract} // Expand) == 0