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

永不漠灭的轨迹-

>> function main()
5 I% i, M7 M5 cclc;clear all;close all;
- Q8 ^+ @+ P4 Etic;                 %程序运行计时
E0=0.001;            %允许误差
$ [* t/ C, [9 x# |/ z( F) Y8 p* jMaxNum=100;           %粒子最大迭代次数
$ M0 G9 X  Q" j7 `narvs=1;            %目标函数自变量个数
9 r1 H6 Z. _$ K( h2 i# Qparticlesize=30;    %粒子群规模
" J! x& ], m+ A/ I* nc1=2;               %每个粒子个体学习因子，也称加速常数
; C3 b/ A& W, e1 x7 ic2=2;               %每个粒子的社会学习因子，也称加速常数
- U8 }7 @* w) l- e' D0 C, l  |& \w=0.6;              %惯性因子
2 m4 z$ ~" T1 l7 k2 P. @vmax=0.8;           %粒子最大飞翔速度
x=-5+10*rand(particlesize,narvs); %粒子所在位置
8 F' N, C) X' ~7 K" Fv=2*rand(particlesize,navrs);        %粒子的飞翔速度
%用inline定义适应度函数以便于将子函数与主程序文件放在一起，
: A" y, P. L3 q4 O. E8 d%目标函数是：y=1+(2.1*(1-x+2*x.^2).*exp(-x.^2/2)）
& O3 w) \+ J1 y$ ]% ~0 g%inline命令定义实应度函数如下
fitness=inline('1/(1+(2.1*(1-x+2*x.^2).*exp(-x.^2/2)))','x');
* f6 K4 I8 K, \$ P%inline定义的实验度函数会使程序运行速度大大降低
for i=1:particlesize
    for j=1:narvs
        f(i)=fitness(x(i,j));
    end
