[求助][原创]共轭梯度法的matlab源程序
<span style="FONT-SIZE: 10pt; FONT-FAMILY: 宋体; mso-ascii-font-family: 'Times New Roman'; mso-hansi-font-family: 'Times New Roman';"><p class="MsoNormal" style="MARGIN: 0cm 0cm 0pt;"><span style="FONT-SIZE: 10pt; FONT-FAMILY: 宋体; mso-ascii-font-family: 'Times New Roman'; mso-hansi-font-family: 'Times New Roman';">以下是同学写的,好像有错,谁能帮忙修改下呀~~或者帮忙写个<span style="FONT-SIZE: 10pt; FONT-FAMILY: 宋体; mso-ascii-font-family: 'Times New Roman'; mso-hansi-font-family: 'Times New Roman';">共轭梯度法编程实现的程序,c,fortran,matlab...都行~~</span></span></p><p class="MsoNormal" style="MARGIN: 0cm 0cm 0pt;"><span style="FONT-SIZE: 10pt; FONT-FAMILY: 宋体; mso-ascii-font-family: 'Times New Roman'; mso-hansi-font-family: 'Times New Roman';">共轭梯度法的</span><span lang="EN-US" style="FONT-SIZE: 10pt;"><font face="Times New Roman">matlab</font></span><span style="FONT-SIZE: 10pt; FONT-FAMILY: 宋体; mso-ascii-font-family: 'Times New Roman'; mso-hansi-font-family: 'Times New Roman';">源程序</span><span lang="EN-US" style="FONT-SIZE: 12pt;"><p></p></span></p><p class="MsoNormal" style="MARGIN: 0cm 0cm 0pt;"><span lang="EN-US" style="FONT-SIZE: 12pt;"><font face="Times New Roman">function g=conjugate_grad_2d(x0,t) <br/>%please input this:conjugate_grad_2d(,0.05) <br/>x=x0; <br/>syms xi yi a <br/>f=xi^2-xi*yi+3*yi^2; <br/>fx=diff(f,xi); <br/>fy=diff(f,yi); <br/>fx=subs(fx,{xi,yi},x0); <br/>fy=subs(fy,{xi,yi},x0); <br/>fi=; <br/>count=0; <br/>while double(sqrt(fx^2+fy^2))>t <br/>s=-fi; <br/>if count<=0 <br/>s=-fi; <br/>else <br/>s=s1; <br/>end <br/>x=x+a*s; <br/>f=subs(f,{xi,yi},x); <br/>f1=diff(f); <br/>f1=solve(f1); <br/>if f1~=0 <br/>ai=double(f1); <br/>else <br/>break <br/>x,f=subs(f,{xi,yi},x),count <br/>end <br/>x=subs(x,a,ai); <br/>f=xi^2-xi*yi+3*yi^2; <br/>fxi=diff(f,xi); <br/>fyi=diff(f,yi); <br/>fxi=subs(fxi,{xi,yi},x); <br/>fyi=subs(fyi,{xi,yi},x); <br/>fii=; <br/>d=(fxi^2+fyi^2)/(fx^2+fy^2); <br/>s1=-fii+d*s; <br/>count=count+1; <br/>fx=fxi; <br/>fy=fyi; <br/>end <br/>x,f=subs(f,{xi,yi},x),count<p></p></font></span></p></span> 想看看 共轭梯度法 求解偏微分方程
页:
[1]