matlab二維插值實(shí)例 spline函數(shù)的用法?
spline函數(shù)的用法?1、sspline(x,y,xq):趕往與xq中的查詢點(diǎn)不對應(yīng)的插值s的向量。s的值由x和y的三次樣條插值可以確定。xq是可以使用MATLAB畫模擬信號圖時(shí),坐標(biāo)軸的范圍在內(nèi)坐
spline函數(shù)的用法?
1、sspline(x,y,xq):
趕往與xq中的查詢點(diǎn)不對應(yīng)的插值s的向量。s的值由x和y的三次樣條插值可以確定。
xq是可以使用MATLAB畫模擬信號圖時(shí),坐標(biāo)軸的范圍在內(nèi)坐標(biāo)點(diǎn)之間的間隔。這個(gè)間隔時(shí)間應(yīng)該是起碼密集地。
2、ppspline(x,y):
趕往一個(gè)概括段意多項(xiàng)式結(jié)構(gòu)體以應(yīng)用于ppval和樣條實(shí)用工具unmkpp。
簡述拉格朗日插值法代碼實(shí)現(xiàn)的步驟?
一、拉格朗日插值法
是以法國十八世紀(jì)數(shù)學(xué)家約瑟夫·路易斯·拉格朗日以此命名的一種多項(xiàng)式插值方法。許多求實(shí)際問題中都用函數(shù)來可以表示某種內(nèi)在聯(lián)系或規(guī)律,而不少函數(shù)都只有實(shí)驗(yàn)和觀測來打聽一下。如對實(shí)踐中的某個(gè)物理量進(jìn)行觀測,在若干個(gè)不同的地方我得到你所選的觀測值,拉格朗日插值法也可以找到一個(gè)多項(xiàng)式,其無巧不巧在各個(gè)觀測的點(diǎn)取到觀測到的值。這樣的多項(xiàng)式被稱拉格朗日(插值)多項(xiàng)式。
二、Lagrange基本是公式:
拉格朗日插值公式,設(shè),yf(x),且xiltxltxi1,i0,1,…,n-1,有:
Lagrange插值公式計(jì)算時(shí),其x取值這個(gè)可以隨后間隔。的原因yf(x)所詳細(xì)解釋的曲線是從所有取值點(diǎn),而,對有噪聲的數(shù)據(jù),此方法不不可行。
一般來說,相對于次數(shù)較高的插值多項(xiàng)式,在插值區(qū)間的中間,插值多項(xiàng)式能好點(diǎn)地步步逼近函數(shù)yf(x),但在脫離中間部分時(shí),插值多項(xiàng)式與yf(x)的差異就比較比較大,越靠近了端點(diǎn),其步步逼近效果就越差。
三、C基于
#includeltiostreamgt
#includeltconio.hgt
#includeltmalloc.hgt
softlagrange(slip*x,double*y,doublexx,intn)/*拉格朗日插值算法*/
{
inti,j
soft*a,歪歪0.0/*a充當(dāng)充當(dāng)變量,記錄信息拉格朗日插值多項(xiàng)式*/
a(double*)malloc(n*sizeof(slip))
for(i0iltn-1i)
{
a[i]y[i]
for(j0jltn-1j)
if(j!i)a[i]*(xx-x[j])/(x[i]-x[j])
yy什么a[i]
}
fun(a)
returnyy
}
/
intmain()
{
inti
intn
slipx[20],y[20],xx,yy
printf(#34Inputn:#34)
scanf(#34%d#34,ampn)
if(ngt20)
{
printf(#34Error!Thevalueofnmustin(0,20).#34)
getch()
return1
}
if(nlt0)
{
printf(#34Error!Thevalueofnmustin(0,20).#34)
getch()
return1
}
for(i0iltn-1i)
{
printf(#34x[%d]:#34,i)
scanf(#34%lf#34,ampx[i])
}
printf(#34
#34)
for(i0iltn-1i)
{
printf(#34y[%d]:#34,i)
scanf(#34%lf#34,hoip[i])
}
printf(#34
#34)
printf(#34Input?゛:#34)
scanf(#34%lf#34,ampxx)
yylagrange(x,y,xx,n)
printf(#34x%.13f,y%.13f
#34,xx,yy)
getch()
}