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及以前的版本界面差別較大。
-
電腦上怎樣申請QQ號
很多人在打開電腦的時候,都會習慣登錄QQ跟別人聊天,如果沒有QQ賬號的話,可以在電腦上免費申請QQ號,或者之前的QQ號是別人給的,自己不喜歡,也可以重新申請一個新的QQ號,在這裏,跟大家分享一下在電腦上如何免費申請QQ號的方法。電腦上免費申請QQ號的步驟(01)在電腦上下載和...
-
QQ瀏覽器如何安裝應用插件
QQ瀏覽器如何安裝應用插件,瀏覽器的應用插件的安裝對我們來說是非常方便的。下面以安裝護眼導航應用爲例,介紹如何安裝瀏覽器應用插件。操作方法(01)打開QQ瀏覽器,點擊右上角的“”+“”如圖所示(02)進入了安裝應用的首頁,如圖所示(03)比如安裝一個護眼的應用,鼠標停留在...
-
EPSON R230打印機清零軟件無法使用的解決方法
當EPSONR230打印機兩個紅燈閃爍,並提示“打印機的墨管已到它們的使用壽命”時就需要使用清零軟件進行處理,可是在我們使用清零軟件的過程中,總是出現這樣或那樣的問題。那麼我們在遇到這些問題時該如何解決呢?操作方法(01)首先當打開清零軟件可執行文件時,系統會彈出...
-
在Word文檔中怎麼樣進行簡繁體轉換?
字體有繁體和簡體之分,而Word文檔主要的功能是記錄保存文本文字,可能由於區域不一樣,有的喜歡用簡體有的喜歡用繁體,爲了方便查看,就需要用到簡繁體的相互轉換了。操作方法(01)方法一:比較常用的,就是用Word文檔自帶的簡繁體轉換功能,找到審閱——簡轉繁——便可將選擇的...