24小时热门版块排行榜    

查看: 1036  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 机械专硕调剂 +3 笨笨兔子 2026-03-12 3/150 2026-03-15 20:02 by 栗子粥?
[考博] 东华理工大学化材专业26届硕士博士申请 +6 zlingli 2026-03-13 6/300 2026-03-15 20:00 by ryzcf
[考研] 0703化学调剂 290分有科研经历,论文在投 +3 腻腻gk 2026-03-14 3/150 2026-03-15 17:28 by 小物理化学
[基金申请] 国自科面上基金字体 +4 iwuli 2026-03-12 5/250 2026-03-15 17:07 by 风云无泪
[考研] 290求调剂 +3 孔志浩 2026-03-12 8/400 2026-03-15 15:30 by 孔志浩
[考研] 材料专硕326求调剂 +4 墨煜姒莘 2026-03-15 4/200 2026-03-15 11:02 by dyw
[考研] 265求调剂 +4 威化饼07 2026-03-12 4/200 2026-03-14 17:23 by userper
[考研] 255求调剂 +3 李嘉慧, 2026-03-12 4/200 2026-03-14 16:58 by 有只狸奴
[考研] 328求调剂 +3 5201314Lsy! 2026-03-13 6/300 2026-03-14 15:31 by hyswxzs
[考研] 307求调剂 +7 超级伊昂大王 2026-03-10 7/350 2026-03-14 00:49 by JourneyLucky
[考研] 0703化学调剂 +4 快乐的香蕉 2026-03-11 4/200 2026-03-13 22:41 by JourneyLucky
[考研] 求材料调剂 +5 隔壁陈先生 2026-03-12 5/250 2026-03-13 22:03 by 星空星月
[考研] 工科,求调剂 +3 我887 2026-03-11 3/150 2026-03-13 21:39 by JourneyLucky
[考研] 【考研调剂求收留】 +3 Ceciilia 2026-03-11 3/150 2026-03-13 20:18 by JourneyLucky
[考研] 求调剂 +5 一定有学上- 2026-03-12 5/250 2026-03-13 18:31 by ms629
[考研] 材料专硕350 求调剂 +4 王金科 2026-03-12 4/200 2026-03-13 16:02 by ruiyingmiao
[考研] 材料调剂,307分 +13 张泳铭1 2026-03-09 17/850 2026-03-13 11:09 by 薛云鹏
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-13 3/150 2026-03-13 10:43 by houyaoxu
[考研] 298求调剂 +3 Vv呀! 2026-03-10 3/150 2026-03-10 22:40 by 剑诗杜康
[考研] 收调剂 +7 调剂的考研学生 2026-03-10 7/350 2026-03-10 17:57 by 麦茶汤圆
信息提示
请填处理意见