 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1778 | 回復(fù): 9 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
zhangluyoa銀蟲 (初入文壇)
|
[求助]
求解擴(kuò)散方程問題
|
||
|
有一個(gè)方程如圖 源碼是 function AA D0=3.55e-19; options=optimset('tolx',1e-9); [D,resnorm,residual,exitflag,output,lambda,jacobian] = lsqnonlin(@fei,D0,[],[],options) function F=fei(D) C0=0.687; d=0.02; Ceq=20; t=60; C=4.136; m=1-(C0-C)/(C0-Ceq); s=[]; k=0:1000; s=sum(8./(pi*(2*k+1)).^2.*exp(-D*(pi*(2*k+1)).^2*t*60/(d^2))); s F=sum((m-s').^2); 而得到的結(jié)果怎么也不對 如下 得到的D值怎么都是跟賦的D0值是一樣的額 s = 0.9998 s = 0.2158 s = Inf s = 0.2447 s = Inf s = 0.2161 Optimization terminated: norm of the current step is less than OPTIONS.TolX. D = 3.5500e-019 resnorm = 0.0010 residual = 0.0318 exitflag = 2 output = firstorderopt: 7.1546e+005 iterations: 2 funcCount: 6 cgiterations: 1 algorithm: 'large-scale: trust-region reflective Newton' message: [1x77 char] lambda = lower: 0 upper: 0 jacobian = (1,1) 2.2480e+007 到底怎么弄得，小弟是菜鳥，希望有高手跟好心人幫我5555 |
回復(fù)此樓» 猜你喜歡
【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專業(yè)調(diào)劑名額
已經(jīng)有7人回復(fù)
-
334分 一志愿武理-080500 材料求調(diào)劑
已經(jīng)有3人回復(fù)
-
0854電子信息求調(diào)劑
已經(jīng)有9人回復(fù)
-
347求調(diào)劑
已經(jīng)有4人回復(fù)
-
0854人工智能方向招收調(diào)劑
已經(jīng)有4人回復(fù)
-
求b區(qū)院校調(diào)劑
已經(jīng)有4人回復(fù)
-
材料與化工考研調(diào)劑
已經(jīng)有9人回復(fù)
-
求調(diào)劑
已經(jīng)有3人回復(fù)
-
總分293求調(diào)劑
已經(jīng)有4人回復(fù)
-
284求調(diào)劑
已經(jīng)有15人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- 求助comsol中各物理場的調(diào)整參數(shù)的意義。
已經(jīng)有5人回復(fù)
- 關(guān)于菲克第二定律中的數(shù)學(xué)求解問題！
已經(jīng)有4人回復(fù)
- 焊接中的擴(kuò)散問題
已經(jīng)有4人回復(fù)
- 對流-擴(kuò)散方程中加入導(dǎo)熱方程的問題
已經(jīng)有3人回復(fù)
- 【500金幣本版活動(dòng)】“經(jīng)驗(yàn)力場分子模擬+第一性原理計(jì)算問題討論互助及資源共享”
已經(jīng)有82人回復(fù)
- 【求助】求助拉格朗日方法追蹤顆粒擴(kuò)散
已經(jīng)有8人回復(fù)
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已經(jīng)有22人回復(fù)
- 【求助】不僅控制方程里含有未知參數(shù)，邊界條件里也有
已經(jīng)有5人回復(fù)
- 【求助】高維PDE的求解
已經(jīng)有11人回復(fù)
- 【求助】求問一個(gè)尋找最優(yōu)解問題，是泛函問題么
已經(jīng)有22人回復(fù)
luyao_zhang
銀蟲 (初入文壇)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 222.2
- 帖子: 11
- 在線: 11小時(shí)
- 蟲號: 2453215
- 注冊: 2013-05-07
- 性別: GG
- 專業(yè): 食品加工學(xué)基礎(chǔ)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
zhangluyoa: 金幣+44, ★★★★★最佳答案, 終于解決了 2013-07-16 17:35:01
xiegangmai: 金幣+2, 謝謝參與！ 2013-07-24 21:03:08
zhangluyoa: 金幣+44, ★★★★★最佳答案, 終于解決了 2013-07-16 17:35:01
xiegangmai: 金幣+2, 謝謝參與！ 2013-07-24 21:03:08
|
精度太小，已經(jīng)小于計(jì)算機(jī)優(yōu)化精度，所以優(yōu)化計(jì)算停止，把D 增大10^8倍，方程中減小10^8倍，可順利求解 D0=3.55; [D,resnorm,residual,exitflag,output,lambda,jacobian] = lsqnonlin(@fei,D0,[],[]) function F=fei(D) C0=0.687; d=0.02; Ceq=20; t=60; C=4.136; m=1-(C0-C)/(C0-Ceq); s=[]; k=0:1000; s=sum(8./(pi*(2*k+1)).^2.*exp(-D*0.00000001*(pi*(2*k+1)).^2*t*60/(d^2))); s F=sum((m-s').^2); 結(jié)果 D = 0.04933156266209 resnorm = 1.200260299201833e-008 D就等于4.93e-10 |
9樓2013-07-16 17:34:30
ya634984094
銅蟲 (小有名氣)
- 應(yīng)助: 39 (小學(xué)生)
- 金幣: 74.7
- 散金: 4
- 紅花: 3
- 帖子: 259
- 在線: 38.5小時(shí)
- 蟲號: 1695458
- 注冊: 2012-03-16
- 專業(yè): 流體力學(xué)
2樓2012-04-03 10:54:59
zhangluyoa
銀蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 40
- 帖子: 41
- 在線: 14.2小時(shí)
- 蟲號: 867634
- 注冊: 2009-10-10
- 性別: GG
- 專業(yè): 食品加工技術(shù)
3樓2012-04-03 11:35:39
ya634984094
銅蟲 (小有名氣)
- 應(yīng)助: 39 (小學(xué)生)
- 金幣: 74.7
- 散金: 4
- 紅花: 3
- 帖子: 259
- 在線: 38.5小時(shí)
- 蟲號: 1695458
- 注冊: 2012-03-16
- 專業(yè): 流體力學(xué)
4樓2012-04-03 12:34:11
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專業(yè)調(diào)劑名額 +4 | 袁奐奐 2026-03-25 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 0854人工智能方向招收調(diào)劑 +4 | 章小魚567 2026-03-24 | 4/200 |
|
|
[考研] 359求調(diào)劑 +3 | 王了個(gè)楠 2026-03-25 | 3/150 |
|
|
|
[考研] B區(qū)考研調(diào)劑 +4 | yqdszhdap－ 2026-03-22 | 5/250 |
|
|
|
[考研] 08工學(xué)調(diào)劑 +12 | 用戶573181 2026-03-20 | 17/850 |
|
|
|
[考研] 一志愿南航材料專317分求調(diào)劑 +5 | 炸呀炸呀炸薯?xiàng)l 2026-03-23 | 5/250 |
|
|
|
[考研] 300求調(diào)劑，材料科學(xué)英一數(shù)二 +5 | leaflight 2026-03-24 | 5/250 |
|
|
|
[考研] 321求調(diào)劑 +4 | Ymlll 2026-03-24 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +7 | 十三加油 2026-03-21 | 7/350 |
|
|
|
[考研] 291求調(diào)劑 +8 | hhhhxn.. 2026-03-23 | 8/400 |
|
|
|
[考研] 工科0856求調(diào)劑 +5 | 沐析汀汀 2026-03-21 | 5/250 |
|
|
|
[考研] 323求調(diào)劑 +6 | 洼小桶 2026-03-18 | 6/300 |
|
|
|
[考研] 293求調(diào)劑 +3 | 濤濤Wjt 2026-03-22 | 5/250 |
|
|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
|
[考研] 求調(diào)劑院校信息 +6 | CX 330 2026-03-21 | 6/300 |
|
|
|
[考研] 296求調(diào)劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
|
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
|
[考研] 321求調(diào)劑 +9 | 何潤采123 2026-03-18 | 11/550 |
|
|
|
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|
|
|
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|