函数的插值方法及matlab程序
4星 · 超过85%的资源 需积分: 5 198 浏览量
2011-12-06
16:05:30
上传
评论 2
收藏 590KB DOC 举报 
6.1 插值问题及其误差
6.1.2 与插值有关的 MATLAB 函数
(一) POLY2SYM 函数
调用格式一：poly2sym (C)
调用格式二：f1=poly2sym(C,'V') 或 f2=poly2sym(C, sym ('V') ),
(二) POLYVAL 函数
调用格式：Y = polyval(P,X)
(三) POLY 函数
调用格式：Y = poly (V)
(四) CONV 函数
调用格式：C =conv (A, B)
例 6.1.2 求三个一次多项式 、 和 的积 .它们的零点
分别依次为 0.4,0.8,1.2.
解 我们可以用两种 MATLAB 程序求之.
方法 1 如输入 MATLAB 程序
>> X1=[0.4,0.8,1.2]; l1=poly(X1), L1=poly2sym (l1)
运行后输出结果为
l1 =
1.0000 -2.4000 1.7600 -0.3840
L1 =
x^3-12/5*x^2+44/25*x-48/125
方法 2 如输入 MATLAB 程序
>> P1=poly(0.4);P2=poly(0.8);P3=poly(1.2);
C =conv (conv (P1, P2), P3) , L1=poly2sym (C)
运行后输出的结果与方法 1 相同.
(五) DECONV 函数
调用格式：[Q,R] =deconv (B,A)
(六) roots(poly(1:n))命令
调用格式：roots(poly(1:n))
(七) det(a*eye(size (A)) - A)命令
调用格式：b=det(a*eye(size (A)) - A)
6.2 拉格朗日（Lagrange）插值及其 MATLAB 程序
6.2.1 线性插值及其 MATLAB 程序
例 6.2.1 已知函数 在 上具有二阶连续导数， ，且满足条件
.求线性插值多项式和函数值 ，并估计其误差.
解 输入程序
>> X=[1,3];Y=[1,2]; l01= poly(X(2))/( X(1)- X(2)), l11=
poly(X(1))/( X(2)- X(1)), l0=poly2sym (l01),l1=poly2sym
(l11), P = l01* Y(1)+ l11* Y(2),
L=poly2sym (P),x=1.5; Y = polyval(P,x)
运行后输出基函数 l
0
和 l
1
及其插值多项式的系数向量 P（略）、插值多项式 L 和插值 Y 为
l0 = l1 = L = Y =
-1/2*x+3/2 1/2*x-1/2 1/2*x+1/2 1.2500
输入程序
>> M=5;R1=M*abs((x-X(1))* (x-X(2)))/2
剩余30页未读，继续阅读
资源评论
LJPZGX19852013-04-09还好 可以参考一下- DeLiang_Y_32014-05-06不错 值得参考学习
- brieflone2013-04-25还可以，作为参考有一定帮助
yinwenbin163
- 粉丝: 0
- 资源: 1