4 v+ x& Z4 D% {8 \; E8 ~2 A, @end
personalbest_x=x;
personalbest_faval=f;
. n2 {( j2 n% r6 _( t[globalbest_faval i]=min(personalbest_faval);
globalbest_x=personalbest_x(i,;
2 q) ]1 c/ g* P  P" [k=1;
- E' Z, K. P) Cwhile k<=MaxNum
    for i=1:particlesize
: H+ b4 o" }% d( K  b7 s        for j=1:narvs;
        f(i)=fitness(x(i,j));
' p( q/ I. m1 w+ t# N$ r        end
    if f(i)<personalbest_faval(i)   %判断当前位置是否历史上最佳位置
( W4 ^, A  C1 L$ u' M0 s. A        personalbest_faval(i)=f(i);
( Q$ e) J7 ?" [# g; p! v* k        personalbest_x(i,=x(i,;
    end
    end
    [globalbest_faval i]=min(personalbest_faval);
    globalbest_x=personalbest_x(i,;
7 S6 V/ b: j0 T; H4 V0 e4 J1 \) L  L    for i=1:particlesize           %更新粒子群里每个个体的最新位置
2 G, c- W8 @! A0 I        v(i,=w*v(i,+c1*rand*(personalbest_x(i,-x(i,)...
% v4 N1 O+ G/ v5 f+ X, d            +c2*rand*(globalbest_x-x(i,);
7 P. d& W# U  `7 B. k         for j=1:narvs             %判断粒子的飞翔速度是否超过了最大飞行速度
             if v(i.j)>vmax;   
                 v(i,j)=vmax;
) ^" b: ~. j; B  k8 k6 g  |$ @             else if v(i,j)<-vmax;
                     v(i,j)=-vmax;
                 end
             end
: q% E! ^2 M0 A- A+ g4 l4 ]4 `             x(i,=x(i,+v(i,;
5 n( Y3 C( ^8 ^5 v         end
% z9 e% X9 w: k  W+ C$ v6 X: }         if abs(globalbest_faval)<E0,break,end
         k=k+1;
/ G1 b: e+ r# z9 {  ~, A- i8 G    end
    Valuel=1/globalbest_faval-1;Valuel=num2str(Valuel);
7 I9 T/ {) {( X. B/ b3 `    %strcat指令可以实现字符的组合输出
    disp(strcat('the maximun valuel','=',Valuel));
  %输出最大值所在的横坐标位置
5 V% @8 r( A  `  S    Value2=globalbest_x; Value2=num2str(Value2);
+ M7 G8 Y6 j; k$ U2 v, s* l9 C% t     disp(strcat('the corresponding coordinate','=',Value2));
     x=-5:0.01:5;
& U0 |& _! {+ k! D" u     y=2.1*(1-x+2*x.^2).*exp(-x.^2/2);
     plot(x,y,'m-','linewidth',3);
     hold on;
; t% c+ T; b, K- d! P     plot(globalbest_x,1/globalbest_faval-1,'kp','linewidth',4);
. c. l' }6 ]; Y! Q. O     legend('目标函数','搜索到的最大值');
+ x" ~/ ?  O# E! R, z     xlabel('x');ylabel('y');grid on;toc;
5 [1 X4 y/ I, v$ k1 {2 l1 @& q
  t: v* C+ a3 D: Z
) u& L% R4 E  w7 l( k6 G: j- P错误提示 ??? function main()    |
Error: Function definitions are not permitted in this context.
( D. o5 c* R$ p0 w3 j9 s5 u
# N; v/ o  G/ h/ @' a% z' m去掉function main()后 提示??? Undefined function or variable 'navrs'.错误。求指点！
* s& o+ s. m1 `- |8 f5 J. Y

平凡之不凡

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

永不漠灭的轨迹-

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

永不漠灭的轨迹-

平凡之不凡 发表于 2014-8-14 09:37
* d. g; U5 u! \7 x代码重新发一下，格式出问题了，你可以调整一下

) S6 y& f5 h& N% c( y$ l
/ {) v% J: z1 P5 C3 s1 R格式是对的  是不是那个自定义函数的问题，我刚刚查了一下，说是要先定义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()
   ( X% x, g* S) Q9 g
2. clc;clear all;close all;
   1 B7 b7 f; W* f+ d$ i( b
3. tic;                              %程序运行计时* Y/ `& n; E' o- _2 Q
   
4. E0=0.001;                        %允许误差
   4 [0 Q% N: w9 M2 l7 h3 }
5. MaxNum=100;                    %粒子最大迭代次数
   9 S2 X; B9 F: o4 ]' }
6. narvs=1;                         %目标函数的自变量个数: F1 x- P# G- n6 ~0 i
   
7. particlesize=30;                    %粒子群规模/ O\" B3 }* d! p9 X5 E: c/ N
   
8. c1=2;                            %每个粒子的个体学习因子，也称为加速常数4 V* h  B* {; E* J9 g) Z/ A% l
   
9. c2=2;                            %每个粒子的社会学习因子，也称为加速常数. e9 _5 J$ X$ k( F7 z\" f  A- n! H
   
10. w=0.6;                           %惯性因子
    1 p) s& t0 ?1 H1 P- }* |' t3 A! u
11. vmax=0.8;                        %粒子的最大飞翔速度& o3 v) v2 t  D: n\" n* o
    
12. x=-5+10*rand(particlesize,narvs);     %粒子所在的位置
    % f9 T3 t+ g* F2 W# G5 A  m2 Y
13. v=2*rand(particlesize,narvs);         %粒子的飞翔速度
    : l* B- L' T' t! ?
14. %用inline定义适应度函数以便将子函数文件与主程序文件放在一起，
    9 n5 k: v4 H' ~' Y# P3 A
15. %目标函数是：y=1+(2.1*(1-x+2*x.^2).*exp(-x.^2/2))6 E/ ~' i: i* p. R( F' H
    
16. %inline命令定义适应度函数如下：
    # u3 a! n3 Y- i- ?1 R: d
17. fitness=inline('1/(1+(2.1*(1-x+2*x.^2).*exp(-x.^2/2)))','x');
    ) F& v( {: z4 }
18. %inline定义的适应度函数会使程序运行速度大大降低
    2 L8 Y2 z9 ^# v# T+ R
19. for i=1:particlesize3 O7 ]7 v1 e6 ^8 Y  M; W
    
20.     for j=1:narvs' @* c- [4 Z\" h' S0 M7 S8 C7 q5 v
    
21.         f(i)=fitness(x(i,j));& l' u& W0 L! f! F
    
22.     end
    1 |1 s( U* ~- V
23. end
    9 {% A, Q3 L1 e8 h7 O
24. personalbest_x=x;
    9 l6 L; l1 o: p- F
25. personalbest_faval=f;
    - t0 Y, @( R7 m+ K9 o\" u% w
26. [globalbest_faval i]=min(personalbest_faval);0 O\" U# b2 P* V' S, U3 W
    
27. globalbest_x=personalbest_x(i,:);
    3 {: s& c0 O$ ?
28. k=1;9 ^4 O, H, Q8 h  `9 Q* }
    
29. while k<=MaxNum
    & u$ S3 f4 Y+ y/ z0 o9 C! W# z% k0 u
30.     for i=1:particlesize2 L5 i) F) w7 J: C+ a) R% w
    
31.         for j=1:narvs
    # M$ s0 c\" ~+ }2 M6 ~% b7 h
32.             f(i)=fitness(x(i,j));9 ]& v% I; }! W0 w( }; h
    
33.         end# C/ q1 S7 E+ N
    
34.         if f(i)<personalbest_faval(i) %判断当前位置是否是历史上最佳位置
    , a3 I  m$ u  J& a/ A
35.             personalbest_faval(i)=f(i);9 d3 \# g  I( ]/ p3 @\" \
    
36.             personalbest_x(i,:)=x(i,:);8 q' v& B( S; Z7 [& M/ A' N1 @( x
    
37.         end
    3 g3 r, ]& w( o# u+ y
38.     end
    , z  V/ M6 d+ V# ~
39.     [globalbest_faval i]=min(personalbest_faval);7 P0 T# c& ]: G5 x* [
    
40.     globalbest_x=personalbest_x(i,:);+ c$ z\" A\" v1 D7 E7 z# X  S- Y
    
41.     for i=1:particlesize %更新粒子群里每个个体的最新位置) x+ z- P- @/ T) G\" e; Z
    
42.         v(i,:)=w*v(i,:)+c1*rand*(personalbest_x(i,:)-x(i,:))...
    ) _+ K5 {3 ~2 o+ M6 y
43.             +c2*rand*(globalbest_x-x(i,:));
    , t3 i2 N: c+ b. M0 A: n4 G
44.         for j=1:narvs    %判断粒子的飞翔速度是否超过了最大飞翔速度  f* l# D\" ~. Z8 U3 W! a. R* e
    
45.             if v(i,j)>vmax;
    + t$ x- O' }- N
46.                 v(i,j)=vmax;5 `\" g4 p0 _$ C9 {
    
47.             elseif v(i,j)<-vmax;5 o2 b' h1 b) C9 J5 \- V
    
48.                 v(i,j)=-vmax;
    ; d# i1 x+ U  z8 p8 `) v) L- W# P
