北京石油化工学院2026年研究生招生接收调剂公告
查看: 2785  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 296求调剂 +6 彼岸t 2026-03-29 6/300 2026-03-29 20:22 by 永字号
[考研] 一志愿郑州大学,080500学硕,总分317分求调剂 +8 举个栗子oi 2026-03-24 9/450 2026-03-29 13:08 by peike
[考研] 315求调剂 +4 akie... 2026-03-28 5/250 2026-03-28 21:05 by zhq0425
[考研] 311(085601)求调剂 +4 liziyeyeye 2026-03-28 4/200 2026-03-28 18:50 by 535743368
[考研] 347求调剂 +3 山顶见α 2026-03-25 3/150 2026-03-28 14:13 by 唐沐儿
[考研] 286求调剂 +4 丢掉懒惰 2026-03-27 7/350 2026-03-28 08:07 by baoball
[考研] 0703化学/290求调剂/本科经历丰富/工科也可 +9 丹青奶盖 2026-03-26 10/500 2026-03-28 07:45 by barnett0632
[考研] 0703化学求调剂,各位老师看看我!!! +5 祁祺祺 2026-03-25 5/250 2026-03-27 21:44 by 东方猪猪
[考研] 化学调剂 +4 爱吃番茄的旭 2026-03-24 5/250 2026-03-27 17:50 by kiokin
[考研] 305求调剂 +5 哇卢卡库 2026-03-26 5/250 2026-03-27 14:01 by laoshidan
[硕博家园] 北京林业大学硕导招生广告 +6 kongweilin 2026-03-26 8/400 2026-03-27 10:18 by FF_16
[考研] 333求调剂 +6 wfh030413@ 2026-03-23 6/300 2026-03-26 22:45 by 学员8dgXkO
[考研] 调剂 +4 柚柚yoyo 2026-03-26 4/200 2026-03-26 20:43 by fmesaito
[考研] 289求调剂 +17 硕星赴 2026-03-23 17/850 2026-03-26 16:18 by 不吃魚的貓
[考研] 一志愿河工大 081700 276求调剂 +4 地球绕着太阳转 2026-03-23 4/200 2026-03-26 14:27 by zzll406
[考研] 化学调剂一志愿上海交通大学336分-本科上海211 +4 小鱼爱有机 2026-03-25 4/200 2026-03-26 10:19 by aa331100
[考研] 0854电子信息求调剂 324 +4 Promise-jyl 2026-03-23 4/200 2026-03-25 11:36 by Sugarlight
[考研] 300分,材料,求调剂,英一数二 +5 超赞的 2026-03-24 5/250 2026-03-24 21:07 by 星空星月
[考研] 材料专硕找调剂 +5 哈哈哈吼吼吼哈 2026-03-23 5/250 2026-03-24 19:07 by 了了了了。。
[考研] 材料专硕331求调剂 +4 鲜当牛 2026-03-24 4/200 2026-03-24 15:58 by JourneyLucky
信息提示
请填处理意见