matlab实验-拉格朗日插值的龙格(Runge)现象

本文通过MATLAB代码详细介绍了如何使用拉格朗日插值法进行函数逼近,展示了从选择插值点到计算插值多项式,再到与原函数对比的全过程。适合数学专业学生和对数值分析感兴趣的读者。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

记一次数学作业
学会用拉格朗日插值法

n=15;  %5,10,15插值点数
syms t;

x = linspace(-5,5,n);
y=f(x);
for i=1:n
L(i)=Lagrange(x,n,i,t);
end
LN=sum(y.*L);
t=-5:0.01:5;
y0=5./(1+t.^2);
 plot(t,y0);
 hold on;
 ezplot(LN,t);
 a=strcat('实验一:',num2str(n),'次插值');
 title(a);
 legend('原函数','插值函数');

function fi=Lagrange(x,n,i,t)

fu=1;
fd=1;

for  j=1:n                      %循环计算 
         if i~=j
             fu=fu.*(t-x(j));
             fd=fd.*(x(i)-x(j));
         end
end
fi=fu/fd;
end
  
function y=f(x)
y=5./(1+x.^2);
end  
评论 8
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值