24小时热门版块排行榜    

查看: 2736  |  回复: 22

zju-hjm

金虫 (初入文坛)

[交流] 【求助】求问一个寻找最优解问题,是泛函问题么 已有6人参与

从一个二阶的偏微分方程出发
dy/dt=d(f(r)*(dy/dr))/dr
其中t与r为自变量;
y=y(t,r),f(r)为r的函数;
dy/dt,dy/dr为一阶导数;
等式右边为f(r)*(dy/dr)对r的导数;
0<=r<=R,R为常数;
现在要使得(dy/dr)在r=R处为定值,怎样确认f(r)达到这种结果,也就是求合适的f(r),这是一个什么问题呢?

怎么最优化f(r),这是泛函问题么?

谢谢帮忙啊

[ Last edited by wuguocheng on 2010-3-31 at 10:12 ]
回复此楼

» 猜你喜欢

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

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

wangwen0811

金虫 (正式写手)


小木虫(金币+0.5):恭喜抢沙发,给个红包
bluesine(金币+0):再次警告! 1-24 14:28
思考思考
2楼2010-01-24 12:19:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bluesine

铁杆木虫 (职业作家)

科苑小木虫


小木虫(金币+0.2):抢了个小板凳,给个红包
zju-hjm(金币+1): 2010-03-31 10:21
感觉楼主的问题听起来有点别扭,呵呵。
直觉是这不是一个泛函问题。
板凳要做十年冷文章不发一个字
3楼2010-01-24 14:39:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bluesine

铁杆木虫 (职业作家)

科苑小木虫

★ ★ ★
Doctorcbw(金币+3):谢谢你的回答 加油 1-24 17:31
zju-hjm(金币+1): 2010-03-31 10:21
结合一下物理意义参考一下:
你的方程就是一个半径为R的球型容器内的温度的变化规律的描述,建议找下学物理的人来帮下忙。了解一下几个符号的意义。dy/dr在r=R处为定值本质上是要求在容器边缘处满足什么条件,(恒定的热对流。。。。具体不懂),反过来在看就明白了
板凳要做十年冷文章不发一个字
4楼2010-01-24 14:53:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zju-hjm

金虫 (初入文坛)

