全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件
1361 0
2010-05-28
我编了个程序,但是不知道为什么运行不了!求高手帮忙。
function f=myfun(x)
F=[acos((x(2)^2+1-(x(3)-x(1))^2)/(2*x(2)));acos(((x(1)+x(3))^2-x(2)^2-1)/(2*x(2)))];
f=min(F(:));
function [c,ceq]=mu(x)
c(1)=45*pi/180-acos((x(2)^2+1-(x(3)-x(1))^2)/(2*x(2)));
c(2)=45*pi/180-acos(((x(1)+x(3))^2-x(2)^2-1)/(2*x(2)));
u=acos((1+x(3)^2-(x(2)+x(1))^2)/(2*x(3)))-acos((1+x(3)^2-(x(2)-x(1))^2)/(2*x(3)));
K=acos((x(2)^2+x(1)^2+cos(u)-1)/(x(2)^2-x(1)^2));
ceq=(pi+K)/(pi-K)-1.1;
A=[1 1 -1;1 -1 -1;1 -1 1];
b=[1;-1;1];
lb=[0;0;0];
ub=[];
option=optimset;option.LargeScale='off';option.Display='off';
[x,f]=fmincon(@myfun,[1;2;1],A,b,[],[],lb,[],@mu,option)
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群