zhangluyoa

銀蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 40
帖子: 41
在線: 14.2小時
蟲號: 867634
注冊: 2009-10-10
性別: GG
專業(yè): 食品加工技術(shù)

[求助] 求解擴散方程問題

有一個方程如圖
源碼是
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ù)此樓

» 猜你喜歡

材料求調(diào)劑已經(jīng)有5人回復(fù)
354求調(diào)劑已經(jīng)有9人回復(fù)
化學(xué)工程321分求調(diào)劑已經(jīng)有22人回復(fù)
考研調(diào)劑已經(jīng)有3人回復(fù)
326求調(diào)劑已經(jīng)有8人回復(fù)
333求調(diào)劑已經(jīng)有7人回復(fù)
一志愿東華大學(xué)控制學(xué)碩320求調(diào)劑已經(jīng)有3人回復(fù)
考研化學(xué)學(xué)碩調(diào)劑，一志愿985 已經(jīng)有7人回復(fù)
0703化學(xué)調(diào)劑已經(jīng)有15人回復(fù)
0703化學(xué)調(diào)劑，六級已過，有科研經(jīng)歷已經(jīng)有14人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

求助comsol中各物理場的調(diào)整參數(shù)的意義。已經(jīng)有5人回復(fù)
關(guān)于菲克第二定律中的數(shù)學(xué)求解問題！已經(jīng)有4人回復(fù)
焊接中的擴散問題已經(jīng)有4人回復(fù)
對流-擴散方程中加入導(dǎo)熱方程的問題已經(jīng)有3人回復(fù)
【500金幣本版活動】“經(jīng)驗力場分子模擬+第一性原理計算問題討論互助及資源共享” 已經(jīng)有82人回復(fù)
【求助】求助拉格朗日方法追蹤顆粒擴散已經(jīng)有8人回復(fù)
【求助】MATLAB 有限差分法（FDM）求解偏微分方程已經(jīng)有22人回復(fù)
【求助】不僅控制方程里含有未知參數(shù)，邊界條件里也有已經(jīng)有5人回復(fù)
【求助】高維PDE的求解已經(jīng)有11人回復(fù)
【求助】求問一個尋找最優(yōu)解問題，是泛函問題么已經(jīng)有22人回復(fù)

1樓 2012-04-02 23:43:14

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

ya634984094

銅蟲 (小有名氣)

應(yīng)助: 39 (小學(xué)生)
金幣: 74.7
散金: 4
紅花: 3
帖子: 259
在線: 38.5小時
蟲號: 1695458
注冊: 2012-03-16
專業(yè): 流體力學(xué)

感謝參與，應(yīng)助指數(shù) +1
臭水溝: 應(yīng)助指數(shù)-1, 非應(yīng)助貼請不要選擇“應(yīng)助回帖”~ 2012-04-03 14:47:18

你這是用uds編的輸出嗎？？

贊一下

回復(fù)此樓

2樓2012-04-03 10:54:59

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhangluyoa

銀蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 40
帖子: 41
在線: 14.2小時
蟲號: 867634
注冊: 2009-10-10
性別: GG
專業(yè): 食品加工技術(shù)

引用回帖:

2樓: Originally posted by ya634984094 at 2012-04-03 10:54:59:
你這是用uds編的輸出嗎？？

什么是ubs啊，用的是matlab，我是菜鳥

贊一下

回復(fù)此樓

3樓2012-04-03 11:35:39

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

ya634984094

銅蟲 (小有名氣)

應(yīng)助: 39 (小學(xué)生)
金幣: 74.7
散金: 4
紅花: 3
帖子: 259
在線: 38.5小時
蟲號: 1695458
注冊: 2012-03-16
專業(yè): 流體力學(xué)

【答案】應(yīng)助回帖

★
臭水溝: 金幣+1, 謝謝交流~~ 2012-04-04 11:15:14

引用回帖:

3樓: Originally posted by zhangluyoa at 2012-04-03 11:35:39:
什么是ubs啊，用的是matlab，我是菜鳥

呃~~matlab我不是很熟，uds是fluent用來求解方程的~~

贊一下

回復(fù)此樓

4樓2012-04-03 12:34:11

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhangluyoa

銀蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 40
帖子: 41
在線: 14.2小時
蟲號: 867634
注冊: 2009-10-10
性別: GG
專業(yè): 食品加工技術(shù)

我就是想求D值，但是不知道為什么給賦D0值是一樣的

贊一下

回復(fù)此樓

5樓2012-04-03 14:58:46

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dbb627

榮譽版主 (著名寫手)

專家經(jīng)驗: +4
仿真EPI: 6
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計算模擬

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
臭水溝: 金幣+3, 謝謝交流~~ 2012-04-04 11:15:22
zhangluyoa: 金幣+4, ★有幫助 2012-04-11 16:35:41

CODE:

D0=3.55e-6;

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:5000;

s=sum(8./(pi*(2*k+1)).^2.*exp(-D*(pi*(2*k+1)).^2[color=red]*t[/color]/(d^2)));

