PSO网络的一个算法，不知道是哪出错了，求指点

永不漠灭的轨迹-

>> function main()
" {0 V" f- f' N/ K5 d9 W, [clc;clear all;close all;
5 E3 v  L& R6 M& f7 d* V$ [tic;                 %程序运行计时
  E! V( S# z# ]4 O4 uE0=0.001;            %允许误差
5 I$ @- d1 L! z9 }# Z) pMaxNum=100;           %粒子最大迭代次数
: |; C8 \- i$ J8 @narvs=1;            %目标函数自变量个数
+ g+ J0 y  A3 q# ^particlesize=30;    %粒子群规模
c1=2;               %每个粒子个体学习因子，也称加速常数
) s+ Y3 Q4 j2 x+ r! u* p6 d) u8 G% ~c2=2;               %每个粒子的社会学习因子，也称加速常数
7 z, N% x# y; O( s% p8 t$ Tw=0.6;              %惯性因子
8 U' F1 Y" j6 ^" o3 {4 O3 \vmax=0.8;           %粒子最大飞翔速度
& M) m6 @4 R# N  n4 Ux=-5+10*rand(particlesize,narvs); %粒子所在位置
7 t! W' A0 R6 i: k; z$ Dv=2*rand(particlesize,navrs);        %粒子的飞翔速度
/ R) K& q- ]  g%用inline定义适应度函数以便于将子函数与主程序文件放在一起，
( Y9 |/ B6 ?  y1 p% V%目标函数是：y=1+(2.1*(1-x+2*x.^2).*exp(-x.^2/2)）
9 [% N6 B3 K1 H8 N%inline命令定义实应度函数如下
fitness=inline('1/(1+(2.1*(1-x+2*x.^2).*exp(-x.^2/2)))','x');
%inline定义的实验度函数会使程序运行速度大大降低
for i=1:particlesize
$ C: I& p  c: U- e! Z" ^    for j=1:narvs
        f(i)=fitness(x(i,j));
    end
end
personalbest_x=x;
, @$ U2 N- q) m8 W4 N! t& {# Kpersonalbest_faval=f;
[globalbest_faval i]=min(personalbest_faval);
globalbest_x=personalbest_x(i,;
9 h$ ~7 i# E( a5 i* D( ^k=1;
while k<=MaxNum
; ?5 C' B" o; z7 n    for i=1:particlesize
$ P8 b4 l- {' z, b; Q' _7 q9 N        for j=1:narvs;
        f(i)=fitness(x(i,j));
" y/ ?) M6 S, m+ g        end
- O% Q( M3 G( v$ j    if f(i)<personalbest_faval(i)   %判断当前位置是否历史上最佳位置
) i7 n, k9 Q' Q# O5 I8 }        personalbest_faval(i)=f(i);
& x; {6 k* O) \, m+ N        personalbest_x(i,=x(i,;
    end
0 I0 y8 Z3 `" _& }    end
    [globalbest_faval i]=min(personalbest_faval);
    globalbest_x=personalbest_x(i,;
; {" }, {: @' m* p" K, v7 K% }    for i=1:particlesize           %更新粒子群里每个个体的最新位置
. m+ R6 d+ q) s$ W) Z        v(i,=w*v(i,+c1*rand*(personalbest_x(i,-x(i,)...
6 W# x; I$ o; S0 v$ `            +c2*rand*(globalbest_x-x(i,);
$ Z: O! Z+ u/ o, m: G         for j=1:narvs             %判断粒子的飞翔速度是否超过了最大飞行速度
             if v(i.j)>vmax;   
0 y% h$ X+ Z& _" b                 v(i,j)=vmax;
% P% p" N8 W3 S/ c; @5 c0 `( w             else if v(i,j)<-vmax;
. d9 W% Z9 i" t/ K! x: h4 O                     v(i,j)=-vmax;
                 end
             end
             x(i,=x(i,+v(i,;
         end
1 g: ~& @" J, J( h; v/ `; H         if abs(globalbest_faval)<E0,break,end
         k=k+1;
    end
    Valuel=1/globalbest_faval-1;Valuel=num2str(Valuel);
% m) ^2 H. q! {' T    %strcat指令可以实现字符的组合输出
. H$ k$ S; m5 q+ h% z    disp(strcat('the maximun valuel','=',Valuel));
  %输出最大值所在的横坐标位置
    Value2=globalbest_x; Value2=num2str(Value2);
     disp(strcat('the corresponding coordinate','=',Value2));
     x=-5:0.01:5;
, P2 V* ^$ x5 a9 h* V, j     y=2.1*(1-x+2*x.^2).*exp(-x.^2/2);
  v! _% _. s4 F# ?8 o     plot(x,y,'m-','linewidth',3);
# P8 s, Z, L2 p( M1 d9 R     hold on;
     plot(globalbest_x,1/globalbest_faval-1,'kp','linewidth',4);
) u. z% J6 e+ ~7 J     legend('目标函数','搜索到的最大值');
) ^/ Q' Q& Q  W$ @     xlabel('x');ylabel('y');grid on;toc;

$ L9 H" y( w$ Q/ O/ y+ y3 _3 Q
% F( J, U6 s9 |( Z错误提示 ??? function main()    |
, D' G4 X; C: S7 `Error: Function definitions are not permitted in this context.

去掉function main()后 提示??? Undefined function or variable 'navrs'.错误。求指点！
/ ?, R& \( }1 c* M$ ]; S6 i; @

平凡之不凡

代码重新发一下，格式出问题了，你可以调整一下

永不漠灭的轨迹-

格式是对的  是不是那个自定义函数的问题，我刚刚查了一下，说是要先定义M文件。

永不漠灭的轨迹-

平凡之不凡 发表于 2014-8-14 09:37
代码重新发一下，格式出问题了，你可以调整一下

格式是对的  是不是那个自定义函数的问题，我刚刚查了一下，说是要先定义M文件。

永不漠灭的轨迹-

function main() clc;clear all;close all; tic;                 %程序运行计时 E0=0.001;            %允许误差 MaxNum=100;           %粒子最大迭代次数 narvs=1;            %目标函数自变量个数 particlesize=30;    %粒子群规模 c1=2;               %每个粒子个体学习因子，也称加速常数 c2=2;               %每个粒子的社会学习因子，也称加速常数 w=0.6;              %惯性因子 vmax=0.8;           %粒子最大飞翔速度 x=-5+10*rand(particlesize,narvs); %粒子所在位置 v=2*rand(particlesize,navrs);        %粒子的飞翔速度 %用inline定义适应度函数以便于将子函数与主程序文件放在一起， %目标函数是：y=1+(2.1*(1-x+2*x.^2).*exp(-x.^2/2)） %inline命令定义实应度函数如下 fitness=inline('1/(1+(2.1*(1-x+2*x.^2).*exp(-x.^2/2)))','x'); %inline定义的实验度函数会使程序运行速度大大降低 for i=1:particlesize     for j=1:narvs         f(i)=fitness(x(i,j));     end end personalbest_x=x; personalbest_faval=f; [globalbest_faval i]=min(personalbest_faval); globalbest_x=personalbest_x(i,; k=1; while k<=MaxNum     for i=1:particlesize         for j=1:narvs;         f(i)=fitness(x(i,j));         end     if f(i)<personalbest_faval(i)   %判断当前位置是否历史上最佳位置         personalbest_faval(i)=f(i);         personalbest_x(i,=x(i,;     end     end     [globalbest_faval i]=min(personalbest_faval);     globalbest_x=personalbest_x(i,;     for i=1:particlesize           %更新粒子群里每个个体的最新位置         v(i,=w*v(i,+c1*rand*(personalbest_x(i,-x(i,)...             +c2*rand*(globalbest_x-x(i,);          for j=1:narvs             %判断粒子的飞翔速度是否超过了最大飞行速度              if v(i.j)>vmax;                      v(i,j)=vmax;              else if v(i,j)<-vmax;                      v(i,j)=-vmax;                  end              end              x(i,=x(i,+v(i,;          end          if abs(globalbest_faval)<E0,break,end          k=k+1;     end     Valuel=1/globalbest_faval-1;Valuel=num2str(Valuel);     %strcat指令可以实现字符的组合输出     disp(strcat('the maximun valuel','=',Valuel));   %输出最大值所在的横坐标位置     Value2=globalbest_x; Value2=num2str(Value2);      disp(strcat('the corresponding coordinate','=',Value2));      x=-5:0.01:5;      y=2.1*(1-x+2*x.^2).*exp(-x.^2/2);      plot(x,y,'m-','linewidth',3);      hold on;      plot(globalbest_x,1/globalbest_faval-1,'kp','linewidth',4);      legend('目标函数','搜索到的最大值');      xlabel('x');ylabel('y');grid on;toc;

平凡之不凡

1. function main()
   6 v% B9 {5 f/ ~$ y2 e5 s/ f\" p5 h
2. clc;clear all;close all;- x\" }( }9 e8 N2 T; i! s
   
3. tic;                              %程序运行计时
   8 s  o- q! T' r+ y' _9 G
4. E0=0.001;                        %允许误差8 h5 M& c! x! ^+ ~
   
5. MaxNum=100;                    %粒子最大迭代次数3 \& F! r: l% U
   
6. narvs=1;                         %目标函数的自变量个数1 ]0 B' Y! j2 f+ ]+ @
   
7. particlesize=30;                    %粒子群规模
   # m: X& V; g5 m1 F9 S1 b
8. c1=2;                            %每个粒子的个体学习因子，也称为加速常数. J. x! ?& m( a8 \5 C
   
9. c2=2;                            %每个粒子的社会学习因子，也称为加速常数\" j7 I: I9 @. t  L3 ?' \\" e
   
10. w=0.6;                           %惯性因子3 }2 {* X$ Q0 U- d/ s
    
11. vmax=0.8;                        %粒子的最大飞翔速度: q  }1 w: X' J0 }6 {8 L
    
12. x=-5+10*rand(particlesize,narvs);     %粒子所在的位置
    2 r* q1 t% s8 U$ h
13. v=2*rand(particlesize,narvs);         %粒子的飞翔速度/ X/ K6 |. u8 \. b! `. i( C. w
    
14. %用inline定义适应度函数以便将子函数文件与主程序文件放在一起，
    , X5 s9 P4 j* R2 X2 ^  E% u
15. %目标函数是：y=1+(2.1*(1-x+2*x.^2).*exp(-x.^2/2))
    : Q( {# a5 x/ K! _. G5 q, z7 P
16. %inline命令定义适应度函数如下：
    * f- z! v$ |1 ?# w$ M! k( M* k! ~+ M
17. fitness=inline('1/(1+(2.1*(1-x+2*x.^2).*exp(-x.^2/2)))','x');5 W# i7 P9 u1 m1 I# Y; Y) J
    
18. %inline定义的适应度函数会使程序运行速度大大降低
    / ^, Y4 x, N9 o4 q( r- q. |; j1 n' h
19. for i=1:particlesize
    6 x2 |+ F: Q( Z
20.     for j=1:narvs1 f7 L7 w) r* L& q9 P
    
21.         f(i)=fitness(x(i,j));
    ; D( K  m  N: p. a# B! H% f
22.     end
    3 |# h+ W( X# z9 E! G$ W
23. end. T! N' B% Q: r2 ?  E8 d! E  M0 c
    
24. personalbest_x=x;
    ; a% W3 r. M4 n/ m, f$ H) b
25. personalbest_faval=f;- F8 q3 M& Q9 t4 }4 v, j
    
26. [globalbest_faval i]=min(personalbest_faval);5 e& j4 Y0 l0 A9 ?: B0 h: e
    
27. globalbest_x=personalbest_x(i,:);; B' ?% \  ~! `- \, H6 z1 Z) e
    
28. k=1;
    : O6 _/ a4 C$ {7 Z$ M4 S
29. while k<=MaxNum& ?- @$ g3 ^# j\" J, |+ A
    
30.     for i=1:particlesize
    + ]+ O7 k6 @, T: J2 P8 q) X
31.         for j=1:narvs: M, e5 X! Y/ {
    
32.             f(i)=fitness(x(i,j));
    2 X% N  p7 a5 b+ ~) O
33.         end
    0 K\" P) ^3 O, |+ |
34.         if f(i)<personalbest_faval(i) %判断当前位置是否是历史上最佳位置) Z\" p: i6 J- j7 v
    
35.             personalbest_faval(i)=f(i);
    2 \8 \+ F- P$ F% q& Y$ [
36.             personalbest_x(i,:)=x(i,:);
    * t+ p7 m9 h! T0 _' b' |; R# w9 i
37.         end
    1 G7 l5 V0 n  S& O; X
38.     end: s  N! T- u) m. r$ ?
    
39.     [globalbest_faval i]=min(personalbest_faval);
    0 p- w5 x8 _\" {2 O! T3 L( c
40.     globalbest_x=personalbest_x(i,:);
    , ?) i\" G1 e\" T2 f% Y/ S4 B( ^4 T
41.     for i=1:particlesize %更新粒子群里每个个体的最新位置! z( A6 A9 s! a! Q% J4 Y
    
42.         v(i,:)=w*v(i,:)+c1*rand*(personalbest_x(i,:)-x(i,:))...
    * m; o# ~5 ?. l$ L( L* {! E
43.             +c2*rand*(globalbest_x-x(i,:));4 E5 Z3 T; g, ?' X: a3 m
    
44.         for j=1:narvs    %判断粒子的飞翔速度是否超过了最大飞翔速度
    + T7 s: [  }( ?\" a
45.             if v(i,j)>vmax;
    # c# V, U8 ^) d2 K
46.                 v(i,j)=vmax;' n7 k5 R7 g  }4 A+ _9 v7 K\" a
    
47.             elseif v(i,j)<-vmax;9 F4 u; w- k5 \; p1 c! X
    
48.                 v(i,j)=-vmax;% x: B+ d. X$ r. @
    
49.             end$ B& K1 R% h& g. {
    
50.         end9 W* n' U) T7 _\" g1 ?
    
51.         x(i,:)=x(i,:)+v(i,:);
    1 l+ |9 j7 v% e3 v/ _6 T$ w
52.     end
    3 {& R: v# l3 }- w. h3 v( W1 {
53.     if abs(globalbest_faval)<E0,break,end
    5 Y% P, H$ N8 S, R
54.     k=k+1;# {& j$ g  z, R\" D\" X
    
55. end
    $ d; h! V: Y6 `1 s7 p: \( D4 l
56. Value1=1/globalbest_faval-1; Value1=num2str(Value1);6 q& l4 L% R4 L0 E\" l\" E
    
57. % strcat指令可以实现字符的组合输出, C- o; z% R( x! k6 P9 H
    
58. disp(strcat('the maximum value','=',Value1));3 h3 ]3 B% b  D& C! p; s
    
59. %输出最大值所在的横坐标位置# S: {1 s5 P2 |( D' A
    
60. Value2=globalbest_x; Value2=num2str(Value2);, J# N! d: W) r; P$ J, h5 y
    
61. disp(strcat('the corresponding coordinate','=',Value2));2 r4 `+ t2 v! p% p
    
62. x=-5:0.01:5;
    \" \: B( K- v/ b1 ~
63. y=2.1*(1-x+2*x.^2).*exp(-x.^2/2);' A7 K% Q; Z) I1 M& x
    
64. plot(x,y,'m-','linewidth',3);
    3 P; K7 J8 k% g, z1 p
65. hold on;6 l* F3 M, N$ o8 U
    
66. plot(globalbest_x,1/globalbest_faval-1,'kp','linewidth',4);
    $ V9 w# w\" E$ R1 u+ U1 @6 P
67. legend('目标函数','搜索到的最大值');xlabel('x');ylabel('y');grid on;toc;
    4 C% P3 h& S. h2 a. F

复制代码

平凡之不凡

这个是调好的

永不漠灭的轨迹-

平凡之不凡 发表于 2014-8-17 08:58

3 o" ^- z  e+ m: T. [" F谢谢老师！