| 24小時(shí)熱門(mén)版塊排行榜 |
| 查看: 2743 | 回復(fù): 22 | |||
zju-hjm金蟲(chóng) (初入文壇)
|
[交流]
【求助】求問(wèn)一個(gè)尋找最優(yōu)解問(wèn)題,是泛函問(wèn)題么 已有6人參與
|
|
從一個(gè)二階的偏微分方程出發(fā) dy/dt=d(f(r)*(dy/dr))/dr 其中t與r為自變量; y=y(t,r),f(r)為r的函數(shù); dy/dt,dy/dr為一階導(dǎo)數(shù); 等式右邊為f(r)*(dy/dr)對(duì)r的導(dǎo)數(shù); 0<=r<=R,R為常數(shù); 現(xiàn)在要使得(dy/dr)在r=R處為定值,怎樣確認(rèn)f(r)達(dá)到這種結(jié)果,也就是求合適的f(r),這是一個(gè)什么問(wèn)題呢? 怎么最優(yōu)化f(r),這是泛函問(wèn)題么? 謝謝幫忙啊 [ Last edited by wuguocheng on 2010-3-31 at 10:12 ] |
金蟲(chóng) (正式寫(xiě)手)
鐵桿木蟲(chóng) (職業(yè)作家)
科苑小木蟲(chóng)

鐵桿木蟲(chóng) (職業(yè)作家)
科苑小木蟲(chóng)

金蟲(chóng) (初入文壇)
鐵桿木蟲(chóng) (職業(yè)作家)
科苑小木蟲(chóng)
|
呵呵,這個(gè)我想繼續(xù)問(wèn)下去。假設(shè)你現(xiàn)在已經(jīng)找到了這個(gè)f(r),就是說(shuō)你能夠讓球形里面的物質(zhì)不斷向外擴(kuò)散出來(lái),那么你這個(gè)球是個(gè)怪物了,它一直向外擴(kuò)散有沒(méi)有耗盡的時(shí)候?或者說(shuō)它是個(gè)黑洞了,可以無(wú)限制向外擴(kuò)散?所以我覺(jué)得你的這個(gè)一直向外以恒定速率擴(kuò)散的要求有點(diǎn)不合理。 |

金蟲(chóng) (初入文壇)
|
球體外部濃度很低,近似為零,球體的物質(zhì)濃度較高,物質(zhì)的擴(kuò)散依靠濃度差,濃度差使物質(zhì)源源不斷地往外擴(kuò)散,可以使得物質(zhì)從球體外表面R釋放出來(lái)的,在一定條件下,在很長(zhǎng)一段時(shí)間內(nèi)能夠?qū)崿F(xiàn)恒速或者接近恒速的,當(dāng)然當(dāng)時(shí)間長(zhǎng)到一定值后釋放速率會(huì)降低的 |
|
樓主的問(wèn)題可以求解。實(shí)際上是個(gè)反問(wèn)題,可以用數(shù)值方法求解。 現(xiàn)在的條件是偏微分方程已知,邊界條件和初始條件:y(0,t)=f1,(dy/dr)在r=R處為定值,初始條件y(r,0)=f2,另外還必須有一個(gè)附加條件,比如y(r,t1)=f3.(這些條件樓主沒(méi)有給全,姑且認(rèn)為是已知的)。 這樣可以通過(guò)Tikhonov正則化和函數(shù)逼近的方法求解。用Matlab編個(gè)程序算一下就搞定了。 Tikhonov正則化方法可以參考相關(guān)資料,主要思想是通過(guò)構(gòu)造一個(gè)正則化算子,把問(wèn)題化為正則泛函的極值問(wèn)題,通過(guò)施加一個(gè)繞動(dòng),對(duì)f(r)進(jìn)行疊加,最后得到得f(r)為基函數(shù)族的線性疊加。在一定精度要求下與理論值能夠很好符合。 |
金蟲(chóng) (初入文壇)
|
還差一個(gè)附加條件,比如y(r,t1)=f(x),即y在某時(shí)刻t1的情況.因?yàn)檫@個(gè)方程肯定有物理背景,y也肯定代表實(shí)際的物理涵義,比如濃度等。某時(shí)刻的y值一般都是可測(cè)的,所以一般可以知道。。你要告訴我一個(gè)這樣的附加條件,我花點(diǎn)時(shí)間把程序?qū)懗鰜?lái)。 關(guān)于Tikhonov正則化方法,這里可能講不清楚,相對(duì)比較復(fù)雜,你可以去查查相關(guān)文獻(xiàn)。 [ Last edited by xinhero on 2010-1-25 at 22:11 ] |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學(xué)336分求調(diào)劑 +4 | zbzihdhd 2026-03-15 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 302求調(diào)劑 +9 | 負(fù)心者當(dāng)誅 2026-03-11 | 9/450 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +3 | 困于星晨 2026-03-17 | 3/150 |
|
|
[考研] 271求調(diào)劑 +12 | 生如夏花… 2026-03-11 | 14/700 |
|
|
[碩博家園] 深圳大學(xué)碩士招生(2026秋,傳感器方向,僅錄取第一志愿) +4 | xujiaoszu 2026-03-11 | 9/450 |
|
|
[考研] 一志愿985,本科211,0817化學(xué)工程與技術(shù)319求調(diào)劑 +5 | Liwangman 2026-03-15 | 5/250 |
|
|
[考研] 311求調(diào)劑 +5 | 26研0 2026-03-15 | 5/250 |
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
[考研]
|
笨笨兔子 2026-03-12 | 3/150 |
|
|
[考博] 東華理工大學(xué)化材專(zhuān)業(yè)26屆碩士博士申請(qǐng) +6 | zlingli 2026-03-13 | 6/300 |
|
|
[考研] 材料工程327求調(diào)劑 +3 | xiaohe12w 2026-03-11 | 3/150 |
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +6 | 邱gl 2026-03-12 | 7/350 |
|
|
[考研] 341求調(diào)劑 +4 | 番茄頭--- 2026-03-10 | 4/200 |
|
|
[考研] 301求調(diào)劑 +6 | Liyouyumairs 2026-03-11 | 6/300 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +7 | 綠豆芹菜湯 2026-03-12 | 7/350 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-12 | 10/500 |
|
|
[考研] 求調(diào)劑 +3 | 程雨杭 2026-03-12 | 3/150 |
|
|
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3 | T123 tt 2026-03-12 | 3/150 |
|