如何用Matlab寫拉格朗日函數(shù)?
網(wǎng)友解答: 謝邀。首先拉格朗日函數(shù)具體公式如下:編寫一個名為lagrange.m的M文件,然后設(shè)n個節(jié)點數(shù)據(jù)以數(shù)組x0, y0輸入(注意Matlab的數(shù)組下標從1開始),m個插值點以數(shù)組
謝邀。首先拉格朗日函數(shù)具體公式如下:
編寫一個名為lagrange.m的M文件,然后設(shè)n個節(jié)點數(shù)據(jù)以數(shù)組x0, y0輸入(注意Matlab的數(shù)組下標從1開始),m個插值點以數(shù)組x輸入,輸出數(shù)組y 為m個插值。
圖片內(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)用編寫的程序,并運行。在Matlab的命令窗口輸入【lagrange (x,y,xh)】按【Enter】鍵即可得到拉格朗日插值函數(shù)計算的插值。
如果你對科學和科技內(nèi)容感興趣,歡迎訂閱我的酷米號。我會在這里發(fā)布所有與科技、科學有關(guān)的有趣文章。偶爾也回答有趣的問題,有問題可隨時在評論區(qū)回復(fù)和討論,看到即回。
(碼字不易,若文章對你幫助可點贊支持~)