F=sum((m-s').^2);

Local minimum possible.

lsqnonlin stopped because the final change in the sum of squares relative to
its initial value is less than the default value of the function tolerance.

D =

  5.0774e-008

resnorm =

  1.1370e-007

residual =

  3.3720e-004

exitflag =

   3

output =

firstorderopt: 39.0251
   iterations: 13
      funcCount: 28
   cgiterations: 0
      algorithm: 'large-scale: trust-region reflective Newton'
      message: [1x458 char]

lambda =

lower: 0
upper: 0

jacobian =

(1,1)    1.1573e+005

贊一下

回復(fù)此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

6樓2012-04-03 17:56:51

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dbb627

榮譽版主 (著名寫手)

專家經(jīng)驗: +4
仿真EPI: 6
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計算模擬

【答案】應(yīng)助回帖

★ ★ ★
臭水溝: 金幣+1, 謝謝交流~~ 2012-04-04 11:15:29
zhangluyoa: 金幣+2, ★有幫助 2012-04-11 16:35:57

s=sum(8./(pi*(2*k+1)).^2.*exp(-D*(pi*(2*k+1)).^2*t*60/(d^2))); 多乘了60？

贊一下

回復(fù)此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

7樓2012-04-03 17:58:17

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhangluyoa

銀蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 40
帖子: 41
在線: 14.2小時
蟲號: 867634
注冊: 2009-10-10
性別: GG
專業(yè): 食品加工技術(shù)

引用回帖:

7樓: Originally posted by dbb627 at 2012-04-03 17:58:17:
s=sum(8./(pi*(2*k+1)).^2.*exp(-D*(pi*(2*k+1)).^2*t*60/(d^2))); 多乘了60？

不是，t是分鐘，乘以60變成秒

贊一下

回復(fù)此樓

8樓2012-04-04 19:08:11

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

luyao_zhang

銀蟲 (初入文壇)

應(yīng)助: 3 (幼兒園)
金幣: 222.2
帖子: 11
在線: 11小時
蟲號: 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

精度太小，已經(jīng)小于計算機優(yōu)化精度，所以優(yōu)化計算停止，把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

贊一下

回復(fù)此樓

熱愛痛苦

9樓2013-07-16 17:34:30

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

1657956529

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 14.6
紅花: 2
帖子: 23
在線: 47.3小時
蟲號: 3806526
注冊: 2015-04-13
專業(yè): 食品科學(xué)

這是解的什么方程，我也想救助一下，fick第二定律怎么求解方程

贊一下

回復(fù)此樓

10樓2016-11-22 18:58:42

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 化學(xué)工程321分求調(diào)劑 +18	大米飯！ 2026-03-15	22/1100	2026-03-21 20:20 by HH領(lǐng)袖
[考研] 313求調(diào)劑 +4	肆叁貳壹22 2026-03-19	4/200	2026-03-21 17:33 by ColorlessPI
[考研] 生物學(xué)一志愿985，分?jǐn)?shù)349求調(diào)劑 +3	zxts12 2026-03-21	3/150	2026-03-21 16:34 by 33來了真來了
[考研] 機械專碩299求調(diào)劑至材料 +3	kkcoco25 2026-03-16	4/200	2026-03-21 03:52 by JourneyLucky
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +12	yangfz 2026-03-17	12/600	2026-03-21 03:30 by JourneyLucky
[考研] 307求調(diào)劑 +3	wyyyqx 2026-03-17	3/150	2026-03-21 03:20 by JourneyLucky
[考研] 271材料工程求調(diào)劑 +8	.6lL 2026-03-18	8/400	2026-03-21 00:58 by JourneyLucky
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境專碩，總分308求調(diào)劑 +3	墨墨漠 2026-03-18	3/150	2026-03-21 00:39 by JourneyLucky
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +9	Ncdx123456 2026-03-19	9/450	2026-03-20 23:41 by lovewei0727
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一） +5	sbdksD 2026-03-19	5/250	2026-03-20 22:10 by luoyongfeng
[考研] 0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章 +22	rare12345 2026-03-18	22/1100	2026-03-20 20:39 by zhukairuo
[考研] 一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑 +3	困困困困坤坤 2026-03-20	3/150	2026-03-20 20:38 by 學(xué)員8dgXkO
[考研] 319求調(diào)劑 +3	小力氣珂珂 2026-03-20	3/150	2026-03-20 19:47 by JourneyLucky
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +9	手機用戶 2026-03-17	9/450	2026-03-20 14:24 by 無懈可擊111
[考研] 材料學(xué)碩318求調(diào)劑 +5	February_Feb 2026-03-19	5/250	2026-03-19 23:51 by 23Postgrad
[考研] 材料與化工求調(diào)劑 +7	為學(xué)666 2026-03-16	7/350	2026-03-19 14:48 by 盡舜堯1
[考研] 【同濟軟件】軟件（085405）考研求調(diào)劑 +3	2026eternal 2026-03-18	3/150	2026-03-18 19:09 by 搏擊518
[考博] 26博士申請 +3	1042136743 2026-03-17	3/150	2026-03-17 23:30 by 輕松不少隨
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級食品/生物方向碩士 +3	1喜春8 2026-03-17	5/250	2026-03-17 17:18 by ber川cool子
[考研] 0856專碩279求調(diào)劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015