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

永不漠灭的轨迹-

>> function main()
clc;clear all;close all;
0 G: l: [. x4 W- D: |7 H7 r' Q, Gtic;                 %程序运行计时
E0=0.001;            %允许误差
5 W' g% R/ `9 UMaxNum=100;           %粒子最大迭代次数
narvs=1;            %目标函数自变量个数
7 Z- y% Q: v# B: E4 s+ I  Dparticlesize=30;    %粒子群规模
& e3 `- [' G1 i, O8 I' ?3 c2 hc1=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)）
: c8 `6 ~5 N8 e! I5 T2 C6 X%inline命令定义实应度函数如下
fitness=inline('1/(1+(2.1*(1-x+2*x.^2).*exp(-x.^2/2)))','x');
%inline定义的实验度函数会使程序运行速度大大降低
for i=1:particlesize
# T# @8 U% U9 N. w% j: l$ p0 T) S    for j=1:narvs
. c5 D4 S# u& v# E        f(i)=fitness(x(i,j));
    end
end
personalbest_x=x;
, B1 N: E/ k  [9 q& e& r/ K* S5 ]personalbest_faval=f;
[globalbest_faval i]=min(personalbest_faval);
/ R0 Y3 C3 n1 b2 l$ r) }globalbest_x=personalbest_x(i,;
k=1;
. [6 U9 C# V1 v% Wwhile k<=MaxNum
! Q4 b9 H/ R4 p    for i=1:particlesize
2 W8 M9 {. _! B* X( `' `' k5 O        for j=1:narvs;
        f(i)=fitness(x(i,j));
        end
    if f(i)<personalbest_faval(i)   %判断当前位置是否历史上最佳位置
        personalbest_faval(i)=f(i);
6 _6 ^0 P0 ^4 Q1 n9 i+ T6 j' ^        personalbest_x(i,=x(i,;
    end
& d' f# I. E7 f: F; J6 b    end
+ d9 |% h9 L) q    [globalbest_faval i]=min(personalbest_faval);
" l" k% ], n# Q! @/ ~0 j, Q    globalbest_x=personalbest_x(i,;
1 C3 t+ X9 \, A    for i=1:particlesize           %更新粒子群里每个个体的最新位置
+ l  k' Q6 r" E" C7 h        v(i,=w*v(i,+c1*rand*(personalbest_x(i,-x(i,)...
* Z& C& g% J( E$ F9 g5 ^6 \. t            +c2*rand*(globalbest_x-x(i,);
% W" j- v2 t: m8 T& s, ^         for j=1:narvs             %判断粒子的飞翔速度是否超过了最大飞行速度
) T" Y/ j: p( o3 }/ u             if v(i.j)>vmax;   
                 v(i,j)=vmax;
             else if v(i,j)<-vmax;
                     v(i,j)=-vmax;
- V. e8 X' Z' P  ~/ Y9 ?                 end
             end
# a/ Y( n4 ?. g  J2 E+ M/ W             x(i,=x(i,+v(i,;
         end
4 J) D7 N2 D* G0 G, T- k' y         if abs(globalbest_faval)<E0,break,end
         k=k+1;
! g0 `% W3 z% L' B    end
    Valuel=1/globalbest_faval-1;Valuel=num2str(Valuel);
    %strcat指令可以实现字符的组合输出
    disp(strcat('the maximun valuel','=',Valuel));
  %输出最大值所在的横坐标位置
    Value2=globalbest_x; Value2=num2str(Value2);
3 }7 G. l4 C0 g9 A     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);
5 t; T, ?) `" X/ r4 w# ~8 f: E     hold on;
  J  U4 z2 U5 }     plot(globalbest_x,1/globalbest_faval-1,'kp','linewidth',4);
& `( O8 @3 B3 t# g4 S- T% ~4 G, B$ I& Y     legend('目标函数','搜索到的最大值');
8 t0 q- c, E& p$ _; M     xlabel('x');ylabel('y');grid on;toc;
# C4 d, P7 p  [2 c0 M
9 k* C; D" _8 d6 ?( [( T
错误提示 ??? function main()    |
Error: Function definitions are not permitted in this context.
" i8 f8 I4 ]) n9 V2 I/ L
去掉function main()后 提示??? Undefined function or variable 'navrs'.错误。求指点！

永不漠灭的轨迹-

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

8 B4 t0 {1 H; B( U' \谢谢老师！

平凡之不凡

这个是调好的

平凡之不凡

1. function main()
     o( v7 G7 _8 W4 I* v
2. clc;clear all;close all;
   6 W% v5 k/ I5 n: E\" @\" j8 }
3. tic;                              %程序运行计时0 _9 S3 n2 D\" q/ f
   
4. E0=0.001;                        %允许误差2 I+ G: ?. ~, H+ B
   
5. MaxNum=100;                    %粒子最大迭代次数
   ( t8 s2 x, S4 y& w4 _7 H$ o% G/ L
6. narvs=1;                         %目标函数的自变量个数  g# d# A/ e5 |- r
   
7. particlesize=30;                    %粒子群规模3 W( [7 O5 _; }# A; P7 {' C& r
   
8. c1=2;                            %每个粒子的个体学习因子，也称为加速常数
   3 S* n  V  P* [
9. c2=2;                            %每个粒子的社会学习因子，也称为加速常数1 v6 z. D  \: \, O
   
10. w=0.6;                           %惯性因子
    0 e# n- o/ _. `\" Y
11. vmax=0.8;                        %粒子的最大飞翔速度
    ' ]6 e% V: Y\" D4 U1 ^# b
12. x=-5+10*rand(particlesize,narvs);     %粒子所在的位置& M5 ~2 k; b7 R! W\" S
    
13. v=2*rand(particlesize,narvs);         %粒子的飞翔速度
    : N0 D6 O7 r. i3 A9 E
14. %用inline定义适应度函数以便将子函数文件与主程序文件放在一起，
    3 Z: s7 G1 V$ o1 j9 ?' N
15. %目标函数是：y=1+(2.1*(1-x+2*x.^2).*exp(-x.^2/2))
    3 O\" Z4 n. S4 ?$ i% }
16. %inline命令定义适应度函数如下：
    0 g1 i4 k5 e$ Z+ L5 X, M! d
17. fitness=inline('1/(1+(2.1*(1-x+2*x.^2).*exp(-x.^2/2)))','x');
    # z  r8 I6 k( b+ g6 L. I$ Q
18. %inline定义的适应度函数会使程序运行速度大大降低8 a' T9 Z. z. e; U
    
19. for i=1:particlesize. i( d4 l( V6 B, ]% ~9 J
    
20.     for j=1:narvs+ o5 Z  J6 B- s( m\" y5 w
    
21.         f(i)=fitness(x(i,j));
    + ^  g\" M* A/ M2 O/ m
22.     end
    * V7 N; L, i1 n% [* Y3 T
23. end2 W8 N0 M6 p& E8 P
    
24. personalbest_x=x;8 C0 G- k4 W2 t\" C& {. S
    
25. personalbest_faval=f;  @/ O+ d4 v+ s2 j4 C. K: G( e0 B% @
    
26. [globalbest_faval i]=min(personalbest_faval);
    7 m' I( U. _4 l\" ^\" Y) k  m, S
27. globalbest_x=personalbest_x(i,:);: T& H2 @# J  y# G
    
28. k=1;
    & L& U) g! }9 `$ G0 p- l% o# T
29. while k<=MaxNum
    \" k1 G  y) w: q. n  I
30.     for i=1:particlesize
    . D\" F1 H1 W8 F, ^7 o0 }
31.         for j=1:narvs- ~+ C  q* y1 F
    
32.             f(i)=fitness(x(i,j));
    2 y7 N  D7 j& b4 x. l1 t# q/ k* @' E
33.         end
      M) K* K6 W% a1 b/ {
34.         if f(i)<personalbest_faval(i) %判断当前位置是否是历史上最佳位置
    $ c/ [  F. K1 j1 n, I6 k4 ~  _
35.             personalbest_faval(i)=f(i);
    7 d: z$ J8 W( t7 E: k$ `; j
36.             personalbest_x(i,:)=x(i,:);( w2 A: L$ U& s: D. d
    
37.         end
    $ C+ ?8 C3 P2 _& L6 g- @
38.     end
    5 S; j. q8 S( u4 S/ ~; u
39.     [globalbest_faval i]=min(personalbest_faval);
    8 b8 d' \8 z+ p& V( w\" `# v
40.     globalbest_x=personalbest_x(i,:);
    ' c3 G. Q2 h! G0 ~1 T4 T5 y/ [
41.     for i=1:particlesize %更新粒子群里每个个体的最新位置% v0 {2 W7 B& y\" P6 m+ `) Z: P
    
42.         v(i,:)=w*v(i,:)+c1*rand*(personalbest_x(i,:)-x(i,:))...
    - [9 M% A1 q4 d0 Q% y, o; L5 E
43.             +c2*rand*(globalbest_x-x(i,:));
    ( i2 \+ G4 i/ ]
44.         for j=1:narvs    %判断粒子的飞翔速度是否超过了最大飞翔速度# I2 X6 L& h8 o4 \
    
45.             if v(i,j)>vmax;
    . }( a0 P$ T, v6 A7 x# ?$ G
46.                 v(i,j)=vmax;
      b\" v, N1 s. k# S9 U
47.             elseif v(i,j)<-vmax;
    6 O: l  \) r  X% Y+ [5 t# Q. p
48.                 v(i,j)=-vmax;
    9 i4 |, F5 m( A1 z
49.             end0 K0 _/ R6 h  {+ B
    
50.         end
    : {- w2 F+ }# q' U/ \
51.         x(i,:)=x(i,:)+v(i,:);* b4 Y2 y3 ?8 r% R1 U0 G- p
    
52.     end
    # `# n' w+ A/ j0 u8 Q9 c
53.     if abs(globalbest_faval)<E0,break,end
    - m3 l: g/ P9 d% P# B5 w9 N, E
54.     k=k+1;
    + L5 C! q$ c. L\" `& p
55. end
    % R\" \- E1 `, g
56. Value1=1/globalbest_faval-1; Value1=num2str(Value1);9 `, T& G* J- f: t' y$ s6 h9 V
    
57. % strcat指令可以实现字符的组合输出. ~9 r: u5 j( @) b. j5 _& R! |! z/ m
    
58. disp(strcat('the maximum value','=',Value1));
    7 }0 J$ L$ _+ K  Z
59. %输出最大值所在的横坐标位置8 e( `, G% \/ W) o
    
60. Value2=globalbest_x; Value2=num2str(Value2);+ i( V6 F2 ]* w+ I3 V0 v
    
61. disp(strcat('the corresponding coordinate','=',Value2));
    . h# r9 V5 o6 P2 o3 u
62. x=-5:0.01:5;, R2 P: y% \2 K/ I
    
63. y=2.1*(1-x+2*x.^2).*exp(-x.^2/2);' [: U4 q1 x( F7 ?( }- k
    
64. plot(x,y,'m-','linewidth',3);0 l: g; G- D8 ~! f; @0 O
    
65. hold on;
    & C& O7 {  C9 y9 s. e
66. plot(globalbest_x,1/globalbest_faval-1,'kp','linewidth',4);
    : d  B% R* Y% N3 z
67. legend('目标函数','搜索到的最大值');xlabel('x');ylabel('y');grid on;toc;
    / n0 c; h6 j2 g) ~' E

复制代码

永不漠灭的轨迹-

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;

永不漠灭的轨迹-

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

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

永不漠灭的轨迹-

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

平凡之不凡

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