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

永不漠灭的轨迹-

>> function main()
clc;clear all;close all;
6 W: a3 b  l9 s/ Y9 F- Dtic;                 %程序运行计时
6 k# A+ u$ ?5 w; J. ?E0=0.001;            %允许误差
9 Z; H3 ]- Y* _' d. a6 tMaxNum=100;           %粒子最大迭代次数
narvs=1;            %目标函数自变量个数
  @% D: d4 t) `8 nparticlesize=30;    %粒子群规模
3 a4 g0 n" D  ^" R- ~7 dc1=2;               %每个粒子个体学习因子，也称加速常数
c2=2;               %每个粒子的社会学习因子，也称加速常数
w=0.6;              %惯性因子
  a+ z, K9 X, M6 Pvmax=0.8;           %粒子最大飞翔速度
* M" s6 H7 i% W! V: K2 J9 V$ ^x=-5+10*rand(particlesize,narvs); %粒子所在位置
v=2*rand(particlesize,navrs);        %粒子的飞翔速度
% V% F  W; ]7 H! ?5 m& F* R: s%用inline定义适应度函数以便于将子函数与主程序文件放在一起，
5 |8 s! A4 a' ~. k+ K0 q3 b' o* k%目标函数是：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定义的实验度函数会使程序运行速度大大降低
' B/ K0 g. _4 Jfor i=1:particlesize
    for j=1:narvs
        f(i)=fitness(x(i,j));
    end
end
personalbest_x=x;
+ ~7 H5 |2 D; I  h7 d+ R2 g8 Y4 _, L5 Gpersonalbest_faval=f;
[globalbest_faval i]=min(personalbest_faval);
globalbest_x=personalbest_x(i,;
k=1;
while k<=MaxNum
    for i=1:particlesize
) {) t1 Q, I6 d6 c7 u2 C5 e        for j=1:narvs;
        f(i)=fitness(x(i,j));
1 P7 J( z6 l/ z# c0 V' r9 ^        end
    if f(i)<personalbest_faval(i)   %判断当前位置是否历史上最佳位置
- G" q- A4 ^$ S2 j/ S, x        personalbest_faval(i)=f(i);
        personalbest_x(i,=x(i,;
9 i/ e1 B% E$ X6 n    end
    end
+ [' e9 ^# k0 v: w  |: H4 z+ T    [globalbest_faval i]=min(personalbest_faval);
    globalbest_x=personalbest_x(i,;
- `0 ^* |4 o# A* C2 K6 V3 P    for i=1:particlesize           %更新粒子群里每个个体的最新位置
        v(i,=w*v(i,+c1*rand*(personalbest_x(i,-x(i,)...
6 y9 t9 V! {# @& T- d0 o5 Y            +c2*rand*(globalbest_x-x(i,);
         for j=1:narvs             %判断粒子的飞翔速度是否超过了最大飞行速度
. H  ]8 ?/ Z2 k             if v(i.j)>vmax;   
  l# x) C; w6 ^9 \4 L' ^                 v(i,j)=vmax;
% p) ~. t  @. e  y0 ~; D9 w1 m             else if v(i,j)<-vmax;
                     v(i,j)=-vmax;
                 end
4 b: ]  V+ L' e9 `- G             end
             x(i,=x(i,+v(i,;
         end
/ g+ n8 f! f9 s! v% z3 D: K7 v         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);
+ o8 A% C3 K) B4 I' W: |     disp(strcat('the corresponding coordinate','=',Value2));
4 g. V" Y. h; K% f9 y2 b     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;
/ z$ h) T, f1 z2 Z( C* Y     plot(globalbest_x,1/globalbest_faval-1,'kp','linewidth',4);
     legend('目标函数','搜索到的最大值');
     xlabel('x');ylabel('y');grid on;toc;

, s! l/ f" G  [( g
错误提示 ??? function main()    |
8 b" T# I" d- MError: Function definitions are not permitted in this context.

1 j4 F) l- g& t# D5 ^: L% s: u1 F去掉function main()后 提示??? Undefined function or variable 'navrs'.错误。求指点！
7 @) i8 y' G6 y

平凡之不凡

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

永不漠灭的轨迹-

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

永不漠灭的轨迹-

平凡之不凡 发表于 2014-8-14 09:37
$ Q9 s+ i, e+ P+ H9 o代码重新发一下，格式出问题了，你可以调整一下

7 O5 @7 a2 _; s: E3 S! m/ w9 v* K格式是对的  是不是那个自定义函数的问题，我刚刚查了一下，说是要先定义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()+ G% ^+ P: g) B2 s+ P; o
   
