matlab数值分析程序--高等数学,数值代数的matlab实现-文字版, matlab电子书, 和matlab 有关的电子书:

6.2.6 拉格朗日插值及其误差估计的MATLAB程序

6.2.6  拉格朗日插值及其误差估计的MATLAB程序

拉格朗日插值及其误差估计的MATLAB主程序

function [y,R]=lagranzi(X,Y,x,M)

n=length(X); m=length(x);

for i=1:m

   z=x(i);s=0.0;

   for k=1:n

       p=1.0; q1=1.0; c1=1.0;

for j=1:n

         if j~=k

p=p*(z-X(j))/(X(k)-X(j));

         end

      q1=abs(q1*(z-X(j)));c1=c1*j;

      end

      s=p*Y(k)+s;

   end

   y(i)=s;

end

R=M*q1/c1;

 

例 6.2.6   已知,用拉格朗日插值及其误差估计的MATLAB主程序求的近似值,并估计其误差.

解  在MATLAB工作窗口输入程序

>> x=2*pi/9; M=1; X=[pi/6 ,pi/4, pi/3];

Y=[0.5,0.7071,0.8660]; [y,R]=lagranzi(X,Y,x,M)

运行后输出插值y及其误差限R为

y =                      R =

0.6434                 8.8610e-004.

欢迎转载,转载请注明来自一手册:http://yishouce.com/book/3/3060206.html
友情链接It题库(ittiku.com)| 版权归yishouce.com所有| 友链等可联系 admin#yishouce.com|粤ICP备16001685号-1