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

永不漠灭的轨迹-

>> function main()
clc;clear all;close all;
) w! r& S  W$ |, B  _tic;                 %程序运行计时
E0=0.001;            %允许误差
MaxNum=100;           %粒子最大迭代次数
narvs=1;            %目标函数自变量个数
! j+ Q* d% U5 G* S0 N  Y; Sparticlesize=30;    %粒子群规模
9 }& f5 q8 _/ L3 H4 t. I7 wc1=2;               %每个粒子个体学习因子，也称加速常数
c2=2;               %每个粒子的社会学习因子，也称加速常数
0 F7 t% D0 j- k5 K8 Nw=0.6;              %惯性因子
vmax=0.8;           %粒子最大飞翔速度
4 `+ E8 r3 Y% W2 Mx=-5+10*rand(particlesize,narvs); %粒子所在位置
! ^6 x9 b" G! T; w) b8 V: N3 m' gv=2*rand(particlesize,navrs);        %粒子的飞翔速度
%用inline定义适应度函数以便于将子函数与主程序文件放在一起，
9 k( y0 b% j3 y6 U' m7 o%目标函数是：y=1+(2.1*(1-x+2*x.^2).*exp(-x.^2/2)）
8 z4 J8 F" E; p2 o+ e%inline命令定义实应度函数如下
: I! `8 f8 A; ?6 `fitness=inline('1/(1+(2.1*(1-x+2*x.^2).*exp(-x.^2/2)))','x');
: T$ A5 X, P( |- y- j%inline定义的实验度函数会使程序运行速度大大降低
9 e6 L, y6 F) ^6 q. _: [$ v  W, t$ hfor i=1:particlesize
+ S* o) U/ }' Z* k    for j=1:narvs
        f(i)=fitness(x(i,j));
    end
end
personalbest_x=x;
personalbest_faval=f;
/ R* H; \# k6 f; n[globalbest_faval i]=min(personalbest_faval);
globalbest_x=personalbest_x(i,;
k=1;
while k<=MaxNum
; @6 `' q& J# J0 s# }# Q1 }    for i=1:particlesize
        for j=1:narvs;
        f(i)=fitness(x(i,j));
        end
