24小时热门版块排行榜    

查看: 1653  |  回复: 7

hudf03

铁杆木虫 (小有名气)

[求助] matlab无约束非线性方程求最小值的问题

我要求一个非线性方程f(x1,x2,x3,...)的最小值以及相应的x1,x2,x3,...的值,试验了matlab中的fminsearch, fminunc和lsqnonlin函数,但是它们给出的最小值比初值代入方程的值还小,这是为什么?有没有更好的函数或算法?
回复此楼

» 猜你喜欢

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

板凳宁坐十年冷,文章不写一句空
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
把原问题以文本格式贴上来看看。
2楼2013-03-12 17:48:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hudf03

铁杆木虫 (小有名气)

引用回帖:
2楼: Originally posted by dingd at 2013-03-12 17:48:56
把原问题以文本格式贴上来看看。

原问题是这样的:

x0 = [2.8;1.5;1.5];
[x,flag,exitflag,output] = fminsearch(@min_search,x0)

% function
function y = min_search(x)

f = 0.305;
phi = 33.8*pi/180;
V = 0.4051;

y = sqrt((tan(phi)*(2*pi*(x(1)-pi/4)*f)-1)^2+(V/(sqrt(cos(phi)/2/pi/(x(1)-pi/4)/f))-1)^2 + ...
    (f*(pi^3*x(3))/((3.6*x(1)-0.5*pi^2*x(2)/V))-1)^2);

return;
板凳宁坐十年冷,文章不写一句空
3楼2013-03-13 09:06:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hudf03

铁杆木虫 (小有名气)

原帖写错了,应该是:
它们给出的最小值比初值代入方程的值还"大"
板凳宁坐十年冷,文章不写一句空
4楼2013-03-13 09:09:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

最小值都是一样的:0.441686294317546

但参数值不唯一,有无穷组,应该和您的函数有关吧(有三角函数,还有x2和x3似乎可以化简为一个):
x1: 1.65080945363066
x3: -1.45001854418982
x2: 1.61353777127811

x1: 1.65080945107221
x3: 2.25348683933654
x2: -1.26157531922869

x1: 1.65080945365762
x3: 0.528027682207919
x2: 0.0779367664781694
5楼2013-03-13 09:38:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hudf03

铁杆木虫 (小有名气)

引用回帖:
5楼: Originally posted by dingd at 2013-03-13 09:38:04
最小值都是一样的:0.441686294317546

但参数值不唯一,有无穷组,应该和您的函数有关吧(有三角函数,还有x2和x3似乎可以化简为一个):
x1: 1.65080945363066
x3: -1.45001854418982
x2: 1.61353777127811
...

您用的什么函数啊? 把x0代入到方程中得到的目标函数值是3.1604,小于0.4417
板凳宁坐十年冷,文章不写一句空
6楼2013-03-13 09:49:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
hudf03: 金币+5, ★★★★★最佳答案 2013-03-13 09:56:50
“3.1604,小于0.4417”,??
7楼2013-03-13 09:50:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hudf03

铁杆木虫 (小有名气)

引用回帖:
7楼: Originally posted by dingd at 2013-03-13 09:50:45
“3.1604,小于0.4417”,??

sorry,是我搞错了,多谢!
板凳宁坐十年冷,文章不写一句空
8楼2013-03-13 09:56:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hudf03 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 u1xkze8ybe 2026-07-30 7/350 2026-07-31 12:10 by syhhffkrym
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 u1xkze8ybe 2026-07-30 3/150 2026-07-31 12:08 by syhhffkrym
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 u1xkze8ybe 2026-07-30 6/300 2026-07-31 12:05 by syhhffkrym
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +6 锦衣卫寒战 2026-07-28 6/300 2026-07-31 11:27 by syl200707
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 小木虫看见有人已经知道结果了 +15 1234567wang 2026-07-25 19/950 2026-07-31 08:30 by ZJTJZ
[基金申请] 2026年国自然面上资助率 +8 布布和一二 2026-07-30 8/400 2026-07-31 01:51 by alongwaytogo
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[基金申请] 微信指数没变化,科研之友没阅读 +14 wangze12014 2026-07-28 17/850 2026-07-30 20:10 by kissu88
[考博] 申请2027年材料/化学类博士 10+3 考研老狗? 2026-07-25 5/250 2026-07-30 19:33 by 白色木鱼
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
[基金申请] 今年的WR进展到哪一步了? +5 wsgjhwz 2026-07-25 9/450 2026-07-30 12:59 by 可淡不可忘
[有机交流] 产物和副产物价值比较 30+4 小汤02 2026-07-27 5/250 2026-07-30 12:44 by czyzsu
[基金申请] 时间戳他又来了 +14 晓晓爱翠翠 2026-07-26 17/850 2026-07-29 16:30 by 呵呵二哥
[有机交流] 这个自发加氧反应的机理是什么? 50+4 YUAN2273 2026-07-25 7/350 2026-07-29 12:27 by xiaobenpu
[基金申请] 没消息就是被刷了呗 +13 lambert2014 2026-07-24 14/700 2026-07-29 10:38 by msh145800
[基金申请] 这种情况还有戏吗 +5 drbart 2026-07-27 11/550 2026-07-29 05:50 by drbart
[基金申请] 同事接到电话了,我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
[基金申请] 准备明年的基金了 +5 Tide man 2026-07-27 5/250 2026-07-28 09:25 by newfuzzy1
[分析] 水蒸气蒸馏法乳化现象 +3 月光娇娇 2026-07-26 3/150 2026-07-26 23:25 by caiyun
信息提示
请填处理意见