查看: 2448  |  回复: 11

citihome

木虫 (正式写手)

[交流] 求不等式方程的解法 已有3人参与

其中既包含值域约束,不等式中也包含导数.请问这样的方程如何求解?
求不等式方程的解法
loss.JPG
求不等式方程的解法-1
inequation.JPG[ Last edited by citihome on 2014-8-25 at 15:50 ]
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : approximate_loss_function.pdf
  • 2014-08-25 15:50:05, 270.4 K

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

peterflyer

木虫之王 (文学泰斗)

peterflyer



小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
10楼: Originally posted by citihome at 2014-08-25 18:29:55
能把你的步骤写写么(这里不仅有y'和y,还有y''.没想到怎么就全部变成y'相关的方程了)
没跟上...

第二个方程:令u=y',
故y''=dy'/dx=dy'/dy*dy/dx=u*du/dy,
代入原方程后,
3/[1+e^(-y)]*u^2-3/2=u*du/dy=1/2*d(u^2)/dy,
整理后得:d(u^2)/dy-6/[1+e^(-y)]*u^2=-3
这是一个u^2的一阶线性非其次常微分方程,代入专门的求解公式后得到:
u^2=(y')^2=C*(1+e^y)^6-3*(1+e^y)/5-3*(1+e^y)^2/4-(1+e^y)^3-3*(1+e^y)^4/2-3*(1+e^y)^5-3*(1+e^y)^6*Ln(1+e^y)-3*y*(1+e^y)^6
Integral{dy/SQRT[C*(1+e^y)^6-3*(1+e^y)/5-3*(1+e^y)^2/4-(1+e^y)^3-3*(1+e^y)^4/2-3*(1+e^y)^5-3*(1+e^y)^6*Ln(1+e^y)-3*y*(1+e^y)^6
]}=A ± x
11楼2014-08-25 19:11:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

Edstrayer

版主 (著名写手)

方寸斗室小天地正气迷漫大世界


小木虫: 金币+0.5, 给个红包,谢谢回帖


则可以构造函数


那么函数p(x)满足下列条件:



青葱岁月圣诞夜,浪漫歌舞迎新年。
2楼2014-08-24 16:02:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็


小木虫: 金币+0.5, 给个红包,谢谢回帖
问题在于, 这类问题跟普通微分方程组不同,不存在简单的"通解";
而且求解的区间有限制时,解的花样可能更多;
所以, 除非给出更详细的问题的背景,否则难以有好的答案.
ส็็็็็็็็็็็็็็็็็็็็
3楼2014-08-24 17:49:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer



小木虫: 金币+0.5, 给个红包,谢谢回帖
因为p∈[0,1],故p≥0
而由题设条件 p*dp/dx≥0
故 dp/dx≥0
因为(dp/dx)^2+p*d^2p/dx^2=d(p*dp/dx)/dx
故有0≤d(p*dp/dx)/dx≤d^2p/dx^2
两边积分后得:
C≤p*dp/dx≤dp/dx + D
F+C*x≤ p^2≤p+D*x + G
C、D、F、G均为积分常数。
由 p^2≥F+C*x,解得:
p(x)≤-SQRT(F+C*x)                                                                                 (1)

p(x)≥SQRT(F+C*x)                                                                                  (2)
由p^2≤p+D*x + G,解得:
1/2*{1-SQRT[1+4*(D*x + G)]}≤p≤1/2*{1+SQRT[1+4*(D*x + G)]}               (3)
联立(1)、(3)或联立(2)、(3)求得p(x)。这两个解的合集即为所要求的解。
4楼2014-08-24 19:48:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer



小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
4楼: Originally posted by peterflyer at 2014-08-24 19:48:25
因为p∈,故p≥0
而由题设条件 p*dp/dx≥0
故 dp/dx≥0
因为(dp/dx)^2+p*d^2p/dx^2=d(p*dp/dx)/dx
故有0≤d(p*dp/dx)/dx≤d^2p/dx^2
两边积分后得:
C≤p*dp/dx≤dp/dx + D
F+C*x≤ p^2≤p+D*x + G
C、D、F、 ...

对了,刚才的求解中忘了将条件p∈[0,1]加上。加上后,(1)的解集和(3)的解集中小于零的部分不合题意被舍去,则得:

因为p∈[0,1],故p≥0
而由题设条件 p*dp/dx≥0
故 dp/dx≥0
因为(dp/dx)^2+p*d^2p/dx^2=d(p*dp/dx)/dx
故有0≤d(p*dp/dx)/dx≤d^2p/dx^2
两边积分后得:
C≤p*dp/dx≤dp/dx + D
F+C*x≤ p^2≤p+D*x + G
C、D、F、G均为积分常数。
由 p^2≥F+C*x,解得:
p(x)≤-SQRT(F+C*x)                                                                                 (1)

p(x)≥SQRT(F+C*x)                                                                                  (2)
由p^2≤p+D*x + G,解得:
1/2*{1-SQRT[1+4*(D*x + G)]}≤p≤1/2*{1+SQRT[1+4*(D*x + G)]}               (3)
因为 p∈[0,1],解集中小于零的部分不合题意被舍去,则得:
令f(x)=Max{SQRT(F+C*x) , 1/2*{1+SQRT[1+4*(D*x + G)]}}≤1
则:
   0≤p(x)≤f(x)≤1  是所求问题的解。
5楼2014-08-24 19:58:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

citihome

木虫 (正式写手)

引用回帖:
5楼: Originally posted by peterflyer at 2014-08-24 19:58:54
对了,刚才的求解中忘了将条件p∈加上。加上后,(1)的解集和(3)的解集中小于零的部分不合题意被舍去,则得:

因为p∈,故p≥0
而由题设条件 p*dp/dx≥0
故 dp/dx≥0
因为(dp/dx)^2+p*d^2p/dx^2=d(p*dp/dx ...

谢谢
如果是这个不等式方程能解么?
见贴图
求不等式方程的解法-2
loss.JPG
6楼2014-08-25 14:51:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

citihome

木虫 (正式写手)

引用回帖:
3楼: Originally posted by cooooldog at 2014-08-24 17:49:11
问题在于, 这类问题跟普通微分方程组不同,不存在简单的"通解";
而且求解的区间有限制时,解的花样可能更多;
所以, 除非给出更详细的问题的背景,否则难以有好的答案.

原始问题是求最优 强凸函数(带约束的最优泛函问题)
尝试了好久,都得不到简单的定解条件(附件).
于是尝试去除"最优",只要求能构造强凸就行
结果又不会解解方程

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : approximatelossfunction.pdf
  • 2014-08-25 15:52:33, 270.4 K
7楼2014-08-25 15:54:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

citihome

木虫 (正式写手)

引用回帖:
2楼: Originally posted by Edstrayer at 2014-08-24 16:02:48

a>0,0<b<\frac{1}{2}e^{-a}
则可以构造函数
p(x)=be^{ax}(x\in)
那么函数p(x)满足下列条件:
Range(p)=
p\frac{dp}{dx}=ab^2e^{2ax}\geqslant 0
0\leqslant\left(\frac{dp}{dx}\right)^2+p\frac ...

谢谢,不知道你怎么想到的
8楼2014-08-25 15:56:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer



小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
6楼: Originally posted by citihome at 2014-08-25 14:51:54
谢谢
如果是这个不等式方程能解么?
见贴图

loss.JPG...

第一个不会做。第二个,令u=y',则会得到一个以u^2为因变量、以y为自变量的一阶线性常微分方程,求解解得(y')^2=C*[y'/(y'-1)]^6-3 ,下面解不出这个以y'为未知数的一元高次方程非线性代数方程,下面的求解也就进行不下去了。
9楼2014-08-25 18:16:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

citihome

木虫 (正式写手)

引用回帖:
9楼: Originally posted by peterflyer at 2014-08-25 18:16:51
第一个不会做。第二个,令u=y',则会得到一个以u^2为因变量、以y为自变量的一阶线性常微分方程,求解解得(y')^2=C*^6-3 ,下面解不出这个以y'为未知数的一元高次方程非线性代数方程,下面的求解也就进行不下去了。...

能把你的步骤写写么(这里不仅有y'和y,还有y''.没想到怎么就全部变成y'相关的方程了)
没跟上
10楼2014-08-25 18:29:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 citihome 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] 化学口download_prp&amp;fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 关于代码变化问题,想知道的进来 +9 且听虎啸 2026-08-07 12/600 2026-08-07 21:53 by 理想三旬2025
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode与中标关系的预测 +3 布布和一二 2026-08-07 3/150 2026-08-07 15:09 by gltch
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[论文投稿] 十年后又回来了,论文投稿求助 +4 哈哈114477 2026-08-01 4/200 2026-08-07 14:39 by jgy194592
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] 关于豆爷回答的JTJC与%2F数量 +4 yang182083 2026-08-06 6/300 2026-08-07 12:43 by 虫友是什么虫
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见