引用回帖:
Originally posted by bluesine at 2010-1-24 14:53:
结合一下物理意义参考一下:
你的方程就是一个半径为R的球型容器内的温度的变化规律的描述,建议找下学物理的人来帮下忙。了解一下几个符号的意义。dy/dr在r=R处为定值本质上是要求在容器边缘处满足什么条件,(恒 ...

我这个是个传质的物理问题,dy/dr是传质推动力,物理模型是一个球形颗粒,里面的物质往外扩散,目标是要让球形里面的物质(可以认为是里面的溶质)以恒定的速率扩散出来,也就是想dy/dr在r=R处为定值。具体的边界条件我可以给,就是
r=R处y的值;
dy/dr在r=0处有限值;
5楼2010-01-24 19:01:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bluesine

铁杆木虫 (职业作家)

科苑小木虫

★ ★ ★
zju-hjm(金币+3): 1-25 10:15
zju-hjm(金币+1): 2010-03-31 10:21
引用回帖:
Originally posted by zju-hjm at 2010-01-24 19:01:03:


我这个是个传质的物理问题,dy/dr是传质推动力,物理模型是一个球形颗粒,里面的物质往外扩散,目标是要让球形里面的物质(可以认为是里面的溶质)以恒定的速率扩散出来,也就是想dy/dr在r=R处为定值。具体的 ...

呵呵,这个我想继续问下去。假设你现在已经找到了这个f(r),就是说你能够让球形里面的物质不断向外扩散出来,那么你这个球是个怪物了,它一直向外扩散有没有耗尽的时候?或者说它是个黑洞了,可以无限制向外扩散?所以我觉得你的这个一直向外以恒定速率扩散的要求有点不合理。
板凳要做十年冷文章不发一个字
6楼2010-01-25 09:47:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zju-hjm

金虫 (初入文坛)

bluesine(金币+0):在一定的时间长度里这是可能的,呵呵。从你的论述看不妨取f(r)为常数试试 1-25 11:33
引用回帖:
Originally posted by bluesine at 2010-01-25 09:47:05:



呵呵,这个我想继续问下去。假设你现在已经找到了这个f(r),就是说你能够让球形里面的物质不断向外扩散出来,那么你这个球是个怪物了,它一直向外扩散有没有耗尽的时候?或者说它是个黑洞了,可以无限制向外 ...

球体外部浓度很低,近似为零,球体的物质浓度较高,物质的扩散依靠浓度差,浓度差使物质源源不断地往外扩散,可以使得物质从球体外表面R释放出来的,在一定条件下,在很长一段时间内能够实现恒速或者接近恒速的,当然当时间长到一定值后释放速率会降低的
7楼2010-01-25 10:23:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xinhero

银虫 (正式写手)

★ ★ ★
zju-hjm(金币+1): 1-25 18:57
zju-hjm(金币+2): 1-25 18:57
zju-hjm(金币+1): 2010-03-31 10:21
楼主的问题可以求解。实际上是个反问题,可以用数值方法求解。
现在的条件是偏微分方程已知,边界条件和初始条件:y(0,t)=f1,(dy/dr)在r=R处为定值,初始条件y(r,0)=f2,另外还必须有一个附加条件,比如y(r,t1)=f3.(这些条件楼主没有给全,姑且认为是已知的)。
这样可以通过Tikhonov正则化和函数逼近的方法求解。用Matlab编个程序算一下就搞定了。
Tikhonov正则化方法可以参考相关资料,主要思想是通过构造一个正则化算子,把问题化为正则泛函的极值问题,通过施加一个绕动,对f(r)进行叠加,最后得到得f(r)为基函数族的线性叠加。在一定精度要求下与理论值能够很好符合。
8楼2010-01-25 16:08:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zju-hjm

金虫 (初入文坛)

引用回帖:
Originally posted by xinhero at 2010-01-25 16:08:49:
楼主的问题可以求解。实际上是个反问题,可以用数值方法求解。
现在的条件是偏微分方程已知,边界条件和初始条件:y(0,t)=f1,(dy/dr)在r=R处为定值,初始条件y(r,0)=f2,另外还必须有一个附加条件,比如y(r,t1)=f ...

谢谢 虫友!!

那 Tikhonov正则化方法与函数逼近法具体怎样在matlab里编程实现呢?

我可以把问题写的清楚些
假定
初始条件 y(r,0)=1;
边界条件 dy/dr=0  在 t>=0,r=0;
         y(r,t)=0  在 t>=0,r=R=1;

如何求反求这个 f(r) ?

[ Last edited by zju-hjm on 2010-1-25 at 19:29 ]
9楼2010-01-25 19:09:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xinhero

银虫 (正式写手)

zju-hjm(金币+1): 2010-03-31 10:21
引用回帖:
Originally posted by zju-hjm at 2010-01-25 19:09:10:

谢谢 虫友!!

那 Tikhonov正则化方法与函数逼近法具体怎样在matlab里编程实现呢?

我可以把问题写的清楚些
假定
初始条件 y(r,0)=1;
边界条件 dy/dr=0  在 t>=0,r=0;
         y(r,t)=0  在 t ...

还差一个附加条件,比如y(r,t1)=f(x),即y在某时刻t1的情况.因为这个方程肯定有物理背景,y也肯定代表实际的物理涵义,比如浓度等。某时刻的y值一般都是可测的,所以一般可以知道。。你要告诉我一个这样的附加条件,我花点时间把程序写出来。
关于Tikhonov正则化方法,这里可能讲不清楚,相对比较复杂,你可以去查查相关文献。

[ Last edited by xinhero on 2010-1-25 at 22:11 ]
10楼2010-01-25 22:06:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zju-hjm 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 欢迎申博同学联系 +3 天道酬勤2026686 2026-03-10 7/350 2026-03-15 19:03 by 天道酬勤2026686
[考研] 304求调剂 +5 小熊joy 2026-03-14 5/250 2026-03-14 21:07 by peike
[考研] 一志愿天津大学,英一数二305分求调剂,四六级已过 +8 小小番的茄 2026-03-09 8/400 2026-03-14 01:53 by JourneyLucky
[考研] 云南财经大学信息学院计算机学硕专硕学位点 +3 zjptai 2026-03-10 5/250 2026-03-14 01:23 by 飞行琦
[考研] 0703求调剂 +7 jtyq001 2026-03-10 7/350 2026-03-14 01:06 by JourneyLucky
[考研] 一志愿安徽大学材料工程专硕313分,求调剂的学校 +8 Yu先生 2026-03-10 10/500 2026-03-14 01:04 by JourneyLucky
[考研] 材料工程专硕,一志愿中国矿业大学,总分314,求调剂 +5 无懈可击的巨人 2026-03-10 5/250 2026-03-14 00:37 by JourneyLucky
[考研] 311求调剂 +5 牛乳糖的卡卡 2026-03-10 5/250 2026-03-14 00:05 by JourneyLucky
[考研] 341求调剂 +4 番茄头--- 2026-03-10 4/200 2026-03-13 23:12 by JourneyLucky
[考研] 求调剂(材料与化工327) +4 爱吃香菜啦 2026-03-11 4/200 2026-03-13 22:11 by JourneyLucky
[考研] 一志愿西南交大,材料专硕317求调剂 +5 lx8568 2026-03-11 5/250 2026-03-13 21:43 by peike
[考研] 26调剂/材料科学与工程/总分295/求收留 +9 2026调剂侠 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
[考研] 301求调剂 +6 Liyouyumairs 2026-03-11 6/300 2026-03-13 20:11 by JourneyLucky
[考研] 310求调剂 +3 【上上签】 2026-03-11 3/150 2026-03-13 16:16 by JourneyLucky
[考研] 302求调剂 +6 负心者当诛 2026-03-11 6/300 2026-03-13 16:11 by JourneyLucky
[考研] 307求调剂 +5 超级伊昂大王 2026-03-12 5/250 2026-03-13 15:56 by 棒棒球手
[考研] 290求调剂 +7 ADT 2026-03-12 7/350 2026-03-13 15:17 by JourneyLucky
[考研] 274求调剂 +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
[考研] 材料301分求调剂 +5 Liyouyumairs 2026-03-12 5/250 2026-03-13 14:42 by JourneyLucky
[考研] 0857环境调剂 +5 熠熠_11 2026-03-10 5/250 2026-03-11 10:59 by wang_dand
信息提示
请填处理意见