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

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

24小時(shí)熱門版塊排行榜

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

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

雨楓林影

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 870
帖子: 133
在線: 34小時(shí)
蟲號(hào): 1444860
注冊(cè): 2011-10-16
性別: MM
專業(yè): 食品科學(xué)基礎(chǔ)

[求助] 求大神幫忙擬合一個(gè)非線性方程，求出模型參數(shù)

我現(xiàn)有一個(gè)動(dòng)力學(xué)模型，數(shù)據(jù)也有，請(qǐng)求大神幫忙進(jìn)行一個(gè)非線性方程的擬合，求出模型中的若干參數(shù)，最好用matlab軟件，實(shí)在不行，別的軟件也可以，只要幫我把參數(shù)算出來(lái)就行了。
應(yīng)用MATLAB軟件根據(jù)發(fā)酵試驗(yàn)數(shù)據(jù)和經(jīng)驗(yàn)獲得模型參數(shù)的初估值，進(jìn)行非線性擬合，采用算法為全局性收斂的Marquardt修正的高斯牛頓法，以誤差平方和最小為目標(biāo)，獲得發(fā)酵動(dòng)力學(xué)模型參數(shù)估計(jì)值，這段話是參考一篇文獻(xiàn)，我和他用同一個(gè)方程，方法，就是試驗(yàn)數(shù)據(jù)不一樣，想重新算出一組參數(shù)。
關(guān)系到本人的畢業(yè)論文，急求，可再追加金幣，如能幫忙，萬(wàn)分感謝。
t    X    S    P
88 6.7 3.17 4.86
92 7.5 2.87 5.17
96 8.4 2.33 5.72
100 9.1 1.87 6.18
104 9.3 1.36 6.73
108 9.1 0.95 7.13
112 8.6 0.52 7.56

360截圖20140321193034293.jpg

360截圖20140321193141731.jpg

回復(fù)此樓

» 猜你喜歡

重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
北工大宋曉艷教授團(tuán)隊(duì)招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！已經(jīng)有1人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有150人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人已經(jīng)有0人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人已經(jīng)有0人回復(fù)
計(jì)算與實(shí)驗(yàn)科學(xué)工程國(guó)際會(huì)議(8月7-11日,香港) 二維材料分會(huì)S5與聲學(xué)換能器分會(huì)S8 已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有4人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來(lái)能源學(xué)院已經(jīng)有26人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有4人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來(lái)能源學(xué)院已經(jīng)有8人回復(fù)

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

【求專家】根據(jù)一組數(shù)據(jù)，求其方程，尋最優(yōu)擬合方程。要求截距為0 已經(jīng)有3人回復(fù)
關(guān)于病態(tài)方程組求解的疑問(wèn) 已經(jīng)有5人回復(fù)
用1stOpt進(jìn)行復(fù)數(shù)曲線擬合時(shí)，結(jié)果出錯(cuò)。急求教！��！已經(jīng)有5人回復(fù)
用1stOpt進(jìn)行多模型多變量非線性擬合得不到結(jié)果已經(jīng)有4人回復(fù)
求用高版本的1stOpt擬合方程已經(jīng)有8人回復(fù)
用matlab 求解參數(shù)方程組已經(jīng)有8人回復(fù)
1stopt擬合問(wèn)題求助已經(jīng)有7人回復(fù)
1stopt能進(jìn)行隱函數(shù)非線性擬合嗎？誰(shuí)能幫擬合一下新型液體的粘度混合計(jì)算公式已經(jīng)有5人回復(fù)
多元變量自定義非線性方程擬合已經(jīng)有13人回復(fù)
誰(shuí)用1stopt這個(gè)軟件啊，離散數(shù)據(jù)擬合方程，方程未知已經(jīng)有5人回復(fù)
求大神幫忙用1stOpt 擬合復(fù)數(shù)數(shù)據(jù) 已經(jīng)有13人回復(fù)
MATLAB微分方程參數(shù)擬合問(wèn)題，求大神已經(jīng)有7人回復(fù)
用1stopt擬合非線性方程結(jié)果與其他軟件擬合結(jié)果差異大已經(jīng)有5人回復(fù)

1樓 2014-03-21 19:34:10

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

月只藍(lán)

主管區(qū)長(zhǎng) (職業(yè)作家)

專家經(jīng)驗(yàn): +1059
計(jì)算強(qiáng)帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68128.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時(shí)
蟲號(hào): 1122189
注冊(cè): 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(qū)

