matlab实验鲍威尔法
鲍威尔法matlab实现
鲍威尔法% 鲍威尔法-计算第2环起始点和搜索方向syms xy1 xy2 fyx m x1 x2 X s1s2 Sf=60-10*x1-4*x2+x1^2+x2^2-x1*x2;disp ' 目标函数:f=60-10*x1-4*x2+x1^2+x2^2-x1*x2.'pretty(f);% 第1环沿e2方向搜索x01=0;x02=0;X0=[x01 x02]; % 第1环初始点x1=5.0000;x2=0;X=[x1 x2]; % 第1环沿e2方向初始点fy1=60-10*x1-4*x2+x1^2+x2^2-x1*x2; % 第1环e2方向初始点函数值e1=[1,0];e2=[0,1]; % 坐标轴单位方向s1=e2(1);s2=e2(2);S=[s1 s2];m=0;% 计算影射点及其函数值xs1=2*xy1-x01;xs2=2*xy2-x02;Xs=[xs1 xs2];fxs=60-10*xs1-4*xs2+xs1^2+xs2^2-xs1*xs2;fx0=60-10*x01-4*x02+x01^2+x02^2-x01*x02; % 第1环初始点函数值disp ' 第1环影射点坐标'disp (Xs)% 判断第1环搜索函数值下降最大方向df01=fx0-fy1; % 第1环沿e1方向初始点函数值-第1环沿e1方向终点函数值df02=fy1-fyx; % 第1环沿e2方向起始点函数值-第1环沿e2方向始点函数值if df01>df02dfm=df01;m=1;elsedfm=df02;m=2;end% 计算POWELL判别式f1=fx0;f2=fyx;f3=fxs;f123=[f1 f2 f3];disp ' 第1环初始点、终点、影射点函数值'disp (f123)fP1=(f1-2*f2+f3)*(f1-f2-dfm)^2;fP2=0.5*dfm*(f1-f3)^2;fP=[fP1 fP2];disp ' 两个POWELL判别式的值'fprintf(1,' POWELL判别式1的值 fP1= %10.4f \n',fP1)fprintf(1,' POWELL判别式1的值 fP2= %10.4f \n',fP2)fprintf(1,' 函数值下降最大方向 m= %2.0f \n',m)if f1>f3 & fP2>fP1sx1=xy1-X0(1); % 同时不满足两个判别式时56428 .产生新方向sx2=xy2-X0(2);sx=[sx1 sx2];if m==1S=[e2,sx]; % 新方向取代上环搜索中函数值下降最大方向m=1elseS=[e1,sx]; % 新方向取代上环搜索中函数值下降最大方向m=2 endelseS=[e1,e2]; % 满足某个判别式时56428 .保留上环搜索方向enddisp ' @@@@ 第2环搜索方向 @@@@'disp (S)% 第1环沿新方向搜索x1=xy1;x2=xy2;X=[x1 x2];s1=xy1-x01;s2=xy2-x02; S=[s1 s2];[xy1,xy2,fyx]=powell(m,x1,x2,X,s1,s2,S);disp '@@@@ 第2环搜索起始点 @@@@'if fyx<fxsx20=[xy1 xy2];disp ' (第1环极小点)'elsex20=[xs1 xs2];disp ' (第1环影射点)'enddisp (x20)% 鲍威尔法-第1环计算% 目标函数中常数项、步长的一次项和二次项系数fxs0=60-10*x1-4*x2+x1^2+x2^2-x1*x2;fxs1=-10*s1-4*s2+2*x1*s1+2*x2*s2-x1*s2-x2*s1;fxs2=s1^2+s2^2-s1*s2;fxs210=[fxs2 fxs1 fxs0];% 计算最优步长af=-fxs1/(2*fxs2);% 计算终点函数值(关于最优步长)fya=fxs0+fxs1*af+fxs2*af^2;% 计算终点坐标值和函数值(关于终点)xy1=x1+s1*af;xy2=x2+s2*af;Xopt=[xy1 xy2];fyx=60-10*xy1-4*xy2+xy1^2+xy2^2-xy1*xy2;if m==0disp ' ******** 计算第1环沿e2方向搜索 ********'elsedisp ' ******** 计算第1环沿新方向搜索 ********'enddisp ' 初始点'disp (X)disp ' 搜索方向'disp (S)disp ' 步长二次项一次项系数常数项'disp (fxs210)disp ' 最优步长'disp (af)if m==0disp ' 终点坐标'disp (Xopt)elsedisp ' 极小点坐标'disp (Xopt)disp ' 极小点函数值(关于最优步长)' disp (fya)disp ' 极小点函数值(关于极小点)' disp (fyx)end。
powell法matlab
powell法matlab
Powell方法是一种用于无约束优化问题的数值优化算法。
它是由Michael J.D. Powell于1964年提出的，是一种直接搜索方法，不需要计算目标函数的梯度。
在MATLAB中，可以使用内置的fminunc函数来实现Powell方法进行优化。
首先，你需要定义一个目标函数，这个函数是你想要优化的目标，比如最小化或最大化的函数。
然后，你可以使用fminunc函数来调用Powell方法进行优化。
fminunc函数的基本语法如下：
matlab.
[x,fval,exitflag,output] = fminunc(fun,x0,options)。
其中，fun是你定义的目标函数，x0是优化的初始点，options 是优化选项。
在fun中，你需要输入目标函数的表达式，并确保它能够接受输入x，并返回一个标量作为目标函数值。
在使用Powell方法时，你需要特别注意初始点的选择，因为初始点的选择可能会影响最终的优化结果。
另外，你也可以通过调整
options来设置一些优化参数，比如迭代次数、容许误差等。
除了使用MATLAB内置的fminunc函数，你还可以自己实现Powell方法的算法，这需要一定的数值计算和优化算法的知识。
你可以参考相关的优化算法书籍或者论文来了解Powell方法的具体实现细节。
总之，Powell方法是一种常用的无约束优化算法，在MATLAB 中可以通过fminunc函数来实现。
希望这些信息对你有所帮助，如果你有其他关于Powell方法或MATLAB优化的问题，也欢迎继续提问。
基于MATLAB的鲍威尔法求极值问题
基于MATLAB的鲍威尔法求极值问题：xxx 学号：xxx（理工大学机械与车辆学院车辆工程，100081）摘要：无约束优化方法主要有七种，按照求导与否把这些方法分为间接法和直接法。
牛顿法的成败与初始点选择有极大关系，其可靠性最差；坐标轮换法、单纯形法和最速下降法对于高维优化问题计算效率很低，有效性差；由于编制变尺度法程序复杂，其简便性不足。
综合考虑后，鲍威尔法、共轭梯度法具有较好的综合性能。
本文首先对鲍威尔法的原理进行阐述，根据其迭代过程给出流程图，并编写MATLAB程序。
最后用此MATLAB程序求解实际的极值问题，并对求解结果进行简要分析。
1.鲍威尔法的基本思想1.1其他优化方法对鲍威尔法形成的影响通过对鲍威尔法的学习，可以很明显看出来其迭代思想中汲取了其他几种优化方法的核心思想。
为了更全面、更深入的学习鲍威尔法，很有必要对其他有影响的优化思想进行学习和梳理。
由最基本的数学基础知识可知，梯度方向是函数增加最快的方向，负梯度方向是函数下降最快的方向，于是，利用这个下降最快方向产生了最速下降法。
每次迭代都沿着负梯度方向进行一维搜索，直到满足精度要求为止。
其特点是相邻两个搜索方向互相正交，所以很明显的一个现象就是刚开始搜索步长比较大，愈靠近极值点其步长愈小，收敛速度愈慢，特别当二维二次目标函数的等值线是较扁的椭圆时，迭代速度更慢。
这时，倘若目标函数是等值线长、短轴都平行于坐标轴的椭圆形，则通过坐标轮换法可以很高效的解决问题。
通过两次分别沿坐标轴进行一维搜索，便可达到极值点。
但对于目标函数的等值线椭圆的长、短轴倾斜于坐标轴时，坐标轮换法的搜索效率也显得极低。
抛开这两种特殊情况，对于一般形态的目标函数，如果在某些明显可以直达最优点的情况下（一般为靠近极值点区域），迭代过程完全可以不沿负梯度方向搜索，取而代之的是找到直达最优点的方向，一步到位。
但这样的直达方向应该如何去找呢？共轭梯度法由此产生。
其基本原理是：任意形式的目标函数在极值点附近的特性都近似一个二次函数，其等值线在极值点附近为近似的同心椭圆簇，而同心椭圆簇有一个特性便是任意两条平行线与椭圆簇切点的连线必通过椭圆的中心。
鲍威尔算法
一阶导数
变尺度法
d k = − Ak g k
A =A
k
k −1
+ ∆A
k −1
一阶导数，使 A k ≈ [G k ]−1
(d j )T ( g k − g k +1 ) = (d j )T G ( x k +1 − x k ) = 0
取 d k = x k +1 − x k 这说只要沿dj 方向分别对函作两次一维搜索， 得到两个极小点xk 和xk+1 ，那么这两点的连线 所给出的方向dk就是与dj一起对G共轭的方向。
k x0k , xnk , xn+1（ xnk+1 = xnk + ( xnk − x0k ) = 2 xnk − x0k ）
分别称为一轮迭代的始点、终点和反射点 。
f i = f ( xik ) (i = 0,1, 2,L, n) 记：
∆ i = f i −1 − f i (i = 1, 2,L, n)
因此 F0 = f 0 , F2 = f n
则在循环中函数下降最多的第m次迭代是
∆ m = max ∆ i = f m−1 − f m
1≤i ≤ n
k 相应的方向为d m 。
为了构成共轭性好的方向组，须遵循下列准则： 在k次循环中，若满足条件：
F3 < F0
和 ( F0 − 2 F2 + F3 )( F0 − F2 − ∆ m ) 2 < 0.5∆ m ( F0 − F3 ) 2
1 d3 方向进行一维搜索得 沿
4 x = 2
2
f ( x 2 ) = −8
检验终止条件
x 2 − x02 = (4 − 3.8) 2 + (2 − 1.7) 2 = 0.36 > ε
鲍威尔算法matlab程序-f
function f=fun(x)f=10*(x(1)+x(2)-5)^2+(x(1)-x(2))^2; function f=fx(x0,alpha,s)x1=x0+alpha*s;f=fun(x1);function f=fsearch(x0,s)%利用进退法确定高低高区间alpha1=0;h=0.1;alpha2=alpha1+h;f1=fx(x0,alpha1,s);f2=fx(x0,alpha2,s);if f1>f2alpha3=alpha2+h;f3=fx(x0,alpha3,s);while f2>f3alpha1=alpha2;alpha2=alpha3;alpha3=alpha3+h;f2=f3;f3=fx(x0,alpha3,s);endelseh=-h;v=alpha1;alpha1=alpha2; alpha2=v;v=f1;f1=f2;f2=v;alpha3=alpha2+h;f3=fx(x0,alpha3,s); while f2>f3alpha1=alpha2; alpha2=alpha3; alpha3=alpha3+h;f2=f3;f3=fx(x0,alpha3,s); endenda=min(alpha1,alpha3); b=max(alpha1,alpha3); %利用黄金分割点法求解alpha1=a+0.382*(b-a); alpha2=a+0.618*(b-a);f1=fx(x0,alpha1,s);f2=fx(x0,alpha2,s); while abs(a-b)>0.001 if f1>f2a=alpha1;alpha1=alpha2;f1=f2;alpha2=a+0.618*(b-a); f2=fx(x0,alpha2,s); elseb=alpha2;alpha2=alpha1;f2=f1;alpha1=a+0.382*(b-a); f1=fx(x0,alpha1,s); endendf=0.5*(a+b);clear%初始点x0=[0;0];%搜索方向e2=[0;1];G0=fun(x0);F0=G0;%第一次迭代%沿着e1alpha1=fsearch(x0,e1);x1=x0+alpha1*e1;F1=fun(x1);delta1=F0-F1;% 沿着方向e2；alpha2=fsearch(x1,e2);x2=x1+alpha2*e2;F2=fun(x2);G2=F2;delta2=F1-F2;deltam=max(delta1,delta2);%映射点x3=2*x2-x0;G3=fun(x3);if G3<G0 & (G0-2*G2+G3)*(G0-G2-deltam)^2<0.5*deltam*(G0-G3) ^2%方向替换e1=e2;e2=s;% 沿着方向s 进行搜索alpha3=fsearch(x2,s);x3=x2+alpha2*s;x0=x3;elseif F2>G3x0=x3;elsex0=x2;endEnd子文件JT,JH进退法程序代码�56555 .function [minx,maxx] = minJT(f,x0,h0,eps) format long;if nargin == 3eps = 1.0e-6;endk = 0;h = h0;while 1x4 = x1 + h;k = k+1;f4 = subs(f, findsym(f),x4); f1 = subs(f, findsym(f),x1); if f4 < f1x2 = x1;x1 = x4;f2 = f1;f1 = f4;h = 2*h;elseif k==1h = -h;x2 = x4;f2 = f4;elsex3 = x2;x2 = x1;break;endendendminx = min(x1,x3);maxx = x1+x3 - minx;format short;黄金分割法程序代码�56555 . function [x,minf] = minHJ(f,a,b,eps) format long;if nargin == 3eps = 1.0e-6;endl = a + 0.382*(b-a);u = a + 0.618*(b-a);k=1;tol = b-a;while tol>eps && k<100000fl = subs(f , findsym(f), l);fu = subs(f , findsym(f), u);a = l;l = u;u = a + 0.618*(b - a);elseb = u;u = l;l = a + 0.382*(b-a);endk = k+1;tol = abs(b - a);endif k == 100000disp('ÕÒ²»μ½215 ?î208 ?¡214 Oμ£¡); x = NaN;minf = NaN;return;endx = (a+b)/2;minf = subs(f, findsym(f),x);format short;。
鲍威尔算法实验报告
一、实验目的1. 理解鲍威尔算法的基本原理和步骤。
2. 掌握鲍威尔算法在求解非线性方程组中的应用。
3. 分析鲍威尔算法的收敛速度和精度。
二、实验原理鲍威尔算法是一种迭代算法，用于求解非线性方程组。
该算法的基本思想是利用相邻迭代的残差向量构造一个线性方程组的系数矩阵，进而求出近似解。
具体步骤如下：1. 初始化：选择初始点 \( x_0 \) 和 \( x_1 \)，计算初始残差向量\( \mathbf{r}_0 = f(x_0) \) 和 \( \mathbf{r}_1 = f(x_1) \)。
2. 构造系数矩阵：根据残差向量 \( \mathbf{r}_0 \) 和 \( \mathbf{r}_1 \) 构造系数矩阵 \( \mathbf{A} \)。
3. 求解线性方程组：求解线性方程组 \( \mathbf{A} \mathbf{x} =\mathbf{r}_1 \)，得到系数 \( \mathbf{x} \)。
4. 更新近似解：根据系数 \( \mathbf{x} \) 更新近似解 \( x_2 = x_1 +\mathbf{x} \)。
5. 检查收敛性：计算新的残差向量 \( \mathbf{r}_2 = f(x_2) \)，如果满足收敛条件，则停止迭代；否则，返回步骤2。
三、实验内容1. 选择非线性方程组：设非线性方程组为\[\begin{cases}f_1(x_1, x_2) = x_1^2 + x_2^2 - 1 = 0 \\f_2(x_1, x_2) = x_1^3 - x_2 - 1 = 0\end{cases}\]2. 选择初始点：取 \( x_0 = (0, 0) \) 和 \( x_1 = (1, 1) \)。
3. 运行鲍威尔算法：根据上述步骤，编写程序求解该非线性方程组。
四、实验结果与分析1. 实验结果：经过多次迭代，鲍威尔算法得到近似解为 \( x_1 \approx 0.5285 \)，\( x_2 \approx 0.8572 \)。
鲍威尔算法matlab程序f.doc
function f=fun(x)f=10*(x(1)+x(2)-5)^2+(x(1)-x(2))^2; function f=fx(x0,alpha,s)x1=x0+alpha*s;f=fun(x1);function f=fsearch(x0,s)%利用进退法确定高低高区间alpha1=0;h=0.1;alpha2=alpha1+h;f1=fx(x0,alpha1,s);f2=fx(x0,alpha2,s);if f1>f2alpha3=alpha2+h;f3=fx(x0,alpha3,s);while f2>f3alpha1=alpha2;alpha2=alpha3;alpha3=alpha3+h;f2=f3;f3=fx(x0,alpha3,s);endelseh=-h;v=alpha1;alpha1=alpha2;alpha2=v;v=f1;f1=f2;f2=v;alpha3=alpha2+h;f3=fx(x0,alpha3,s);while f2>f3alpha1=alpha2;alpha2=alpha3;alpha3=alpha3+h;f2=f3;f3=fx(x0,alpha3,s);endenda=min(alpha1,alpha3);b=max(alpha1,alpha3);%利用黄金分割点法求解alpha1=a+0.382*(b-a);alpha2=a+0.618*(b-a);f1=fx(x0,alpha1,s);f2=fx(x0,alpha2,s);while abs(a-b)>0.001if f1>f2a=alpha1;alpha1=alpha2;f1=f2;alpha2=a+0.618*(b-a);f2=fx(x0,alpha2,s);elseb=alpha2;alpha2=alpha1;f2=f1;alpha1=a+0.382*(b-a);f1=fx(x0,alpha1,s);endendf=0.5*(a+b);clear%初始点x0=[0;0];%搜索方向e1=[1;0];e2=[0;1];G0=fun(x0);F0=G0;%第一次迭代%沿着e1alpha1=fsearch(x0,e1);x1=x0+alpha1*e1;F1=fun(x1);delta1=F0-F1;% 沿着方向e2；alpha2=fsearch(x1,e2);x2=x1+alpha2*e2;F2=fun(x2);G2=F2;delta2=F1-F2;deltam=max(delta1,delta2);%映射点x3=2*x2-x0;G3=fun(x3);if G3<G0 & (G0-2*G2+G3)*(G0-G2-deltam)^2<0.5*deltam*(G0-G3)^2s=x2-x0;%方向替换e1=e2;e2=s;% 沿着方向s 进行搜索alpha3=fsearch(x2,s);x3=x2+alpha2*s;x0=x3;elseif F2>G3x0=x3;elsex0=x2;endEnd子文件JT,JH进退法程序代码56555 .function [minx,maxx] = minJT(f,x0,h0,eps) format long;if nargin == 3eps = 1.0e-6;endx1 = x0;k = 0;h = h0;while 1x4 = x1 + h;k = k+1;f4 = subs(f, findsym(f),x4);f1 = subs(f, findsym(f),x1);if f4 < f1x2 = x1;x1 = x4;f2 = f1;f1 = f4;h = 2*h;elseif k==1h = -h;x2 = x4;f2 = f4;elsex3 = x2;x2 = x1;x1 = x4;break;endendendminx = min(x1,x3);maxx = x1+x3 - minx;format short;黄金分割法程序代码56555 .function [x,minf] = minHJ(f,a,b,eps) format long;if nargin == 3eps = 1.0e-6;endl = a + 0.382*(b-a);u = a + 0.618*(b-a);k=1;tol = b-a;while tol>eps && k<100000fl = subs(f , findsym(f), l);fu = subs(f , findsym(f), u);if fl > fua = l;l = u;u = a + 0.618*(b - a);elseb = u;u = l;l = a + 0.382*(b-a);endk = k+1;tol = abs(b - a);endif k == 100000disp('ÕÒ²»μ½215 ?î208 ?¡214 Oμ£¡); x = NaN;minf = NaN;return;endx = (a+b)/2;minf = subs(f, findsym(f),x);format short;。
鲍威尔法实验报告
多维优化设计实验报告——修正鲍威尔法求极值一、 实验内容用多维优化设计方法解决以下问题： 求无约束优化问题min f(x)=60－10x1－4x2+x21+x22－xx 21的最优点[]x x XT*2*1,*= ；给定初始点[]0,0)0()1(0TXX==计算精度ε=0.0001。
利用修正鲍威尔共轭方向法，求得满足计算精度的最小极值点的近似解。
二、程序流程图修正鲍威尔法的计算框图如下所示：三、程序源代码#include "stdio.h" #include "stdlib.h" #include "math.h"double objf(double x[]){double ff;ff=x[0]*x[0]+x[1]*x[1]-x[0]*x[1]-10*x[0]-4*x[1]+60;return(ff);}void jtf(double x0[],double h0,double s[],int n,double a[],double b[]) {int i;double *x[3],h,f1,f2,f3;for(i=0;i<3;i++)x[i]=(double *)malloc(n*sizeof(double));h=h0;for(i=0;i<n;i++)*(x[0]+i)=x0[i];f1=objf(x[0]);for(i=0;i<n;i++)*(x[1]+i)=*(x[0]+i)+h*s[i];f2=objf(x[1]);if(f2>=f1){h=-h0;for(i=0;i<n;i++)*(x[2]+i)=*(x[0]+i);f3=f1;for(i=0;i<n;i++){*(x[0]+i)=*(x[1]+i);*(x[1]+i)=*(x[2]+i);}f1=f2;f2=f3;}for(;;){h=2*h;for(i=0;i<n;i++)*(x[2]+i)=*(x[1]+i)+h*s[i]; f3=objf(x[2]);if(f2<f3) break;else{ for(i=0;i<n;i++){*(x[0]+i)=*(x[1]+i);*(x[1]+i)=*(x[2]+i);}f1=f2;f2=f3;}}if(h<0)for(i=0;i<n;i++){a[i]=*(x[2]+i);b[i]=*(x[0]+i);}elsefor(i=0;i<n;i++){a[i]=*(x[0]+i);b[i]=*(x[2]+i);}for(i=0;i<3;i++)free(x[i]);}double gold(double a[],double b[],double eps,int n,double xx[]) {int i;double f1,f2,*x[2],ff,q,w;for(i=0;i<2;i++)x[i]=(double *)malloc(n*sizeof(double));for(i=0;i<n;i++){*(x[0]+i)=a[i]+0.618*(b[i]-a[i]);*(x[1]+i)=a[i]+0.382*(b[i]-a[i]);}f1=objf(x[0]);f2=objf(x[1]);do{if(f1>f2){for(i=0;i<n;i++){b[i]=*(x[0]+i);*(x[0]+i)=*(x[1]+i);}f1=f2;for(i=0;i<n;i++)*(x[1]+i)=a[i]+0.382*(b[i]-a[i]);f2=objf(x[1]);}else{ for(i=0;i<n;i++){a[i]=*(x[1]+i);*(x[1]+i)=*(x[0]+i);}f2=f1;for(i=0;i<n;i++)*(x[0]+i)=a[i]+0.618*(b[i]-a[i]);f1=objf(x[0]);}q=0;for(i=0;i<n;i++)q=q+(b[i]-a[i])*(b[i]-a[i]);w=sqrt(q);}while(w>eps);for(i=0;i<n;i++)xx[i]=0.5*(a[i]+b[i]);ff=objf(xx);for(i=0;i<2;i++)free(x[i]);return(ff);}double oneoptim(double x0[],double s[],double h0,double epsg,int n,double x[]){double *a,*b,ff;a=(double *)malloc(n*sizeof(double));b=(double *)malloc(n*sizeof(double));jtf(x0,h0,s,n,a,b);ff=gold(a,b,epsg,n,x);free(a);free(b);return (ff);}double powell(double p[],double h0,double eps,double epsg,int n,double x[]){int i,j,m;double *xx[4],*ss,*s;double f,f0,f1,f2,f3,fx,dlt,df,sdx,q,d;ss=(double *)malloc(n*(n+1)*sizeof(double));s=(double *)malloc(n*sizeof(double));for(i=0;i<n;i++){for(j=0;j<=n;j++)*(ss+i*(n+1)+j)=0;*(ss+i*(n+1)+i)=1;}for(i=0;i<4;i++)xx[i]=(double *)malloc(n*sizeof(double));for(i=0;i<n;i++)*(xx[0]+i)=p[i];for(;;){for(i=0;i<n;i++){*(xx[1]+i)=*(xx[0]+i);x[i]=*(xx[1]+i);}f0=f1=objf(x);dlt=-1;for(j=0;j<n;j++){for(i=0;i<n;i++){*(xx[0]+i)=x[i];*(s+i)=*(ss+i*(n+1)+j);}f=oneoptim(xx[0],s,h0,epsg,n,x);df=f0-f;if(df>dlt){dlt=df;m=j;}}sdx=0;for(i=0;i<n;i++)sdx=sdx+fabs(x[i]-(*(xx[1]+i)));if(sdx<eps){free(ss);free(s);for(i=0;i<4;i++)free(xx[i]);return(f);}for(i=0;i<n;i++)*(xx[2]+i)=x[i];f2=f;for(i=0;i<n;i++){*(xx[3]+i)=2*(*(xx[2]+i)-(*(xx[1]+i))); x[i]=*(xx[3]+i);}fx=objf(x);f3=fx;q=(f1-2*f2+f3)*(f1-f2-dlt)*(f1-f2-dlt); d=0.5*dlt*(f1-f3)*(f1-f3);if((f3<f1)||(q<d)){if(f2<=f3)for(i=0;i<n;i++)*(xx[0]+i)=*(xx[2]+i);elsefor(i=0;i<n;i++)*(xx[0]+i)=*(xx[3]+i);}else{for(i=0;i<n;i++){*(ss+(i+1)*(n+1))=x[i]-(*(xx[1]+i));*(s+i)=*(ss+(i+1)*(n+1));}f=oneoptim(xx[0],s,h0,epsg,n,x);for(i=0;i<n;i++)*(xx[0]+i)=x[i];for(j=m+1;j<=n;j++)for(i=0;i<n;i++)*(ss+i*(n+1)+j-1)=*(ss+i*(n+1)+j);}}}void main(){double p[]={1,2};double ff,x[2];ff=powell(p,0.3,0.001,0.0001,2,x);printf("x[0]=%f,x[1]=%f,ff=%f\n",x[0],x[1],ff); getchar();}四、程序运行结果。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验报告
实验名称：鲍威尔法
院（系）：机电学院
专业班级：机械制造及其自动化
姓名：
学号：
2013年5 月13 日
实验一：鲍威尔法实验日期：2013年5 月
13 日
一、实验目的
了解MATLAB的基本运用
了解MATLB在优化中的使用
二、实验原理
鲍威尔法也是一种共轭法，利用函数值来构造共轭方向，同时引入坐标轮换的概念，利用搜索前后两个点之间的连线形成新的共轭方向，替换旧的共轭方向。
三、实验内容
鲍威尔法程序：
x0=[12;10];
xk=x0;
ie=10^(-7);
ae=1;
%初始化搜索方向
d=zeros(2,2);
d(:,1)=[1;0];
d(:,2)=[0;1];
Inc=zeros(2,1);
k=0;
MLN=100;
%迭代求解
while (ae>ie&&k<MLN)
syms x1
syms x2
xktemp=xk;
fun1=fun(x1,x2);
fun1=inline(fun1);
f0=feval(fun1,xk(1),xk(2));
F0=f0;
if k>0
F0=eval(F0);
end
%沿d1方向进行一维搜索
syms a
syms x1;
syms x2;
xk1=xk+a*d(:,1);
x1=xk1(1);
x2=xk1(2);
fun1=fun(x1,x2);
fxa=diff(fun1,'a');
a=solve(fxa);
xk1=inline(xk1);
xk1=feval(xk1,a);
xk1(1)=eval(xk1(1));
xk1(2)=eval(xk1(2));
syms x1;
syms x2;
fun1=fun(x1,x2);
fun1=inline(fun1);
f1=feval(fun1,xk1(1),xk1(2)); f1=eval(f1);
Inc(1)=f0-f1;
%沿d2方向进行搜索
syms a;
syms x1;
syms x2;
xk2=xk1+a*d(:,2);
x1=xk2(1);
x2=xk2(2);
fun1=fun(x1,x2);
fxa=diff(fun1,'a');
a=solve(fxa);
xk2=inline(xk2);
xk2=feval(xk2,a);
xk2(1)=eval(xk2(1));
xk2(2)=eval(xk2(2));
syms x1;
syms x2;
fun1=fun(x1,x2);
fun1=inline(fun1);
f2=feval(fun1,xk2(1),xk2(2));
f2=eval(f2);
F2=f2;
Inc(2)=f1-f2;
[Incm,row]=max(Inc);
x3=2*xk2-xk;%计算反射点
syms x1;
syms x2;
fun1=fun(x1,x2);
fun1=inline(fun1);
f3=feval(fun1,x3(1),x3(2));
f3=eval(f3);
F3=f3;
temp1=(F0-2*F2+F3)*(F0-F2-Incm)^2; temp2=0.5*Incm*(F0-F3)^2;
%判断是否更换搜索方向
if (F3<F0&&temp1<temp2)
syms a;
syms x1;
syms x2;
d(:,row)=xk2-xk;
xk=xk2+a*d(:,row);
x1=xk(1);
x2=xk(2);
fun1=fun(x1,x2);
fxa=diff(fun1,'a');
a=solve(fxa);
xk=inline(xk);
xk=feval(xk,a);
%不更换搜索方向else if F2<F3
xk=xk2;
else
xk=x3;
end
end
xkerror=eval(xk2-xktemp); ae=norm(xkerror);
k=k+1;
end
x=eval(xk)
函数程序：
function [f]=fun(x1,x2)
f=2*x1^2+4*x1*x2+x2^2
执行结果：
x =
四、实验小结
通过本实验了解了了matlab的基本操作方法，了解鲍威尔法的原理与基本运用。