版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

feifei_ney

木蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1673.5
散金: 20
紅花: 2
帖子: 228
在線: 113.1小時
蟲號: 847106
注冊: 2009-09-13
性別: MM
專業(yè): 等離子體物理

[求助] 用Matlab求解方程

有個遷移模型方程

已知MF,∞=9.18mg,Lp=0.009,
MF,t與t的一組相關(guān)實驗數(shù)據(jù)如下所示
t MF,t
7200 3.28265
86400 5.57149
345600 6.84803
864000 7.47618
7.776E6 9.17744
想求解得到D值。
請大家?guī)兔η蠼庖幌�，不勝感激，希望能夠提供詳�?xì)的求解步驟。

回復(fù)此樓

» 猜你喜歡

307求調(diào)劑已經(jīng)有20人回復(fù)
285求調(diào)劑已經(jīng)有8人回復(fù)
材料與化工（0856）304求B區(qū)調(diào)劑已經(jīng)有8人回復(fù)
370求調(diào)劑已經(jīng)有3人回復(fù)
070300化學(xué)279求調(diào)劑已經(jīng)有5人回復(fù)
福建理工大學(xué)材料學(xué)院先進(jìn)合金團(tuán)隊招收考研調(diào)劑學(xué)生已經(jīng)有4人回復(fù)
英一數(shù)一408，總分284，二戰(zhàn)真誠求調(diào)劑已經(jīng)有5人回復(fù)
083000環(huán)境科學(xué)與工程調(diào)劑，總分281 已經(jīng)有4人回復(fù)
求調(diào)劑已經(jīng)有3人回復(fù)
食品工程專碩一志愿中海洋309求調(diào)劑已經(jīng)有10人回復(fù)

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

Matlab求解偏微分方程組已經(jīng)有22人回復(fù)
matlab求解非線性方程組已經(jīng)有16人回復(fù)
求matlab求解一個方程組的問題已經(jīng)有6人回復(fù)
求高人指點(diǎn)用matlab求解非線性方程組，解決了追加100金幣；已經(jīng)有11人回復(fù)
matlab的fsove 命令求解非線性方程組已經(jīng)有6人回復(fù)
【求助】matlab 求解微分方程中的未知參數(shù) 已經(jīng)有20人回復(fù)
【求助】向各位大俠求助matlab求解微分方程組遇到的一個問題已經(jīng)有21人回復(fù)
【求助】MATLAB 有限差分法（FDM）求解偏微分方程已經(jīng)有22人回復(fù)
【求助】積分微分方程matlab求解已經(jīng)有6人回復(fù)
【求助】matlab怎么求解偏微分方程組啊，先謝謝了已經(jīng)有13人回復(fù)
【資源】用matlab求解超越方程已經(jīng)有13人回復(fù)
【求助】求MAtlab求解程序�。�！已經(jīng)有20人回復(fù)

1樓 2011-11-08 08:41:40

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

wjx0162

銀蟲 (著名寫手)

應(yīng)助: 3 (幼兒園)
金幣: 468.9
紅花: 1
帖子: 1152
在線: 68.4小時
蟲號: 1187390
注冊: 2011-01-11
專業(yè): 機(jī)械動力學(xué)

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

★
臭水溝(金幣+1): 謝謝應(yīng)助~~ 2011-11-08 11:22:41

兩邊取對數(shù)，把指數(shù)方程變?yōu)榇鷶?shù)方程，直接求代數(shù)方程的根。

贊一下(1人)

回復(fù)此樓

2樓2011-11-08 09:03:25

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

feifei_ney

木蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1673.5
散金: 20
紅花: 2
帖子: 228
在線: 113.1小時
蟲號: 847106
注冊: 2009-09-13
性別: MM
專業(yè): 等離子體物理

引用回帖:

2樓: Originally posted by wjx0162 at 2011-11-08 09:03:25:
兩邊取對數(shù)，把指數(shù)方程變?yōu)榇鷶?shù)方程，直接求代數(shù)方程的根。

這是一個求和的公式，我給出的只是其中的幾個實驗參數(shù)，我希望能通過曲線擬合得到更準(zhǔn)確的結(jié)果。

贊一下

回復(fù)此樓

3樓2011-11-08 12:52:05

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

wjx0162

銀蟲 (著名寫手)

應(yīng)助: 3 (幼兒園)
金幣: 468.9
紅花: 1
帖子: 1152
在線: 68.4小時
蟲號: 1187390
注冊: 2011-01-11
專業(yè): 機(jī)械動力學(xué)

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

★
xiegangmai(金幣+1): 謝謝應(yīng)助！ 2011-11-08 16:23:29

