共軛梯度法--MATLAB程序.docx_第1頁(yè)
共軛梯度法--MATLAB程序.docx_第2頁(yè)
共軛梯度法--MATLAB程序.docx_第3頁(yè)
共軛梯度法--MATLAB程序.docx_第4頁(yè)
免費(fèi)預(yù)覽已結(jié)束，剩余1頁(yè)可下載查看

下載本文檔

版權(quán)說明：本文檔由用戶提供并上傳，收益歸屬內(nèi)容提供方，若內(nèi)容存在侵權(quán)，請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)

文檔簡(jiǎn)介

1、%conjugate gradient methods%method:FR,PRP,HS,DY,CD,WYL,LS%精確線搜索 , 梯度終止準(zhǔn)則function m,k,d,a,X,g1,fv = conjgradme( G,b,c,X,e,method)ifnargin=ek(i-1)=(m/m1)2;d(:,i)=-g1(:,i)+k(i-1)*d(:,i-1);a(i)=-(d(:,i)*g1(:,i)/(d(:,i)*G*d(:,i);%a1(i)=-(X(:,i)*G*d(:,i)+b*d(:,i)/(d(:,i)*G*d(:, i);a(i)=g1(:,i)*g1(:,i)/(d(

2、:,i)*G*d(:,i);X(:,i+1)=X(:,i)+a(i)*d(:,i);g1=g1 subs(subs(g,x1,X(1,i+1),x2,X(2,i+1);m1=m;m=norm(g1(:,i+1);i=i+1;endcasePRPwhilem=ek(i-1)=g1(:,i)*(g1(:,i)-g1(:,i-1)/(norm(g1(:,i-1)2;d(:,i)=-g1(:,i)+k(i-1)*d(:,i-1);a(i)=-(d(:,i)*g1(:,i)/(d(:,i)*G*d(:,i);X(:,i+1)=X(:,i)+a(i)*d(:,i);g1=g1 subs(subs(g,x1

3、,X(1,i+1),x2,X(2,i+1);m=norm(g1(:,i+1);i=i+1;endcaseHSwhilem=ek(i-1)=g1(:,i)*(g1(:,i)-g1(:,i-1)/(d(:,i-1)*(g1(:,i)-g1(:,i-1);d(:,i)=-g1(:,i)+k(i-1)*d(:,i-1);a(i)=-(d(:,i)*g1(:,i)/(d(:,i)*G*d(:,i);X(:,i+1)=X(:,i)+a(i)*d(:,i);g1=g1 subs(subs(g,x1,X(1,i+1),x2,X(2,i+1);m=norm(g1(:,i+1);i=i+1;endcaseDYwh

4、ilem=ek(i-1)=g1(:,i)*g1(:,i)/(d(:,i-1)*(g1(:,i)-g1(:,i-1);d(:,i)=-g1(:,i)+k(i-1)*d(:,i-1);a(i)=-(d(:,i)*g1(:,i)/(d(:,i)*G*d(:,i);X(:,i+1)=X(:,i)+a(i)*d(:,i);g1=g1 subs(subs(g,x1,X(1,i+1),x2,X(2,i+1);m=norm(g1(:,i+1);i=i+1;endcaseLSwhilem=ek(i-1)=g1(:,i)*(g1(:,i)-g1(:,i-1)/(d(:,i-1)*(-g1(:,i-1);d(:,i

5、)=-g1(:,i)+k(i-1)*d(:,i-1);a(i)=-(d(:,i)*g1(:,i)/(d(:,i)*G*d(:,i); %a(i)=-(X(:,i)*G*d( :,i)+b*d(:,i)/(d(:,i)*G*d(:,i);X(:,i+1)=X(:,i)+a(i)*d(:,i);g1=g1 subs(subs(g,x1,X(1,i+1),x2,X(2,i+1);m=norm(g1(:,i+1);i=i+1;endcaseCDwhilem=ek(i-1)=g1(:,i)*g1(:,i)/(d(:,i-1)*(-g1(:,i-1);d(:,i)=-g1(:,i)+k(i-1)*d(:,

6、i-1);a(i)=-(d(:,i)*g1(:,i)/(d(:,i)*G*d(:,i);X(:,i+1)=X(:,i)+a(i)*d(:,i);g1=g1 subs(subs(g,x1,X(1,i+1),x2,X(2,i+1);m=norm(g1(:,i+1);i=i+1;endcaseWYLwhilem=ek(i-1)=g1(:,i)*(g1(:,i)-(m/m1)*g1(:,i-1)/(m12);d(:,i)=-g1(:,i)+k(i-1)*d(:,i-1);a(i)=-(d(:,i)*g1(:,i)/(d(:,i)*G*d(:,i);%a(i)=-(X(:,i)*G*d(:,i)+b*d

7、(:,i)/(d(:,i)*G*d(:,i);X(:,i+1)=X(:,i)+a(i)*d(:,i);g1=g1 subs(subs(g,x1,X(1,i+1),x2,X(2,i+1);m1=m;m=norm(g1(:,i+1);i=i+1;endendfv=subs(subs(f,x1,X(1,i),x2,X(2,i);endl1=X(1,i);l2=X(2,i);w1=X(1,1);w2=X(2,1);v1=min(l1,w1)-abs(l1-w1)/10:abs(l1-w1)/10:max(l1,w1)+abs(l1-w1)/10;v2=min(l2,w2)-abs(l2-w2)/10:abs(l2-w2)/10:max(l2,w2)+abs(l2-w2)/10;x,y=meshgrid(v1,v2);s=size(x);z=zeros(size(x);fori=1:s(1)forj=1:s(2)z(i,j)=1/2*x(i,j),y(i,j)*G*x(i,j);y(i,j)+b*x(i,j);y(i,j)+c;endendpx,py = gradient(z,.2,.2);contour(v1,v2,z), holdon , qui

溫馨提示

  • 1. 本站所有資源如無特殊說明，都需要本地電腦安裝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ù)覽，若沒有圖紙預(yù)覽就沒有圖紙。
  • 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ì)自己和他人造成任何形式的傷害或損失。

最新文檔

評(píng)論

0/150

提交評(píng)論