坐标轮换的matlab程序,基于matlab的坐标轮换法程序

《基于matlab的坐标轮换法程序》由会员分享，可在线阅读，更多相关《基于matlab的坐标轮换法程序(5页珍藏版)》请在人人文库网上搜索。

1、function y=ff(x1,x2)%y=10*(x1+x2-5)2+(x1-x2)2; y=x12+x22-x1*x2-10*x1-4*x2+60;function alpha=goldmethod2(x01,x02,d,h0)% 输入分割比gamagama=0.618; % 确定搜索区间a,ba,b=search2(x01,x02,d,h0); %计算a1,a2和函数值y1,y2, 确定最优步长 alphaa1=b-gama*(b-a);y1=ff(x01+d(1)*a1,x02+d(2)*a1);a2=a+gama*(b-a);y2=ff(x01+d(1)*a2,x02+d(2)*a。

2、2);for n=1:100if y1=y2x1(n)=a1;x2(n)=a2;yp1(n)=y1;yp2(n)=y2;a=a1;a1=a2;y1=y2;a2=a+gama*(b-a);y2=ff(x01+d(1)*a2,x02+d(2)*a2);elsex1(n)=a1;x2(n)=a2;yp1(n)=y1;yp2(n)=y2;b=a2;a2=a1;y2=y1;a1=b-gama*(b-a);y1=ff(x01+d(1)*a1,x02+d(2)*a1);endaa(n)=(a+b)/2; %输出极值点y(n)=ff(x01+d(1)*aa(n),x02+d(2)*aa(n); %输出极值点函。

3、数值e(n)=abs(b-a);alpha=(a+b)/2;if abs(b-a)=Delta_m/2*(0-F3)2x01=x(n,1);x02=x(n,2);alpha=goldmethod2(x01,x02,d,h0); % 确定最优步长 alphaalphax(n+1,:)=x(n,:)+alpha*d;e(j,:)=;e(n,:)=d;x0=x(n+1,:);else if F2y1 %反向搜索h=-h;a3=a1;y3=y1;a1=a2;y1=y2;a2=a3;y2=y3;enda3=a2+h;y3=ff(x01+d(1)*a3,x02+d(2)*a3);% 继续搜索 %while y30a=a1;b=a3;elsea=a3;b=a1;end% a,b。