24小时热门版块排行榜    

查看: 2448  |  回复: 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的回帖

loujing

铁杆木虫 (正式写手)

引用回帖:
13楼: Originally posted by xzczd at 2016-01-08 14:44:20
刚不放心所以又带具体的数算了一遍,如果各个式子的定义真的如你所写,那么分母各项的指数毫无疑问就是2j而不是j+1。...

您好,感谢您的回复。w_j最终的闭解,是原书上给出的。所以我想利用软件验证一下。

我在网上搜了一下,李航老师在他新浪博客上给出了勘误表,http://blog.sina.com.cn/s/blog_7ad48fee01017dpi.html,如截图。
勘误表上,原错误w_j的结果分子上是x_i,而我的书上是x_i^j。看来,我这本书是已经修订过的了(购于2015年6月,2015年3月印刷的版本)。
关于用Maple或Mathematica实现符号多项式拟合-3

您能否将您的代码给出,或将.nb以附件形式上传,让我研究一下,非常感谢。
14楼2016-01-08 17:03:58
已阅   回复此楼   关注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实现符号多项式拟合-4



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