版塊導航: 正在加載中...

登錄注冊

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

24小時熱門版塊排行榜

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

返回列表

當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖

yupeifeng

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 960.3
散金: 34
紅花: 1
帖子: 220
在線: 78.4小時
蟲號: 1421266
注冊: 2011-09-28
專業(yè): 無機材料化學

[交流] matlab求解方程中的參數(shù) 已有4人參與

搞了這么久還是沒有解決這個問題�，F(xiàn)在求助于高手給予指點。謝謝了！！
實驗數(shù)據(jù)為：（t,c）=（0,0.69）（2,0.645）（4,0.635）（8,0.62）（24,0.61）（48,0.61）.其中t為時間，c為某離子的濃度。
動力學方程模型為：-dc/dt=k*(c0-c)^(1/3)*(c-c~).
其中c0為初始濃度可以取0.7，c~為平衡濃度取0.61.
怎么樣才能擬合出參數(shù)k的值呢？謝謝大家給出程序代碼，再次感謝

回復(fù)此樓

» 猜你喜歡

計算與實驗科學工程國際會議(8月7-11日,香港) 二維材料分會S5與聲學換能器分會S8 已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有4人回復(fù)
物理學I論文潤色/翻譯怎么收費? 已經(jīng)有294人回復(fù)
歡迎調(diào)劑到贛南師范大學智能制造與未來能源學院已經(jīng)有28人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有4人回復(fù)
歡迎調(diào)劑到贛南師范大學智能制造與未來能源學院已經(jīng)有26人回復(fù)
重慶交大26年碩士生招生擬調(diào)劑通知已出！歡迎加入光子學微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
UJN物理學專業(yè)調(diào)劑已經(jīng)有3人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有0人回復(fù)

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

matlab有限差分法求解復(fù)雜的偏微分方程組已經(jīng)有19人回復(fù)
matlab微分方程組參數(shù)擬合，以周為單位求解，匯總后以年為單位進行數(shù)值比較已經(jīng)有4人回復(fù)
Matlab求解偏微分方程組已經(jīng)有22人回復(fù)
matlab求解非線性方程組已經(jīng)有16人回復(fù)
求回歸方程參數(shù)估計的最優(yōu)化算法matlab代碼已經(jīng)有9人回復(fù)
求matlab求解一個方程組的問題已經(jīng)有6人回復(fù)
用Matlab求解方程已經(jīng)有6人回復(fù)
求高人指點用matlab求解非線性方程組，解決了追加100金幣；已經(jīng)有11人回復(fù)
matlab的fsove 命令求解非線性方程組已經(jīng)有6人回復(fù)
【求助】matlab 求解微分方程中的未知參數(shù) 已經(jīng)有20人回復(fù)
【求助】MATLAB 有限差分法（FDM）求解偏微分方程已經(jīng)有22人回復(fù)
【求助】積分微分方程matlab求解已經(jīng)有6人回復(fù)
【討論】CFD軟件是求解PDE的，matlab中也可以求解PDE 已經(jīng)有6人回復(fù)
【求助】matlab怎么求解偏微分方程組啊，先謝謝了已經(jīng)有13人回復(fù)
【求助】求MAtlab求解程序！��！已經(jīng)有20人回復(fù)
【求助】求教matlab解非線性方程組已經(jīng)有9人回復(fù)

憧憬未來

1樓 2012-07-02 16:58:15

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

dbb627

榮譽版主 (著名寫手)

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

★ ★ ★ ★
fegg7502: 金幣+4, 應(yīng)助指數(shù)+1, 專家考核, 多謝交流 2012-07-09 08:07:38

結(jié)果如下
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.

Estimated Parameters by Lsqnonlin():
k = 1.0138 ± 0.3000
The sum of the residual squares is: 1.5e-004

贊一下

