



全文預(yù)覽已結(jié)束
下載本文檔
版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進行舉報或認(rèn)領(lǐng)
文檔簡介
應(yīng)用共軛梯度法求解方程組的根。初始值分析:將方程組轉(zhuǎn)化為優(yōu)化問題中的極值問題然后應(yīng)用共軛梯度法進行求解。令,只需求得,使得取得最小值。則有:與比較則有:回到題目問題中,則對應(yīng)有也就是應(yīng)用共軛梯度法求點使取得最小值。程序清單:#include #includedouble a,b,s2; /*全局變量*/double E=1e-6;double F(double x1,double x2)double y;y=2*x1*x1+2*x2*x2-12*x1-4*x2+20;return(y);void qujian(double x1,double x2) /*用前進后退法求a的探索區(qū)間*/double a0=0,h=1,a1,a2,f1,f2,X1,X2,X3,X4;a1=a0;a2=a0+h;X1=x1+a1*s0;X2=x2+a1*s1;f1=F(X1,X2);X3=x1+a2*s0;X4=x2+a2*s1;f2=F(X3,X4);if (f1f2)while(1)h=h*2;a2=a2+h;f1=f2;X3=x1+a2*s0;X4=x2+a2*s1;f2=F(X3,X4);if (f1f2) a1=a2-h;elsea=a1;b=a2;break;elseh=-h/4;while(1)a1=a1+h;f2=f1;X1=x1+a1*s0;X2=x2+a1*s1;f1=F(X1,X2);if(f1f2) a2=a1-h;h=2*h;elsea=a1;b=a2;break;double MIN(double x1,double x2) /*二次插值法求a的最小值*/double x01,x02,x03,x0,xmin;double f1,f2,f3,f;double X1,X2,X3,X4,X5,X6,X7,X8;double h=(b-a)/2;double k1,k2;x01=a;x02=a+h;x03=b;X1=x1+x01*s0;X2=x2+x01*s1;X3=x1+x02*s0;X4=x2+x02*s1;X5=x1+x03*s0;X6=x2+x03*s1;f1=F(X1,X2);f2=F(X3,X4);f3=F(X5,X6);while (1)k1=(f3-f1)/(x03-x01);k2=(f2-f1)/(x02-x01)-k1)/(x02-x03);x0=(x01+x03-k1/k2)/2;X7=x1+x0*s0;X8=x2+x0*s1;f=F(X7,X8);if(fabs(x0-x02)=E)xmin=x0;break;else if(x02x0)if(f2f)x03=x0;f3=f;else x01=x02;f1=f2;x02=x0;f2=f;elseif(f2f)x01=x0;f1=f;elsex03=x02;f3=f2;x02=x0;f2=f;return(xmin);void main()double g3,m1,m2,b;int n=2,k=0;double x1=0,x2=0; /*起始點*/g1=4*x1-12;g2=4*x2-4;m1=g1*g1+g2*g2;s0=-g1;s1=-g2;while(1)qujian(x1,x2);a=MIN(x1,x2); /*求a的探索區(qū)間*/x1=x1+a*s0;x2=x2+a*s1;g1=4*x1-12;m2=g1*g1+g2*g2;if(sqrt(g1*g1+g2*g2)=E)break;if(k+1=n)g1=4*x1-12;g2=4*x2-4;s0=-g1;s1=-g2;else b=m2/m1;s0=-g1+b*s0;s1=-g2
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 檔案行政管理辦法規(guī)定
- 地理知識梳理與綜合能力提升策略
- 北京護送車輛管理辦法
- 村民務(wù)工補貼管理辦法
- 因果復(fù)句的歷史演變與語言學(xué)分析
- 廢舊農(nóng)膜回收與處置制度困境與完善路徑探究
- 公共住房資產(chǎn)管理辦法
- 決策咨詢工作管理辦法
- 銀行金融產(chǎn)品的精準(zhǔn)營銷策略
- 內(nèi)部孵化項目管理辦法
- 鍋爐澆注料施工方案
- GB/T 17394.1-2014金屬材料里氏硬度試驗第1部分:試驗方法
- GB/T 1606-2008工業(yè)碳酸氫鈉
- 葛的栽培技術(shù)
- 《綠色建筑概論》整套教學(xué)課件
- 山東中醫(yī)藥大學(xué)2020-2021學(xué)年內(nèi)科護理學(xué)試題及答案2
- 2022年綿陽江油市社區(qū)工作者招聘考試模擬試題及答案解析
- 初中道德與法治學(xué)科教學(xué)經(jīng)驗交流
- 工程測量、定位放線控制點復(fù)核記錄表
- 申辦出入境證件的函
- 安全評估收費指導(dǎo)意見
評論
0/150
提交評論