版塊導(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)擊這里查看本話題的所有回帖

yupeifeng

金蟲 (小有名氣)

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

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

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

回復(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ù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有162人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院已經(jīng)有28人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有4人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院已經(jīng)有26人回復(fù)
重慶交大26年碩士生招生擬調(diào)劑通知已出！歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
UJN物理學(xué)專業(yè)調(diào)劑已經(jīng)有3人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有0人回復(fù)

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

matlab有限差分法求解復(fù)雜的偏微分方程組已經(jīng)有19人回復(fù)
matlab微分方程組參數(shù)擬合，以周為單位求解，匯總后以年為單位進(jìn)行數(shù)值比較已經(jīng)有4人回復(fù)
Matlab求解偏微分方程組已經(jīng)有22人回復(fù)
matlab求解非線性方程組已經(jīng)有16人回復(fù)
求回歸方程參數(shù)估計(jì)的最優(yōu)化算法matlab代碼已經(jīng)有9人回復(fù)
求matlab求解一個(gè)方程組的問題已經(jīng)有6人回復(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 有限差分法（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的回帖

yupeifeng

金蟲 (小有名氣)

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

引用回帖:

5樓: Originally posted by dingd at 2012-07-03 08:52:49
1stOpt求解：

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
...

我把你的程序復(fù)制后，提示沒有定義參數(shù)，這是怎么回事啊

贊一下

回復(fù)此樓

憧憬未來

9樓2012-07-03 10:16:30

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

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

dbb627

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

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

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

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

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的回帖

dbb627

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

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

★ ★ ★ ★
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的回帖

yupeifeng

金蟲 (小有名氣)

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

引用回帖:

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.

< ...

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

贊一下

回復(fù)此樓

憧憬未來

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

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

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

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

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

[考研] 材料工程322分 +7	哈哈哈吼吼吼哈 2026-04-01	7/350	2026-04-02 07:05 by yjolah
[考研] 08工科275分求調(diào)劑 +13	AaAa7420 2026-03-31	13/650	2026-04-02 06:19 by 286640313
[考研] 291求調(diào)劑 +20	Y-cap 2026-03-29	25/1250	2026-04-01 23:49 by 欣喜777
[考研] 085801 總分275 本科新能源求調(diào)劑 +10	bradoner 2026-04-01	11/550	2026-04-01 22:57 by 來看流星雨10
[考研] 材料專碩322分 +10	哈哈哈吼吼吼哈 2026-04-01	10/500	2026-04-01 22:19 by Dyhoer
[考研] 339求調(diào)劑，想調(diào)回江蘇 +7	烤麥芽 2026-03-27	10/500	2026-04-01 21:35 by 495374996
[考研] 086502化學(xué)工程342求調(diào)劑 +7	阿姨復(fù)古不過 2026-03-27	7/350	2026-04-01 16:14 by yanflower7133
[考研] 材料專碩306英一數(shù)二 +7	z1z2z3879 2026-03-31	7/350	2026-04-01 14:50 by ZXlzxl0425
[考研] 材料專業(yè)求調(diào)劑 +7	月月鳥木 2026-04-01	7/350	2026-04-01 11:36 by wangjy2002
[考研] 求調(diào)劑，一志愿北林食品與營(yíng)養(yǎng)095500，301分，已過六級(jí)，有科研經(jīng)歷 +4	快樂儲(chǔ)蓄罐 2026-03-31	4/200	2026-04-01 09:26 by JourneyLucky
[考研] 080500-315分復(fù)試調(diào)劑 +9	上岸3821 2026-03-31	9/450	2026-03-31 17:29 by 唐沐兒
[考研] 0856 335分 +9	cccchenso 2026-03-29	9/450	2026-03-31 16:37 by lishahe
[考研] 材料工程專碩求調(diào)劑 +10	hyl3153942 2026-03-29	10/500	2026-03-31 16:31 by hypershenger
[考研] 調(diào)劑求院校招收 +7	鶴鯨鴿 2026-03-28	7/350	2026-03-31 11:21 by oooqiao
[考研] 085600，專業(yè)課化工原理，321分求調(diào)劑 +5	大饞小子 2026-03-28	5/250	2026-03-29 08:56 by qingfeng258
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +6	崔wj 2026-03-26	6/300	2026-03-29 01:11 by hanserlol
[考研] 數(shù)一英一271專碩（085401）求調(diào)劑，可跨 +7	前行必有光 2026-03-28	8/400	2026-03-28 23:22 by 小木蟲tim
[考研] 316求調(diào)劑 +7	江辭666 2026-03-26	7/350	2026-03-28 21:28 by sanrepian
[考研] 308求調(diào)劑 +7	墨墨漠 2026-03-27	7/350	2026-03-28 07:43 by 熱情沙漠
[考研] 復(fù)試調(diào)劑，一志愿南農(nóng)083200食品科學(xué)與工程 +5	XQTJZ 2026-03-26	5/250	2026-03-27 14:49 by 狂炫麥當(dāng)當(dāng)

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

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

yupeifeng

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

» 猜你喜歡

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

yupeifeng

dbb627

dbb627

yupeifeng

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

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