24小时热门版块排行榜

返回列表

lynettecc

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 379.5
帖子: 13
在线: 15.9小时
虫号: 1759069
注册: 2012-04-16
专业: 运筹学

[求助] matlab程序出错调试已有2人参与

我的程序出现问题，但是不知道原因请求各位高手帮忙，不胜感激！

function f = statefunction(xx)
global mu lam10 Tmp TmE
global n
x=xx(1:n);
x=reshape(x,n,1);
u=xx(end-1:end);
%x1=x(1),x2=x(2),x3=x(3),x4=x(4),x5=x(5),x6=x(6),x7=x(7),x8=x(8),x9=x(9),x10=x(10),x11=x(11),x12=x(12)
% lamb7'=x(13),lamb8'=x(14),lamb9'=x(15),lamb11'=x(16),lamb12'=x(17)
%控制变量u1,u2
u1=u(1);
u2=u(2);

%控制变量u3,u4
uu1=controlE1(x);
u3=uu1(1);
u4=uu1(2);

function y=controlE1(x)

%x1=x(1),x2=x(2),x3=x(3),x4=x(4),x5=x(5),x6=x(6),x7=x(7),x8=x(8),x9=x(9),x10=x(10),x11=x(11),x12=x(12)
% lamb7'=x(13),lamb8'=x(14),lamb9'=x(15),lamb11'=x(16),lamb12'=x(17)

%x1=x(1),x2=x(2),x3=x(3),x4=x(4),x5=x(5),x6=x(6),x7=x(7),x8=x(8),x9=x(9),x10=x(10),x11=x(11),x12=x(12)
%lamb1=x(13),lamb2=x(14),lamb3=x(15),lamb4=x(16),lamb5=x(17),lamb6=x(18),...
%lamb7=x(19),lamb8=x(20),lamb9=x(21),lamb10=x(22),lamb11=x(23),lamb12=x(24);
options = optimset('TolFun',1e-6,'TolCon',1e-8);
options = optimset(options,'Algorithm','sqp');
uuE = fmincon(@(xx) myfun1(xx,x),[0;0],[],[],[],[],[-pi;-pi],[pi;pi],@(xx) cons1(xx,x),options);
u3=uuE(1);
u4=uuE(2);

y=zeros(2,1);
y(1)=u3;
y(2)=u4;

function q=myfun1(p,x)
global TmE

u3=p(1);
u4=p(2);
a=TmE*cos(u4)/x(8)*(x(15)*cos(u3)-x(8)*x(14)*sin(u3));
b=-TmE*x(14)*sin(u4)*cos(u3)+TmE/x(8)*(x(17)*cos(u4)/cos(x(9))-x(15)*sin(u3)*sin(u4));
q=a^2+b^2;

function [c ceq]=cons1(p,x)
global TmE

u3=p(1);
u4=p(2);
H11=-TmE*cos(u4)/x(8)*(x(8)*x(14)*cos(u3)+x(15)*sin(u3));
H12=TmE*sin(u4)/x(8)*(x(8)*x(14)*sin(u3)-x(15)*cos(u3));
H22=-TmE*(1/x(8)*(x(15)*sin(u3)*cos(u4)+x(17)*sin(u4)/cos(x(9)))+x(14)*cos(u3)*cos(u4));
c=zeros(2,1);
c(1)=H11+H22;
c(2)=-(H11*H22-H12*H12);
ceq=[];

如上程序运行时报错，请问哪位高手知道哪里出错吗？
??? Error using ==> sqpLineSearch at 20
Objective function is undefined at initial point. Fmincon cannot continue.

Error in ==> fmincon at 816
[X,FVAL,EXITFLAG,OUTPUT,LAMBDA,GRAD,HESSIAN] =
sqpLineSearch(funfcn,X,full(A),full(B),full(Aeq),full(Beq), ...

Error in ==> controlE1 at 17
uuE = fmincon(@(xx) myfun1(xx,x),[0;0],[],[],[],[],[-pi;-pi],[pi;pi],@(xx) cons1(xx,x),options);

Error in ==> statefunction at 14
uu1=controlE1(x);

回复此楼

» 猜你喜欢

有谁可曾问过你过的还好吗？已经有21人回复
E0414, 我的本子有没有希望？已经有7人回复
一篇论文同时出现在两个期刊，一模一样，这算不算学术不端，请各位老师斧正。已经有12人回复
希望面上有个好结果已经有7人回复
三区计算机方向期刊推荐已经有5人回复
sci论文二审求助已经有5人回复
函评已经有7人回复

» 本主题相关价值贴推荐，对您同样有帮助:

【求助】求大神，matlab程序运行出错已经有6人回复
浓度一时间比法确定反应级数和计算速率常数附带MATLAB程序出错，求大神教我已经有0人回复
浓度一时间比法确定反应级数和计算速率常数 MATLAB程序出错已经有0人回复
matlab程序出错已经有0人回复
MATLAB程序总是出错? 已经有1人回复
Matlab程序求助已经有6人回复
matlab程序运行出错，哪位大神能帮忙解决下，谢谢了已经有1人回复
Matlab程序出错已经有3人回复
实在是看不出这段matlab程序里哪里出错了，请大家帮忙看看！已经有10人回复
帮我看一段Matlab程序哪里出错了已经有1人回复
对流扩散方程matlab程序出错已经有0人回复
matlab 问题求助，一段程序出错，求纠错已经有7人回复
大家爱帮我看看这段MATLAB程序为啥出错了已经有5人回复
matlab程序问题，总是出错，麻烦大家看看问题在那里？已经有6人回复
matlab程序运行出错，求指点修改已经有4人回复
求助一简单matlab程序的出错问题，谢谢已经有3人回复
【求助】那位高手帮看下MATLAB程序出错了已经有4人回复

1楼 2012-11-06 13:14:18

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军

程序强帖: 40
应助: 817 (博后)
金币: 12914.1
红花: 64
帖子: 2238
在线: 287.3小时
虫号: 696514
注册: 2009-02-05
专业: 计算机软件

【答案】应助回帖

★
感谢参与，应助指数 +1
jjdg: 金币+1, 感谢参与 2012-11-07 10:24:22

??? Error using ==> sqpLineSearch at 20
Objective function is undefined at initial point. Fmincon cannot continue.

sqpLineSearch函数中,目标函数在初始点未定义.sqpLineSearch 这个函数体呢?

赞一下(1人)

回复此楼

» 本帖已获得的红花（最新10朵）

lynettecc

matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com

2楼2012-11-06 16:43:04

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

lynettecc

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 379.5
帖子: 13
在线: 15.9小时
虫号: 1759069
注册: 2012-04-16
专业: 运筹学

送鲜花一朵

引用回帖:

2楼: Originally posted by libralibra at 2012-11-06 16:43:04
??? Error using ==> sqpLineSearch at 20
Objective function is undefined at initial point. Fmincon cannot continue.

sqpLineSearch函数中,目标函数在初始点未定义.sqpLineSearch 这个函数体呢?

我没有编辑sqpLineSearch 这个函数体，他是fmincon函数体的一部分，这是matlab工具箱内自带的啊！

我只是调用fmincon函数，其内部怎么算我不了解！并且我的程序能运行一段时间，是中间突然报错的，百思不得其解

赞一下

回复此楼

3楼2012-11-07 09:57:07

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