查看: 612  |  回复: 6
【悬赏金币】回答本帖问题,作者CPYHH将赠送您 100 个金币

CPYHH

金虫 (著名写手)

[求助] 求Matlab代码的改进已有1人参与

请问这个matlab代码有什么需要改进的地方?

计算lagrange插值多项式的代码,输入x=[1 2 3 4 ],y=[1 4 9 16] 和 x0=1.5以及 x=[1 2 3 4 ],y=[1 4 9 16]和x0=2.5输出都是y0=30,请问怎么改进?
代码如下:
function y0 = lagrange2(x,y,x0 )
n=length(x);
l=ones(1,n);
for k=1:n
for j=1:n
if j~=k
l(k)=l(k)*(x0-x(j))/(x(k)-x(j));
end
end
end
y0=sum(y.*1);
end

求Matlab代码的改进


求Matlab代码的改进-1


求Matlab代码的改进-2


发自小木虫Android客户端
回复此楼

» 猜你喜欢

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

sspa9999

至尊木虫 (著名写手)

你算了一通 l(k) 后面完全没有用到,y0 y

发自小木虫Android客户端
2楼2020-10-02 21:30:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sspa9999

至尊木虫 (著名写手)

y0 只与y相关,y不变y0就不会变,这不是很正常吗?

发自小木虫Android客户端
3楼2020-10-02 21:32:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sspa9999

至尊木虫 (著名写手)

另外这个不叫改进,程序完全是错的,只能叫修错

发自小木虫Android客户端
4楼2020-10-02 21:34:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sspa9999

至尊木虫 (著名写手)

y0那个式子里,我猜是y与L相乘,而不是与1相乘

发自小木虫Android客户端
5楼2020-10-03 02:02:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hzlhm

至尊木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
题主给出计算Lagrange插值多项式的代码,经计算当x0=1.5和x0=2.5输出都是y0=30,从结果来看明显有错。出错的原因有:
1、没有完全理解拉格朗日插值公式的结构,该公式是两部分组成,即求和部分,求积部分
2、虽然题主考虑了用双循环语句,但没有明确计算初值,即求和的初值以及求积的初值
3、在双循环语句里,求和应采用累加的形式求解,在这里使用用sum()函数不是太合理,即为
s=p*y0(k)+s; %s为求和的值
4、完善代码,运行可以得到如下结果。
求Matlab代码的改进-3
123.jpg

QQ:2120156492
6楼2020-10-03 12:44:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

CPYHH

金虫 (著名写手)

引用回帖:
5楼: Originally posted by sspa9999 at 2020-10-03 02:02:05
y0那个式子里,我猜是y与L相乘,而不是与1相乘

你是对的

发自小木虫Android客户端
7楼2020-10-03 14:12:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 CPYHH 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 有的审稿人太无底线了 +26 流浪的YANG 2022-06-29 31/1550 2022-07-05 12:49 by cfkybfq
[论文投稿] under review 45天还没找到审稿人正常吗? +12 20食品 2022-07-04 21/1050 2022-07-05 09:25 by absu
[文学芳草园] 望果节的青稞‖藏地之行 +5 幽兰星空 2022-07-02 14/700 2022-07-05 09:03 by ahsoarli
[考博] 硕士毕业,六级已过,一篇一作cell reports,一篇二作advanced science +3 dong10 2022-07-04 4/200 2022-07-05 00:31 by 宁宝胖虎
[论文投稿] ACS AMI 20+3 DDDDD小白 2022-06-29 10/500 2022-07-04 17:18 by DDDDD小白
[访问学者] 2022年访学结果出来了,没有通过 +21 lanpin 2022-06-30 32/1600 2022-07-04 17:10 by wnyq3767
[基金申请] 博士后特别资助基金 +16 linxiyuan68 2022-06-30 25/1250 2022-07-04 15:29 by 怅然若失的
[基金申请] 武汉大学2003级博士研究生抄袭,现任图书馆纪委书记 +9 苏东坡二世 2022-07-03 9/450 2022-07-04 11:18 by shl2112501
[论文投稿] 审稿人审稿慢 +11 明月清风糖 2022-06-30 14/700 2022-07-04 09:37 by 心诚则灵3
[教师之家] 北京的女大学老师征婚 (EPI+-1)(金币-50) +9 木兰2022 2022-07-03 11/550 2022-07-04 08:49 by mydandelion
[教师之家] 不用辞职的读博,现在没有了吗? +8 protrans 2022-06-28 9/450 2022-07-04 02:41 by protrans
[有机交流] 求购甲基磺酰氯 +4 笑看人生1993 2022-06-29 7/350 2022-07-04 00:04 by 笑看人生1993
[考博] 23申博什么时候联系导师合适 +13 ylwwww 2022-07-03 15/750 2022-07-03 21:47 by Joey*
[考博] 2023考博求助 +3 一季繁华 2022-07-03 8/400 2022-07-03 18:43 by 一季繁华
[基金申请] 工材面上项目上会名单出来了吗? +4 mrkang 2022-07-01 5/250 2022-07-02 17:22 by jxwsytx
[论文投稿] 编辑这是大修还是小修的意思? +8 liao8540 2022-07-01 8/400 2022-07-02 13:38 by Nano2016
[考博] 马上研三了有点迷茫 +5 宁静36 2022-07-01 10/500 2022-07-01 23:19 by 复旦彭超
[功能材料] 光照自由基测试 +4 zuocuiping 2022-06-29 5/250 2022-07-01 19:04 by hgxt3368
[论文投稿] 论文大修 +9 明月清风糖 2022-06-30 13/650 2022-07-01 16:06 by QZ1205
[教师之家] 材料纳米论文漫天飞,不是为了材料研究和应用,只是为了影响因子和帽子。 +4 有余12 2022-06-30 4/200 2022-07-01 03:41 by 小木脑虫
信息提示
请填处理意见