




版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡介
1、精品文檔精品文檔西京學(xué)院數(shù)學(xué)軟件實(shí)驗(yàn)任務(wù)書課程名稱數(shù)學(xué)軟件實(shí)驗(yàn)班級數(shù) 0901學(xué)號0912020119姓名王震實(shí)驗(yàn)課題非線性方程的弦割法,拋物線法實(shí)驗(yàn)?zāi)康氖煜し蔷€性方程的弦割法,拋物線法實(shí)驗(yàn)要求運(yùn)用 Matlab/C/C+/Java/Maple/Mathematica 等具中一種語言完成實(shí)驗(yàn)內(nèi)容非線性方程的弦割法非線性方程的拋物線法成績教師【實(shí)驗(yàn)課題】非線性方程的弦割法,拋物線法【實(shí)驗(yàn)?zāi)康摹渴煜し蔷€性方程的弦割法,拋物線法【實(shí)驗(yàn)內(nèi)容】1、弦割法弦割法是函數(shù)逼近法的一種,基本思想是用用區(qū)間Xk山Xk上的割線近似代替目標(biāo)函數(shù)的導(dǎo)函數(shù)的曲線。并用割線與橫軸交點(diǎn)的橫坐標(biāo)作為方程式的根的近似。對于牛頓迭
2、代式f XkxXkq k= 0,1,2,f xk由于上式存在微分,為了避免計(jì)算導(dǎo)數(shù)值,我們用差商來代替導(dǎo)數(shù)。即:Xk 1 = Xkf(Xk)f (Xk) - f(Xk4)(Xk -Xk)算法過程:1、定義函數(shù)f(X),輸入控制精度總和迭代初始值Xo,Xi ,并求出 f (Xo), f (Xi)2、進(jìn)行循環(huán)迭代,若相鄰的兩個(gè)數(shù)值之差的絕對值小于所給定的則跳出循環(huán),輸出Xk2、拋物線法拋物線法是在函數(shù)的最優(yōu)點(diǎn)附近取三個(gè)構(gòu)造點(diǎn),然后用這三個(gè)點(diǎn)構(gòu)造一條拋物線,把這條拋物線的極值點(diǎn)作為函數(shù)的極 值點(diǎn)的近似。每次構(gòu)造一條拋物線后,拋物線的極值點(diǎn)就可作 為一個(gè)新的構(gòu)造點(diǎn),新的構(gòu)造點(diǎn)與原來的三個(gè)構(gòu)造點(diǎn)經(jīng)過某種
3、 算法,得到下一步拋物線逼近的三個(gè)構(gòu)造點(diǎn),這就是拋物線法 的算法過程。算法過程:選定初始值Xo,Xi,X2,并計(jì)算f (%), f (Xi), f(X2)和差分f(X2”) JX2)- f(Xl)X2 - Xif(Xi,X。)J(Xl)-f(X0)Xi -Xof(X2,Xi) - f(Xi,Xo)f(X2,Xi,Xo) =-LX2 - Xo注意,Xo, Xi, X2三點(diǎn)不能共線。2、用牛頓插值法對三點(diǎn)(X0,f(X0)、(x1, f(Xi) (x2, f(X2)進(jìn)行插值得到一條拋物線,存在兩個(gè)根:-B _、B2 -4ACX3 2C其中,A=f(X2), C = f(X2, x, Xo)B =
4、f (X2,Xi) f (X2,Xi,Xo)(X2 -Xi)兩個(gè)根中取其中靠近X2的那個(gè)根。3、重復(fù)以上步驟,并控制精度。精品文檔精品文檔【程序】%function x1,n=gexian()x0=input(請輸入初始值x0=);x1=input(請輸入初始值x1=);err=input( 請輸入誤差極限err=);N=input( 請輸入最大迭代次數(shù)N=);x2=x1-fun(x1)*(x1-x0)/(fun(x1)-fun(x0);n=0;while abs(x1-x0)errx0=x1;x1=x2;x2=x1-fun(x1)*(x1-x0)/(fun(x1)-fun(x0);n=n+1
5、;if nNdisp( 迭代次數(shù)過大!);returnendendend%function root=paowuxian()x0=input(請輸入初值x0=);x1=input(請輸入初值x1=);x2=input(請輸入初值x2=);N=input( 請輸入最大迭代次數(shù)N=);f1=fun(x0);f2=fun(x1);if f1=0root=x0;endif f2=0root=x1;endif f1*f20disp( 兩點(diǎn)之積大于0! );returnelseerr=input( 請輸入誤差極限err=);f1=fun(x0);f2=fun(x1);f3=fun(x2);d1=(f2-f
6、1)/(x1-x0);d2=(f3-f2)/(x2-x1);d3=(d2-d1)/(x2-x0);B=d2+d3*(x2-x1);root=x2-2*f3/(B+sign(B)*sqrt(BA24f3*d3);t(1)=x0;t(2)=x1;t(3)=x2;n=0;tol=abs(root-t(3);while tolerrn=n+1;if nNdisp(迭代次數(shù)過大!);returnendt(1)=t(2);t(2)=t(3);t(3)=root;f1=fun(t(1);f2=fun(t(2);f3=fun(t(3);d1=(f2-f1)/(t(2)-t(1);d2=(f3-f2)/(t(3)-t(2);d3=(d2-d1)/(
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲(chǔ)空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 2025至2030年中國老鼠浴沙數(shù)據(jù)監(jiān)測研究報(bào)告
- 2025至2031年中國醬鹵香辛料包行業(yè)投資前景及策略咨詢研究報(bào)告
- 2025至2031年中國電動(dòng)對開多葉調(diào)節(jié)器行業(yè)投資前景及策略咨詢研究報(bào)告
- 2025至2031年中國油壓鉆床行業(yè)投資前景及策略咨詢研究報(bào)告
- 2025至2031年中國可拆卸水射器行業(yè)投資前景及策略咨詢研究報(bào)告
- 2025至2031年中國全粘膠無放布行業(yè)投資前景及策略咨詢研究報(bào)告
- 2025至2031年中國中段排氣管行業(yè)投資前景及策略咨詢研究報(bào)告
- 2025至2031年中國CO2激光標(biāo)記機(jī)行業(yè)投資前景及策略咨詢研究報(bào)告
- 科技環(huán)保趨勢下的可持續(xù)發(fā)展
- 網(wǎng)絡(luò)改造施工合同范本
- 2025年黑龍江農(nóng)業(yè)工程職業(yè)學(xué)院單招職業(yè)適應(yīng)性測試題庫及答案1套
- 《勞動(dòng)法常識(第3版)》中職全套教學(xué)課件
- 2025年勞動(dòng)合同延期補(bǔ)充協(xié)議模板
- 2025年日歷表(含農(nóng)歷、節(jié)假日、記事、A4打印版)
- 《反家庭暴力》課件
- 二零二五年度房地產(chǎn)預(yù)售合同協(xié)議4篇
- 2025-2030年中國天線行業(yè)市場需求狀況規(guī)劃研究報(bào)告
- 2024年南京旅游職業(yè)學(xué)院高職單招職業(yè)技能測驗(yàn)歷年參考題庫(頻考版)含答案解析
- 如何提升自我管理能力
- 2025年潛江市城市建設(shè)發(fā)展集團(tuán)招聘工作人員【52人】高頻重點(diǎn)提升(共500題)附帶答案詳解
- 人教版(新)九年級下冊化學(xué)全冊教案教學(xué)設(shè)計(jì)及教學(xué)反思
評論
0/150
提交評論