| 查看: 2741 | 回復(fù): 22 | |||
zju-hjm金蟲 (初入文壇)
|
[交流]
【求助】求問一個(gè)尋找最優(yōu)解問題,是泛函問題么 已有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è)什么問題呢? 怎么最優(yōu)化f(r),這是泛函問題么? 謝謝幫忙啊 [ Last edited by wuguocheng on 2010-3-31 at 10:12 ] |
金蟲 (正式寫手)
鐵桿木蟲 (職業(yè)作家)
科苑小木蟲

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

金蟲 (初入文壇)
鐵桿木蟲 (職業(yè)作家)
科苑小木蟲

金蟲 (初入文壇)
|
球體外部濃度很低,近似為零,球體的物質(zhì)濃度較高,物質(zhì)的擴(kuò)散依靠濃度差,濃度差使物質(zhì)源源不斷地往外擴(kuò)散,可以使得物質(zhì)從球體外表面R釋放出來的,在一定條件下,在很長(zhǎng)一段時(shí)間內(nèi)能夠?qū)崿F(xiàn)恒速或者接近恒速的,當(dāng)然當(dāng)時(shí)間長(zhǎng)到一定值后釋放速率會(huì)降低的 |
|
樓主的問題可以求解。實(shí)際上是個(gè)反問題,可以用數(shù)值方法求解。 現(xiàn)在的條件是偏微分方程已知,邊界條件和初始條件:y(0,t)=f1,(dy/dr)在r=R處為定值,初始條件y(r,0)=f2,另外還必須有一個(gè)附加條件,比如y(r,t1)=f3.(這些條件樓主沒有給全,姑且認(rèn)為是已知的)。 這樣可以通過Tikhonov正則化和函數(shù)逼近的方法求解。用Matlab編個(gè)程序算一下就搞定了。 Tikhonov正則化方法可以參考相關(guān)資料,主要思想是通過構(gòu)造一個(gè)正則化算子,把問題化為正則泛函的極值問題,通過施加一個(gè)繞動(dòng),對(duì)f(r)進(jìn)行疊加,最后得到得f(r)為基函數(shù)族的線性疊加。在一定精度要求下與理論值能夠很好符合。 |
金蟲 (初入文壇)
|
還差一個(gè)附加條件,比如y(r,t1)=f(x),即y在某時(shí)刻t1的情況.因?yàn)檫@個(gè)方程肯定有物理背景,y也肯定代表實(shí)際的物理涵義,比如濃度等。某時(shí)刻的y值一般都是可測(cè)的,所以一般可以知道。。你要告訴我一個(gè)這樣的附加條件,我花點(diǎn)時(shí)間把程序?qū)懗鰜怼?br /> 關(guān)于Tikhonov正則化方法,這里可能講不清楚,相對(duì)比較復(fù)雜,你可以去查查相關(guān)文獻(xiàn)。 [ Last edited by xinhero on 2010-1-25 at 22:11 ] |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 293求調(diào)劑 +7 | zjl的號(hào) 2026-03-16 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 312求調(diào)劑 +4 | 陌宸希 2026-03-16 | 5/250 |
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +8 | 一鴨鴨喲 2026-03-14 | 10/500 |
|
|
[考研] 一志愿蘇州大學(xué)材料工程(085601)專碩有科研經(jīng)歷三項(xiàng)國(guó)獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
[考研] 283求調(diào)劑 +3 | 聽風(fēng)就是雨; 2026-03-16 | 3/150 |
|
|
[考研] 286求調(diào)劑 +3 | lemonzzn 2026-03-16 | 5/250 |
|
|
[考研] 環(huán)境工程調(diào)劑 +6 | 大可digkids 2026-03-16 | 6/300 |
|
|
[基金申請(qǐng)]
今年的國(guó)基金是打分制嗎?
50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
[考研] 321求調(diào)劑 +5 | 大米飯! 2026-03-15 | 5/250 |
|
|
[考研] 085600材料與化工 求調(diào)劑 +13 | enenenhui 2026-03-13 | 14/700 |
|
|
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +7 | Ncdx123456 2026-03-13 | 8/400 |
|
|
[考研] 0703 物理化學(xué)調(diào)劑 +3 | 我可以上岸的對(duì)?/a> 2026-03-13 | 5/250 |
|
|
[考研] 326求調(diào)劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
[基金申請(qǐng)]
有必要更換申報(bào)口嗎
20+3
|
fannyamoy 2026-03-11 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 快樂的香蕉 2026-03-11 | 4/200 |
|
|
[考研] 0856材料與化工301求調(diào)劑 +5 | 奕束光 2026-03-13 | 5/250 |
|
|
[考研] 0703化學(xué)一志愿211 總分320求調(diào)劑 +5 | 瑪卡巴卡啊哈 2026-03-11 | 5/250 |
|
|
[考研] 四川大學(xué)085601材料工程專碩 初試294求調(diào)劑 +4 | 祝我們好在冬天 2026-03-11 | 4/200 |
|
|
[考研] 301求調(diào)劑 +6 | Liyouyumairs 2026-03-11 | 6/300 |
|
|
[考研] 290求調(diào)劑 +3 | ADT 2026-03-13 | 3/150 |
|