Matlab程式設計實現拉格朗日插值法計算插值
在數值分析中,拉格朗日插值法是以法國十八世紀數學家約瑟夫·拉格朗日命名的一種多項式插值方法。許多實際問題中都用函式來表示某種內在聯絡或規律,而不少函式都只能通過實驗和觀測來了解。如對實踐中的某個物理量進行觀測,在若干個不同的地方得到相應的觀測值,拉格朗日插值法可以找到一個多項式,其恰好在各個觀測的點取到觀測到的值。
操作方法
(01)理解多項式插值的內涵,明確拉格朗日插值多項式的性質和計算方法。
(02)開啟Matlab軟體。
(03)新建M檔案。1:在Matlab中點選左上角的New加號按鈕。2:在彈出的下拉選單中選擇Script。(Ps:可直接使用快捷鍵【Ctrl+N)】軟體就會彈出M函式編輯介面。)
(04)編寫拉格朗日插值函式程式:function yh=lagrange (x,y,xh)n = length(x);m = length(xh);yh = zeros(1,m);c1 = ones(n-1,1);c2 = ones(1,m);for i=1:nxp = x([1:i-1 i+1:n]);yh = yh + y(i)*prod((c1*xh-xp'*c2)./(x(i)-xp'*c2));end注:該程式只可一次計算實現一個插值計算。可實現多個插值計算的程式如下:function yh=lagrange(x,y,xh)n = length(x);m = length(xh);x = x(:);y = y(:);xh = xh(:);yh = zeros(m,1);c1 = ones(1,n-1);c2 = ones(m,1);for i=1:n,xp = x([1:i-1 i+1:n]);yh = yh + y(i) * prod((xh*c1-c2*xp')./(c2*(x(i)*c1-xp')),2);end
(05)儲存檔案,需要特別注意的是檔名應與函式名相同。
(06)在Matlab中輸入給定資料x和y,以及需要計算插值點的值xh。如圖輸入後按【Enter】鍵,載入資料。
(07)呼叫編寫的程式,並執行。在Matlab的命令視窗輸入【lagrange (x,y,xh)】按【Enter】鍵即可得到拉格朗日插值函式計算的插值。
特別提示
本文筆者所用的Matkab版本為2014a,其他版本均適用,只是2012a及以前的版本介面差別較大。
-
如何利用ps去水印
當你在網路上看到一張圖片的時候很想儲存下來的時候,卻無奈的發現這張圖片有水印,通常我們使用的方法就是截圖,用裁剪的方式去掉有水印的地方,但是如此就會破壞整張圖的佈局美感,下面就教你如何用ps的填充功能快速去除水印。操作方法(01)相信大家對於如何安裝ps已經很...
-
怎麼在Excel列A資料中篩選出列B資料?
已知excel列A資料和列B資料有部分內容相同,怎樣才能篩選出哪些混在列A中的列B資料呢?最笨的方法就是複製列B的資料查詢。這裡,我給大家分享一種方法,既準確又快捷的方法。不信,您嘗試一下。操作方法(01)開啟需要編輯的Excel檔案,觀察列A和列B,列A列B有部分內容重複,現需...
-
win10任務計劃在哪裡 如何開啟計劃任務程式
在使用Windows10的過程中,有時候我們需要禁用任務計劃程式裡的某個計劃,這時就需要開啟任務計劃程式,那麼我們怎麼樣開啟計劃任務程式呢?下面就與大家分享兩個方法吧。方法/步驟1(01)在Windows10系統桌面,依次點選“開始/Windows系統/控制面板”選單項。(02)在開啟的Win...
-
如何將圖片轉換成word文字
在工作生活中,我們經常為圖片文字識別的問題而煩惱,要知道圖片中的文字是不可編輯的,我們如果想將圖片中的文字提取出來,就只能使用識別工具對圖片進行識別,將其轉換成Word文字。可是,我們要選擇哪一款識別工具呢?大家如何將圖片轉換成Word文字呢?下面我們就一起來看看...