查看: 3269|回复: 7

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

字体大小: 正常放大

6 主题	15 听众	148 积分

升级 24%

TA的每日心情

	开心 2016-1-28 16:30

签到天数: 52 天

[LV.5]常住居民I

自我介绍: 第一次参加国赛大家多多照顾！

群组: 数学建摸协会

群组: Matlab讨论组

电梯直达

1^#

发表于 2014-8-14 09:21 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

>> 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;

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

去掉function main()后提示??? Undefined function or variable 'navrs'.错误。求指点！

zan

转播0 淘帖0 分享0 收藏0 支持0 反对0 微信

使用道具举报

平凡之不凡

90 主题	105 听众	9482 积分

TA的每日心情

	怒 2017-1-30 23:35

签到天数: 318 天

[LV.8]以坛为家I

2010挑战赛参赛者

2011挑战赛参赛者

2012挑战赛参赛者

2013挑战赛参赛者

群组: 物联网工程师培训

群组: Matlab讨论组

群组: 2013认证赛C题讨论群组

群组: 数学建模

群组: 2013认证赛A题讨论群组

2^#

发表于 2014-8-14 09:37 |只看该作者 |招呼Ta 关注Ta

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

点评

永不漠灭的轨迹- 格式是对的是不是那个自定义函数的问题，我刚刚查了一下，说是要先定义M文件。详情回复发表于 2014-8-14 10:21

我也说一句

发表

使用道具举报

永不漠灭的轨迹-

6 主题	15 听众	148 积分

升级 24%

TA的每日心情

	开心 2016-1-28 16:30

签到天数: 52 天

[LV.5]常住居民I

自我介绍: 第一次参加国赛大家多多照顾！

群组: 数学建摸协会

群组: Matlab讨论组

3^#

发表于 2014-8-14 10:20 |只看该作者 |招呼Ta 关注Ta

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

我也说一句

发表

使用道具举报

永不漠灭的轨迹-

6 主题	15 听众	148 积分

升级 24%

TA的每日心情

	开心 2016-1-28 16:30

签到天数: 52 天

[LV.5]常住居民I

自我介绍: 第一次参加国赛大家多多照顾！

群组: 数学建摸协会

群组: Matlab讨论组

4^#

发表于 2014-8-14 10:21 |只看该作者 |招呼Ta 关注Ta

平凡之不凡发表于 2014-8-14 09:37
+ h4 A5 y9 X. n9 l8 |: w3 N代码重新发一下，格式出问题了，你可以调整一下

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

收起(5)

平凡之不凡 :嗯，对的，可能是这样的问题

2014-08-14 11:12 回复

永不漠灭的轨迹- :回复平凡之不凡 :求老师指点，昨天在我同学机子上一样的程序可以运行，然后我拷贝过来，有出现上面相同的问题，不知到要怎么弄，希望老师解答！

2014-08-15 09:54 回复

永不漠灭的轨迹- :回复平凡之不凡 :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-08-15 10:05 回复

平凡之不凡 :你调的怎么样了，我改了以后，发现小错误好多

2014-08-15 17:14 回复

永不漠灭的轨迹- :回复平凡之不凡 :啊，这样啊我不会啊，希望老师指点！

2014-08-16 10:04 回复

我也说一句

发表

使用道具举报

永不漠灭的轨迹-

6 主题	15 听众	148 积分

升级 24%

TA的每日心情

	开心 2016-1-28 16:30

签到天数: 52 天

[LV.5]常住居民I

自我介绍: 第一次参加国赛大家多多照顾！

群组: 数学建摸协会

群组: Matlab讨论组

5^#

发表于 2014-8-15 10:00 |只看该作者 |招呼Ta 关注Ta

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;

我也说一句

发表

使用道具举报

平凡之不凡

90 主题	105 听众	9482 积分

TA的每日心情

	怒 2017-1-30 23:35

签到天数: 318 天

[LV.8]以坛为家I

2010挑战赛参赛者

2011挑战赛参赛者

2012挑战赛参赛者

2013挑战赛参赛者

群组: 物联网工程师培训

群组: Matlab讨论组

群组: 2013认证赛C题讨论群组

群组: 数学建模

群组: 2013认证赛A题讨论群组

6^#

发表于 2014-8-17 08:58 |只看该作者 |招呼Ta 关注Ta

