並行插值演算法
1. 插值的計算方法是什麼
計算方法:假設與A1對應的數據是B1,與A2對應的數據是B2,現在已知與A對應的數據是B,A介於A1和A2之間,則可以按照(A1-A)/(A1-A2)=(B1-B)/(B1-B2)計算得出A的數值,其中A1、A2、B1、B2、B都是已知數據。
根據(A1-A)/(A1-A2)=(B1-B)/(B1-B2)可知:(A1-A)=(B1-B)/(B1-B2)×(A1-A2)
A=A1-(B1-B)/(B1-B2)×(A1-A2)=A1+(B1-B)/(B1-B2)×(A2-A1)
插值法又稱「內插法」,是利用函數f (x)在某區間中已知的若干點的函數值,作出適當的特定函數,在區間的其他點上用這特定函數的值作為函數f (x)的近似值,這種方法稱為插值法。如果這特定函數是多項式,就稱它為插值多項式。
如果只需要求出某一個x所對應的函數值,可以用「圖解內插」。它利用實驗數據提供要畫的簡單曲線的形狀,然後調整它,使得盡量靠近這些點。
如果還要求出因變數p(x)的表達式,這就要用「表格內插」。通常把近似函數p(x)取為多項式(p(x)稱為插值多項式),最簡單的是取p(x)為一次式,即線性插值法。在表格內插時,使用差分法或待定系數法(此時可以利用拉格朗日公式)。在數學、天文學中,插值法都有廣泛的應用。
2. 請列一下插值法的計算公式,並舉個例子。
舉個例子。
2008年1月1日甲公司購入乙公司當日發行的面值600 000元、期限3年、票面利率8%、每年年末付息且到期還本的債券作為可供出售金融資產核算,實際支付的購買價款為620 000元。
則甲公司2008年12月31日因該可供出售金融資產應確認的投資收益是()元。(已知PVA(7%,3)=2.2463,PVA(6%,3)=2.673,PV(7%,3)=0.8163,PV(6%,3)=0.8396)
題目未給出實際利率,需要先計算出實際利率。600 000×PV(r,3)+600 000×8%×PVA(r,3)=620 000,採用內插法計算,得出r=6.35%。甲公司2008年12月31日因該可供出售金融資產應確認的投資收益=620 000×6.35%=39 370(元)。
插值法計算過程如下:
已知PVA(7%,3)=2.2463,PVA(6%,3)=2.673,PV(7%,3)=0.8163,PV(6%,3)=0.8396)
600 000×PV(r,3)+600 000×8%×PVA(r,3)=620 000
R=6%時
600000*0.8396+600000*8%*2.673=503760+128304=632064
R=7%時
600000*0.8163+600000*8%*2.2463=489780+107823=597603
6% 632064
r 620000
7% 597603
(6%-7%)/(6%-R)=(632064-597603)/(632064-620000)
解得R=6.35%
注意上面的式子的數字順序可以變的,但一定要對應。如可以為
(R-7%)/(7%-6%)=(620000-597603)/(597603-632064)也是可以的,當然還有其他的順序。"
(2)並行插值演算法擴展閱讀:
若函數f(x)在自變數x一些離散值所對應的函數值為已知,則可以作一個適當的特定函數p(x),使得p(x)在這些離散值所取的函數值,就是f(x)的已知值。從而可以用p(x)來估計f(x)在這些離散值之間的自變數所對應的函數值,這種方法稱為插值法。
如果只需要求出某一個x所對應的函數值,可以用「圖解內插」。它利用實驗數據提供要畫的簡單曲線的形狀,然後調整它,使得盡量靠近這些點。
如果還要求出因變數p(x)的表達式,這就要用「表格內插」。通常把近似函數p(x)取為多項式(p(x)稱為插值多項式),最簡單的是取p(x)為一次式,即線性插值法。
在表格內插時,使用差分法或待定系數法(此時可以利用拉格朗日公式)。在數學、天文學中,插值法都有廣泛的應用。
3. 插值法計算公式是什麼
公式就是:Y=Y1+(Y2-Y1)×(X-X1)/(X2-X1)。
通俗地講,線性內插法就是利用相似三角形的原理,來計算內插點的數據。
內插法又稱插值法。根據未知函數f(x)在某區間內若干點的函數值,作出在該若干點的函數值與f(x)值相等的特定函數來近似原函數f(x),進而可用此特定函數算出該區間內其他各點的原函數f(x)的近似值,這種方法,稱為內插法。
按特定函數的性質分,有線性內插、非線性內插等;按引數(自變數)個數分,有單內插、雙內插和三內插等。
介紹:
線性插值是指插值函數為一次多項式的插值方式,其在插值節點上的插值誤差為零。線性插值相比其他插值方式,如拋物線插值,具有簡單、方便的特點。
線性插值的幾何意義即為概述圖中利用過A點和B點的直線來近似表示原函數。線性插值可以用來近似代替原函數,也可以用來計算得到查表過程中表中沒有的數值。
4. 數字圖像處理——圖像插值
網上有很多介紹插值演算法的,但感覺收獲都不大
介紹三種圖像插值演算法:最近鄰內插,雙線性內插,雙三次內插(雙立方內插)
三次插值即用三階多項式擬合原函數(也應該有其他用途)。假設三次擬合函數為
在matlab中,圖像被定義為一個三維向量,若不考慮圖像的通道數,可以將圖像看作一個二維矩陣處理。matlab圖像矩陣中坐標值映射到二維坐標系中,每一個像素塊對應的是一個點,但實際的像素塊是有一定尺寸的。
在進行雙線性插值和雙三次插值時,需要用坐標值擬合函數,為了簡化計算,總是選取 作為局部坐標系原點,其中 為待插值坐標。
當出現這些情況時,補充這些像素的灰度值為圖像內最相鄰像素塊的灰度值。
進行坐標變換後,選取與內插點 歐式距離最近的像素值進行插值。在程序中,使用將 按照四捨五入的舍入方式選取最近鄰的像素塊。
雙線性內插是線性內插的二維實現,在x維度先進行線性插值,再由得到的值對y維度進行插值。在局部坐標系中,選取 相鄰的四個像素進行雙線性內插。由在數學原理中的推導可知
雙三次內插是三次插值的二維實現。選取與 相鄰的16個像素進行雙三次內插,局部坐標系中x與y坐標范圍均為 。由數學原理中的推到可知
最近鄰插值法的優點是計算量很小,運算速度較快。但它僅使用離待測采樣點最近的像素的灰度值作為該采樣點的灰度值,而沒考慮其他相鄰像素點的影響,因而重新采樣後灰度值有明顯的不連續性,會產生明顯的馬賽克和鋸齒現象。
雙線性插值法效果要好於最近鄰插值,計算量較大。縮放後圖像質量高,基本克服了最近鄰插值灰度值不連續的特點,因為它考慮了待測采樣點周圍四個直接鄰點對該采樣點的相關性影響。但是,此方法未考慮到各鄰點間灰度值變化率的影響, 具有低通濾波器的性質, 從而導致縮放後圖像的高頻分量受到損失, 圖像邊緣在一定程度上變得較為模糊,丟失了一些細節信息。
雙立方插值計算量最大,運算速度慢。雙立方插值用三階函數逼近,不僅考慮到周圍四個直接相鄰像素點灰度值的影響,還考慮到它們灰度值變化率的影響,能夠產生比雙線性插值更為平滑的邊緣,計算精度很高,處理後的圖像細節損失最少,效果最佳。
5. 線性插值法如何計算
線性插值是數學、計算機圖形學等領域廣泛使用的一種簡單插值方法。 常用計算方法如下:假設我們已知坐標 (x0,y0)與 (x1,y1),要得到 [x0,x1]區間內某一位置x在直線上的值。 我們可以得到(y-y0) (x-x0)/ (y1-y0) (x1-x0) 假設方程兩邊的值為α,那麼這個值就是插值系數—從x0到x的距離與從x0到x1距離的比值。 由於x值已知,所以可以從公式得到α的值 α= (x-x0)/ (x1-x0) 同樣,α= (y-y0)/ (y1-y0) 這樣,在代數上就可以表示成為: y = (1- α)y0 + αy1 或者, y = y0 + α (y1 - y0) 這樣通過α就可以直接得到 y。
6. 線性插值法計算公式是什麼
線性插值法計算公式:Y=Y1+(Y2-Y1)×(X-X1)/(X2-X1)。其中Y2>Y1,X2>X>X1。線性插值是指插值函數為一次多項式的插值方式,其在插值節點上的插值誤差為零。線性插值相比其他插值方式,如拋物線插值,具有簡單、方便的特點。線性插值可以用來近似代替原函數,也可以用來計算得到查表過程中表中沒有的數值。
線性插值使用的原因
目前,線性插值演算法使用比較廣泛。在很多場合我們都可以使用線性插值。其中,最具代表性的使用方法是變數之間的對應關系沒有明確的對應關系,無法使用公式來描述兩個變數之間的對應關系,在這種情況下使用線性插值是比較好的解決辦法。可以在變數的變化區間上取若干個離散的點,以及對應的輸出值,然後將對應關系分成若干段,當計算某個輸入對應的輸出時,可以進行分段線性插值。
7. 什麼是插值演算法
插值法又稱「內插法」,是利用函數f (x)在某區間中插入若干點的函數值,作出適當的特定函數,在這些點上取已知值,在區間的其他點上用這特定函數的值作為函數f (x)的近似值,這種方法稱為插值法。如果這特定函數是多項式,就稱它為插值多項式。
1、Lagrange插值:
Lagrange插值是n次多項式插值,其成功地用構造插值基函數的 方法解決了求n次多項式插值函數問題;
★基本思想將待求的n次多項式插值函數pn(x)改寫成另一種表示方式,再利 用插值條件⑴確定其中的待定函數,從而求出插值多項式。
2、Newton插值:
Newton插值也是n次多項式插值,它提出另一種構造插值多項式的方法,與Lagrange插值相比,具有承襲性和易於變動節點的特點;
★基本思想將待求的n次插值多項式Pn(x)改寫為具有承襲性的形式,然後利用插值條件⑴確定Pn(x)的待定系數,以求出所要的插值函數。
3、Hermite插值:
Hermite插值是利用未知函數f(x)在插值節點上的函數值及導數值來構造插值多項式的,其提法為:給定n+1個互異的節點x0,x1,……,xn上的函數值和導數值
求一個2n+1次多項式H2n+1(x)滿足插值條件
H2n+1(xk)=yk
H'2n+1(xk)=y'k k=0,1,2,……,n ⒀
如上求出的H2n+1(x)稱為2n+1次Hermite插值函數,它與被插函數
一般有更好的密合度;
★基本思想
利用Lagrange插值函數的構造方法,先設定函數形式,再利
用插值條件⒀求出插值函數.
4、分段插值:
插值多項式余項公式說明插值節點越多,誤差越小,函數逐近越好,但後來人們發現,事實並非如此,例如:取被插函數,在[-5,5]上的n+1個等距節點:計算出f(xk)後得到Lagrange插值多項式Ln(x),考慮[-5,5]上的一點x=5-5/n,分別取n=2,6,10,14,18計算f(x),Ln(x)及對應的誤差Rn(x),得下表
從表中可知,隨節點個數n的增加,誤差lRn(x)l不但沒減小,反而不斷的增大.這個例子最早是由Runge研究,後來人們把這種節點加密但誤差增大的現象稱為Runge現象.出現Runge現象的原因主要是當節點n較大時,對應
的是高次插值多項式,此差得積累"淹沒"了增加節點減少的精度.Runge現象否定了用高次插值公式提高逼近精度的想法,本節的分段插值就是克服Runge現象引入的一種插值方法.
分段多項式插值的定義為
定義2: a=x0<x1<…<xn=b: 取[a,b]上n+1個節點 並給定在這些節點 上的函數值f(xR)=yR R=0,1,…,n
如果函數Φ(x)滿足條件
i) Φ(x)在[a,b]上連續
ii) Φ(xr)=yR,R =0,1,…,n
iii) Φ(x)zai 每個小區間[xR,xR+1]是m次多項式,
R=0,1,…,n-1則稱Φ(x)為f(x)在[a,b]上的分段m次插值多項式
實用中,常用次數不超過5的底次分段插值多項式,本節只介紹分段線性插值和分段三次Hermite插值,其中分段三次Hermite插值還額外要求分段插值函數Φ(x)
在節點上與被插值函數f(x)有相同的導數值,即
★基本思想將被插值函數f〔x〕的插值節點 由小到大 排序,然後每對相鄰的兩個節點為端點的區間上用m 次多項式去近似f〔x〕.
例題
例1 已知f(x)=ln(x)的函數表為:
試用線性插值和拋物線插值分別計算f(3.27)的近似值並估計相應的誤差。
解:線性插值需要兩個節點,內插比外插好因為3.27 (3.2,3.3),故選x0=3.2,x1=3.3,由n=1的lagrange插值公式,有
所以有,為保證內插對拋物線插值,選取三個節點為x0=3.2,x1=3.3,x2=3.4,由n=2的lagrange插值公式有
故有
所以線性插值計算ln3.27的誤差估計為
故拋物線插值計算ln3.27的誤差估計為:
顯然拋物線插值比線性插值精確;
5、樣條插值:
樣條插值是一種改進的分段插值。
定義 若函數在區間〖a,b〗上給定節點a=x0<x1<;…<xn=b及其函數值yj,若函數S(x)滿足
⒈ S(xj)=yj,j=0,1,2,…,n;
插值法主要用於道路橋梁,機械設計,電子信息工程等 很多工科領域的優化方法。
8. 一行在插值演算法方面有哪些貢獻
今天常用的牛頓插值公式,其不等間距的形式比等間距的形式要復雜得多。天算史界有一種流行的看法,認為在中國古代,唐朝天文學家、數學家一行在其《大衍歷》中發明了二次不等間距插值法,且一行還有意識地應用了三次差內插法近似公式。因此,一行在插值法方面的貢獻備受中外天算史研究者的關注。
中國古代非線性插值法,是劉焯在其《皇極歷》(604年)中考慮到太陽運動不均勻性為計算太陽行度改正值時首創的。有關中國古代插值法的算理研究的新成果表明,劉焯二次等間距插值法的造術原理建立在源於《九章算術》描述勻變速運動的模型基礎之上,認為太陽每日的運行速度之值構成一等差數列。換言之,所用數學方法就是構造一等差數列並求其前若干項和。
一行的插值法並沒有人們所想像那樣的推廣意義。就插值演算法本身,一行演算法與劉焯演算法實質完全相同。所不同的是,《皇極歷》是在以平氣為間隔的日躔表基礎上插值。而《大衍歷》是在以定氣為間隔的日躔表上插值。
《太初歷》以後,各歷都以平分一回歸年365.25日為24等份而得每節氣長15.22日,這樣規定的二十四氣稱為「常氣」,或叫「平氣」。張子信指出「日行春分後則遲,秋分後則速」,於是劉焯造《皇極歷》時體會到二十四氣皆應有「定日」,他說:「春、秋分定日去冬至各八十八日有奇,去夏至各九十三日有奇。」但劉焯並沒有搞清楚太陽速度的加減和季節的關系,他的日躔表是把秋分定日後到春分定日前平均分為12段,每氣14.54日;春分定日後到秋分定日前也平分為12段,每氣15.45日。這顯然不是「定氣」。
一行認為,太陽在一回歸年365.2444日中共行365.2444度,每氣行15.2185度。冬至附近日行速度最急,故二氣間所需運行時間最短,夏至附近日行速度最緩,故二氣間的時間最長。
實際上,《大衍歷》這里首先提出了平分黃道為24等份,以太陽實際走完每個等份的時間長度為各節氣長度,這就是通常所稱的「定氣」概念。一行提出正確的定氣概念以後,在計算太陽改正時自然就以定氣為插值間隔。至於插值法本身則完全是沿用劉焯的方法。
值得一提的是,劉焯在日躔表中規定太陽視運動一年內的變化規律是:冬至最快,冬至後漸慢,到立春時開始加快,春分時又達到最快,冬至到春分這段時間內日速比平均速度快。春分後太陽視運動的速度突變為最慢,之後逐漸加快,到立夏時又開始減慢,夏至達到最慢。春分到夏至時段內比平均速度慢。夏至以後的變化情況以夏至處為鏡面對稱。
《大衍歷》盈縮分一年內的變化趨勢將盈縮分在冬至附近最大,以後逐漸變小,夏至時最小,之後又逐漸增大。這相當於把冬至作為太陽視運動的近日點,夏至為遠日點。這種認識是正確的,而《皇極歷》的規定是不符合實際的。
說一行有意識地應用了三次差內插法的近似公式,是指《大衍歷》的月亮極黃緯演算法和五星中心差改正演算法中所用的插值法。當對中國古代歷法中的插值法的構造原理有了深入的認識之後,研究者進一步通過將這兩處插值法的有關術文與劉焯二次等間距插值法的術文進行對比研究,證明兩者在實質上也是相同的。
人們之所以會認為《大衍歷》使用了三次差插值法,是因為《大衍歷》在上述兩種演算法的插值法中都引入了「中差」概念的緣故。
但實際上一行引入「中差」的原因在於,劉焯日躔表中的各氣陟降率之差是相等的,而《大衍歷》月亮極黃緯等數表相鄰兩欄的差一般不等。這種現象的出現,正是一行受命改歷時作了大量天文觀測的結果。若仍照搬《皇極歷》的做法,就會出現同一點處有可能得到兩個不同的值的現象,這就迫使一行必須在計算方法上進行一點細節上的調整。
一行作為科學家,在中國科技史上具有重要的地位,作為佛教高僧,一行傳承胎藏和金剛兩大部密法,在密宗史上的作用,不只系統組織密教的教義教規,也把兩大部融合起來。集科學家與高僧於一身這個特殊身份本身,也說明佛法和科學技術在一定條件下的相融性。