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瀏覽器如何安裝應用插件,瀏覽器的應用插件的安裝對我們來說是非常方便的。下面以安裝護眼導航應用爲例,介紹如何安裝瀏覽器應用插件。操作方法(01)開啟QQ瀏覽器,點擊右上角的“”+“”如圖所示(02)進入了安裝應用的首頁,如圖所示(03)比如安裝一個護眼的應用,鼠標停留在...
-
單片機案例點亮一顆發光二極管
大家好,今天和大家分享的是單片機案例——點亮一顆led燈。這個可以說是每個人在初步接觸單片機是都會學到的一個程序,我們來看看怎麼做吧操作方法(01)這裏我就直接進入程序了,前面的建立工程呀,設定參數那些我就不說了,在我其他的經驗中也有講到了,不懂的可以看看我以...
-
怎樣製作電子版的賀卡?
在這個互聯網如此發達的時代,每逢親朋好友過生日,各種節日的到來之際,我們總想借助互聯網送去自己的祝福。往往想給親朋好友郵寄一張賀卡,然而網絡搜尋到的賀卡有時不盡人意,或沒有創意,或沒法表達自己的感情。那麼,爲什麼不嘗試親手製作一張賀卡呢?步驟/方法(01)設定賀...
-
電腦工具欄到了電腦的右邊,怎麼恢復到底部?
一般電腦的工具欄是在電腦的底部,可是有些時候就會跑到電腦兩側,今天我要講的就是恢復方法。操作方法(01)工具欄跑到了電腦的右側(02)在工具欄處點擊鼠標右鍵,找到“屬性”點擊,出現任務欄對話框(03)在鎖定任務欄前打√,任務欄在屏幕上的位置選擇底部,點擊確定即可...