function main()
. D. T/ E' S( n+ ^
clc;clear all;close all;6 k- X/ O9 E1 |) Z6 g l% u1 |
tic; %程序运行计时
2 ]( z\" z1 h/ W
E0=0.001; %允许误差
% ]! X4 O% d* P% R' a4 R! D3 r
MaxNum=100; %粒子最大迭代次数
+ a+ S1 N- h, d; L, d
narvs=1; %目标函数的自变量个数, U) o+ b0 a! L/ a' ]% b! h. O
particlesize=30; %粒子群规模& x- z: |7 m$ O\" \$ j! W. u& o
c1=2; %每个粒子的个体学习因子，也称为加速常数
v' U% `; ^5 [
c2=2; %每个粒子的社会学习因子，也称为加速常数$ y4 b6 u2 t3 Y( i$ D6 S& c! [
w=0.6; %惯性因子
/ S2 i& B& \5 I$ P! B* J
vmax=0.8; %粒子的最大飞翔速度
e, `& S3 R U- M
x=-5+10*rand(particlesize,narvs); %粒子所在的位置
( O\" I, }\" Q& r' s. E3 V! U2 V: C7 u
v=2*rand(particlesize,narvs); %粒子的飞翔速度- @; z: A- O: \- c! z9 }$ y5 n7 o
%用inline定义适应度函数以便将子函数文件与主程序文件放在一起，9 g. B! I6 |& x8 C
%目标函数是：y=1+(2.1*(1-x+2*x.^2).*exp(-x.^2/2)), l& C: G1 k8 z& J. J, Q2 E7 C2 L% \
%inline命令定义适应度函数如下：\" O# g: l+ }+ _/ ]4 o* ~
fitness=inline('1/(1+(2.1*(1-x+2*x.^2).*exp(-x.^2/2)))','x');8 o$ B) s/ a: a4 m! N1 {* N2 _
%inline定义的适应度函数会使程序运行速度大大降低
8 y* ^+ L) I\" a
for i=1:particlesize$ Y# p. V, w9 C9 q. Y
for j=1:narvs
9 M! v7 B. F2 N; L. V7 L- b4 H3 m' V
f(i)=fitness(x(i,j));$ L* d) y+ [, l
end
) E/ O6 z% m, w W
end. Y7 r% R/ R$ a
personalbest_x=x;7 K; O$ B6 S v% Z' D1 \, r
personalbest_faval=f;
* H: _9 `$ o9 e( |$ E0 n$ a
[globalbest_faval i]=min(personalbest_faval);( J3 |# c. W% V4 J
globalbest_x=personalbest_x(i,:);' P9 D4 I1 J3 f
k=1;) y Y& ?. Q+ K) s( M
while k<=MaxNum
7 X; p8 n8 A9 D- G( R$ h# J
for i=1:particlesize
for j=1:narvs
6 l* S4 m' A$ W! F# t# I
f(i)=fitness(x(i,j));
6 g0 {- ?. [$ p2 }& _, c5 { Z& A
end3 B; m' i\" H# j ?
if f(i)<personalbest_faval(i) %判断当前位置是否是历史上最佳位置
/ X3 H. X0 A- Z8 o
personalbest_faval(i)=f(i);) o4 o3 ^( i5 [! M, U8 R3 v0 d, R
personalbest_x(i,:)=x(i,:);
& |2 i N2 Z: @3 ^ _\" i
end3 U, q& u& u# m. N( D9 W. B
end
/ B- q4 g. A- C! I
[globalbest_faval i]=min(personalbest_faval);
& S! z, [# X1 g; e6 Q4 L
globalbest_x=personalbest_x(i,:);
2 u& F7 ?% h5 w) C9 o K8 n+ m' T
for i=1:particlesize %更新粒子群里每个个体的最新位置
) Y- H6 n+ F p/ _
v(i,:)=w*v(i,:)+c1*rand*(personalbest_x(i,:)-x(i,:))...
0 R# \/ n+ E# v4 n
+c2*rand*(globalbest_x-x(i,:));8 ]6 C- Y& o$ T3 R
for j=1:narvs %判断粒子的飞翔速度是否超过了最大飞翔速度
; T- y% U2 E* Z. o1 K8 m\" P3 _, s
if v(i,j)>vmax;4 B# V* x# p1 u( f8 l2 m
v(i,j)=vmax;
\" Y$ [' e' o1 I e$ R) l4 O7 D: ~
elseif v(i,j)<-vmax;
& w: w' i; }3 g0 y$ T
v(i,j)=-vmax;
& U3 ^9 |8 f/ D
end1 C0 L) s; \8 q. r# M9 b\" t7 _
end4 G' b6 B* P, n\" x9 X
x(i,:)=x(i,:)+v(i,:);
8 u) u\" Z9 D5 q& E9 }: W3 g G! H
end
- I8 ^# T% L' B- m6 ~* u# X
if abs(globalbest_faval)<E0,break,end
: W7 O\" z* ^8 U6 C: a8 V% E1 n
k=k+1;* r5 [/ g! x0 E7 ~: P1 P
end7 Z K d& r# H2 ~% A; I
Value1=1/globalbest_faval-1; Value1=num2str(Value1);8 q3 F2 P0 K& k3 b* [
% strcat指令可以实现字符的组合输出# u3 U- r( c* B! d [# |& o
disp(strcat('the maximum value','=',Value1));8 s8 H- k, r7 H+ \
%输出最大值所在的横坐标位置# s1 Z% f2 W+ F9 u
Value2=globalbest_x; Value2=num2str(Value2);
9 X# { s# K2 f& o) u
disp(strcat('the corresponding coordinate','=',Value2));9 g2 l( M$ ~; G1 r3 e! o# P
x=-5:0.01:5;
% P! b& w0 [6 z0 b, g9 N
y=2.1*(1-x+2*x.^2).*exp(-x.^2/2);' H8 W+ d- V6 L
plot(x,y,'m-','linewidth',3);
' p, F$ i! X2 T- Y* T; _- W
hold on;
. x3 L9 c* w! E) r* S) w4 h0 G6 J
plot(globalbest_x,1/globalbest_faval-1,'kp','linewidth',4);
1 h K6 N9 v' Y* W( z7 R; v/ T; ~5 }
legend('目标函数','搜索到的最大值');xlabel('x');ylabel('y');grid on;toc;
& j% E' O5 G/ D0 T' s' C