 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 286 | 回復(fù): 0 | |||
wodaifei銀蟲 (小有名氣)
|
[交流]
【求助】各位師哥師姐幫小弟檢查一下復(fù)雜參數(shù)回歸的程序吧��！小弟萬分感謝�。�！
|
|
各位師哥師姐，小弟求助： dy1/dt=k1y1y2-k2y3 dy2/dt=k1y1y2-k2y3 dy3/dt=k2y3 y1(0)=10,y2(0)=8,y3(0)=0 小弟想用四階-龍閣庫塔積分和非線性最小二乘估計(jì)以上的動(dòng)力學(xué)參數(shù)k1,k2,k3,數(shù)據(jù)如下：t y1 y2 y3 0 10 8 0 20 8 7 5 40 6 6 7 60 5 5 8 120 3 3 10 180 2 2 12 300 1 1 14 以下是小弟的程序： k0 = [1,2,3]; %設(shè)參數(shù)初值 y0 = [10,8,0]; tspan = [0,20,40,60,120,180,300]; %反應(yīng)時(shí)間 y1 = [10,8,6,5,3,2,1]';%測得的第一組數(shù)據(jù) y2 = [8,7,6,5,3,2,1]';%測得的第二組數(shù)據(jù) y3 = [0,5,7,8,10,12,14]';%測得的第三組數(shù)據(jù) yexp = [y1;y2;y3]; function dy = KineticEqs(t,y,k)%動(dòng)力學(xué)參數(shù)方程 dy = zeros(3,1); dy(1)=k1*y(1)*y(2)-k2*y(3); dy(2)=k1*y(1)*y(2)-k2*y(3); dy(3)=k2*y(3); function f =OptObjFunc(k,y0,yexp) %目標(biāo)函數(shù) [t Y] = ode45(@KineticEqs,tspan,y,[],k); f1 = y(1) - y1; f2 = y(2) - y2; f3 = y(3) - y3; f = [f1; f2; f3]; [k,resnorm,resid,exitflag,output,lambda,jacobian] = ...%最小二乘法 lsqnonlin(@OptObjFunc,k0,[0 0],[],[],tspan,y0) ci = nlparci(k,resid,jacobian)%置信區(qū)間 fprintf('Estimated Parameters:\n') fprintf('\tk1 = %.4f ± %.4f\n',k1) fprintf('\tk2 = %.4f ± %.4f\n',k2) 可是運(yùn)行不出來，小弟是新手，不知道錯(cuò)了哪里，希望各位高手幫我一下，小弟真的是急死了�。。。∽詈媚芙o出正確的程序，小弟剛?cè)胝搲痪�，金幣不多，但小弟愿拿出一半給大家一表自己的心意。謝謝了��！ [ Last edited by wodaifei on 2011-3-26 at 10:07 ] |
回復(fù)此樓» 猜你喜歡
284求調(diào)劑
已經(jīng)有10人回復(fù)
-
一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑
已經(jīng)有4人回復(fù)
-
07化學(xué)280分求調(diào)劑
已經(jīng)有4人回復(fù)
-
298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑
已經(jīng)有12人回復(fù)
-
求材料，環(huán)境專業(yè)調(diào)劑
已經(jīng)有3人回復(fù)
-
335求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿吉大化學(xué)322求調(diào)劑
已經(jīng)有4人回復(fù)
-
環(huán)境學(xué)碩288求調(diào)劑
已經(jīng)有8人回復(fù)
-
341求調(diào)劑(一志愿湖南大學(xué)070300)
已經(jīng)有6人回復(fù)
找到一些相關(guān)的精華帖子，希望有用哦~
SFDA認(rèn)證中心2010版藥品GMP認(rèn)證檢查結(jié)果判定程序
已經(jīng)有14人回復(fù)
- 求助格氏反應(yīng)，小弟感謝！��！
已經(jīng)有9人回復(fù)
- 小弟初學(xué)生物信息學(xué)，完全摸不著頭腦，求啟蒙
已經(jīng)有7人回復(fù)
- 關(guān)于XPS的幾個(gè)問題，小弟先再次謝過了
已經(jīng)有8人回復(fù)
- 求助申請文部省獎(jiǎng)學(xué)金的師兄師姐們幫小弟解解惑呀
已經(jīng)有11人回復(fù)
- 請質(zhì)譜高人幫忙解析質(zhì)譜，萬分感謝！
已經(jīng)有14人回復(fù)
- 請教各位師兄師姐們了~~不勝感激
已經(jīng)有5人回復(fù)
- 求助，各位出國讀博的師兄師姐，請給我點(diǎn)建議，小弟在此拜謝！
已經(jīng)有8人回復(fù)
- 【求助】各位師哥師姐幫小弟檢查一下復(fù)雜參數(shù)回歸的程序吧��！小弟萬分感謝！��！
已經(jīng)有3人回復(fù)
- 【求助】師兄師姐們啊 wingslib快搞死我了 幫幫小弟吧
已經(jīng)有4人回復(fù)
點(diǎn)擊這里搜索更多相關(guān)資源科研從小木蟲開始，人人為我，我為人人
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 wodaifei 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料專碩找調(diào)劑 +3 | 哈哈哈吼吼吼哈 2026-03-23 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +5 | 葵梓衛(wèi)隊(duì) 2026-03-18 | 7/350 |
|
|
[考研] 328求調(diào)劑 +4 | LHHL66 2026-03-23 | 4/200 |
|
|
|
[考研] 北科281學(xué)碩材料求調(diào)劑 +8 | tcxiaoxx 2026-03-20 | 9/450 |
|
|
|
[考研] 323求調(diào)劑 +6 | 洼小桶 2026-03-18 | 6/300 |
|
|
|
[考研] 280分求調(diào)劑 一志愿085802 +4 | PUMPT 2026-03-22 | 7/350 |
|
|
|
[考研] 求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項(xiàng)目，四級已過 +6 | 幸運(yùn)哩哩 2026-03-22 | 10/500 |
|
|
|
[考研] 308求調(diào)劑 +3 | 墨墨漠 2026-03-21 | 3/150 |
|
|
|
[考研] 求調(diào)劑院校信息 +6 | CX 330 2026-03-21 | 6/300 |
|
|
|
[考研] 286分人工智能專業(yè)請求調(diào)劑愿意跨考！ +4 | lemonzzn 2026-03-17 | 8/400 |
|
|
|
[考研] 材料工程專碩 348分求調(diào)劑 +3 | 冬辭. 2026-03-17 | 5/250 |
|
|
|
[考研] 299求調(diào)劑 +5 | shxchem 2026-03-20 | 7/350 |
|
|
|
[基金申請]
學(xué)校已經(jīng)提交到NSFC，還能修改嗎？
 40+4
|
babangida 2026-03-19 | 9/450 |
|
|
|
[考研] 332求調(diào)劑 +3 | 鳳凰院丁真 2026-03-20 | 3/150 |
|
|
|
[考研] 303求調(diào)劑 +5 | 睿08 2026-03-17 | 7/350 |
|
|
|
[考研] 299求調(diào)劑 +6 | △小透明* 2026-03-17 | 6/300 |
|
|
|
[考研] 一志愿華南師大 070300（化學(xué)）304分求調(diào)劑 +3 | 0703武芊慧雪304 2026-03-18 | 3/150 |
|
|
|
[考研] A區(qū)線材料學(xué)調(diào)劑 +5 | 周周無極 2026-03-20 | 5/250 |
|
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
|
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|