那也一樣，兩邊取對數(shù)，先把指數(shù)方程變?yōu)榇鷶?shù)方程，然后用數(shù)據(jù)點(diǎn)擬合待定系數(shù)D的代數(shù)方程，最后求出D

贊一下(1人)

回復(fù)此樓

4樓2011-11-08 15:10:03

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

feifei_ney

木蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1673.5
散金: 20
紅花: 2
帖子: 228
在線: 113.1小時
蟲號: 847106
注冊: 2009-09-13
性別: MM
專業(yè): 等離子體物理

引用回帖:

4樓: Originally posted by wjx0162 at 2011-11-08 15:10:03:
那也一樣，兩邊取對數(shù)，先把指數(shù)方程變?yōu)榇鷶?shù)方程，然后用數(shù)據(jù)點(diǎn)擬合待定系數(shù)D的代數(shù)方程，最后求出D

能否麻煩您幫我解一下？我實在不知道怎么處理。

贊一下

回復(fù)此樓

5樓2011-11-08 17:25:47

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

dbb627

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

專家經(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)助回帖

★ ★
feifei_ney(金幣+10): 非常感謝您，希望能夠繼續(xù)幫我解答。我對Matlab軟件不熟悉，能否告訴我每步所輸入的函數(shù)的意義。還有個疑問，D0的值是自己定義的嗎？ 2011-11-08 20:14:31
xiegangmai(金幣+2): 謝謝應(yīng)助！ 2011-11-09 09:28:43

引用回帖:

1樓: Originally posted by feifei_ney at 2011-11-08 08:41:40:
有個遷移模型方程

已知MF,∞=9.18mg,Lp=0.009,
MF,t與t的一組相關(guān)實驗數(shù)據(jù)如下所示
t MF,t
7200 3.28265
86400 5.57149
345600 6.84803
864 ...

用最小二乘法解了下
不知道你看結(jié)果如何

CODE:

function AA

D0=3.55e-15;

%D0=3.55e-19;

% options = optimset('Algorithm','interior-point');

% [D,fval,flag] = fmincon(@fei,D0,[],[],[],[],0,[],[],options)

[D,resnorm,residual,exitflag,output,lambda,jacobian] = lsqnonlin(@fei,D0,0,[],[])

function F=fei(D)

Lp=0.009;

s=[];

A=[7200        3.28265

86400        5.57149

345600        6.84803

864000        7.47618

7.776e6       9.17744];

t=A(:,1);

MFt=A(:,2);

MF=9.18;

Lp=0.009;

y=1-MFt/MF;

k=0:100000;

for i=1:length(t)

s(i)=sum(8./(pi*(2*k+1)).^2.*exp(-D*(pi*(2*k+1)).^2*t(i)^2/(4*Lp^2)));

end

s

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

整個文件復(fù)制到editor中保存為m文件后運(yùn)行

D0區(qū)不同初值的結(jié)果
D0=3.55e-15
s =

0.9462 0.3616 0.0000 0.0000       0

s =

   0    0    0    0    0

Initial point is a local minimum.

Optimization completed because the size of the gradient at the initial point
is less than the default value of the function tolerance.

D =

  3.5500e-015

resnorm =

0.0370

residual =

0.1923

exitflag =

   1

output =

firstorderopt: 2.1708e-008
   iterations: 0
      funcCount: 2
   cgiterations: 0
      algorithm: 'large-scale: trust-region reflective Newton'
      message: [1x500 char]

lambda =

lower: 6.1148e+006
upper: 0

D0=3.55e-19
s =

0.9995 0.9935 0.9742 0.9355 0.4218

s =

   0    0    0    0    0

s =

   0    0    0    0    0

s =

   0    0    0    0    0

Local minimum found.

Optimization completed because the size of the gradient is less than
the default value of the function tolerance.

D =

  2.4089e-008

resnorm =

0.4438

residual =

0.6662

exitflag =

   1

output =

firstorderopt: 0
   iterations: 1
      funcCount: 4
   cgiterations: 0
      algorithm: 'large-scale: trust-region reflective Newton'
      message: [1x427 char]

lambda =

lower: 0
upper: 0

jacobian =

All zero sparse: 1-by-1

贊一下(2人)