. N. N: {, ~6 |2 K9 z+ m/ C    if f(i)<personalbest_faval(i)   %判断当前位置是否历史上最佳位置
1 h% D! N0 |- n. P& C9 a        personalbest_faval(i)=f(i);
        personalbest_x(i,=x(i,;
    end
0 c4 {4 M* {# |* o4 q) K    end
    [globalbest_faval i]=min(personalbest_faval);
2 S9 {  H4 g$ C5 u    globalbest_x=personalbest_x(i,;
    for i=1:particlesize           %更新粒子群里每个个体的最新位置
% U% e9 s. m6 S3 v) |) m        v(i,=w*v(i,+c1*rand*(personalbest_x(i,-x(i,)...
0 ?# X2 q& b0 q* D            +c2*rand*(globalbest_x-x(i,);
         for j=1:narvs             %判断粒子的飞翔速度是否超过了最大飞行速度
             if v(i.j)>vmax;   
  j2 `. U; d0 a/ L4 t( e                 v(i,j)=vmax;
6 K, S$ K0 h+ R7 R6 l' _5 V- \             else if v(i,j)<-vmax;
                     v(i,j)=-vmax;
                 end
             end
- S  l( t+ W& \( \# J8 @% Y             x(i,=x(i,+v(i,;
         end
         if abs(globalbest_faval)<E0,break,end
         k=k+1;
% f5 G4 C, I3 ?7 r5 i( G5 i3 g) \    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;
  i4 n) \; A1 V7 j3 c/ c5 Q/ Z     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('目标函数','搜索到的最大值');
3 }3 i6 J& A" x7 R     xlabel('x');ylabel('y');grid on;toc;
! v9 g: x. Z7 \0 i3 K
5 d/ l# Q( h. g
错误提示 ??? function main()    |
Error: Function definitions are not permitted in this context.
5 [, d; _' F: z, L$ [; i# ?
$ X& G' S8 ^& z) o' A. w; |4 e去掉function main()后 提示??? Undefined function or variable 'navrs'.错误。求指点！

平凡之不凡

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

永不漠灭的轨迹-

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

永不漠灭的轨迹-

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

  u) k7 n7 e" Q, @  s/ Z; b
+ b* K) S+ \2 @1 U+ u5 r1 o4 d格式是对的  是不是那个自定义函数的问题，我刚刚查了一下，说是要先定义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()) y# G5 s% S7 j/ C- a, g
   
2. clc;clear all;close all;) B7 V5 w9 [- c& M/ L6 ~4 ^
   
3. tic;                              %程序运行计时7 r+ M; U# A/ n: w8 f\" d4 y, m% N\" ^
   
4. E0=0.001;                        %允许误差
     ?5 I\" X) T5 ~( h8 }
5. MaxNum=100;                    %粒子最大迭代次数+ {\" a1 U7 s5 }. a8 x4 H
   
6. narvs=1;                         %目标函数的自变量个数
   * |; C2 C6 v* i2 i; x
7. particlesize=30;                    %粒子群规模
   & M5 m+ I; x6 Q# g
8. c1=2;                            %每个粒子的个体学习因子，也称为加速常数
   2 D' m& ^3 e; j* ~) K
9. c2=2;                            %每个粒子的社会学习因子，也称为加速常数
   * h3 e4 F  O9 h( d# p! S& Z9 G- l# `
10. w=0.6;                           %惯性因子7 z# i$ I% ?5 ]: b\" _
    
11. vmax=0.8;                        %粒子的最大飞翔速度
    * t/ \* m7 u+ o; R5 n) C
12. x=-5+10*rand(particlesize,narvs);     %粒子所在的位置
    . U# x2 ^. m+ b! l4 @. Q
13. v=2*rand(particlesize,narvs);         %粒子的飞翔速度( U' g. W\" h( _) t9 A* I
    
14. %用inline定义适应度函数以便将子函数文件与主程序文件放在一起，
    3 l! l- ~  c7 M\" u
15. %目标函数是：y=1+(2.1*(1-x+2*x.^2).*exp(-x.^2/2))
    ( l+ b% ]' k* ^  M
16. %inline命令定义适应度函数如下：
    3 M4 n  X# ~+ U' K
17. fitness=inline('1/(1+(2.1*(1-x+2*x.^2).*exp(-x.^2/2)))','x');( a! w& r0 y& A1 v, _8 d\" ?3 n
    
18. %inline定义的适应度函数会使程序运行速度大大降低, M/ @5 F0 C1 j- g
    
19. for i=1:particlesize! G( w) W8 C) t( E
    
20.     for j=1:narvs0 B0 n7 D* G( ]+ v3 Q+ C' t8 {
    
21.         f(i)=fitness(x(i,j));  g$ K1 L, e, P' A4 Z$ Z7 Y  G' ]
    
22.     end
    3 m1 D: y3 z5 Z0 l- k3 c
23. end
    3 v; r: i1 g1 b& j4 F5 m9 j
24. personalbest_x=x;6 Y3 j5 p3 t$ U8 C8 {1 x: ?: A. P
    
25. personalbest_faval=f;\" V' Q2 @1 ?9 a2 o3 U# m7 l
    
26. [globalbest_faval i]=min(personalbest_faval);
    3 Y% W% ?) a/ @
27. globalbest_x=personalbest_x(i,:);1 r. k, g3 t9 q1 H9 h. k# {
    
28. k=1;8 l: s! f$ a# G3 D& ?/ m
    
29. while k<=MaxNum
    ; B/ i, C2 [# T, o: d
30.     for i=1:particlesize
    . ~7 B: T/ o/ F0 ^: r
31.         for j=1:narvs
    0 Y$ R+ L* U  L! |$ r\" S
32.             f(i)=fitness(x(i,j));& O  Z* O7 `) C( ^5 T0 N& o% _, ^
    
33.         end9 C\" z\" o, {4 V. T\" P
    
34.         if f(i)<personalbest_faval(i) %判断当前位置是否是历史上最佳位置% T! q, X, T3 R! L( Z  \! |
    
35.             personalbest_faval(i)=f(i);) u, B; N9 p3 a7 o# Q1 J0 p2 d1 t
    
36.             personalbest_x(i,:)=x(i,:);
    4 |, g0 ~2 [0 c4 d$ o% m
37.         end
    3 ~; y/ M0 w! B& S! G2 Y, V
38.     end
    ' U+ D9 b4 O0 F1 ]
39.     [globalbest_faval i]=min(personalbest_faval);
    $ U  S4 B2 B! N
40.     globalbest_x=personalbest_x(i,:);; u: }4 v7 k' w# q# y, F
    
41.     for i=1:particlesize %更新粒子群里每个个体的最新位置
    4 e) |2 r/ u7 h
42.         v(i,:)=w*v(i,:)+c1*rand*(personalbest_x(i,:)-x(i,:))..., Y2 X: e+ `$ I
    
43.             +c2*rand*(globalbest_x-x(i,:));4 O. j) k4 k0 H0 ?
    
44.         for j=1:narvs    %判断粒子的飞翔速度是否超过了最大飞翔速度- a1 _, L4 \* D# q2 ~, M
    
45.             if v(i,j)>vmax;: X+ p+ q5 ^$ r. T* x& f
    
46.                 v(i,j)=vmax;
    ! G- X1 n: r9 z. v! F
47.             elseif v(i,j)<-vmax;
    4 v. E. {\" `9 i9 p' ^; t
48.                 v(i,j)=-vmax;3 |9 g& B* I3 p7 Y
    
49.             end
    ( D7 h5 J$ s% I7 `, h7 j8 u
50.         end\" v% I$ j7 P8 M% @9 e
    
51.         x(i,:)=x(i,:)+v(i,:);
    \" j$ {$ {% {- k2 Q$ O! v/ o
52.     end
    , J* E: H  c3 r; M+ y
53.     if abs(globalbest_faval)<E0,break,end6 I% E2 [' O) r: M% q5 w% L
    
54.     k=k+1;% `7 ~2 R. o' u
    
55. end
    4 h' W3 J5 i7 c& s9 l/ Z
56. Value1=1/globalbest_faval-1; Value1=num2str(Value1);
    4 N- n0 J$ {; \5 d2 ]2 W! w$ J
57. % strcat指令可以实现字符的组合输出1 T7 f& }% N9 k  ^6 J8 \9 W
    
58. disp(strcat('the maximum value','=',Value1));$ [8 O9 ^1 G- ?
    
59. %输出最大值所在的横坐标位置
    % q% j; z. V0 m
60. Value2=globalbest_x; Value2=num2str(Value2);! I7 ]! F, G/ `+ J. Q: y
    
61. disp(strcat('the corresponding coordinate','=',Value2));/ _7 W, R/ ^8 `
    
62. x=-5:0.01:5;
    2 i4 u: i$ |6 m5 @) R/ j7 C4 Y- x
63. y=2.1*(1-x+2*x.^2).*exp(-x.^2/2);8 }2 H- m& l- q: l# x8 K8 P
    
64. plot(x,y,'m-','linewidth',3);
    + C9 [) _3 \9 g2 D
65. hold on;$ \0 H% i; n/ Z' F5 K7 e: }
    
66. plot(globalbest_x,1/globalbest_faval-1,'kp','linewidth',4);8 K. o: G4 f1 f
    
67. legend('目标函数','搜索到的最大值');xlabel('x');ylabel('y');grid on;toc;
    8 I6 T( L2 \& J/ `5 d+ q

复制代码

平凡之不凡

这个是调好的

永不漠灭的轨迹-

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

谢谢老师！