24小时热门版块排行榜    

查看: 2440  |  回复: 22
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

loujing

铁杆木虫 (正式写手)

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

最近在看李航老师的统计学习方法,想用Maple或Mathematica验证一下第一章中的多项式拟合。
关于用Maple或Mathematica实现符号多项式拟合
关于用Maple或Mathematica实现符号多项式拟合-1
关于用Maple或Mathematica实现符号多项式拟合-2

看起来结果都不太对。
这是我异想天开了(这问题用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

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

xzczd

木虫 (小有名气)

【答案】应助回帖

引用回帖:
21楼: Originally posted by xzczd at 2016-01-08 21:36:58
说结束也算结束了,说没结束也没结束,4楼的推导在这之后又多了一步,是把之前的结果写成了矩阵乘法的形式,这个其实是矩阵点乘规则的直接应用,,但是如果你想把这步也交给软件推导的话,我可以告诉你这个规则又是 ...

啊,第一段代码漏贴了eqn的定义:
CODE:
eqn = (D[l /. rulej1[w@j], w[j]] /. {rulej2, rulej3} //. {ruleExpand, ruleDistribute,
      ruleExtract} // Expand) == 0

小木虫Mathematica版块已毁(当然原本也不咋的),建议大家前往百度贴吧或Stackexchange。
22楼2016-01-08 21:45:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 23 个回答

赵梦92

木虫 (著名写手)

大前提错了吧。我的书上和你的不一样呢

发自小木虫Android客户端
2楼2016-01-08 12:31:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

loujing

铁杆木虫 (正式写手)

引用回帖:
2楼: Originally posted by 赵梦92 at 2016-01-08 12:31:47
大前提错了吧。我的书上和你的不一样呢

谢谢答复,损失函数L(w)是没问题的,所以我就想用Maple或者Mathematica来验证一下w_j是否是这个结果。
3楼2016-01-08 12:34:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

赵梦92

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
loujing: 金币+20, 有帮助 2016-01-08 13:09:37
我的是这个,你看一下吧
关于用Maple或Mathematica实现符号多项式拟合-3



发自小木虫Android客户端
4楼2016-01-08 12:35:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见