| 查看: 3576 | 回复: 1 | |||
[求助]
求一个可运行的二次插值法（抛物线法）的matlab程序
|
|
谢谢大神了 发自小木虫IOS客户端 |
回复此楼» 猜你喜欢
哪位老哥知道今年的国自然具体哪一天放榜？
已经有3人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有37人回复
-
好奇怪的filecode
已经有7人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有43人回复
-
filecode
已经有9人回复
-
请教兼职经验
已经有3人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有5人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有13人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有11人回复
-
读博的好处
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【转帖】关于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
