如何用Matlab寫拉格朗日函數(shù)?
網(wǎng)友解答: 謝邀。首先拉格朗日函數(shù)具體公式如下:編寫一個(gè)名為lagrange.m的M文件,然后設(shè)n個(gè)節(jié)點(diǎn)數(shù)據(jù)以數(shù)組x0, y0輸入(注意Matlab的數(shù)組下標(biāo)從1開始),m個(gè)插值點(diǎn)以數(shù)組
謝邀。首先拉格朗日函數(shù)具體公式如下:
編寫一個(gè)名為lagrange.m的M文件,然后設(shè)n個(gè)節(jié)點(diǎn)數(shù)據(jù)以數(shù)組x0, y0輸入(注意Matlab的數(shù)組下標(biāo)從1開始),m個(gè)插值點(diǎn)以數(shù)組x輸入,輸出數(shù)組y 為m個(gè)插值。
圖片內(nèi)容如下:
純文本內(nèi)容如下(可直接復(fù)制使用):
function y=lagrange(x0,y0,x)
n=length(x0);m=length(x);
for i=1:m
z=x(i);
s=0.0;
for k=1:n
p=1.0;
for j=1:n
if j~=k
p=p*(z-x0(j))/(x0(k)-x0(j));
end
end
s=p*y0(k)+s;
end
y(i)=s;
end
保存后調(diào)用編寫的程序,并運(yùn)行。在Matlab的命令窗口輸入【lagrange (x,y,xh)】按【Enter】鍵即可得到拉格朗日插值函數(shù)計(jì)算的插值。
如果你對(duì)科學(xué)和科技內(nèi)容感興趣,歡迎訂閱我的酷米號(hào)。我會(huì)在這里發(fā)布所有與科技、科學(xué)有關(guān)的有趣文章。偶爾也回答有趣的問題,有問題可隨時(shí)在評(píng)論區(qū)回復(fù)和討論,看到即回。
(碼字不易,若文章對(duì)你幫助可點(diǎn)贊支持~)