声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1205|回复: 4

紧急求助!(MATLAB辅助优化设计)

[复制链接]
发表于 2006-2-28 14:30 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
<P>  本人编制优化程序时遇到一问题,在M文件中用积分函数运行不起来!为什么?请高手指点一下,谢谢了!<BR>目标函数M file文件:<BR>function f=myfun1(x)<BR>global t;<BR>v=0.8*1.256*t-3.1415+acos((x(1)^2+x(4)^2-2*x(1)*x(4)*cos(1.256*t+2.669)+x(3)^2-x(2)^2)/(2*x(3)*sqrt(x(1)^2+x(4)^2-2*x(1)*x(4)*cos(1.256*t+2.669))))+acos((x(1)^2+x(4)^2-2*x(1)*x(4)*cos(1.256*t+2.669)+x(4)^2-x(1)^2)/(2*x(4)*sqrt(x(1)^2+x(4)^2-2*x(1)*x(4)*cos(1.256*t+2.669))));<BR>f=int(v,t,0,2.5);<BR>约束条件M file 文件:<BR>function [c,ceq]=mycon1(x)<BR>c(1)=x(1)^2+x(3)^2-(x(4)-x(1))^2-2*x(2)*x(3)*cos(3.1415*40/180);<BR>c(2)=-x(2)^2-x(3)^2+(x(4)+x(1))^2-2*x(2)*x(3)*cos(3.1415*40/180);<BR>ceq=0;<BR><BR>MATLAB窗口:<BR>&gt;&gt; A=[1 -1 -1 1;1 -1 1 -1;1 1 -1 -1];<BR>&gt;&gt; b=[0;0;0];<BR>&gt;&gt; x0=[10;51;28.5;53.5];<BR>&gt;&gt; lb=zeros(4,1);<BR>&gt;&gt; [x,fval,exitflag,output,lambda]=fmincon(@myfun1,x0,A,b,[],[],lb,[],@mycon1)<BR><BR>就是目标函数M文件中我想把求积分当作目标函数,加了个f=int(v,t,0,2.5),运行时就出现下面的错误:<BR>??? Error using ==&gt; fmincon<BR>FMINCON cannot continue because user supplied objective function<BR>failed with the following error:</P>
<P>No appropriate methods for function int.<BR><BR>大家帮忙解决一下吧!</P>
回复
分享到:

使用道具 举报

发表于 2006-2-28 15:01 | 显示全部楼层

回复:(monster421)紧急求助!(MATLAB辅助优化设计...

function f=myfun1(x)<BR>global t;<BR>v=0.8*1.256*t-3.1415+acos((x(1)^2+x(4)^2-2*x(1)*x(4)*cos(1.256*t+2.669)+x(3)^2-x(2)^2)/(2*x(3)*sqrt(x(1)^2+x(4)^2-2*x(1)*x(4)*cos(1.256*t+2.669))))+acos((x(1)^2+x(4)^2-2*x(1)*x(4)*cos(1.256*t+2.669)+x(4)^2-x(1)^2)/(2*x(4)*sqrt(x(1)^2+x(4)^2-2*x(1)*x(4)*cos(1.256*t+2.669))));<BR>f=int(v,t,0,2.5);<BR><BR><BR>int用错了,你这样用v就变成数值了<BR>而对于int这个函数而言v应该是表达式
 楼主| 发表于 2006-2-28 19:38 | 显示全部楼层
能具体给个做法吗? v是关于x和t的表达式啊!我还是不理解?
发表于 2006-3-2 09:18 | 显示全部楼层

回复:(monster421)紧急求助!(MATLAB辅助优化设计...

doc int看一下这个函数应该怎么用吧
发表于 2006-3-2 14:29 | 显示全部楼层
help sym/int.m 有例子说明
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-20 23:23 , Processed in 0.059211 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表