此外，β哪里有用到？

贊一下

回復(fù)此樓

MATLAB、MS小問(wèn)題、普通問(wèn)題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無(wú)償私信求助。

4樓2014-03-22 15:04:20

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

查看全部 16 個(gè)回答

雨楓林影

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 870
帖子: 133
在線: 34小時(shí)
蟲號(hào): 1444860
注冊(cè): 2011-10-16
性別: MM
專業(yè): 食品科學(xué)基礎(chǔ)

拜托，大神快出來(lái)

回復(fù)此樓

2樓2014-03-22 11:04:14

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

月只藍(lán)

主管區(qū)長(zhǎng) (職業(yè)作家)

專家經(jīng)驗(yàn): +1059
計(jì)算強(qiáng)帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68128.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時(shí)
蟲號(hào): 1122189
注冊(cè): 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(qū)

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

★
感謝參與，應(yīng)助指數(shù) +1
fegg7502: 金幣+1, 鼓勵(lì)交流 2014-03-24 08:31:54

時(shí)間t=0時(shí)，X P S的初值呢？

贊一下

回復(fù)此樓

MATLAB、MS小問(wèn)題、普通問(wèn)題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無(wú)償私信求助。

3樓2014-03-22 14:21:27

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

月只藍(lán)

主管區(qū)長(zhǎng) (職業(yè)作家)

專家經(jīng)驗(yàn): +1059
計(jì)算強(qiáng)帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68128.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時(shí)
蟲號(hào): 1122189
注冊(cè): 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(qū)

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
雨楓林影: 金幣+100, ★★★很有幫助 2014-03-23 16:42:18
fegg7502: 金幣+3, 應(yīng)助指數(shù)+1, 3ks 2014-03-24 08:32:10

1. 8個(gè)參數(shù)的初值，你自己再調(diào)一下；
2. 原問(wèn)題沒(méi)有給出t=0時(shí)的初值，代碼編寫時(shí)以t=88為零時(shí)刻。

代碼和計(jì)算結(jié)果如下：

CODE:

function ode_fit

format long

clear all

clc

x0 = [6.7 4.86 3.17];

k0 = [1.5 1.5 1.5 6.8 2/3 0.1 0.1 0.01];   

lb = -[1 1 1 1 1 1 1 1 ]*1e9;

ub = [1 1 1 1 1 1 1 1 ]*1e9;

data=...

    [

88        6.7        3.17        4.86

92        7.5        2.87        5.17                                   

96        8.4        2.33        5.72                 

100        9.1        1.87        6.18

104        9.3        1.36        6.73               

108        9.1        0.95        7.13

112        8.6        0.52        7.56

];

