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

永不漠灭的轨迹-

>> function main()
clc;clear all;close all;
/ u2 q9 z# ^( s3 q7 B5 e% H; ftic;                 %程序运行计时
E0=0.001;            %允许误差
) g2 e. A& |8 t; K6 AMaxNum=100;           %粒子最大迭代次数
narvs=1;            %目标函数自变量个数
particlesize=30;    %粒子群规模
c1=2;               %每个粒子个体学习因子，也称加速常数
( H' Z" j- a: _( F) }9 tc2=2;               %每个粒子的社会学习因子，也称加速常数
w=0.6;              %惯性因子
+ A! y' z1 n; |" [vmax=0.8;           %粒子最大飞翔速度
x=-5+10*rand(particlesize,narvs); %粒子所在位置
v=2*rand(particlesize,navrs);        %粒子的飞翔速度
9 b: ?+ N. }1 M" [- r: b3 V* f+ @%用inline定义适应度函数以便于将子函数与主程序文件放在一起，
) l2 g: G, j6 c%目标函数是：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
+ F, @$ G, A) w    for j=1:narvs
        f(i)=fitness(x(i,j));
    end
$ Z9 X: V. [/ ~$ l# gend
personalbest_x=x;
; \& ?$ a9 Y0 N9 d* E- ipersonalbest_faval=f;
[globalbest_faval i]=min(personalbest_faval);
globalbest_x=personalbest_x(i,;
7 t+ z  i& i9 c' Hk=1;
( }5 p, m; R) x, }; Cwhile k<=MaxNum
    for i=1:particlesize
  J+ }" b" }  e, ?/ t) [5 g6 q/ X        for j=1:narvs;
9 x) t# B1 u6 X3 l$ f* \        f(i)=fitness(x(i,j));
5 S% z$ g4 ^: M# f; E- r1 D2 {        end
% M6 ]0 j" F; s" |5 E; r8 ?" ]' o1 X    if f(i)<personalbest_faval(i)   %判断当前位置是否历史上最佳位置
        personalbest_faval(i)=f(i);
' A) m; z/ ^: [! [, C        personalbest_x(i,=x(i,;
    end
& k3 {; n. o, f! m( R6 k# u    end
+ W1 L5 R9 G$ S! ^) F    [globalbest_faval i]=min(personalbest_faval);
; O0 ~, _. T) c8 `    globalbest_x=personalbest_x(i,;
/ D/ m) G: M: }3 p    for i=1:particlesize           %更新粒子群里每个个体的最新位置
        v(i,=w*v(i,+c1*rand*(personalbest_x(i,-x(i,)...
/ V8 C3 ?' B$ {: a8 }' @            +c2*rand*(globalbest_x-x(i,);
         for j=1:narvs             %判断粒子的飞翔速度是否超过了最大飞行速度
% c! U: F( g' d) i5 s             if v(i.j)>vmax;   
4 q$ i# j* [3 O% A                 v(i,j)=vmax;
' u% ?! ]" g2 ~' H' C4 a* S( q             else if v(i,j)<-vmax;
                     v(i,j)=-vmax;
; S  W. |  {: {( D                 end
3 _2 Y* I+ U+ M' o) G! X             end
             x(i,=x(i,+v(i,;
; }: a( N" R+ u0 \, T         end
' x9 I9 ]' d! k: v! n$ _         if abs(globalbest_faval)<E0,break,end
         k=k+1;
& [( @3 v* ?3 O  i: Q1 P) \    end
    Valuel=1/globalbest_faval-1;Valuel=num2str(Valuel);
2 y( ^) ^* b6 P" m- z    %strcat指令可以实现字符的组合输出
    disp(strcat('the maximun valuel','=',Valuel));
  %输出最大值所在的横坐标位置
    Value2=globalbest_x; Value2=num2str(Value2);
     disp(strcat('the corresponding coordinate','=',Value2));
/ B  E' u8 }8 A/ y# \5 c& n     x=-5:0.01:5;
4 P+ @5 [) I% O) y. C4 R     y=2.1*(1-x+2*x.^2).*exp(-x.^2/2);
     plot(x,y,'m-','linewidth',3);
     hold on;
9 t1 W* C6 Y% j5 y3 {8 e     plot(globalbest_x,1/globalbest_faval-1,'kp','linewidth',4);
     legend('目标函数','搜索到的最大值');
* i7 {5 n- ~5 ^9 Z     xlabel('x');ylabel('y');grid on;toc;


4 s5 v8 u  i% r9 P/ }% @% E错误提示 ??? function main()    |
Error: Function definitions are not permitted in this context.
5 ^0 O' x. W7 m# T% A
去掉function main()后 提示??? Undefined function or variable 'navrs'.错误。求指点！

平凡之不凡

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

永不漠灭的轨迹-

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

永不漠灭的轨迹-

平凡之不凡 发表于 2014-8-14 09:37
1 C; V* m- H0 P  {4 `代码重新发一下，格式出问题了，你可以调整一下

  f9 ?; Y2 o3 ?& N+ H# p格式是对的  是不是那个自定义函数的问题，我刚刚查了一下，说是要先定义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()) m+ p& K8 i, y
   
