我编的2维寻优powell算法

ggyyll8683 · 发表于 2005-5-19 07:17:27

<

>同大家交流一下：

<

>function f=powellfun(x) f=x(1)^2+2*x(2)^2-2*x(1)*x(2)-4*x(1);

<

>close all; clear;
e=1e-5;n=2; %%%%%% d=[1 0;0 1]; x=[];x(1,

=[1 1];f=[];a=[];m=0;k=1;ap=0; f(1)=powellfun(x(1,

); r=f(1); for p=1:1:n for q=1:1:n for i=0:0.01:3 x((k+1),

=x(k,:)+i*d(k,:); b=powellfun(x((k+1),:)); if b<r f(k+1)=b; r=b; a(k)=i; end end x((k+1),:)=x(k,:)+a(k)*d(k,:); if k==n if abs(x(k+1,:)-x(1,:))<=e lastx=x(k+1,:); lastf=powellfun(lastx); hold else MAX=0; MIN=0; b=[]; for i=1:1:k b(i)=f(i)-f(i+1); end MIN=b(1); MAX=b(1); for i=2:1:k if MIN>b(i) MIN=b(i); end if MAX<b(i) MAX=b(i); end end for i=1:1:k if MAX==b(i) m=i-1; end end fl=powellfun((2*x(k+1,:)-x(1,:))); if (fl>=f(1))|(((f(1)-2*f(k+1)+fl)*(f(1)-f(k+1)-MAX)^2)>((f(1)-fl)^2*MAX)/2) f(1)=f(k+1); x(1,:)=x(k,:); k=1; else for i=1:1:m d(i,:)=d(i,:); end for i=(m+1):1:n-1 d(i,:)=d((i+1),:); end d(k,:)=(x(k+1,:)-x(1,:)); fp=f(k+1); for i=0:0.01:3 xp=x(k+1,:)+i*d(k,:); s=powellfun(xp); if s<fp fp=s; s=fp; ap=i; end end x(1,:)=x(k+1,:)+ap*d(k,:); f(1)=powellfun(x(1,:)); k=1; end end else k=k+1; end end end lastx=x(1,:) lastf=f(1)

要做多维的改进一下就可以了，在多维的情况下速度快！
而且可做无目标函数寻优！！

wangkw · 发表于 2005-11-22 05:10:42

<

>你好，请问你有DFP算法的程序吗，我用了POWELL方法，解的结果很不理想，不知道是不是我写错了，帮忙看看，谢谢，我的邮箱是<a href="mailtwangkwmail@163.com" target="_blank" >wangkwmail@163.com</A>可以与我联系，谢谢，期待你的回答
<

>不知道你的一维搜索采用的是什么方法，能求多极值的问题吗
<

>function x=test global x; global d; global e; dlt=1; x0=[15.1,5.3,0.1,0.1]; d=eye(4,4); %d1=eye(4,4); e=zeros(1,4); epnal=0.0000001; alfa=zeros(1,4); alfa5=0; k=1; %i=1; x=x0; x00=x; %while feixianxing(x)>=0.05 g=1; while g>=0.000001 for i=1:1:4 alfa(i)=yiweisousuo(x,d,i); x=x00+alfa(i)*d(i,

; x00=x; end r=norm(x0-x,'fro'); e=(x-x0)/r; alfa5=yiweisousuop(x,e) x=x+alfa5*e alfamax=max(alfa) %delt=alfamax/r; for i=1:1:4 if alfa(i)==max(alfa); s=i end end if alfamax*dlt/r>=epnal for i=1:1:4 if i==s d(i,

=e else d(i,

=d(i,:) end end dlt=alfamax*dlt/r; end g=feixianxing(x) x0=x; if g>0.001 break; end end
%end

wangkw · 发表于 2005-11-22 05:11:57

<

>你好，请问你有DFP算法的程序吗，我用了POWELL方法，解的结果很不理想，不知道是不是我写错了，帮忙看看，谢谢，我的邮箱是<a href="mailtwangkwmail@163.com" target="_blank" >wangkwmail@163.com</A>可以与我联系，谢谢，期待你的回答
<

>不知道你的一维搜索采用的是什么方法，能求多极值的问题吗
<

>function x=test global x; global d; global e; dlt=1; x0=[15.1,5.3,0.1,0.1]; d=eye(4,4); %d1=eye(4,4); e=zeros(1,4); epnal=0.0000001; alfa=zeros(1,4); alfa5=0; k=1; %i=1; x=x0; x00=x; %while feixianxing(x)>=0.05 g=1; while g>=0.000001 for i=1:1:4 alfa(i)=yiweisousuo(x,d,i); x=x00+alfa(i)*d(i,

; x00=x; end r=norm(x0-x,''fro''); e=(x-x0)/r; alfa5=yiweisousuop(x,e) x=x+alfa5*e alfamax=max(alfa) %delt=alfamax/r; for i=1:1:4 if alfa(i)==max(alfa); s=i end end if alfamax*dlt/r>=epnal for i=1:1:4 if i==s d(i,

=e else d(i,

=d(i,:) end end dlt=alfamax*dlt/r; end g=feixianxing(x) x0=x; if g>0.001 break; end end
%end

		自动登录	找回密码
密码			注-册-帐-号