tspan = [data(:,1)'-88];

yexp = [data(2:end,2) data(2:end,4) data(2:end,3)];

[k,resnorm,residual,exitflag,output,lambda,jacobian] = ...

    lsqnonlin(@ObjFunc,k0,lb,ub,[],tspan,x0,yexp);       

ci = nlparci(k,residual,jacobian);

fprintf('\n\n使用函數(shù)lsqnonlin()估計(jì)得到的參數(shù)值為:\n')

fprintf('\tk1 = %.9f \n',k(1))

fprintf('\tk2 = %.9f \n',k(2))

fprintf('\tk3 = %.9f \n',k(3))

fprintf('\tk4 = %.9f \n',k(4))

fprintf('\tk5 = %.9f \n',k(5))

fprintf('\tk6 = %.9f \n',k(6))

fprintf('\tk7 = %.9f \n',k(7))

fprintf('\tk8 = %.9f \n',k(8))

figure(1)

ts=0:((max(tspan)-min(tspan))/100):max(tspan);

[ts ys] = ode45(@KineticsEqs,ts,x0,[],k);

yy = [data(:,2) data(:,4) data(:,3)];

figure(1)

plot(ts,ys(:,1),'b',tspan,yy(:,1),'bo',ts,ys(:,2),'r',tspan,yy(:,2),'ro',ts,ys(:,3),'k',tspan,yy(:,3),'ko'),

legend('X的計(jì)算值','X的實(shí)驗(yàn)值','P的計(jì)算值','P的實(shí)驗(yàn)值','S的計(jì)算值','S的實(shí)驗(yàn)值','Location','best');

function f = ObjFunc(k,tspan,x0,yexp)           % 目標(biāo)函數(shù)

[t Xsim] = ode45(@KineticsEqs,tspan,x0,[],k);

Xsim1=Xsim(:,1);

Xsim2=Xsim(:,2);

Xsim3=Xsim(:,3);

ysim(:,1) = Xsim1(2:end);

ysim(:,2) = Xsim2(2:end);

ysim(:,3) = Xsim3(2:end);

size(ysim(:,1));

size(ysim(:,2));

size(yexp(:,1));

size(yexp(:,2));

f = [(ysim(:,1)-yexp(:,1)) (ysim(:,2)-yexp(:,2)) (ysim(:,3)-yexp(:,3))];

function dCdt = KineticsEqs(t,C,k)              % ODE模型方程

umax=k(1);Ks=k(2);Ki=k(3);Pmax=k(4);a=k(5);Yxs=k(6);Yps=k(7);m=k(8);

X=C(1);P=C(2);S=C(3);

dXdt=umax*(S./(S+Ks)).*(1./(1+S/Ki)).*(1-P/Pmax).*X;

dPdt=a*dXdt;

dSdt=-Yxs*dXdt-Yps*dPdt-m*X;

dCdt = [dXdt; dPdt;dSdt];

附圖1.jpg

贊一下(1人)

回復(fù)此樓

MATLAB、MS小問(wèn)題、普通問(wèn)題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無(wú)償私信求助。

5樓2014-03-22 15:36:01

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

查看全部 16 個(gè)回答

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

[考研] 285求調(diào)劑 +6	AZMK 2026-03-30	10/500	2026-03-31 11:41 by AZMK
[考研] 南京大學(xué)化學(xué)調(diào)劑 +11	景隨風(fēng) 2026-03-29	16/800	2026-03-31 10:14 by herarysara
[考研] 求調(diào)劑生物學(xué) 377分 +4	zzll03 2026-03-31	4/200	2026-03-31 10:07 by zhyzzh
[考研] 物理學(xué)調(diào)劑 +3	小羊36 2026-03-30	3/150	2026-03-31 09:20 by jp9609
[考研] 308求調(diào)劑 +8	墨墨漠 2026-03-25	8/400	2026-03-31 09:15 by dgd2013
[考研] 083000環(huán)境科學(xué)與工程調(diào)劑，總分281 +4	橙子（勝意） 2026-03-30	4/200	2026-03-31 00:44 by Linzejun
[考研] 085602化工求調(diào)劑（331分） +8	111@127 2026-03-30	8/400	2026-03-30 21:23 by 研究僧導(dǎo)導(dǎo)
[考研] 材料工程專碩求調(diào)劑 +8	hyl3153942 2026-03-29	8/400	2026-03-30 20:45 by dophin1985
[考研] 289求調(diào)劑 +16	新時(shí)代材料 2026-03-27	16/800	2026-03-30 19:04 by Wang200018
[考研] 105500藥學(xué)求調(diào)劑，一志愿山東大學(xué)藥學(xué)，348分 +3	gr哈哈哈 2026-03-28	3/150	2026-03-30 18:56 by 源_2020
[考研] 286求調(diào)劑 +5	Faune 2026-03-30	5/250	2026-03-30 17:15 by 無(wú)際的草原
[考研] 324求調(diào)劑 +9	hanamiko 2026-03-26	11/550	2026-03-30 14:27 by JourneyLucky
[考研] 一志愿北京工業(yè)大學(xué)，324分求調(diào)劑 +6	零八# 2026-03-28	6/300	2026-03-29 21:20 by nanaliuyun
[考研] 求調(diào)劑 +4	QiMing7 2026-03-25	5/250	2026-03-29 21:10 by 唐沐兒
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3	@taotao 2026-03-28	3/150	2026-03-28 10:26 by JourneyLucky
[考研] 材料求調(diào)劑一志愿哈工大324 +7	閆旭東 2026-03-28	9/450	2026-03-28 08:51 by Xu de nuo
[考研] 化學(xué)調(diào)劑 +4	愛(ài)吃番茄的旭 2026-03-24	5/250	2026-03-27 17:50 by kiokin
[考研] 085600，材料與化工321分調(diào)劑 +4	大饞小子 2026-03-27	6/300	2026-03-27 14:11 by 松花缸1201
[考研] 求調(diào)劑 +6	林之夕 2026-03-24	6/300	2026-03-27 08:38 by hypershenger
[考研] 302求調(diào)劑 +4	錦衣衛(wèi)藤椒 2026-03-25	4/200	2026-03-25 16:29 by 功夫瘋狂