matlab中拉格朗日插值的多項式和係數怎麼求?急急急

時間 2021-08-11 17:29:59

1樓:匿名使用者

function yy=lagrange(x1,y1,xx)%本程式為lagrange1插值,其中x1,y1%為插值節點和節點上的函式值,輸出為插值點xx的函式值,%xx可以是向量。

syms x

n=length(x1);

for i=1:n

t=x1;t(i)=;l(i)=prod((x-t)./(x1(i)-t));% l向量用來存放插值基函式

endu=sum(l.*y1);

p=simplify(u) % p是簡化後的lagrange插值函式(字串)

yy=subs(p,x,xx);

clfplot(x1,y1,'ro',xx,yy,'*')

2樓:匿名使用者

假設已知n個插值資料(x1,y1),(x2,y2),(x3,y3),……(xn,yn),則拉格朗日插值多項式的形式為:

代入x值即可求出相應的y值

3樓:匿名使用者

多看看拉格朗日插值公式,多看幾遍,觀察下係數的關係

4樓:軒轅汐玥

%lagr 返回所求拉格朗日多項式的係數

% x為每個點的橫座標,y為對應的縱座標

% 輸入n個點,求n-1次拉格朗日多項式

n=length(x);

%輸入點的個數

syms x;

%引數x

sum=0;

for i=1:n

%n次求和

p=y(i);

%子項初始值為y的值

for j=1:n

if j~=i

p=p*((x-x(j))/(x(i)-x(j)));

%連乘 1 2 3 4 5 6

endend

sum=p+sum;

%求和end% sum=polyfit(x,sum,n);

result=coeffs(sum,x);

%提取多項式的係數

result=vpa(result);

%將係數化為小數end

matlab求拉格朗日插值多項式並畫圖

5樓:匿名使用者

你需要把上面的部分儲存為lagrangenew。m檔案。

檔案頭應加  function s=lagrangenew(x,y,x0)

後面是nx=length(x);

ny=length(y);

。。。。。。

然後,可

回以命令視窗執行下答面部分的命令

結果如下

拉格朗日乘數法中可以為零嗎,拉格朗日乘數法系數 可不可以為0

叫俄小博 拉格朗日乘數的數值是按照實際演算獲取的,不排除為0的可能性。根據推導過程可知,是不可以等於0的。如果等於0,f對x求導,就是原函式對x求導 f對y求導,就是原函式對y求導 上面兩個式子一般是不可能解出來的 由拉格朗日乘數法的推導過程可以看出,0,否則駐點 x0,y0 滿足的式子就變成了 f...

怎樣理解流體力學中的拉格朗日描述和尤拉描述

滿腹經綸公子 lagrange描述和euler描述是描述物體運動的兩種方法 拉格朗日法用來描述一個質點的運動,用初始時刻的座標來標記質點,記錄這個質點每時每刻所在的位置。用數學來表達就是r a,b,c,t 這裡a,b,c就是初始時刻質點的座標。拉格朗日描述其實就是理論力學裡的方法。尤拉法描述固定的空...

高等數學利用拉格朗日證明不等式的問題

你好!你理解的非常正確,那個點 或者可能有不止一個 是依存與函式f和區間 a,b 而客觀存在的,如果直接人為指定那個點的值,那是絕對錯誤的!但是我們仍然可以運用拉格朗日中值定理來證明不等式,原因並不在於我們可以指定任意一點c的值,而是在於我們可以找出f c 的範圍,因為c是在區間 a,b 上的,所以...