回復(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.

3樓2012-07-02 17:39:27

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

查看全部 22 個回答

dbb627

榮譽版主 (著名寫手)

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

★ ★ ★ ★ ★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
fegg7502: 金幣+4, 多謝交流 2012-07-09 08:07:15

這個其實論壇里有很多例子，參考就能寫出來。
給你寫了個

CODE:

function parafit

clear all;

t=[0 2 4 8 24 48];

y=[0.69 0.645 0.635 0.62 0.61 0.61];

y0=0.69;

% Nonlinear least square estimate using lsqnonlin()

beta0=0.5;

lb=[0];ub=[inf];

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

    lsqnonlin(@Func,beta0,lb,ub,[],t,y,y0);         

ci = nlparci(beta,residual,jacobian);

beta;

% result

fprintf('\n Estimated Parameters by Lsqnonlin():\n')

fprintf('\t k = %.4f ± %.4f\n',beta(1),ci(1,2)-beta(1))

fprintf('  The sum of the residual squares is: %.1e\n\n',sum(residual.^2))

% plot of fit results

tspan = [0  max(t)];

[tt yc] = ode45(@ModelEqs,tspan,y0,[],beta);

tc=linspace(0,max(t),200);

yca = spline(tt,yc,tc);

plot(t,y,'ro',tc,yca,'r-');

hold on

xlabel('Time');

ylabel('Concentration');

hold off

% =======================================

function f1 = Func(beta,t,y,y0)        % Define objective function

tspan =t;

[tt yy] = ode45(@ModelEqs,tspan,y0,[],beta);

yc= spline(tt,yy,t);

f1=y-yc;

% ==================================

function dydt = ModelEqs(t,y,beta)          % Model equations

dydt = -beta*(0.7-y).^(1/3)*(y-0.61);

贊一下(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.

2樓2012-07-02 17:38:01

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

yupeifeng

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 960.3
散金: 34
紅花: 1
帖子: 220
在線: 78.4小時
蟲號: 1421266
注冊: 2011-09-28
專業(yè): 無機材料化學

引用回帖:

3樓: Originally posted by dbb627 at 2012-07-02 17:39:27
結(jié)果如下
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.

< ...

首先感謝這位學長啊，我有點看不懂啊，你幫我分析下思路，步驟是什么？就是先做設(shè)么計算，

贊一下

回復(fù)此樓

憧憬未來

4樓2012-07-02 22:01:21

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

dingd

鐵桿木蟲 (職業(yè)作家)

計算強帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

★ ★ ★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
fegg7502: 金幣+2, 多謝交流 2012-07-09 08:07:55

1stOpt求解：

CODE:

Constant c0=0.7,cp=0.61;

Variable t,c;

ODEFunction c'=-k*(c0-c)^(1/3)*(c-cp)

Data;

t,c

0,0.69

2,0.645

4,0.635

8,0.62

24,0.61

48,0.61

均方差(RMSE): 0.00535058522113925
殘差平方和(SSE): 0.000143143811043369
相關(guān)系數(shù)(R): 0.955733024504782
相關(guān)系數(shù)之平方(R^2): 0.913425614129058
決定系數(shù)(DC): 0.852429060780033

參數(shù) 最佳估算
-------------------- -------------
k 1.02159845408008

贊一下

回復(fù)此樓

5樓2012-07-03 08:52:49

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

查看全部 22 個回答

普通表情龍兔虎貓高級回復(fù) (可上傳附件)

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

[考研] 279求調(diào)劑 +7	莫xiao 2026-04-01	7/350	2026-04-01 22:05 by 客爾美德
[考研] 0703化學338求調(diào)劑！ +8	Zuhui0306 2026-03-26	9/450	2026-04-01 21:59 by xhai2011
[考研] 0856，材料與化工321分求調(diào)劑 +13	大饞小子 2026-03-27	14/700	2026-04-01 15:58 by yanflower7133
[考研] 求調(diào)劑 +4	圖鑒212 2026-03-30	5/250	2026-04-01 15:32 by 圖鑒212
[考研] 302求調(diào)劑一志愿北航070300，本科鄭大化學 +8	圣日耳曼條 2026-04-01	10/500	2026-04-01 14:03 by ZXlzxl0425
[考研] 085600，321分求調(diào)劑 +13	大饞小子 2026-03-31	13/650	2026-04-01 12:35 by chemdavid
[考研] 一志愿北交材料工程總分358 +5	cs0106 2026-04-01	7/350	2026-04-01 11:45 by wangjy2002
[考研] 358求調(diào)劑 +3	王向陽花 2026-03-31	3/150	2026-04-01 09:56 by zzchen2000
[考研] 254材料與化工求調(diào)劑 +3	翰冬林楠 2026-03-30	4/200	2026-03-31 17:53 by yishunmin
[考研] 求調(diào)劑 +9	akdhjs 2026-03-31	10/500	2026-03-31 17:14 by zhyzzh
[考研] 322求調(diào)劑：一志愿湖南大學材料與化工（085600），已過六級。 +10	XX小鄧 2026-03-29	10/500	2026-03-31 16:46 by 不吃魚的貓
[考研] 求化學調(diào)劑 +12	wulanna 2026-03-28	12/600	2026-03-31 16:38 by 690616278
[考研] 吉大生物學326分求調(diào)劑 +3	sunnyupup 2026-03-31	3/150	2026-03-31 09:28 by longlotian
[考研] 一志愿大連理工大學材料求調(diào)劑 +6	Gymno 2026-03-30	6/300	2026-03-31 07:26 by 無際的草原
[考研] 303求調(diào)劑 +7	DLkz1314. 2026-03-30	7/350	2026-03-30 21:07 by peike
[考研] 一志愿廈門大學材料工程專碩354找調(diào)劑�。。� +5	貝唄鋇鋇 2026-03-30	5/250	2026-03-30 18:16 by 無際的草原
[考研] 332求調(diào)劑 +6	@MZB382400 2026-03-28	6/300	2026-03-30 16:57 by 無際的草原
[考研] 數(shù)一英一271專碩（085401）求調(diào)劑，可跨 +7	前行必有光 2026-03-28	8/400	2026-03-28 23:22 by 小木蟲tim
[考研] 本科新能源科學與工程，一志愿華理能動285求調(diào)劑 +7	AZMK 2026-03-28	11/550	2026-03-28 21:01 by xxxsssccc
[考研] 藥學105500求調(diào)劑 +3	Ssun。。 2026-03-28	3/150	2026-03-28 11:24 by lxf170613

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

yupeifeng

[交流] matlab求解方程中的參數(shù) 已有4人參與

» 猜你喜歡

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

dbb627

dbb627

yupeifeng

dingd

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

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