24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2013  |  回复: 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 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见