第04章.zip

  • 烟花太阳
    了解作者
  • matlab
    开发工具
  • 7KB
    文件大小
  • zip
    文件格式
  • 0
    收藏次数
  • 1 积分
    下载积分
  • 1
    下载次数
  • 2021-04-22 15:03
    上传日期
MATLAB2019素材,共同学习,共同进步
第04章.zip
  • 4.2
  • banana.m
    124B
  • qiujizhi.m
    276B
  • F.m
    127B
  • fx3.m
    169B
  • fx1.m
    176B
  • lagrangechazhi.m
    1.2KB
  • eulerff2.m
    113B
  • fsolvefun.m
    152B
  • dfx3.m
    185B
  • dfx2.m
    169B
  • threepoints.m
    344B
  • dfx1.m
    186B
  • eulerff.m
    107B
  • jacobi.m
    255B
  • SOR.m
    268B
  • staticiterate.m
    182B
  • diffqiudao.m
    270B
  • euler.m
    275B
  • newtoniterate.m
    275B
  • fx.m
    138B
  • broydeniterate.m
    360B
  • fx2.m
    156B
  • gauss.m
    252B
  • twopoint.m
    298B
  • eulerff3.m
    111B
  • ff1.m
    112B
内容介绍
% Lagrange插值 clear clc %----------------------------- n=10; %结点个数 lb=-1; %下界 ub=1; %上界 step=0.01; %作图点步长 %----------------------------- % 原始函数图形 x0=lb:step:ub; y0=1./(1+25*x0.^2); plot(x0,y0,'r-'); hold on %----------------------------- % 插值函数 for i=1:n+1 xi(i)=lb+(ub-lb)*(i-1)/n; yi(i)=1/(1+25*xi(i)^2); end %------------------------------ count=1; for x=lb:step:ub fl=0; %-------------------------- %求出pn(xk) for k=1:n+1 up=1; dn=1; %---------------------- %求出f(xk) for i=1:n+1 if k~=i up=up*(x-xi(i)); dn=dn*(xi(k)-xi(i)); end end %---------------------- fl=fl+yi(k)*up/dn; end pn(count)=fl; %-------------------------- fi(count)=1/(1+25*x^2);%求原函数的值 count=count+1; end %------------------------------ % L插值函数图 x=lb:step:ub; plot(x,pn,'g--') %------------------------------ num=(ub-lb)/step+1; for i=1:num p_f(i)=pn(i)-fi(i); end center=fix(num/2); scale=fix(num/10); a=center-scale; b=center+scale; disp ' pn(i)-fi(i) 的值为:' p_f(a:b)
评论
    相关推荐