24小时热门版块排行榜    

查看: 1037  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 308求调剂 +4 是Lupa啊 2026-03-09 4/200 2026-03-16 01:21 by Xttdmn
[基金申请] NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗 20+5 NSFC2026我来了 2026-03-10 14/700 2026-03-15 23:53 by 不负韶华的虎
[考研] 东南大学364求调剂 +3 JasonYuiui 2026-03-15 3/150 2026-03-15 18:57 by 无际的草原
[考研] 321求调剂 +3 大米饭! 2026-03-15 3/150 2026-03-15 17:48 by 哈哈哈哈嘿嘿嘿
[考研] 290求调剂 +3 孔志浩 2026-03-12 8/400 2026-03-15 15:30 by 孔志浩
[考研] 309求调剂 +4 花与叶@ 2026-03-10 4/200 2026-03-14 21:26 by a不易
[考研] 290求调剂 +4 @将就将就看 2026-03-10 8/400 2026-03-14 14:23 by 千千运气
[考研] 0703化学求调剂 +5 很老实人 2026-03-09 5/250 2026-03-14 02:57 by JourneyLucky
[考研] 环境调剂 +6 晓看天暮看云 2026-03-09 6/300 2026-03-14 01:16 by JourneyLucky
[考研] 材料工程专硕,一志愿中国矿业大学,总分314,求调剂 +5 无懈可击的巨人 2026-03-10 5/250 2026-03-14 00:37 by JourneyLucky
[考研] 332求调剂 +3 zjy101327 2026-03-11 6/300 2026-03-13 22:48 by JourneyLucky
[考研] 泣血叩求调剂恩,愿以丹心报师恩 +6 Iuruoh 2026-03-11 6/300 2026-03-13 22:06 by JourneyLucky
[考研] 304求调剂 +7 7712b 2026-03-13 7/350 2026-03-13 21:42 by peike
[考研] 【考研调剂求收留】 +3 Ceciilia 2026-03-11 3/150 2026-03-13 20:18 by JourneyLucky
[考研] 工科278分求调剂 +5 周慢热啊 2026-03-12 7/350 2026-03-13 15:49 by JourneyLucky
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂! +3 不想理你 2026-03-12 3/150 2026-03-13 14:18 by JourneyLucky
[考研] 26考研求调剂 +5 丶宏Sir 2026-03-13 5/250 2026-03-13 13:05 by JourneyLucky
[考研] 268求调剂 +4 好运连绵不绝 2026-03-12 4/200 2026-03-13 10:45 by hyswxzs
[考研] 290求调剂 +3 柯淮然 2026-03-10 8/400 2026-03-11 13:48 by 柯淮然
[考研] 一志愿:武汉理工,材料工程,英二数二 总分314 +3 2202020125 2026-03-10 4/200 2026-03-10 13:54 by xiongyaxuan
信息提示
请填处理意见