| 查看: 3586 | 回复: 1 | ||
[求助]
求一个可运行的二次插值法（抛物线法）的matlab程序
|
|
谢谢大神了 发自小木虫IOS客户端 |
回复此楼» 猜你喜欢
建议基金发布提前给出明确的时间点
已经有8人回复
-
看来今天不会放榜了？
已经有11人回复
-
时间戳又变了
已经有19人回复
-
filecode，4个jtjc了
已经有14人回复
-
让我中一个面上吧！
已经有13人回复
-
人气不行了
已经有7人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有18人回复
-
感觉是下周放榜了
已经有12人回复
-
我面上完蛋了
已经有8人回复
-
只有每年这种时候来逛逛小木虫
已经有25人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【转帖】关于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
40