24小时热门版块排行榜    

查看: 1037  |  回复: 5

研究僧战士

木虫 (正式写手)

[求助] 求解L最大值,生生把工艺问题变成了数学题,插值来看有最大值 已有3人参与

回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer


【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
研究僧战士: 金币+20, ★★★很有帮助 2017-04-19 09:00:49
这是一个条件极值问题,只不过极值的条件是个隐函数而已。
建议楼主试试拉格朗日乘子法进行求解,也就是再造一个函数F=k+Cosx+λ*{sqrt[1-k^2]+k*ArcCos(-k)-Sinx-k*x}。然后在0<k<1、π<x<3*π的区域内求解F的极值。其中,λ为拉格朗日乘子。
也可以试试用罚函数法求函数F=k+Cosx-1000*{sqrt[1-k^2]+k*ArcCos(-k)-Sinx-k*x}^2的极大值,而0<k<1、π<x<3*π。1000是取的罚函数因子,是个比原函数的取值范围大很多的数。
2楼2017-04-18 11:16:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

研究僧战士

木虫 (正式写手)

引用回帖:
2楼: Originally posted by peterflyer at 2017-04-18 11:16:09
这是一个条件极值问题,只不过极值的条件是个隐函数而已。
建议楼主试试拉格朗日乘子法进行求解,也就是再造一个函数F=k+Cosx+λ*{sqrt+k*ArcCos(-k)-Sinx-k*x}。然后在0<k<1、π<x<3*π的区域内求解F ...

非常感谢!
拉格朗日乘子法进行求解后
得到一个关于k, x 的超越方程组。不太会解这个东西,据说matlab可以解决是吗?
怎么写程序呢?
求解L最大值,生生把工艺问题变成了数学题,插值来看有最大值-1
webwxgetmsgimg.jpg

3楼2017-04-19 08:01:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer


【答案】应助回帖

引用回帖:
3楼: Originally posted by 研究僧战士 at 2017-04-19 08:01:03
非常感谢!
拉格朗日乘子法进行求解后
得到一个关于k, x 的超越方程组。不太会解这个东西,据说matlab可以解决是吗?
怎么写程序呢?

webwxgetmsgimg.jpg
...

看一下“数值分析”课程中的“非线性方程组的迭代解法”相关章节。一般用牛顿-拉斐逊方法。在给定的初值X=X0附近给个扰动ΔX,将非线性方程组全部在初值处线性化了,得到:F(X0)+dF/dX*ΔX=0。其中,X0为初始向量,ΔX为未知向量的改变。然后以求解出未知数ΔX。再令X(i+1))=X0(i))+ΔX(i),反复迭代下去,直至向量ΔX(i)的范数与向量X(i)的范数之比小于某个小正数如0.000001,迭代收敛。最后得到欲求的数值。
4楼2017-04-19 09:39:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

supervb

铁杆木虫 (文坛精英)

【答案】应助回帖

感谢参与,应助指数 +1
看了回帖,发现本来简单的问题被你们越弄越复杂了!
约束方程才一个,直接对k和x的关系进行隐函数绘图(我刚画过图了),会发现 x在[Pi,2.5Pi]内是随着k的增大单调递减,考虑到余弦函数周期性且最大值为1,以及x位于[2.5Pi,3Pi]内时k值恒小于0.2的情况,可以知道优化参数空间需要在k的[0,1]和x的[Pi,2.5Pi]内即可;
然后将L=k+cos(x)变为k=L-cos(x)并代入原来的隐函数中,得到关于L和x的隐函数,且已知x的有效定义域为[Pi,2.5Pi],由于x取值小程度波动对L值影响极小。
经过隐函数绘图法,粗略求解得到x=6时有L的最大值为1.26
咱是做非线性随机动力学的哟
5楼2017-04-19 10:53:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

my768188

金虫 (正式写手)

【答案】应助回帖

楼主书写要规范啊
最大值1.26,具体见图
求解L最大值,生生把工艺问题变成了数学题,插值来看有最大值-2
捕获2.PNG

6楼2017-05-07 21:52:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 研究僧战士 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 288求调剂 +15 于海海海海 2026-03-19 15/750 2026-03-19 22:41 by 学员8dgXkO
[考研] 一志愿南京理工大学085701资源与环境302分求调剂 +3 葵梓卫队 2026-03-18 5/250 2026-03-19 19:35 by 给你你注意休息
[考研] 288求调剂,一志愿华南理工大学071005 +5 ioodiiij 2026-03-17 5/250 2026-03-19 18:22 by zcl123
[考研] 本人考085602 化学工程 专硕 +17 不知道叫什么! 2026-03-15 19/950 2026-03-19 15:06 by 尽舜尧1
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +11 yangfz 2026-03-17 11/550 2026-03-19 15:06 by houyaoxu
[考研] 324分 085600材料化工求调剂 +3 llllkkkhh 2026-03-18 3/150 2026-03-19 14:22 by houyaoxu
[考研] 311求调剂 +4 冬十三 2026-03-18 4/200 2026-03-18 21:47 by 尽舜尧1
[考研] 化学工程321分求调剂 +15 大米饭! 2026-03-15 18/900 2026-03-18 14:52 by haxia
[考研] 08工科 320总分 求调剂 +5 梨花珞晚风 2026-03-17 5/250 2026-03-18 14:49 by haxia
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[考研] 293求调剂 +11 zjl的号 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
[考研] 302求调剂 +4 小贾同学123 2026-03-15 8/400 2026-03-17 10:33 by 小贾同学123
[考研] 283求调剂 +3 听风就是雨; 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] 11408 一志愿西电,277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 290求调剂 +3 ADT 2026-03-13 3/150 2026-03-13 10:19 by peike
信息提示
请填处理意见