24小时热门版块排行榜    

查看: 1023  |  回复: 1

hustyhwang

新虫 (初入文坛)

[求助] 含约束的解方程问题

%%
%Constraints约束
%  (sum(1./(da-p5)) / ( sum(dt./(da-p5))-((t(n)-t(2))*(n-2))/(au(n-1)-au(1)-(n-2)*p5)))>0
%  (au(n-1)-au(1)-(n-2)*p5)/(t(n)-t(2))/(sum(1./(da-p5)) / ( sum(dt./(da-p5))-((t(n)-t(2))*(n-2))/(au(n-1)-au(1)-(n-2)*p5)))>0
%%
% Solution of p5 when F=0  (others parameters are known.)

F=sum(dt.*log(da-p5))-...
    (t(n)-t(2))*log((au(n-1)-au(1)-(n-2)*p5)/(t(n)-t(2))/(sum(1./(da-p5)) /( sum(dt./(da-p5))-((t(n)-t(2))*(n-2))/(au(n-1)-au(1)-(n-2)*p5))))-...
    sum(   psi((   sum(1./(da-p5))/(sum(dt./(da-p5))-((t(n)-t(2))*(n-2))/(au(n-1)-au(1)-(n-2)*p5))  ).*dt)  .*dt    );

想在满足上述两个约束的p5范围内搜索F=0的解,什么matlab命令可以实现?求指点,感激。。

[ Last edited by hustyhwang on 2013-6-21 at 13:05 ]
回复此楼

» 猜你喜欢

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

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

hustyhwang

新虫 (初入文坛)

上述问题源自一个三元方程组求解,将p3p4消元后等到p5的方程求解,约束是由于p3和p4的范围,不加约束直接用fsolve和fminsearch等都会提示psi参数不能为负的错误!

若能指点以下方程组求解,问题同样解决:
(au(n-1)-au(1)-(n-2)*p5-p3.*p4.*(t(n)-t(2)))=0……(1)
(sum( (1-p3.*(t2-t1))./(a2-a1-p5)  )+(n-2)./p4)=0……(2)
temp=p3.*(t2-t1);%这个temp必须大于0,否则psi(temp)无解会报错
(sum(  (t2-t1) .* log(a2-a1-p5)  )-(t(n)-t(2))*log(p4)-sum(psi(temp).*(t2-t1)))=0……(3)

解两天未果,特来请教。
2楼2013-06-21 13:12:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hustyhwang 的主题更新
信息提示
请填处理意见