实验1帮助文件插值和拟合案例.docVIP

1.1 插值问题及其误差 与插值有关的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 例 1 求三个一次多项式、和的积.它们的零点分别依次为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 1.2 拉格朗日（Lagrange）插值及其MATLAB程序 一．线性插值及其MATLAB程序 例2 已知函数在上具有二阶连续导数，，且满足条件.求线性插值多项式和函数值，并估计其误差. 解 输入程序 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 运行后输出基函数l0和l1及其插值多项式的系数向量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 运行后输出误差限为 R1 1.8750 例3 求函数e在上线性插值多项式，并估计其误差. 解 输入程序 X [0,1]; Y exp -X , 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 , 运行后输出基函数l0和l1及其插值多项式的系数向量P和插值多项式L为 l0 l1 P -x+1 x -0.6321 1.0000 L -1座机电话号码/座机电话号码1座机电话号码*x+1 输入程序 M 1;x 0:0.001:1; R1 M*max abs x-X 1 .* x-X 2 ./2 运行后输出误差限为 R1 0.1250. 二．抛物线插值及其MATLAB程序 例4 求将区间 [0, π/2] 分成等份，用产生个节点，然后分别作线性插值函数和抛物线插值函数.用它们分别计算cos π/6 取四位有效数字 ，并估计其误差. 解 1 线性插值 输入程序 X [0,pi/2]; Y cos X , 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 pi/6; Y polyval P,x 运行后输出基函数l0和l1及其插值多项式的系数向量P、插值多项式和插值为 l0 -座机电话号码/座机电话号码*x+1 l1 座机电话号码/座机电话号码*x P -0.6366 1.0000 L -座机电话号码/座机电话号码*x+1 Y 0.6667 输入程序 M 1;x pi/6; R1 M*abs x-X 1 * x-X 2