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

永不漠灭的轨迹-

>> function main()
3 h2 H  h7 R" N7 X+ W, M- yclc;clear all;close all;
tic;                 %程序运行计时
8 o, R. v% v6 k+ ]  \: E/ Z4 y# oE0=0.001;            %允许误差
MaxNum=100;           %粒子最大迭代次数
narvs=1;            %目标函数自变量个数
3 d+ [; c# f+ `* a8 P: tparticlesize=30;    %粒子群规模
: K$ p; Y3 B2 E0 s' a) |c1=2;               %每个粒子个体学习因子，也称加速常数
9 B0 ]& @0 k0 I# R. Cc2=2;               %每个粒子的社会学习因子，也称加速常数
w=0.6;              %惯性因子
2 C& e' s8 p, {. F1 ^$ O! yvmax=0.8;           %粒子最大飞翔速度
x=-5+10*rand(particlesize,narvs); %粒子所在位置
' X5 ]; w+ s2 \) I+ Hv=2*rand(particlesize,navrs);        %粒子的飞翔速度
. j  g9 ]& p3 c) w%用inline定义适应度函数以便于将子函数与主程序文件放在一起，
%目标函数是：y=1+(2.1*(1-x+2*x.^2).*exp(-x.^2/2)）
, y7 I0 h# V6 c/ w- N) N%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));
: a/ F2 p& r/ g' F0 _    end
end
personalbest_x=x;
personalbest_faval=f;
) }/ b+ Y9 t" ]% o+ G0 ][globalbest_faval i]=min(personalbest_faval);
( B8 B/ }2 M# _' Aglobalbest_x=personalbest_x(i,;
k=1;
4 p9 p9 x1 q6 c! [/ o, Z/ cwhile k<=MaxNum
0 A( X5 i* G, ^$ S- y. ~4 l5 z    for i=1:particlesize
- p% Z$ n" h" `  o4 T' a        for j=1:narvs;
        f(i)=fitness(x(i,j));
7 q5 K1 W( r; g! Y1 d3 y) o& X        end
0 W0 p  I$ u5 n% m" b    if f(i)<personalbest_faval(i)   %判断当前位置是否历史上最佳位置
" u0 W2 Y5 Z( j0 j% V# \8 A* v        personalbest_faval(i)=f(i);
7 E1 U1 z* B; d  S        personalbest_x(i,=x(i,;
    end
    end
. G6 }  K  U; {( D1 w: }    [globalbest_faval i]=min(personalbest_faval);
0 f3 D# `5 f6 K$ F9 F! z, G    globalbest_x=personalbest_x(i,;
( v3 o9 j5 b2 b7 l: P; S    for i=1:particlesize           %更新粒子群里每个个体的最新位置
        v(i,=w*v(i,+c1*rand*(personalbest_x(i,-x(i,)...
6 G! b  c' v& A7 U+ U& M5 S5 q, h' b            +c2*rand*(globalbest_x-x(i,);
7 w1 Q1 O/ [, x# Z4 T0 F" Y         for j=1:narvs             %判断粒子的飞翔速度是否超过了最大飞行速度
             if v(i.j)>vmax;   
4 L  r5 Z* G# N                 v(i,j)=vmax;
% S: \0 [4 V9 e' W9 B             else if v(i,j)<-vmax;
- l1 M$ }% M$ y7 m& ?# ~                     v(i,j)=-vmax;
                 end
             end
0 ~5 B$ w$ d; ]5 w( L) v$ ]& J             x(i,=x(i,+v(i,;
         end
. \" z6 C( ]2 u1 u8 S: L& h         if abs(globalbest_faval)<E0,break,end
         k=k+1;
    end
    Valuel=1/globalbest_faval-1;Valuel=num2str(Valuel);
    %strcat指令可以实现字符的组合输出
1 _! Q; A- u/ f: U$ r$ k    disp(strcat('the maximun valuel','=',Valuel));
  %输出最大值所在的横坐标位置
    Value2=globalbest_x; Value2=num2str(Value2);
     disp(strcat('the corresponding coordinate','=',Value2));
2 o3 i( P0 @6 U7 N     x=-5:0.01:5;
6 r5 s& T9 p9 ~6 {! V5 p; A! K4 n     y=2.1*(1-x+2*x.^2).*exp(-x.^2/2);
     plot(x,y,'m-','linewidth',3);
6 I4 Y  h0 \$ P! z& i0 g# S: z& K     hold on;
4 `. X8 h! n, _( z3 D     plot(globalbest_x,1/globalbest_faval-1,'kp','linewidth',4);
     legend('目标函数','搜索到的最大值');
     xlabel('x');ylabel('y');grid on;toc;
+ ]  S% ]% y3 x
, d9 l4 ]. m2 h' }' i6 ?% z
错误提示 ??? function main()    |
  o9 m% [! b. N3 ^3 W) T* jError: Function definitions are not permitted in this context.

去掉function main()后 提示??? Undefined function or variable 'navrs'.错误。求指点！
9 [8 Y8 I( E. j/ H3 ?

平凡之不凡

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

永不漠灭的轨迹-

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

永不漠灭的轨迹-

平凡之不凡 发表于 2014-8-14 09:37
7 i, s2 E* ^' {$ U. k  J4 W代码重新发一下，格式出问题了，你可以调整一下

5 d' V7 R5 g* q
" l! S: G: u3 Z  _8 C, ~, u格式是对的  是不是那个自定义函数的问题，我刚刚查了一下，说是要先定义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()
   . \/ V0 i5 ~9 h
2. clc;clear all;close all;
   $ u0 V$ \. D; v8 W2 c' D
3. tic;                              %程序运行计时$ |8 \4 l* T1 T% U: I
   
4. E0=0.001;                        %允许误差
   ' E. W7 b9 z* g) Z/ z% M
5. MaxNum=100;                    %粒子最大迭代次数5 e- e+ M6 C+ f\" x( m! p
   
6. narvs=1;                         %目标函数的自变量个数
   1 A7 S4 ?. D1 z
7. particlesize=30;                    %粒子群规模
   : \7 n/ e2 E3 J7 G) ?  C0 |1 k
8. c1=2;                            %每个粒子的个体学习因子，也称为加速常数
   ! m  u0 a- B9 ]; d# z$ \( N
9. c2=2;                            %每个粒子的社会学习因子，也称为加速常数7 V  X! a8 M, C9 o4 _+ `
   
10. w=0.6;                           %惯性因子: b$ X- F- p8 C\" @9 X5 M; r
    
11. vmax=0.8;                        %粒子的最大飞翔速度) y1 M/ J# ^' z. a# V; g
    
12. x=-5+10*rand(particlesize,narvs);     %粒子所在的位置
    % p1 e+ w4 n# ~+ J5 R9 y6 i
13. v=2*rand(particlesize,narvs);         %粒子的飞翔速度. o1 e% l7 _5 [0 _& Z+ k- F/ q
    
14. %用inline定义适应度函数以便将子函数文件与主程序文件放在一起，) U0 a8 M' G) n# v5 z. X9 E
    
15. %目标函数是：y=1+(2.1*(1-x+2*x.^2).*exp(-x.^2/2)): h# R$ e( ?  w\" m# ?8 R
    
16. %inline命令定义适应度函数如下：' Z  t, r$ q0 h$ W
    
17. fitness=inline('1/(1+(2.1*(1-x+2*x.^2).*exp(-x.^2/2)))','x');5 O* j# S7 g2 l7 f
    
18. %inline定义的适应度函数会使程序运行速度大大降低
    1 O( `: A9 p& h- {% b3 t
19. for i=1:particlesize
    8 z0 t\" X; f; W0 s( E# c: ~9 J
20.     for j=1:narvs/ N7 r$ K7 R' ]* `3 ?5 S  a
    
21.         f(i)=fitness(x(i,j));, W9 d. b* t* G3 q* Y* D( R2 h
    
22.     end
    . C7 `$ ~8 E\" O\" w8 D4 K/ x
23. end% N1 w; x, ]9 b5 E# m4 k: i2 r: }
    
24. personalbest_x=x;
    - y# H+ {+ A9 ]: l2 c& ?
25. personalbest_faval=f;
    # E# i$ n$ w4 z' |' W) j
26. [globalbest_faval i]=min(personalbest_faval);
    4 V! W4 }7 V5 }  h& @' g& \\" q! }
27. globalbest_x=personalbest_x(i,:);
    * p- `  ^4 D# c& ~+ d- ]
28. k=1;
    ) Q4 S7 Y5 y3 E: h
29. while k<=MaxNum- S+ @' t& C9 b7 i1 S4 r( C: s
    
30.     for i=1:particlesize
    - j0 y5 @$ w: _
31.         for j=1:narvs. ]' @* |- `, j0 o
    
32.             f(i)=fitness(x(i,j));
    1 [1 v6 T8 Z# r, I
33.         end
    ; ]3 s* r8 N\" A! z
34.         if f(i)<personalbest_faval(i) %判断当前位置是否是历史上最佳位置  B0 c' v1 i' L% m& `0 |* L# g
    
35.             personalbest_faval(i)=f(i);
    2 t  J6 E  }7 a
36.             personalbest_x(i,:)=x(i,:);: h! s+ r6 P# Q. P- G
    
37.         end2 @' E2 @! K0 \* b* V( @2 E3 \
    
38.     end& {* A8 ~0 l: {, f) r\" Z  _
    
39.     [globalbest_faval i]=min(personalbest_faval);9 A& S8 `$ t2 r9 n
    
40.     globalbest_x=personalbest_x(i,:);
    * f: b, f# Y* S: U6 k; ?3 C0 G
41.     for i=1:particlesize %更新粒子群里每个个体的最新位置8 P7 k( b* z  j( c! G% N4 ^7 g# V
    
42.         v(i,:)=w*v(i,:)+c1*rand*(personalbest_x(i,:)-x(i,:))...: _# Q, c& w5 I7 q+ O; J
    
43.             +c2*rand*(globalbest_x-x(i,:));
    ; G; x  ]. V+ N$ X
44.         for j=1:narvs    %判断粒子的飞翔速度是否超过了最大飞翔速度
    1 ^% `. P) x- ~1 z2 x
45.             if v(i,j)>vmax;7 K# V+ A0 w\" z2 y, c$ r$ [
    
46.                 v(i,j)=vmax;5 w% h: z# H8 s8 l0 ^/ m
    
47.             elseif v(i,j)<-vmax;
    # y' ?9 X8 D% J. R% A
48.                 v(i,j)=-vmax;
    - O# k8 e7 Z2 a' b% i0 P9 g( T
49.             end
    3 r3 d0 c  s9 @2 f' r  P9 F\" h
50.         end
    - J0 v  J- y% ~\" L
51.         x(i,:)=x(i,:)+v(i,:);/ _) \9 v4 j8 o
    
52.     end
    8 i* D* ]4 W2 E+ r% k8 X
53.     if abs(globalbest_faval)<E0,break,end
    0 m* P: ~0 W- C! ~
54.     k=k+1;
    ; @$ e$ k; _8 {( x# @2 t/ x0 z
55. end  K* _/ K/ a/ ~9 U
    
56. Value1=1/globalbest_faval-1; Value1=num2str(Value1);1 ^7 F3 v+ H) y6 \\" u. X
    
57. % strcat指令可以实现字符的组合输出/ r: I/ c% j) n8 N- \3 A3 k  ^
    
58. disp(strcat('the maximum value','=',Value1));
    6 _+ \6 X1 h( N\" b
59. %输出最大值所在的横坐标位置3 R. q6 A( G4 u1 f5 ?
    
60. Value2=globalbest_x; Value2=num2str(Value2);
    4 b4 e- R; @5 Q$ \\" D6 d+ ]
61. disp(strcat('the corresponding coordinate','=',Value2));
    * P; V7 X9 q9 r4 R' I; R
62. x=-5:0.01:5;5 A8 h) X  O+ ^\" w3 R0 K
    
63. y=2.1*(1-x+2*x.^2).*exp(-x.^2/2);
    9 u. ?% ~. w$ m* u/ }8 x
64. plot(x,y,'m-','linewidth',3);1 e+ c. g. w! i/ `0 l1 J
    
65. hold on;8 U: s0 w$ o3 M( V
    
66. plot(globalbest_x,1/globalbest_faval-1,'kp','linewidth',4);
    * r\" A7 n7 J7 C2 t& M/ A7 s
67. legend('目标函数','搜索到的最大值');xlabel('x');ylabel('y');grid on;toc;. ]7 ~! p5 Q$ c$ k0 R0 }7 P+ {  [
    

复制代码

平凡之不凡

这个是调好的

永不漠灭的轨迹-

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

# y+ X1 x! ?- g7 t/ O& X谢谢老师！