2. clc;clear all;close all;
   9 \: c# M0 Q. P- D6 {5 F1 n; l
3. tic;                              %程序运行计时
   / u! R% S2 [' N1 l. ?5 K
4. E0=0.001;                        %允许误差/ |2 v8 x9 d# r% Z4 a
   
5. MaxNum=100;                    %粒子最大迭代次数# S0 E! l* ?4 P% q# R
   
6. narvs=1;                         %目标函数的自变量个数
   0 O' y% _+ T% f9 ]3 X5 H
7. particlesize=30;                    %粒子群规模! B# s- M% B5 C0 ], p! K
   
8. c1=2;                            %每个粒子的个体学习因子，也称为加速常数
   + u% G5 j) D3 @7 O
9. c2=2;                            %每个粒子的社会学习因子，也称为加速常数+ k2 X2 f- ]& u) F, b6 Q9 v
   
10. w=0.6;                           %惯性因子5 b8 D8 O  z3 ^+ {, R8 E2 p  C
    
11. vmax=0.8;                        %粒子的最大飞翔速度
    ' \7 @9 Y! s8 ?. U! \. M1 C
12. x=-5+10*rand(particlesize,narvs);     %粒子所在的位置- t\" ~: ~3 x/ t1 ]  i# U5 O
    
13. v=2*rand(particlesize,narvs);         %粒子的飞翔速度
    + r0 |& a/ y; O' }# [
14. %用inline定义适应度函数以便将子函数文件与主程序文件放在一起，4 r# F0 G$ m( q+ q% h
    
15. %目标函数是：y=1+(2.1*(1-x+2*x.^2).*exp(-x.^2/2))
    * m  x% w( t2 U
16. %inline命令定义适应度函数如下：$ ]% R- l9 ~0 A) ]
    
17. fitness=inline('1/(1+(2.1*(1-x+2*x.^2).*exp(-x.^2/2)))','x');
    5 W! u% L, i6 n$ W  d8 D\" j! R
18. %inline定义的适应度函数会使程序运行速度大大降低
    / i4 X, t3 s  E+ x
19. for i=1:particlesize8 f) y3 ?$ G\" j( n
    
20.     for j=1:narvs
    8 s5 n4 H1 ~. F: v* p
21.         f(i)=fitness(x(i,j));
    ( V+ c\" I9 P7 y/ x& {7 y4 z, n
22.     end0 F) ~8 A2 M2 f0 k
    
23. end. Z; K! o4 Y. y( v\" z
    
24. personalbest_x=x;
    3 }1 W) c$ g) _, S# }) {2 n
25. personalbest_faval=f;
    9 \9 @- b; U; W( n6 e* C* h
26. [globalbest_faval i]=min(personalbest_faval);
    ; \  n0 H& I8 K3 z) u
27. globalbest_x=personalbest_x(i,:);
    1 V: u% x. A9 T, c( Q/ ]  ~
28. k=1;1 o8 G7 s$ ~' I6 ?4 e2 y
    
