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

永不漠灭的轨迹-

>> function main()
* Q* d- ]5 v( Nclc;clear all;close all;
+ h6 \+ f0 l, _3 l' L1 ftic;                 %程序运行计时
# U( m- E" A9 q! ]1 z% P; wE0=0.001;            %允许误差
MaxNum=100;           %粒子最大迭代次数
narvs=1;            %目标函数自变量个数
particlesize=30;    %粒子群规模
9 R5 B" u/ S" n& \1 Ac1=2;               %每个粒子个体学习因子，也称加速常数
c2=2;               %每个粒子的社会学习因子，也称加速常数
w=0.6;              %惯性因子
. w# T1 `9 e+ n1 L# k& Z3 V2 M7 ^vmax=0.8;           %粒子最大飞翔速度
( p' f6 I9 V! B4 Tx=-5+10*rand(particlesize,narvs); %粒子所在位置
v=2*rand(particlesize,navrs);        %粒子的飞翔速度
% t+ l, c- T- i9 _; [4 \# n6 a%用inline定义适应度函数以便于将子函数与主程序文件放在一起，
! H* P. ~8 y1 n' ]% n" G%目标函数是：y=1+(2.1*(1-x+2*x.^2).*exp(-x.^2/2)）
4 i2 D1 K" B+ Z' C( Q%inline命令定义实应度函数如下
fitness=inline('1/(1+(2.1*(1-x+2*x.^2).*exp(-x.^2/2)))','x');
%inline定义的实验度函数会使程序运行速度大大降低
for i=1:particlesize
2 K7 P' i$ J- U/ ?, J    for j=1:narvs
( m+ a# y3 H, {% R) f9 ~! {        f(i)=fitness(x(i,j));
    end
end
& V0 G1 r5 u6 B) Gpersonalbest_x=x;
personalbest_faval=f;
* k$ U! \8 }% R# q* Z" C. f[globalbest_faval i]=min(personalbest_faval);
/ w8 c0 F  ~1 p7 Q0 Eglobalbest_x=personalbest_x(i,;
k=1;
while k<=MaxNum
    for i=1:particlesize
( W, n7 Q* O, J" ?, x* a) o4 I        for j=1:narvs;
        f(i)=fitness(x(i,j));
        end
    if f(i)<personalbest_faval(i)   %判断当前位置是否历史上最佳位置
. O+ s/ M7 F' B: I, O3 p8 m: ?; ]4 c        personalbest_faval(i)=f(i);
* R& M6 g1 Y6 u9 C( {) V  s        personalbest_x(i,=x(i,;
    end
    end
: Y; u/ m$ z. O    [globalbest_faval i]=min(personalbest_faval);
9 f: U5 U) w9 N# {4 \5 L) l1 |* ]* `    globalbest_x=personalbest_x(i,;
% [% q5 C* i4 U    for i=1:particlesize           %更新粒子群里每个个体的最新位置
: Z" B/ \1 T5 E6 @6 a+ n        v(i,=w*v(i,+c1*rand*(personalbest_x(i,-x(i,)...
            +c2*rand*(globalbest_x-x(i,);
         for j=1:narvs             %判断粒子的飞翔速度是否超过了最大飞行速度
, s- u- l2 ~5 ]/ C# W3 o0 J* ~             if v(i.j)>vmax;   
                 v(i,j)=vmax;
             else if v(i,j)<-vmax;
                     v(i,j)=-vmax;
8 i; i* A4 a9 t: B                 end
9 m  c: ~# u  r             end
             x(i,=x(i,+v(i,;
6 l) l% m- B% m& A3 w3 Q! `) s' W         end
         if abs(globalbest_faval)<E0,break,end
5 T$ g+ ?% [, W* H: D         k=k+1;
    end
/ M( J) Z8 Y0 ]: j' ^/ ~    Valuel=1/globalbest_faval-1;Valuel=num2str(Valuel);
0 S) l- F1 v" b/ O    %strcat指令可以实现字符的组合输出
    disp(strcat('the maximun valuel','=',Valuel));
4 S0 t* h6 j9 S- `$ w/ y  x  %输出最大值所在的横坐标位置
    Value2=globalbest_x; Value2=num2str(Value2);
3 M) y; x! |0 O! R0 S/ E2 M     disp(strcat('the corresponding coordinate','=',Value2));
) _- O" N2 N/ r; c) A& [     x=-5:0.01:5;
! J1 h9 D4 g3 H7 r3 t: B+ j     y=2.1*(1-x+2*x.^2).*exp(-x.^2/2);
3 t- U( ?% e; W( Q     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;


错误提示 ??? function main()    |
Error: Function definitions are not permitted in this context.

去掉function main()后 提示??? Undefined function or variable 'navrs'.错误。求指点！
6 C" X8 z3 s* D" C

平凡之不凡

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

永不漠灭的轨迹-

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

永不漠灭的轨迹-

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

7 M- @* Y3 a* T, j# ^
格式是对的  是不是那个自定义函数的问题，我刚刚查了一下，说是要先定义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()/ [  t' c; r2 M! T9 g) ^2 w& R
   
2. clc;clear all;close all;
   . b0 }; O3 w5 I/ t1 |/ o* p
3. tic;                              %程序运行计时
     q, l3 ?5 i6 n$ p, v0 n3 j, M
4. E0=0.001;                        %允许误差7 o2 G5 H; o! H
   
5. MaxNum=100;                    %粒子最大迭代次数
   . L! z; A0 q2 A3 ^7 h
6. narvs=1;                         %目标函数的自变量个数2 D) t! U' Y5 I
   
7. particlesize=30;                    %粒子群规模; `; Z8 o' h1 u) U% {! l% q
   
8. c1=2;                            %每个粒子的个体学习因子，也称为加速常数( H7 Z1 a. n, J8 s2 ?2 b7 F
   
9. c2=2;                            %每个粒子的社会学习因子，也称为加速常数; u. @2 T+ V- f1 r
   
10. w=0.6;                           %惯性因子% y/ \\" E4 U) q2 t, ~$ a6 ?! \
    
11. vmax=0.8;                        %粒子的最大飞翔速度
    1 H( Q\" q! O+ A& S& N9 K5 ?
12. x=-5+10*rand(particlesize,narvs);     %粒子所在的位置8 }: n0 _2 y) u) O' {5 h% h
    
13. v=2*rand(particlesize,narvs);         %粒子的飞翔速度
    : @2 s! Z: i6 f' b8 }' G\" A4 K\" _
14. %用inline定义适应度函数以便将子函数文件与主程序文件放在一起，
    7 r6 K: S: @$ y& c# A
15. %目标函数是：y=1+(2.1*(1-x+2*x.^2).*exp(-x.^2/2))
    9 f8 Z; C! i* f9 u
16. %inline命令定义适应度函数如下：
      A9 F/ G( ^- ~( ~+ X  z
17. fitness=inline('1/(1+(2.1*(1-x+2*x.^2).*exp(-x.^2/2)))','x');
    , [7 }& l2 ^3 `
18. %inline定义的适应度函数会使程序运行速度大大降低
    ; k) K7 [2 ^& c3 K! p\" S
19. for i=1:particlesize
    - F% O4 @' s2 n
20.     for j=1:narvs
    $ m! n0 x, A. x\" _3 b2 T
21.         f(i)=fitness(x(i,j));\" n1 p+ [$ {3 ?
    
22.     end
    ; F$ k/ L( c8 U* [\" l! P, w; l+ {3 h
23. end% z' C- J5 S5 v9 Q$ x5 Z+ F
    
24. personalbest_x=x;
    ) D  Q# W1 Y! L) e
25. personalbest_faval=f;  {; Q( d1 S3 G8 C( g+ r- P
    
26. [globalbest_faval i]=min(personalbest_faval);
    6 G/ F, [% u# _0 @' g' ]: T  ^
27. globalbest_x=personalbest_x(i,:);
    9 j2 R' N4 V4 O$ R) J9 J
28. k=1;9 E5 r1 ]) n! c( g  t
    
29. while k<=MaxNum
    9 i2 k! ]; g% k8 _6 U5 @2 b) s/ [
30.     for i=1:particlesize' z. |+ W) G1 U  ~7 ~8 F
    
31.         for j=1:narvs\" L$ z+ n; O/ d* z: R* P
    
32.             f(i)=fitness(x(i,j));
    - Q' t+ P) ~2 r% }1 [# w5 [  ]& r
33.         end
      _! `. Q4 d- m5 s$ q3 A
34.         if f(i)<personalbest_faval(i) %判断当前位置是否是历史上最佳位置7 I- b- m& o\" h5 v( ~
    
35.             personalbest_faval(i)=f(i);- v\" P& e\" V- {# q
    
36.             personalbest_x(i,:)=x(i,:);
    7 d6 }/ e( W7 B: n  j) N
37.         end3 w  o; R- K- h  x) p- G
    
38.     end  Q  e% @; ]' e( `. ^7 X; C, R3 d
    
39.     [globalbest_faval i]=min(personalbest_faval);
    2 B) o# ]3 |, F0 D) Q: b
40.     globalbest_x=personalbest_x(i,:);
    ! o* P7 J) J& w0 K2 P  t& ~
41.     for i=1:particlesize %更新粒子群里每个个体的最新位置
    / ?2 e6 i# L' a% M* {
42.         v(i,:)=w*v(i,:)+c1*rand*(personalbest_x(i,:)-x(i,:))...
    / X% Q9 R; A$ h: T
43.             +c2*rand*(globalbest_x-x(i,:));
    5 I! F. _! i0 G! W: y
44.         for j=1:narvs    %判断粒子的飞翔速度是否超过了最大飞翔速度
    + J2 e  u/ f$ f
45.             if v(i,j)>vmax;
    / t- D- P; }\" s6 x
46.                 v(i,j)=vmax;3 ?( }# E4 V* K# H1 K, A) F2 c
    
47.             elseif v(i,j)<-vmax;
    0 f. ?' j' ]- q( }' I. C\" p
48.                 v(i,j)=-vmax;
    \" R5 L! j, @  n9 [
49.             end
    4 X, d6 P, }7 y  c! P7 M
50.         end0 {& o. S( b) [
    
51.         x(i,:)=x(i,:)+v(i,:);
    , H$ v; X+ T8 E; B
52.     end
    . m/ B' x0 w2 c( `: x0 v
53.     if abs(globalbest_faval)<E0,break,end
    / {# Z8 ~% G- F4 Q! _% V
54.     k=k+1;: ]* C4 |* }% N4 }5 \) Z
    
55. end
    9 _$ e# z! ~4 M. B# J' H
56. Value1=1/globalbest_faval-1; Value1=num2str(Value1);
    ( A/ K& L* L; U\" R4 S
57. % strcat指令可以实现字符的组合输出
    9 ]) \5 l/ M6 j
58. disp(strcat('the maximum value','=',Value1));
    % H5 Z! n5 N; E; }7 p  G3 ?% E
59. %输出最大值所在的横坐标位置
    / F4 S# ?; a$ o9 J/ l) x
60. Value2=globalbest_x; Value2=num2str(Value2);/ {& ^2 u: x: T/ W# u
    
61. disp(strcat('the corresponding coordinate','=',Value2));! E1 k2 x8 Q. e; S: u; _; m
    
62. x=-5:0.01:5;
    : B( y1 s' i\" V0 ]# G7 R6 N0 Y
63. y=2.1*(1-x+2*x.^2).*exp(-x.^2/2);- r# N\" V- D( K6 J+ A\" ]
    
64. plot(x,y,'m-','linewidth',3);
    ( w8 m( p\" e! N. a# M
65. hold on;# D/ F1 x\" u/ B# S. T( N
    
66. plot(globalbest_x,1/globalbest_faval-1,'kp','linewidth',4);9 W! g7 F7 U3 r3 M& Q: J6 L
    
67. legend('目标函数','搜索到的最大值');xlabel('x');ylabel('y');grid on;toc;
    3 x0 o' {1 r# ^4 |5 O# t

复制代码

平凡之不凡

这个是调好的

永不漠灭的轨迹-

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

6 i0 d  q" P# n6 e7 o0 U) {谢谢老师！