zhangluyoa

銀蟲 (初入文壇)

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

[求助] 求解擴(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é)果怎么也不對(duì) 如下
得到的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)有4人回復(fù)
085600材料與化工306 已經(jīng)有4人回復(fù)
286求調(diào)劑已經(jīng)有10人回復(fù)
328求調(diào)劑，英語六級(jí)551，有科研經(jīng)歷已經(jīng)有9人回復(fù)
一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑已經(jīng)有4人回復(fù)
286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考！已經(jīng)有8人回復(fù)
資源與環(huán)境調(diào)劑申請(qǐng)(333分) 已經(jīng)有5人回復(fù)
280求調(diào)劑已經(jīng)有12人回復(fù)
269專碩求調(diào)劑已經(jīng)有5人回復(fù)
材料學(xué)碩301分求調(diào)劑已經(jīng)有7人回復(fù)

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

求助comsol中各物理場(chǎng)的調(diào)整參數(shù)的意義。已經(jīng)有5人回復(fù)
關(guān)于菲克第二定律中的數(shù)學(xué)求解問題！已經(jīng)有4人回復(fù)
焊接中的擴(kuò)散問題已經(jīng)有4人回復(fù)
對(duì)流-擴(kuò)散方程中加入導(dǎo)熱方程的問題已經(jīng)有3人回復(fù)
【500金幣本版活動(dòng)】“經(jīng)驗(yàn)力場(chǎng)分子模擬+第一性原理計(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ù)

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小時(shí)
蟲號(hào): 1695458
注冊(cè): 2012-03-16
專業(yè): 流體力學(xué)

感謝參與，應(yīng)助指數(shù) +1
臭水溝: 應(yīng)助指數(shù)-1, 非應(yīng)助貼請(qǐ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小時(shí)
蟲號(hào): 867634
注冊(cè): 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小時(shí)
蟲號(hào): 1695458
注冊(cè): 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小時(shí)
蟲號(hào): 867634
注冊(cè): 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

榮譽(yù)版主 (著名寫手)

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

【答案】應(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

榮譽(yù)版主 (著名寫手)

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

【答案】應(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小時(shí)
蟲號(hào): 867634
注冊(cè): 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小時(shí)
蟲號(hào): 2453215
注冊(cè): 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)小于計(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

贊一下

回復(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小時(shí)
蟲號(hào): 3806526
注冊(cè): 2015-04-13
專業(yè): 食品科學(xué)

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

贊一下

回復(fù)此樓

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

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

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

[考研] 求調(diào)劑院校信息 +4	CX 330 2026-03-21	4/200	2026-03-21 23:48 by ms629
[考研] 354求調(diào)劑 +6	Tyoumou 2026-03-18	9/450	2026-03-21 20:47 by lbsjt
[考研] 一志愿東華大學(xué)控制學(xué)碩320求調(diào)劑 +3	Grand777 2026-03-21	3/150	2026-03-21 19:23 by 簡(jiǎn)之-
[基金申請(qǐng)] 山東省面上項(xiàng)目限額評(píng)審 +3	石瑞0426 2026-03-19	3/150	2026-03-21 18:59 by RRRKKK
[考研] 317求調(diào)劑 +9	申子申申 2026-03-19	15/750	2026-03-21 17:31 by 學(xué)員8dgXkO
[考研] 生物學(xué)一志愿985，分?jǐn)?shù)349求調(diào)劑 +3	zxts12 2026-03-21	3/150	2026-03-21 16:34 by 33來了真來了
[考研] 279求調(diào)劑 +5	紅衣隱官 2026-03-21	5/250	2026-03-21 14:59 by lature00
[考研] 08工科 320總分求調(diào)劑 +6	梨花珞晚風(fēng) 2026-03-17	6/300	2026-03-21 03:40 by JourneyLucky
[考研] 265求調(diào)劑 +3	Jack?k?y 2026-03-17	3/150	2026-03-21 03:17 by JourneyLucky
[考研] 一志愿西南交大，求調(diào)劑 +5	材化逐夢(mèng)人 2026-03-18	5/250	2026-03-21 00:26 by JourneyLucky
[考研] 296求調(diào)劑 +6	www_q 2026-03-18	10/500	2026-03-20 23:56 by JourneyLucky
[考研] 材料專碩英一數(shù)二306 +7	z1z2z3879 2026-03-18	7/350	2026-03-20 23:48 by JourneyLucky
[考研] 287求調(diào)劑 +7	晨昏線與星海 2026-03-19	8/400	2026-03-20 22:19 by JourneyLucky
[考研] 290求調(diào)劑 +7	^O^乜 2026-03-19	7/350	2026-03-20 21:43 by JourneyLucky
[考研] 材料學(xué)求調(diào)劑 +4	Stella_Yao 2026-03-20	4/200	2026-03-20 20:28 by ms629
[考研] 一志愿南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +5	@taotao 2026-03-20	5/250	2026-03-20 20:16 by JourneyLucky
[考研] 085600材料與化工調(diào)劑 324分 +10	llllkkkhh 2026-03-18	12/600	2026-03-19 14:33 by llllkkkhh
[考研] 286求調(diào)劑 +6	lemonzzn 2026-03-16	10/500	2026-03-19 14:31 by lemonzzn
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +5	share_joy 2026-03-16	5/250	2026-03-18 14:48 by haxia
[考研] 材料專碩326求調(diào)劑 +6	墨煜姒莘 2026-03-15	7/350	2026-03-17 17:10 by ruiyingmiao