29. while k<=MaxNum
    . n- R) }: h/ L( f
30.     for i=1:particlesize: O, r6 G' P0 s6 Q+ i3 f, q) [% F3 P
    
31.         for j=1:narvs! Q: `7 p\" p. C; r  w: o/ w  K
    
32.             f(i)=fitness(x(i,j));
      H+ L  V: W1 c; G- A- ^8 W+ B6 I
33.         end# U) [' N+ e# r
    
34.         if f(i)<personalbest_faval(i) %判断当前位置是否是历史上最佳位置8 c% J- m1 _\" ^& ^+ A, X* o# Z) A
    
35.             personalbest_faval(i)=f(i);
    & {8 M2 ?6 `! y
36.             personalbest_x(i,:)=x(i,:);
    ; |9 [& E$ ^9 ^
37.         end' R+ w/ o2 J0 Y' `5 }# F6 @% j
    
38.     end5 B  e; Q. z  R9 x/ l) @  `\" |7 D
    
39.     [globalbest_faval i]=min(personalbest_faval);
    , \3 p. K( A' }3 }+ o: P( Y8 G$ ~
40.     globalbest_x=personalbest_x(i,:);( ~$ w3 g% T, w6 }& p4 Y
    
41.     for i=1:particlesize %更新粒子群里每个个体的最新位置3 q. y/ z7 {8 E: H# f% a; J5 x
    
42.         v(i,:)=w*v(i,:)+c1*rand*(personalbest_x(i,:)-x(i,:))...
    - q8 r' H5 v9 T( F. X: Y& z
43.             +c2*rand*(globalbest_x-x(i,:));, U* z2 X' P) j+ T5 Z- U
    
44.         for j=1:narvs    %判断粒子的飞翔速度是否超过了最大飞翔速度0 q- n' y8 A* f/ ]/ H  b# g5 k, q
    
45.             if v(i,j)>vmax;; ~1 O, K: W: C) S/ w. O# j
    
46.                 v(i,j)=vmax;4 P. d7 X2 v$ ~: T: `
    
47.             elseif v(i,j)<-vmax;* u  b+ z9 a* u; q5 Q! p' \( W# h
    
48.                 v(i,j)=-vmax;
    ; k: P2 H  v: C
49.             end
    + j+ [9 F' j( J/ |( M- L. K; |
50.         end
    7 L& m+ T  J6 R+ `' k, @\" }% Z
51.         x(i,:)=x(i,:)+v(i,:);8 K; o9 U4 a( n5 h8 G\" C
    
52.     end) \6 w2 ~4 b6 R  @- d7 P+ K% |
    
53.     if abs(globalbest_faval)<E0,break,end
    0 ?% [2 x# _/ J* B
54.     k=k+1;
    5 l( b+ f- T- J  `
55. end
    5 U( w\" m\" X7 L  a
56. Value1=1/globalbest_faval-1; Value1=num2str(Value1);
    $ a. t( |7 r0 \  p; X; P6 Z
57. % strcat指令可以实现字符的组合输出+ j\" F8 |# ~9 J: G
    
58. disp(strcat('the maximum value','=',Value1));
    0 B6 J5 P) z! @2 f, ]2 o\" ~
59. %输出最大值所在的横坐标位置
    4 p& m; u$ e1 m1 q
60. Value2=globalbest_x; Value2=num2str(Value2);
    2 K/ M  ]2 k, _. s/ z  K' g( X
61. disp(strcat('the corresponding coordinate','=',Value2));, Y7 n, I9 \& ~* ?6 b  [# A
    
62. x=-5:0.01:5;
    \" m. p0 g, Z  q7 r$ ]1 r; J
63. y=2.1*(1-x+2*x.^2).*exp(-x.^2/2);
    4 n( a8 _& A6 H2 I
64. plot(x,y,'m-','linewidth',3);0 r. m% u8 Q' {) J/ \6 D
    
65. hold on;* ]4 ?% {: q8 w6 m* V& [
    
66. plot(globalbest_x,1/globalbest_faval-1,'kp','linewidth',4);0 g; [# e5 P  U. j4 X7 c, M2 C
    
67. legend('目标函数','搜索到的最大值');xlabel('x');ylabel('y');grid on;toc;2 I; n; O' K! Q5 V5 X$ w
    

复制代码

平凡之不凡

这个是调好的

永不漠灭的轨迹-

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

谢谢老师！