| 查看: 3491 | 回复: 1 | ||
[求助]
求一个可运行的二次插值法（抛物线法）的matlab程序
|
|
谢谢大神了 发自小木虫IOS客户端 |
回复此楼» 猜你喜欢
参与限项
已经有3人回复
-
假如你的研究生提出不合理要求
已经有7人回复
-
实验室接单子
已经有4人回复
-
全日制（定向）博士
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
不自信的我
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【转帖】关于NP，NP-hard，P，NPC等相关问题的讨论
已经有13人回复
- 【软件】超强数学公式编辑器MathType 能够自动编号
已经有23人回复
- Mathtype使用技巧（转帖）
已经有1人回复
mlypplover
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 5823.4
- 散金: 50
- 红花: 3
- 帖子: 219
- 在线: 21.6小时
- 虫号: 3688030
- 注册: 2015-02-13
- 性别: GG
- 专业: 控制理论与方法
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
jsjs0827: 金币+10, ★有帮助 2016-01-03 10:15:10
jsjs0827: 金币+10, ★有帮助 2016-01-03 10:15:10
|
举个例子吧 已知插值点(-2.00,17.00), (0.00,1.00), (1.00,2.00), (2.00,17.00),求三次插值多项式，并计算 . 其Matlab程序如下： function yy=M10(x,y,xx) m=length(x);n=length(y); if m~=n, error('向量x与y的长度必须一致');end s=0; for i=1:n t=ones(1,length(xx)); for j=1:n if j~=i t=t.*(xx-x(j))/(x(i)-x(j)); end end s=s+t*y(i); end yy=s; 在Matlab窗口中执行： x=[-2.00 0.00 1.00 2.00]; y=[17.00 1.00 2.00 17.00]; xx=0.6; M10(x,y,xx) 结果如下： ans = 0.25600000000000 |
2楼2015-12-07 16:04:03