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

n次拉格朗日函数插值 程序： function LagrangesNs %用于求过n点的拉格朗日n-1次插值多项式 options Name of data file ; title Lagranges_points; lineNo 2; def Lagranges.dat ; 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 0.8,0.8,func ; end %创建一个Vandermonde行列式 function XX Vandermonde x,m %创建x的列向量 if m 1 n length x ; XX zeros n ; for i 1:n XX :,i x.^ i-1 ; end else n length x ; XX zeros n,m ; for i 1:m XX :,i x.^ i-1 ; end end end 运行情况： 按“run”运行时，弹出窗口 注：在Lagranges 1.dat文件中数据为： 0 1 1 0.3679 2 0.1353 点击图框中的“OK”，在“command window”中输出结果为： 将图框中的相关数据更改为：