 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 470 | 回復(fù): 2 | |||
wangjunlitylg金蟲 (小有名氣)
|
[求助]
參數(shù)擬合：能運(yùn)行但初值不進(jìn)行迭代，輸入什么值結(jié)果還是什么值
|
|
能運(yùn)行的程序如下： function wenxian1 global k0 T s1 data=xlsread('D:文獻(xiàn).xlsx');T = data(:,1)';v= data(:,2); E0=232000;k0=8e12;kam=7130;s1=10/6;v0=0; beta0=[E0 kam]; [beta,resnorm,resid,exitflag,output,lambda] = lsqnonlin(@OptObjFunc,beta0,[],[],[],T,v0,v) vg=KineticsEqs(beta); figure,plot(T,v,'b-',T,vg,'k-') legend('Exp','Model') function f = OptObjFunc(beta,T,v0,v) vg = KineticsEqs(beta); f =v - vg'; function vg = KineticsEqs(beta) global a b c T s1 k0 d a=2*beta(2)^2; b=-k0*8.314/s1; c = 1/(beta(2)*sqrt(2*pi)); d=beta(1) ss=beta(1)+200*beta(2); xx=0; E=xx:ss; vg= quadv(@func,xx,ss); function y=func(E) global k0 T s1 a b c d y=k0./s1.*exp(-E./(8.314.*T)+b.*T.^2./E.*exp(-E./(8.314.*T))).*(c.*exp(-(E-d).^2./a)); 運(yùn)行結(jié)果：還是原來附的初值，且模擬和實(shí)驗(yàn)差別很大，即好像不進(jìn)行擬合， 我很努力了但就解決不了，謝謝大家，如果解決很感謝！ |
回復(fù)此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : 文獻(xiàn).xlsx
2016-04-02 11:20:32, 10.81 K
» 猜你喜歡
招收調(diào)劑碩士
已經(jīng)有12人回復(fù)
-
296求調(diào)劑
已經(jīng)有8人回復(fù)
-
招收博士1-2人
已經(jīng)有3人回復(fù)
-
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有11人回復(fù)
-
279分求調(diào)劑 一志愿211
已經(jīng)有9人回復(fù)
-
274求調(diào)劑
已經(jīng)有8人回復(fù)
-
290求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑
已經(jīng)有4人回復(fù)
-
295復(fù)試調(diào)劑
已經(jīng)有6人回復(fù)
-
工科材料085601 279求調(diào)劑
已經(jīng)有8人回復(fù)
月只藍(lán)
主管區(qū)長 (職業(yè)作家)
-
專家經(jīng)驗(yàn): +1059 - 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68127.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時(shí)
- 蟲號(hào): 1122189
- 注冊(cè): 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計(jì)算模擬區(qū)
2樓2016-04-02 11:32:26
wangjunlitylg
金蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 1048
- 紅花: 2
- 帖子: 91
- 在線: 47.7小時(shí)
- 蟲號(hào): 337417
- 注冊(cè): 2007-04-02
- 性別: MM
- 專業(yè): 傳遞過程
3樓2016-04-03 16:44:30
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 wangjunlitylg 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|
|---|---|---|---|---|
|
|
[考研] 材料080500調(diào)劑求收留 +6 | 一顆meteor 2026-03-13 | 6/300 |
|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +8 | 手機(jī)用戶 2026-03-17 | 8/400 |
|
|
|
[考研] 081700化工學(xué)碩調(diào)劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|
|
[考研] 288求調(diào)劑 +15 | 于海海海海 2026-03-19 | 15/750 |
|
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
|
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|
|
[考研] 344求調(diào)劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
|
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
|
[考研] 331求調(diào)劑（0703有機(jī)化學(xué) +7 | ZY-05 2026-03-13 | 8/400 |
|
|
|
[考研] 304求調(diào)劑 +12 | 小熊joy 2026-03-14 | 13/650 |
|
|
|
[考研] 303求調(diào)劑 +4 | 睿08 2026-03-17 | 6/300 |
|
|
|
[基金申請(qǐng)] 被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
|
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15 | 7/350 |
|
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
|
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼？~+123456 2026-03-14 | 3/150 |
|