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

永不漠灭的轨迹-

>> function main()
clc;clear all;close all;
tic;                 %程序运行计时
E0=0.001;            %允许误差
MaxNum=100;           %粒子最大迭代次数
narvs=1;            %目标函数自变量个数
6 w" n/ b5 ]+ E6 {, p5 Mparticlesize=30;    %粒子群规模
c1=2;               %每个粒子个体学习因子，也称加速常数
. `/ l  t/ B8 P0 |) [c2=2;               %每个粒子的社会学习因子，也称加速常数
w=0.6;              %惯性因子
vmax=0.8;           %粒子最大飞翔速度
1 e: Q! Z! w  ^+ t: zx=-5+10*rand(particlesize,narvs); %粒子所在位置
1 c  Z1 F/ _. M6 cv=2*rand(particlesize,navrs);        %粒子的飞翔速度
6 d* o7 Z/ \, i) A0 L4 h% A+ g%用inline定义适应度函数以便于将子函数与主程序文件放在一起，
0 E- \( k& V* J5 U. T1 @; ]%目标函数是：y=1+(2.1*(1-x+2*x.^2).*exp(-x.^2/2)）
! W; t+ \) Z* V) P" w6 ~: l%inline命令定义实应度函数如下
  U& G$ r3 @$ |3 jfitness=inline('1/(1+(2.1*(1-x+2*x.^2).*exp(-x.^2/2)))','x');
%inline定义的实验度函数会使程序运行速度大大降低
for i=1:particlesize
    for j=1:narvs
% g! l$ C5 z5 Y        f(i)=fitness(x(i,j));
6 d9 r1 i5 u; ]" ^& Q0 i! Z    end
end
personalbest_x=x;
personalbest_faval=f;
% E; Z& |" B4 o0 D% g[globalbest_faval i]=min(personalbest_faval);
; @1 x& Z% K8 ^; g2 @  Vglobalbest_x=personalbest_x(i,;
k=1;
) }' e' w7 ]4 Qwhile k<=MaxNum
    for i=1:particlesize
        for j=1:narvs;
% x* S+ G' j4 h. o2 N9 i        f(i)=fitness(x(i,j));
" r0 y% e, t6 R. w/ j0 V        end
: F( t: U* x! ~    if f(i)<personalbest_faval(i)   %判断当前位置是否历史上最佳位置
        personalbest_faval(i)=f(i);
        personalbest_x(i,=x(i,;
  C8 n0 f; t* v4 B# S    end
: B/ C6 \3 f0 c# C8 s    end
8 F& z- X$ w9 y; r& y    [globalbest_faval i]=min(personalbest_faval);
. e0 e! k: F1 @2 l" q    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             %判断粒子的飞翔速度是否超过了最大飞行速度
: F0 J7 M1 o- a1 x             if v(i.j)>vmax;   
                 v(i,j)=vmax;
             else if v(i,j)<-vmax;
, B/ q. u1 d9 I$ _/ }% E6 y6 L                     v(i,j)=-vmax;
" B  U1 f% U. N* N3 g! h6 I. `( A                 end
( t6 p$ |( Y* Q' G             end
             x(i,=x(i,+v(i,;
         end
         if abs(globalbest_faval)<E0,break,end
+ ~$ [% @0 J8 ?0 q         k=k+1;
    end
+ k  C* {, h' H( G7 y" d6 W2 C    Valuel=1/globalbest_faval-1;Valuel=num2str(Valuel);
( X. o7 p- s8 x, `    %strcat指令可以实现字符的组合输出
    disp(strcat('the maximun valuel','=',Valuel));
  %输出最大值所在的横坐标位置
    Value2=globalbest_x; Value2=num2str(Value2);
$ b% L. S' V! ?+ Q     disp(strcat('the corresponding coordinate','=',Value2));
- {8 P9 ]2 [3 C7 ^/ b" r- K( [     x=-5:0.01:5;
/ ]2 F: r9 x- b9 G. k1 n& l     y=2.1*(1-x+2*x.^2).*exp(-x.^2/2);
$ |' z2 G& w9 t* o     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;
/ `% O' K) f- ^/ L8 K; J( {
8 t0 w, u3 J+ `! s
错误提示 ??? function main()    |
: y8 V! n6 b8 D; p6 ^- YError: Function definitions are not permitted in this context.
4 `$ u8 g- s6 y& Q8 X1 `
8 O, r0 H7 h* A9 R( l" n! M去掉function main()后 提示??? Undefined function or variable 'navrs'.错误。求指点！
& [! q/ Z" y, F% _

平凡之不凡

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

永不漠灭的轨迹-

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

永不漠灭的轨迹-

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

4 l* B1 [6 Q# f- i0 T
, I2 B  F" w: B+ z- e格式是对的  是不是那个自定义函数的问题，我刚刚查了一下，说是要先定义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()1 e- K! ~/ `/ w- S7 i
   
2. clc;clear all;close all;% a: p0 K+ \0 k
   
3. tic;                              %程序运行计时5 f3 p% ]' J8 N5 [! k
   
4. E0=0.001;                        %允许误差
   * A* j. a\" f% L& |7 ~' @# u
5. MaxNum=100;                    %粒子最大迭代次数
   & ^5 W1 F9 {' k. w) D/ {2 q
6. narvs=1;                         %目标函数的自变量个数\" g* @1 z7 K& A3 H
   
7. particlesize=30;                    %粒子群规模
   \" d! t- v) p2 p2 e5 I  }1 ~
8. c1=2;                            %每个粒子的个体学习因子，也称为加速常数
   3 z2 Y' h, X) l; u) X: s7 y9 X/ z; O
9. c2=2;                            %每个粒子的社会学习因子，也称为加速常数8 N7 |5 `0 z5 v1 n' p7 r2 b
   
10. w=0.6;                           %惯性因子' h7 f% J, ]8 p+ g; W
    
11. vmax=0.8;                        %粒子的最大飞翔速度: Z/ j4 A* e* s5 N/ X* j
    
12. x=-5+10*rand(particlesize,narvs);     %粒子所在的位置
    8 t  H% h7 A9 _
13. v=2*rand(particlesize,narvs);         %粒子的飞翔速度
    9 _+ ?! E$ a- ?5 V
14. %用inline定义适应度函数以便将子函数文件与主程序文件放在一起，4 g2 _& I! {+ b5 b4 P) \! J\" [
    
15. %目标函数是：y=1+(2.1*(1-x+2*x.^2).*exp(-x.^2/2))
    - T. D4 n( w- E; s
16. %inline命令定义适应度函数如下：& ]7 E  \% p- H& R3 ^. C
    
17. fitness=inline('1/(1+(2.1*(1-x+2*x.^2).*exp(-x.^2/2)))','x');
    6 I3 k' |7 G4 O! E' S
18. %inline定义的适应度函数会使程序运行速度大大降低
    5 X' Q; o  I: Z7 W' x
19. for i=1:particlesize% t. [3 B9 w1 |0 J
    
20.     for j=1:narvs
    . s- Q! j( l2 R- l$ b
21.         f(i)=fitness(x(i,j));/ a: s: V% B/ J; p- g$ _
    
22.     end
    ! J1 F/ P! [& ~) f8 e& V
23. end
    ' r- z; d, F$ A% q/ Z: j' p  ~
24. personalbest_x=x;0 z5 R/ u9 U& }! R) L
    
25. personalbest_faval=f;
    4 D+ O- h$ Q& B1 S1 |  S/ q
26. [globalbest_faval i]=min(personalbest_faval);) v2 N* t$ u5 R* E: a; C
    
27. globalbest_x=personalbest_x(i,:);
    : O( ~8 Y1 X- M
28. k=1;
    9 `  W: F; Q4 v0 Y+ o& r/ Q+ d* ?- ]
29. while k<=MaxNum* s7 p3 f' C3 H4 @
    
30.     for i=1:particlesize
    % O% `! H& C9 f5 L  T
31.         for j=1:narvs
    ( u9 k9 g4 ^6 E* D5 n
32.             f(i)=fitness(x(i,j));5 T9 T5 \/ i: h* b
    
33.         end
    5 e6 v* H  x0 x0 _/ s& t
34.         if f(i)<personalbest_faval(i) %判断当前位置是否是历史上最佳位置
    # m/ C- P2 E+ i# J& ]
35.             personalbest_faval(i)=f(i);
    , \' ^, q1 R5 U, a\" b9 p
36.             personalbest_x(i,:)=x(i,:);) t6 h+ V/ E' v2 A. B) Z
    
37.         end! l$ b2 u3 |) @) A: C/ ^
    
38.     end
    % f4 \* @: A7 U. j/ }' d! v3 G
39.     [globalbest_faval i]=min(personalbest_faval);- z: p# P6 O! ?# U
    
40.     globalbest_x=personalbest_x(i,:);. [  }6 |9 {$ j# Y' ?8 @4 F
    
41.     for i=1:particlesize %更新粒子群里每个个体的最新位置
    6 Z/ V) N* ^$ E, _
42.         v(i,:)=w*v(i,:)+c1*rand*(personalbest_x(i,:)-x(i,:))...
    $ @5 u2 h) t+ T; S; u7 T* q
43.             +c2*rand*(globalbest_x-x(i,:));
    ! _; i  Y' M+ D: V. V4 L
44.         for j=1:narvs    %判断粒子的飞翔速度是否超过了最大飞翔速度6 z) Q& C8 \3 ?
    
45.             if v(i,j)>vmax;9 m3 n+ r! l, U- w* M( r
    
46.                 v(i,j)=vmax;
    4 b: p( L2 a; n7 o  G5 W
47.             elseif v(i,j)<-vmax;2 b6 l) q% S0 S7 f
    
48.                 v(i,j)=-vmax;
    3 I) Z# V3 p, r8 `' m8 q
49.             end# R; l$ n4 o9 ~
    
50.         end0 d9 X3 ^: t& P1 B, J4 {) J- G
    
51.         x(i,:)=x(i,:)+v(i,:);8 p8 P6 V  v9 }: A1 N- s
    
52.     end
    , ?) ^0 z\" {& x/ J
53.     if abs(globalbest_faval)<E0,break,end
    4 s, I1 {, V3 p* [5 k/ v( `, O
54.     k=k+1;
    9 j% j. N\" P2 b9 E; d- Q! P& Y0 U
55. end
    , b; V% i' g) l6 R
56. Value1=1/globalbest_faval-1; Value1=num2str(Value1);
      {) E+ w# E3 m, P\" Z
57. % strcat指令可以实现字符的组合输出) X  G$ o  a% x* b0 S
    
58. disp(strcat('the maximum value','=',Value1));
    , ]5 \( s( G4 u+ N+ A: |
59. %输出最大值所在的横坐标位置
    ; p* Y: y( {* k/ ]5 ?' X7 W' ~$ x
60. Value2=globalbest_x; Value2=num2str(Value2);
    \" L$ E9 V; Y5 O9 q( I
61. disp(strcat('the corresponding coordinate','=',Value2));
    + ?' b0 M$ Y. |
62. x=-5:0.01:5;
    8 G( |& `4 L/ i, @5 [
63. y=2.1*(1-x+2*x.^2).*exp(-x.^2/2);
    # c2 H2 G1 E8 Z7 v
64. plot(x,y,'m-','linewidth',3);
    # F0 u1 P6 Q) i9 I
65. hold on;: D6 @! g) A6 \/ K
    
66. plot(globalbest_x,1/globalbest_faval-1,'kp','linewidth',4);
    ( Z: O5 J; ^* A) r
67. legend('目标函数','搜索到的最大值');xlabel('x');ylabel('y');grid on;toc;
    3 E\" m0 r; J* y) `7 @0 H$ N' w

复制代码

平凡之不凡

这个是调好的

永不漠灭的轨迹-

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

谢谢老师！