下载此文档

MATLAB编辑n次拉格朗日函数插值法的程序.doc


文档分类:IT计算机 | 页数:约3页 举报非法文档有奖
1/3
下载提示
  • 1.该资料是网友上传的,本站提供全文预览,预览什么样,下载就什么样。
  • 2.下载该文档所得收入归上传者、原创者。
  • 3.下载的文档,不会出现我们的网址水印。
1/3 下载此文档
文档列表 文档介绍
n次拉格朗日函数插值
程序:
function LagrangesNs() %用于求过n点的拉格朗日n-1次插值多项式
options={'Name of data file'};
title='Lagranges_points';
lineNo=2;
def={''};
outval=inputdlg(options,title,lineNo,def);
if isempty(outval)==1,return,end
filename=outval{1};
data=load(filename);
x=data(:,1);
y=data(:,2);
lagrangesN(x,y);
end
function lagrangesN(x,y)
%画出已知n个点的位置
plot(x,y,'*');
hold on
%n次拉格朗日多项式为 y=a0+a1*x+a2*x^2+…+a(n-1)*x^(n-1)
%其中a0 a1 a2…a(n-1)为待求系数
n=length(x);
X=Vandermonde(x,1);
A=X\y;
%绘制插值函数图象
x1=linspace(0,max(x));
x2=Vandermonde(x1',n);
y1=x2*A;
plot(x1',y1);
hold on
%显示公式
func=['y= ',num2str(A(1))];
for i=2:n;
b=['+ ',num2str(A(i)),'*x^',num2str(i-1)];
func=[func,b];
end
text(,,func);
end
%创建一个Vandermonde行列式

MATLAB编辑n次拉格朗日函数插值法的程序 来自淘豆网www.taodocs.com转载请标明出处.

非法内容举报中心
文档信息
  • 页数3
  • 收藏数0 收藏
  • 顶次数0
  • 上传人笨小孩
  • 文件大小0 KB
  • 时间2011-08-29