24小时热门版块排行榜    

查看: 1035  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 梁成伟老师课题组欢迎你的加入 +6 一鸭鸭哟 2026-03-14 7/350 2026-03-15 22:12 by Winj1e
[考研] 中科院材料273求调剂 +3 yzydy 2026-03-15 3/150 2026-03-15 21:15 by ms629
[考研] 326求调剂 +3 上岸的小葡 2026-03-15 4/200 2026-03-15 18:50 by 无际的草原
[基金申请] NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗 20+5 NSFC2026我来了 2026-03-10 13/650 2026-03-15 12:20 by NSFC2026我来了
[考研] 085600材料与化工 求调剂 +11 enenenhui 2026-03-13 12/600 2026-03-14 21:48 by laoshidan
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
[考研] 328求调剂 +3 5201314Lsy! 2026-03-13 6/300 2026-03-14 15:31 by hyswxzs
[考研] 学硕285求调剂 +13 Wisjxn 2026-03-12 46/2300 2026-03-14 10:33 by JourneyLucky
[考研] 265求调剂 +9 小木虫085600 2026-03-09 12/600 2026-03-14 01:11 by JourneyLucky
[考研] 312求调剂 +6 陌宸希 2026-03-10 6/300 2026-03-14 00:40 by JourneyLucky
[考研] 26考研调剂 +3 ying123. 2026-03-10 3/150 2026-03-14 00:18 by JourneyLucky
[考研] 279求调剂 +3 抓着星星的女孩 2026-03-10 3/150 2026-03-13 23:47 by userper
[考研] 279求调剂 +3 Dizzy123@ 2026-03-10 3/150 2026-03-13 23:02 by JourneyLucky
[考研] 308求调剂 +5 是Lupa啊 2026-03-11 5/250 2026-03-13 22:13 by JourneyLucky
[考研] 【考研调剂求收留】 +3 Ceciilia 2026-03-11 3/150 2026-03-13 20:18 by JourneyLucky
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂! +3 不想理你 2026-03-12 3/150 2026-03-13 14:18 by JourneyLucky
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-13 3/150 2026-03-13 10:43 by houyaoxu
[考研] 333求调剂 +3 152697 2026-03-12 4/200 2026-03-13 07:08 by Iveryant
[考研] 327分求调剂086 +4 西红柿?小帅 2026-03-09 7/350 2026-03-10 14:47 by ruiyingmiao
[考研] 337求调剂 +3 睡醒,。 2026-03-09 3/150 2026-03-09 10:02 by 求调剂zz
信息提示
请填处理意见