24小时热门版块排行榜    

查看: 1681  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 时间戳变了,能看出什么问题? +17 基诺咪客 2026-08-17 22/1100 2026-08-20 16:47 by 基诺咪客
[基金申请] 时间戳又变了 +11 wuchongjun 2026-08-20 14/700 2026-08-20 16:26 by 虫友是什么虫
[基金申请] 科研孤儿太难了 +4 我4大白菜 2026-08-20 4/200 2026-08-20 15:58 by 陈晨晨陈啊
[基金申请] 言之凿凿,明天官网会通知放榜? +14 医学老男孩 2026-08-20 16/800 2026-08-20 15:47 by LNP@mRNA
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +16 Ldrop2023 2026-08-13 20/1000 2026-08-20 11:07 by xskun
[基金申请] 今天基金会出结果吗?20260819 +14 kkkl_v 2026-08-19 15/750 2026-08-20 08:21 by sunzitan
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
请填处理意见