word版可編輯-機(jī)械優(yōu)化設(shè)計(jì)matlab程序精心整理docxdocx-文庫(kù)吧資料
2025-07-21 12:24本頁(yè)面
【正文】 迭代次數(shù)while norm(x0x02)=esp k=k+1。esp=1e12。x02=[0。d02=[0。 用鮑威爾法在MATLAB的M文件編輯器中編寫(xiě)的M文件，如下：function [x,fx,k]=bowell(A,b,c,x0)%鮑威爾法d01=[1。0],0,[2。運(yùn)行結(jié)果如下：》 [x,fx,k]=bianchidufa([3 1。*A*x1+b39。 break endend x=x1。 s0=a0*d0。 d0=d1。*G0*d1)。 a1=d139。 k=k+1。*H0)/(s039。H0*y0*s039。*H0*y0/(s039。*H0*y0)。*y0)H0*y0*y039。 % H1=H0+s0*s039。 y0=A*a0*d0。 x1=x0+a0*d0。*G0*d0)。 a0=d039。 k=0。 g0=A*x0+b。x1=x0+a0*d0。y0=A*a0*d0。s0=a0*d0。*g0/(d039。d0=H0*g0。H0=eye(2)。g0=A*x0+b。*x+c。1])y = 1x = k = 1 (2) 用變尺度法在MATLAB的M文件編輯器中編寫(xiě)的M文件，如下：function [x,fx,k]=bianchidufa(A,b,c,x0)%用變尺度法求fx=*x39。1 1],[2。*x+c。y=*x39。 k=k+1。 r0=r1。 r2=A*x1+b。*A*p1)。 alpha=(r139。*r0)。while norm(r1)eps beta=(r139。r1=A*x1+b。*A*p0)。a=r039。 k=0。*A*x+b39。if norm(r0)=eps x=x0。.*x0+c。%迭代停機(jī)原則%fx=*x039。*A*X+b39。運(yùn)行結(jié)果如下：[x,fx,k]=zuniniudunfa([1。endx=x0。 G1=subs(G,{x1,x2},{x0(1,1),x0(2,1)})。 x0=x0+a0*p。*g1/(p39。while(norm(a0*p)epson) p=G1\g1。*G1*p)。a0=p39。k=0。%停機(jī)原則g1=subs(df,{x1,x2},{x0(1,1),x0(2,1)})。G=jacobian(df,v)。df=jacobian(f,v)。fx=0。運(yùn)行結(jié)果如下： [x,fx,k]=niudunfa([1。endx=x0。 G1=subs(G,{x1,x2},{x0(1,1),x0(2,1)})。 x0=x0+p。x0=x0+p。k=0。g1=subs(df,{x1,x2},{x0(1,1),x0(2,1)})。G=jacobian(df,v)。df=jacobian(f,v)。fx=0。3. 用牛頓法、阻尼牛頓法及變尺度法求函數(shù)的極小點(diǎn)。函數(shù)值f0=%\n39。區(qū)間中點(diǎn)t=%\n39。函數(shù)值f2=%\n39。試點(diǎn)2坐標(biāo)值t2=%\n39。迭代區(qū)間右端b=%\n39。函數(shù)值f1=%\n39。試點(diǎn)1坐標(biāo)值t1=%\n39。迭代區(qū)間左端a=%\n39。迭代次數(shù)k=% \n39。 f0=f(t)。 end t=*(b+a)。 t2=a+*(ba)。t1=t2。f1=f(t1)。f2=f1。while abs(ba)=epsilon if f1f2 b=t2。f2=f(t2)。f1=f(t1)。epsilon=。a=0。t39。t^(2/3)(t^2+1)^(1/3)39。,f)運(yùn)行結(jié)果如下：迭代計(jì)算k= 7極小點(diǎn)坐標(biāo)t= 2函數(shù)值f=2. 用黃金分割法求函數(shù)的極小值，精度e=。,t)fprintf(1,39。,k)fprin
教學(xué)教案相關(guān)推薦
鄂ICP備17016276號(hào)-1