回復(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樓2011-11-08 19:23:30

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

dbb627

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

專家經(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)助回帖

★ ★
xiegangmai(金幣+2): 謝謝應(yīng)助！ 2011-11-09 09:29:01

[quote]6樓: Originally posted by dbb627 at 2011-11-08 19:23:30:
用最小二乘法解了下
不知道你看結(jié)果如何

CODE:

function AA

D0=3.55e-15;

%D0=3.55e-19;

% options = optimset('Algorithm','interior-point');

% [D,fval,flag] = fmincon(@fei,D0,[],[],[],[],0,[],[] ... [/quote]

[code]function AA

D0=3.55e-15;%給需要計算的D一個迭代計算的初值

%D0=3.55e-19;%給需要計算的D一個其他迭代計算的初值

% options = optimset('Algorithm','interior-point');

% [D,fval,flag] = fmincon(@fei,D0,[],[],[],[],0,[],[],options)%去掉這三行的%將用另外算法fmincon求D

[D,resnorm,residual,exitflag,output,lambda,jacobian] = lsqnonlin(@fei,D0,0,[],[])%采用lsqnonlin算法求D

%====編寫函數(shù)關(guān)于D的方程======

function F=fei(D)

%===========賦值====

Lp=0.009;

s=[];

A=[7200        3.28265

86400        5.57149

345600        6.84803

864000        7.47618

7.776e6       9.17744];

t=A(:,1);

MFt=A(:,2);

MF=9.18;

Lp=0.009;

y=1-MFt/MF;

%=========計算無窮級數(shù)（用前100000項逼近）===

k=0:100000;

for i=1:length(t)

s(i)=sum(8./(pi*(2*k+1)).^2.*exp(-D*(pi*(2*k+1)).^2*t(i)^2/(4*Lp^2)));

end

%=======顯示無窮級數(shù)計算結(jié)果====

s

%====構(gòu)造需要優(yōu)化的函數(shù)（1-MFt/MF與無窮級數(shù)的差的平方和）=======

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

贊一下(1人)

回復(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樓2011-11-09 09:26:31

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 feifei_ney 的主題更新

返回列表

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

[考研] 材料科學(xué)與工程求調(diào)劑 +9	深V宿舍吧 2026-03-29	9/450	2026-03-31 00:15 by jp9609
[考研] 303求調(diào)劑 +7	DLkz1314. 2026-03-30	7/350	2026-03-30 21:07 by peike
[考研] 083000學(xué)碩274求調(diào)劑 +9	Li李魚 2026-03-26	9/450	2026-03-30 21:06 by dophin1985
[考研] 調(diào)劑 +4	GK72 2026-03-30	4/200	2026-03-30 20:32 by dick_runner
[考研] 309求調(diào)劑 +15	誰不是少年 2026-03-29	15/750	2026-03-30 19:28 by 2026材料調(diào)劑
[考研] 0703化學(xué)321分求調(diào)劑 +10	三dd. 2026-03-30	11/550	2026-03-30 19:24 by markhwc
[考研] 085600，專業(yè)課化工原理，320分求調(diào)劑 +5	大饞小子 2026-03-29	5/250	2026-03-30 19:02 by 源_2020
[考研] 26考研-291分-廈門大學(xué)（085601）-柔性電子學(xué)院材料工程專業(yè)求調(diào)劑 +5	min3 2026-03-24	6/300	2026-03-30 18:42 by 544594351
[考研] 化學(xué)308分調(diào)劑 +7	你好明天你好 2026-03-30	8/400	2026-03-30 13:02 by 你好明天你好
[考研] 一志愿中南大學(xué)化學(xué)0703總分337求調(diào)劑 +6	niko- 2026-03-27	6/300	2026-03-30 10:25 by herarysara
[考研] 337求調(diào)劑 +6	《樹》 2026-03-29	6/300	2026-03-30 10:15 by herarysara
[考研] 311求調(diào)劑 +6	冬十三 2026-03-24	6/300	2026-03-29 20:45 by 無際的草原
[考研] 考研調(diào)劑 +7	小蠟新筆 2026-03-29	7/350	2026-03-29 19:00 by 學(xué)員8dgXkO
[考研] 0856材料化工調(diào)劑總分330 +14	zhubinhao 2026-03-27	14/700	2026-03-29 10:01 by Sjndkwm
[考研] 266求調(diào)劑 +11	陽陽哇塞 2026-03-27	12/600	2026-03-27 17:56 by yu221
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +4	Ripcord上岸 2026-03-27	4/200	2026-03-27 17:00 by zhyzzh
[考研] 調(diào)劑 +4	柚柚yoyo 2026-03-26	4/200	2026-03-26 20:43 by fmesaito
[考研] 085602 289分求調(diào)劑 +8	WWW西西弗斯 2026-03-24	8/400	2026-03-26 16:33 by 不吃魚的貓
[考研] 07化學(xué)303求調(diào)劑 +5	睿08 2026-03-25	5/250	2026-03-25 22:46 by 418490947
[考研] 285求調(diào)劑 +3	AZMK 2026-03-24	3/150	2026-03-25 12:23 by userper

24小時熱門版塊排行榜

[求助] 用Matlab求解方程

» 猜你喜歡

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

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

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

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

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

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