2. clc;clear all;close all;
   / o5 |. }2 N) _4 X\" j# A
3. tic;                              %程序运行计时* W. h1 S' a# Z
   
4. E0=0.001;                        %允许误差
   + {) m) r4 I3 d- b
5. MaxNum=100;                    %粒子最大迭代次数
   8 G1 S  o7 R. x6 c
6. narvs=1;                         %目标函数的自变量个数) L* {2 e9 d' U6 A
   
7. particlesize=30;                    %粒子群规模
   - [5 O9 X3 |% _\" I
8. c1=2;                            %每个粒子的个体学习因子，也称为加速常数
   * L- r7 [4 L) X2 P# \
9. c2=2;                            %每个粒子的社会学习因子，也称为加速常数
   1 O5 E! g  s  k+ {% j' s% d
10. w=0.6;                           %惯性因子  \  J* i6 a! z0 Y, {4 Z
    
11. vmax=0.8;                        %粒子的最大飞翔速度
    2 h\" C* k* [\" x) y0 {+ S
12. x=-5+10*rand(particlesize,narvs);     %粒子所在的位置
      h  z1 |4 x( Q% t! I7 n
13. v=2*rand(particlesize,narvs);         %粒子的飞翔速度
    : u! s/ E3 F8 w' B5 h4 S( p
14. %用inline定义适应度函数以便将子函数文件与主程序文件放在一起，, T% E5 ^, _8 J+ k0 s1 Z/ L* \/ L. ^- L
    
15. %目标函数是：y=1+(2.1*(1-x+2*x.^2).*exp(-x.^2/2))3 K\" G3 d. W! ?: J
    
16. %inline命令定义适应度函数如下：
    ; B1 B% ^7 \. Z  y
17. fitness=inline('1/(1+(2.1*(1-x+2*x.^2).*exp(-x.^2/2)))','x');, J) a% b- w\" Z/ ^$ ^
    
18. %inline定义的适应度函数会使程序运行速度大大降低
    4 ^/ O4 c$ g2 y+ d7 N8 z8 f
19. for i=1:particlesize/ i& i5 ]. Z  T: q* f
    
20.     for j=1:narvs& f7 [) `) W1 X2 N2 L4 g0 x; q\" U
    
21.         f(i)=fitness(x(i,j));, r) ~. J. P, j9 Q7 a/ E
    
22.     end& v9 E\" L/ F' r0 j0 Z
    
23. end\" }% J2 W* \5 D
    
24. personalbest_x=x;
    1 k\" h0 @1 S8 K2 V. J' d
25. personalbest_faval=f;# a6 f1 @; i$ k
    
26. [globalbest_faval i]=min(personalbest_faval);
    - a  s& z. r( ~$ ?2 l+ H8 c
27. globalbest_x=personalbest_x(i,:);1 q; c8 T8 _% Y
    
28. k=1;$ _7 h# @+ H  C
    
29. while k<=MaxNum2 U. F8 Z9 [  k7 I* `: ~
    
30.     for i=1:particlesize
      G) v* u6 Z/ i: n
31.         for j=1:narvs2 ~2 f; O5 \* J$ U+ ]6 _
    
32.             f(i)=fitness(x(i,j));
    $ D2 H, t! z( G0 ?9 M) K
33.         end
    ; k2 I3 N( p% _9 O# d
34.         if f(i)<personalbest_faval(i) %判断当前位置是否是历史上最佳位置4 P% F* N$ Z& @$ X
    
35.             personalbest_faval(i)=f(i);8 Y! d4 V0 n$ g2 u& D
    
36.             personalbest_x(i,:)=x(i,:);& b% U) l2 H- k
    
37.         end
    + B4 X+ P* i. u
38.     end
    5 T/ d( e  [) |! E: e
39.     [globalbest_faval i]=min(personalbest_faval);
    4 K' y7 Q0 P3 r/ D- D! H
40.     globalbest_x=personalbest_x(i,:);+ u- K9 T4 y9 `2 }- ]' l; Z
    
41.     for i=1:particlesize %更新粒子群里每个个体的最新位置( E( B& u+ ]. U/ b* v\" @9 a
    
42.         v(i,:)=w*v(i,:)+c1*rand*(personalbest_x(i,:)-x(i,:))...
    ; y8 s# w( Z7 t
43.             +c2*rand*(globalbest_x-x(i,:));
    4 J1 X7 ^. ^/ |
44.         for j=1:narvs    %判断粒子的飞翔速度是否超过了最大飞翔速度* Y+ N+ D3 ?% [. ~0 v2 U
    
45.             if v(i,j)>vmax;
    ) ^- A0 b6 q% N+ S8 n# R4 k
46.                 v(i,j)=vmax;3 r8 v# D\" H9 x
    
47.             elseif v(i,j)<-vmax;4 `9 p4 y1 l/ n: n5 A' f- v/ D
    
48.                 v(i,j)=-vmax;* P! v# E9 D) w& _; ~, G
    
49.             end5 z$ ]4 ~0 r( s
    
50.         end2 n4 I' ^3 r) m, S( r6 c+ D9 C
    
51.         x(i,:)=x(i,:)+v(i,:);
    6 a2 ~\" Q* g/ `. ]5 J/ p
52.     end2 e* d6 o+ Z; G6 N8 d  ]
    
53.     if abs(globalbest_faval)<E0,break,end' m- T\" O1 W/ d, }) ^# G
    
54.     k=k+1;6 B/ V0 i0 f7 f3 t* Q( e
    
55. end/ a; [+ p0 {9 o  b# u2 l7 _
    
56. Value1=1/globalbest_faval-1; Value1=num2str(Value1);
    \" i\" m* [. ?% V0 }: M7 K: X: @1 b
57. % strcat指令可以实现字符的组合输出3 k* A, |% H/ d\" ?4 F
    
58. disp(strcat('the maximum value','=',Value1));3 O. @\" \+ w; y# a8 Y
    
59. %输出最大值所在的横坐标位置, B# h& D) \% i0 f
    
60. Value2=globalbest_x; Value2=num2str(Value2);) q% Q$ }& j  {8 e3 o9 R2 j. K
    
61. disp(strcat('the corresponding coordinate','=',Value2));
    0 y* B\" I  t$ D6 K! q
62. x=-5:0.01:5;
    3 R; e; |\" d( b
63. y=2.1*(1-x+2*x.^2).*exp(-x.^2/2);
    + l1 [2 K% J5 u% J
64. plot(x,y,'m-','linewidth',3);# P, n; }6 P) i% o0 [, E* Z
    
65. hold on;( d, J, D% [7 k2 R
    
66. plot(globalbest_x,1/globalbest_faval-1,'kp','linewidth',4);
    6 }1 |/ C& q\" P
67. legend('目标函数','搜索到的最大值');xlabel('x');ylabel('y');grid on;toc;) A. ^8 U  R$ S7 Y7 f3 l: [
    

复制代码

平凡之不凡

这个是调好的

永不漠灭的轨迹-

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

谢谢老师！