




版權(quán)說(shuō)明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
1、庶爾牘n理2厚賣 驗(yàn) 報(bào) 告課程名稱數(shù)值分析實(shí)驗(yàn)項(xiàng)目名稱插值實(shí)驗(yàn)實(shí)驗(yàn)類型上機(jī)實(shí)驗(yàn)學(xué)吋2班級(jí)20111131學(xué)號(hào)2011113130姓名張振指導(dǎo)教師沈艷實(shí)驗(yàn)室名稱理學(xué)樓407實(shí)驗(yàn)時(shí)間2013.10.11實(shí)驗(yàn)成績(jī)預(yù)習(xí)部分實(shí)驗(yàn)過(guò)程 表現(xiàn)實(shí)驗(yàn)報(bào)告 部分總成績(jī)教師簽字日期哈爾濱工程大學(xué)教務(wù)處制實(shí)驗(yàn)一插值實(shí)驗(yàn)一.插值法的基本思想1 .lagrange插值多項(xiàng)式的基本思想及公式:lagrange插值多項(xiàng)式的基本思想是把多項(xiàng)式的構(gòu)造問(wèn)題轉(zhuǎn)化成n+1個(gè)基函 數(shù).(x)的構(gòu)造。其公式如下:l”(x)二外仏.其中 ljx)=(7)(f)(f m)用0,1,n.(耳_兀0)(忑("一忑+j(耳-斗)2.ne
2、wton插值多項(xiàng)式基本思想及公式:newton插值多項(xiàng)式是lagrange插值多項(xiàng)式的改進(jìn),是由基函數(shù)1, (xx°),,(xxo).(xx“_)逐次地推得到的多項(xiàng)式。其公式如下:p“(x) = ao+a】(xxo)+ art (x-x0).(x-xw_1). 其中 = fxo,x,xr ,k=0,l,n.3.hermite插值的基木思想及公式:hermite插值不僅滿足各插值點(diǎn)上函數(shù)值相等,而且各節(jié)點(diǎn)上導(dǎo)數(shù)值也想等, 甚至高階導(dǎo)數(shù)值也相等。在這里只給出兩種典型的hermite插值: p(x) = f(x()+fx0,x)(x-x()+ fx(),xl,x2(x-x0)(x-xi)+
3、a(xx0)(x-x!)(x-x2) 其屮系數(shù) a二 f'(xjfxo,x|fxo,x|,x2(x|xo)八 八(xj -x 0)(x, -x2 )余項(xiàng) r(x)=lf(4)( © ) (xx()(xx j2(xx2).4!其中(1 + 2x-xi *xoj-州丿/、/ 、21 + 2 兀一 x。i ®x)l坷-“0丿h°(x)h (兀)=/ 、2【x- xxho(兀)=(兀一兀o)衛(wèi)0 一州丿h (%) = hq(x)yq + /?j (x)y +(兀)幾+h(x)y、2w1(x) = (x-x1)余項(xiàng)r3(x) = f(x)-h3(x)=(x-x0)2
4、(x-x()2 兵(a,b)4 三次樣條函數(shù)s(x)wl a,b,且在每個(gè)小區(qū)間抹沖鬥上是三次多項(xiàng)式,其中a =x0<x1<.< xzj = b是給定節(jié)點(diǎn),則稱s(x)是節(jié)點(diǎn)x()x.x“上的三次樣條函數(shù)。若在節(jié)點(diǎn)上給定函數(shù)yj二 f(x7.).(j=o, 1,.,n),并成立 s(xy. )=y . .(j= 0, 1,., n),則稱 s(x)為三次樣條插值函 數(shù)。二實(shí)驗(yàn)題目及實(shí)驗(yàn)?zāi)康?. (第2章習(xí)題20)給定數(shù)據(jù)表如下:x/0.250.300.390.450.53y;0.50000.54770.62450.67080.7280試求三次樣條插值s (x),并滿足條件:(
5、1) s (0.25) =1.0000, s (0.53) =0.6868;(2) s" (0.25) =s (0.53) =0.2. (第2章計(jì)算實(shí)習(xí)題3)下列數(shù)據(jù)點(diǎn)的插值x01491625364964y012345678可以得到平方根函數(shù)的近似,在區(qū)間0,64上作圖。(1)用這9個(gè)點(diǎn)作8次多項(xiàng)式插值l* (x).(2)用三次樣條(第一邊界條件)程序求s (x).從得到結(jié)果看在0,64上,哪個(gè)插值更精確;在區(qū)間0,1上,兩種插值哪個(gè)更精 確?三.實(shí)驗(yàn)手段:指操作環(huán)境和平臺(tái):'¥ir)7系統(tǒng)b* matlab r2009a程序語(yǔ)言:一種類似c語(yǔ)言的程序語(yǔ)言,但比c語(yǔ)
6、言要寬松得多,非常方便。四程序1.三次樣條插值程序:function y = scyt(x0, yo, x)n = length(xo); m = length(x);y = zeros(1# m);p(n) = 0; q(l) = 1; d(l) = zl.oooo; d(n) = 0.6868;for k = 2:n-1,h(1) = xo(k) - xo(k-1) ; h(2) = xo(k+1) - xo(k);p(k) = h(l) / (h(l) + h(2);q(k) = h(2) / (h(l) + h(2);d(k) = 6* ( (yo (k+l) - yo (k) ) /
7、h(2) - (yo (k) - yo (k-1) ) / h(l) ) / (h(l) + h (2); endb(l) = q(l) / 2;for k = 2 zn-l,b(k) = q(k) / (2 - p(k)*b(k-1);endtemp(1) = d(1) / 2;for k = 2 zn-l,temp (k) = (d (k.) p (k) *temp (k-1) )/ (2 - p (k) *b (k-1);endm(l) = 0;m(n) = 0;for k =1:-1:2,m(k) = temp(k) - b(k)*m(k+l); endfor i = 1:m#for
8、k = 1:n-l,if (x(i) >= xo(k) && (x(i) <= xo(k+1),h = xo(k+1)- xo(k);y(i)= m(k)*(xo(k+1)- x(i)3 / (6*h)x(i) )/h)-xo(k)/h);+ m (k+1)*(x(i)- xo(k)a3 / (6*h)+ (yo(k) m(k)*h“2 / 6)*(xo(k+1) + (yo(k+1) - m(k+1)*h2 / 6)*(x(i) endendend 控制臺(tái)輸入代碼:x0= 0.25 0.30 0.39 0.45 0.53;y0= 0.5000 0.5477 0.6
9、245 0.6708 07280; x= 0.25:0.02:0.53;y=scyt(xo z yo z x)2. 拉格朗日插值程序:function y=lagrange(xo,yo,x) n=length(xo);m=length(x); for i=l:mz=x(i);s=00;for k=l:np=1.0;for j =1:nif j =kp=p*(z-xo(j)/(xo(k)-xo(j );endends=p*yo(k)+s;endy (i)=s;end三次樣條插值程序:function y = scyt(xo, yo, x)n = length(xo); m = length(x)
10、;y = zeros(1, m);p (n) = 0 ; q (1) = 0 ; d (1) = 0 ; d (n) =0 ;for k = 2h(1) = xo(k)- xo(k-1); h(2) = xo(k+1)- xo(k);p(k) = h(l) / (h(l) + h(2);q(k) = h(2) / (h(l) + h(2);d(k) = 6*(y0(k+l) - y0(k)/h(2) - (yo(k) yo(kl) / h(l) ) / (h(l) + h (2);endb(l) = q(l) / 2;for k = 2:n-1,b(k) = q(k) / (2 p(k)*b(
11、k-l);end temp(1) = d(1) / 2 ;for k = 2:n,temp(k) = (d(k)- p(k)*temp(k-1)/ (2 - p(k)*b(k-1);endb(k)*m(k+1);m(n) = temp(n); for k = n-1:-1:1,m (k) = temp(k) end for i = 1:m,for k = 1:n-l,if (x(i)>= xo(k) && (x(i) <= xo(k+1),h = xo(k+1)- xo(k);y(i) = m(k)*(xo(k+1)- x(i)a3 / (6*h)+ m(k+1)*
12、(x(i)- xo(k)人3/(6*h)+ (yo(k)- m(k)*hx2 / 6)*(xo(k+1)- x(i)/h)+ (yo (k+1) puk+lphy / 6 ) * ( (x (i)- xo (k) ) /h);endendend 控制臺(tái)輸入代碼:x0= 0 1 4 9 16 25 36 49 64;y0= 0 1234567 8;x= 0:4 : 64;yl=sqrt(xo)plot(xo,yl)hold ony2 = lagrange(xo z yo z x)plot(x,y2,1r1)hold ony3= scyt(xo, yo, x)plot(x,y3,'g*)五.實(shí)驗(yàn)結(jié)果比較與分析1.圖像:分析與結(jié)果:從圖像上看,用三次樣條插值得到的擬合曲線比較完美,無(wú)“龍格”現(xiàn)彖出現(xiàn)。2.圖像比較:0,1區(qū)間上圖像放大結(jié)果
溫馨提示
- 1. 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 2025服裝店鋪面房屋租賃合同
- 定制工廠出品合同范本
- 預(yù)防醫(yī)學(xué)(安徽中醫(yī)藥大學(xué))知到課后答案智慧樹(shù)章節(jié)測(cè)試答案2025年春安徽中醫(yī)藥大學(xué)
- 伐木機(jī)械租賃合同范本
- 2025標(biāo)準(zhǔn)委托設(shè)計(jì)合同
- 2025商業(yè)綜合體中央空調(diào)系統(tǒng)投資合作合同
- 2024年四川阿壩州人民醫(yī)院招聘緊缺衛(wèi)生專業(yè)技術(shù)人員真題
- 新版?zhèn)€人租房合同范本
- 2024年佳木斯市郊區(qū)招聘公益性崗位人員真題
- 2024年北京協(xié)和醫(yī)院后勤保障處宿舍管理人員招聘筆試真題
- 《檔案編研工作》課件
- 《山水林田湖草生態(tài)保護(hù)修復(fù)工程指南(試行)》
- 初中英語(yǔ)牛津深圳版單詞表(按單元順序)七年級(jí)至九年級(jí)
- 槍支安全及使用指南
- 《肝衰竭診治指南(2024版)》解讀
- 國(guó)省道公路標(biāo)志標(biāo)線維護(hù)方案投標(biāo)文件(技術(shù)方案)
- 【MOOC】科技英語(yǔ)寫作-西安電子科技大學(xué) 中國(guó)大學(xué)慕課MOOC答案
- 電動(dòng)汽車課件
- 原始點(diǎn)醫(yī)學(xué)(201904第15版)
- 網(wǎng)絡(luò)安全應(yīng)急處置工作預(yù)案
- 住宅物業(yè)消防安全管理 XF1283-2015知識(shí)培訓(xùn)
評(píng)論
0/150
提交評(píng)論