49.             end
    * Y5 O: Z- T5 N3 C5 n
50.         end# l\" P\" ~$ U( |0 e! V( _
    
51.         x(i,:)=x(i,:)+v(i,:);) t) k* @' Z\" P9 M/ z2 }& P+ T1 m
    
52.     end) L3 `+ F/ }6 _8 x- r0 d# ~6 S; O$ |. O
    
53.     if abs(globalbest_faval)<E0,break,end* s+ H6 w\" ~# f2 V0 l\" `
    
54.     k=k+1;
    2 k( F$ B9 Z: s, \
55. end
    9 j2 G; }) o% F
56. Value1=1/globalbest_faval-1; Value1=num2str(Value1);6 z/ _3 m\" Q  i* b\" j% H6 D$ Y+ b\" z
    
57. % strcat指令可以实现字符的组合输出. P2 w6 J8 V$ |7 B. {9 W2 e
    
58. disp(strcat('the maximum value','=',Value1));9 b' g! t\" w  d
    
59. %输出最大值所在的横坐标位置
    2 N2 J9 R0 u, I9 Q+ n# c3 `
60. Value2=globalbest_x; Value2=num2str(Value2);\" T\" T7 D  G+ \9 E$ N3 ~
    
61. disp(strcat('the corresponding coordinate','=',Value2));- ^0 |- M. q+ {\" a! N
    
62. x=-5:0.01:5;& W' l* l; u  u- k4 Z; V  Y) E
    
63. y=2.1*(1-x+2*x.^2).*exp(-x.^2/2);
    % ?- d\" D\" i( K- b
64. plot(x,y,'m-','linewidth',3);
    / {# D: I, Q6 L6 v
65. hold on;% Q5 W5 f4 W6 k2 k) S3 H
    
66. plot(globalbest_x,1/globalbest_faval-1,'kp','linewidth',4);
    ; D# l' E: l( G7 X. ~' @
67. legend('目标函数','搜索到的最大值');xlabel('x');ylabel('y');grid on;toc;$ U, \. y( E( ]; {& R7 E* e7 F
    

复制代码

平凡之不凡

这个是调好的

永不漠灭的轨迹-

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

谢谢老师！