heshuangrong

應(yīng)助: 2 (幼兒園)
金幣: 784.2
散金: 25
帖子: 96
在線: 93.5小時(shí)
蟲(chóng)號(hào): 1498593
注冊(cè): 2011-11-18
專業(yè): 化學(xué)反應(yīng)工程

[求助] MATLAB擬合動(dòng)力學(xué)問(wèn)題已有1人參與

程序是書(shū)上的，店Debug--Run沒(méi)法運(yùn)行，求教各位大俠

function KineticDataFit
clear all
clc

% 動(dòng)力學(xué)數(shù)據(jù)
t = [0 20  40  60  120  180  300];
CA = [10  8  6  5  3  2  1];

% 用最小二乘樣條擬合法計(jì)算微分dCA/dt--使用不經(jīng)過(guò)實(shí)驗(yàn)點(diǎn)的B樣條插值函數(shù)
knots = 3;
K = 3;    % 三次B樣條
sp = spap2(knots,K,t,CA)
sp = spap2(newknt(sp),K,t,CA);
pp = fnder(sp)       % 計(jì)算B樣條函數(shù)的導(dǎo)函數(shù)
dCAdt = fnval(pp,t)    % 計(jì)算t處的導(dǎo)函數(shù)值

% 繪制圖形
ti = linspace(t(1),t(end),200);
CAi = fnval(sp,ti)
plot(t,CA,'ro',ti,CAi,'b-')
xlabel('t')
ylabel('C_A')
figure
fnplt(pp)
% dCAdti = fnval(pp,ti)
% plot(ti,dCAdti,'-')
xlabel('t')
ylabel('dC/dt')

% 線性擬合
rA = dCAdt;
y = log(-rA);
x = log(CA);
p = polyfit(x,y,1);
k = exp(p(2))
n = p(1)

回復(fù)此樓

» 猜你喜歡

中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博士后入站、德國(guó)兩年全職已經(jīng)有0人回復(fù)
面向全球申請(qǐng)人的中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博后入站德國(guó)兩年全職已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有90人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(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é)理學(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)有1人回復(fù)

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

向大神求教一個(gè)用matlab擬合動(dòng)力學(xué)參數(shù)的問(wèn)題，求指導(dǎo) 已經(jīng)有5人回復(fù)
給我大神求教一個(gè)用matlab擬合動(dòng)力學(xué)參數(shù)的問(wèn)題，求指導(dǎo) 已經(jīng)有0人回復(fù)
Matlab擬合動(dòng)力學(xué)參數(shù)問(wèn)題，求指導(dǎo)！已經(jīng)有9人回復(fù)
matlab 擬合多參數(shù)動(dòng)力學(xué)方程求程序已經(jīng)有4人回復(fù)
MATLAB擬合動(dòng)力學(xué)問(wèn)題已經(jīng)有11人回復(fù)
matlab擬合動(dòng)力學(xué)參數(shù)，非常謝謝~ 已經(jīng)有0人回復(fù)
求助matlab擬合動(dòng)力學(xué)K值已經(jīng)有4人回復(fù)
matlab擬合吸附動(dòng)力學(xué) 已經(jīng)有0人回復(fù)
使用MATLAB擬合動(dòng)力學(xué)參數(shù)問(wèn)題已經(jīng)有18人回復(fù)
拜求！matlab擬合動(dòng)力學(xué)參數(shù) 已經(jīng)有14人回復(fù)
求助用matlab擬合動(dòng)力學(xué)方程已經(jīng)有9人回復(fù)
求助用matlab擬合動(dòng)力學(xué)方程已經(jīng)有13人回復(fù)
matlab 擬合反應(yīng)動(dòng)力學(xué)參數(shù) 已經(jīng)有8人回復(fù)
matlab擬合酶促反應(yīng)機(jī)制已經(jīng)有7人回復(fù)
求大家?guī)蛶兔纯磎atlab擬合動(dòng)力學(xué)方程的問(wèn)題，已經(jīng)走投無(wú)路了，求大家?guī)蛶兔Π。?/a> 已經(jīng)有19人回復(fù)

matlab 擬合反應(yīng)動(dòng)力學(xué)參數(shù)結(jié)果很差。大家?guī)兔匆幌?/a> 已經(jīng)有14人回復(fù)

用matlab 擬合動(dòng)力學(xué)級(jí)數(shù) 已經(jīng)有2人回復(fù)
matlab擬合擬均相動(dòng)力學(xué)參數(shù) 已經(jīng)有6人回復(fù)
matlab擬合動(dòng)力學(xué)參數(shù) 已經(jīng)有4人回復(fù)
【求助】使用Matlab擬合反應(yīng)動(dòng)力學(xué)方程問(wèn)題已經(jīng)有7人回復(fù)

1樓 2015-09-01 14:23:46

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

lvgao116

鐵桿木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 4 (幼兒園)
金幣: 20624.4
紅花: 13
沙發(fā): 1
帖子: 1691
在線: 153小時(shí)
蟲(chóng)號(hào): 907551
注冊(cè): 2009-11-19
性別: GG
專業(yè): 應(yīng)用地球物理學(xué)

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

感謝參與，應(yīng)助指數(shù) +1

不要輸入帶有function語(yǔ)句的那一行�；蛘哂�%注釋。我已經(jīng)試過(guò)沒(méi)問(wèn)題。

贊一下(1人)

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

heshuangrong

lovelinan

2樓2015-09-01 15:21:47

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

heshuangrong

新蟲(chóng) (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 784.2
散金: 25
帖子: 96
在線: 93.5小時(shí)
蟲(chóng)號(hào): 1498593
注冊(cè): 2011-11-18
專業(yè): 化學(xué)反應(yīng)工程

送紅花一朵

引用回帖:

2樓: Originally posted by lvgao116 at 2015-09-01 15:21:47
不要輸入帶有function語(yǔ)句的那一行�；蛘哂�%注釋。我已經(jīng)試過(guò)沒(méi)問(wèn)題。

為什么我點(diǎn)Debug-run后，在命令窗口紅色顯示undefined function or variable 'spap2'，第一行我都把它刪了

發(fā)自小木蟲(chóng)Android客戶端

贊一下

回復(fù)此樓

3樓2015-09-01 15:49:15

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

lvgao116

鐵桿木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 4 (幼兒園)
金幣: 20624.4
紅花: 13
沙發(fā): 1
帖子: 1691
在線: 153小時(shí)
蟲(chóng)號(hào): 907551
注冊(cè): 2009-11-19
性別: GG
專業(yè): 應(yīng)用地球物理學(xué)

引用回帖:

3樓: Originally posted by heshuangrong at 2015-09-01 15:49:15
為什么我點(diǎn)Debug-run后，在命令窗口紅色顯示undefined function or variable 'spap2'，第一行我都把它刪了
...

spap2是一個(gè)matlab自帶函數(shù)，你運(yùn)行后matlab給你的提示是沒(méi)有定義該函數(shù)或變量，說(shuō)明你的matlab版本中沒(méi)有這個(gè)函數(shù)。ps：我的matlab版本是R2008a，7.6.0.

贊一下

回復(fù)此樓

lovelinan

4樓2015-09-01 17:01:02

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

lvgao116

鐵桿木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 4 (幼兒園)
金幣: 20624.4
紅花: 13
沙發(fā): 1
帖子: 1691
在線: 153小時(shí)
蟲(chóng)號(hào): 907551
注冊(cè): 2009-11-19
性別: GG
專業(yè): 應(yīng)用地球物理學(xué)

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

引用回帖:

sp =

   form: 'B-'
   knots: [0 0 0 50 120 300 300 300]
   coefs: [10.0579 6.9123 3.5802 1.3462 0.9975]
number: 5
   order: 3
   dim: 1

pp =

   form: 'B-'
   knots: [0 0 78.8178 182.1889 300 300]
   coefs: [-0.1187 -0.0312 -0.0104 -0.0061]
number: 4
   order: 2
   dim: 1

dCAdt =

-0.1187 -0.0965 -0.0743 -0.0521 -0.0229 -0.0109 -0.0061

CAi =

  Columns 1 through 11

10.0307 9.8531 9.6780 9.5054 9.3353 9.1678 9.0027 8.8402 8.6803 8.5228 8.3679

  Columns 12 through 22

8.2155 8.0656 7.9182 7.7734 7.6310 7.4913 7.3540 7.2192 7.0870 6.9573 6.8301

  Columns 23 through 33

6.7054 6.5833 6.4637 6.3466 6.2320 6.1200 6.0105 5.9035 5.7990 5.6970 5.5976

  Columns 34 through 44

5.5007 5.4063 5.3144 5.2251 5.1382 5.0539 4.9722 4.8929 4.8162 4.7420 4.6703

  Columns 45 through 55

4.6011 4.5345 4.4703 4.4088 4.3497 4.2931 4.2391 4.1876 4.1386 4.0916 4.0452

  Columns 56 through 66

3.9992 3.9536 3.9085 3.8639 3.8197 3.7760 3.7328 3.6900 3.6476 3.6057 3.5643

  Columns 67 through 77

3.5233 3.4828 3.4427 3.4031 3.3640 3.3253 3.2870 3.2492 3.2119 3.1750 3.1386

  Columns 78 through 88

3.1027 3.0672 3.0321 2.9975 2.9634 2.9297 2.8965 2.8638 2.8315 2.7996 2.7682

  Columns 89 through 99

2.7373 2.7068 2.6768 2.6472 2.6181 2.5895 2.5613 2.5335 2.5062 2.4794 2.4530

  Columns 100 through 110

2.4271 2.4017 2.3767 2.3521 2.3280 2.3044 2.2812 2.2585 2.2363 2.2144 2.1931

  Columns 111 through 121

2.1722 2.1518 2.1318 2.1123 2.0932 2.0746 2.0564 2.0387 2.0215 2.0047 1.9884

  Columns 122 through 132

1.9725 1.9568 1.9413 1.9258 1.9104 1.8950 1.8798 1.8646 1.8496 1.8346 1.8197

  Columns 133 through 143

1.8048 1.7901 1.7754 1.7608 1.7463 1.7319 1.7176 1.7033 1.6892 1.6751 1.6611

  Columns 144 through 154

1.6472 1.6333 1.6196 1.6059 1.5923 1.5788 1.5654 1.5521 1.5388 1.5256 1.5126

  Columns 155 through 165

1.4996 1.4866 1.4738 1.4610 1.4484 1.4358 1.4233 1.4108 1.3985 1.3862 1.3740

  Columns 166 through 176

1.3619 1.3499 1.3380 1.3262 1.3144 1.3027 1.2911 1.2796 1.2682 1.2568 1.2456

  Columns 177 through 187

1.2344 1.2233 1.2122 1.2013 1.1905 1.1797 1.1690 1.1584 1.1479 1.1374 1.1271

  Columns 188 through 198

1.1168 1.1066 1.0965 1.0865 1.0765 1.0667 1.0569 1.0472 1.0376 1.0281 1.0186

  Columns 199 through 200

1.0093 1.0000

k =

0.0053

n =

1.3854

贊一下

回復(fù)此樓

lovelinan

5樓2015-09-01 17:01:44

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

heshuangrong

新蟲(chóng) (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 784.2
散金: 25
帖子: 96
在線: 93.5小時(shí)
蟲(chóng)號(hào): 1498593
注冊(cè): 2011-11-18
專業(yè): 化學(xué)反應(yīng)工程

引用回帖:

4樓: Originally posted by lvgao116 at 2015-09-01 17:01:02
spap2是一個(gè)matlab自帶函數(shù)，你運(yùn)行后matlab給你的提示是沒(méi)有定義該函數(shù)或變量，說(shuō)明你的matlab版本中沒(méi)有這個(gè)函數(shù)。ps：我的matlab版本是R2008a，7.6.0....

我的是6.5,去下個(gè)高點(diǎn)的版本看看

發(fā)自小木蟲(chóng)Android客戶端

贊一下

回復(fù)此樓

6樓2015-09-01 17:28:16

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

lvgao116

鐵桿木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 4 (幼兒園)
金幣: 20624.4
紅花: 13
沙發(fā): 1
帖子: 1691
在線: 153小時(shí)
蟲(chóng)號(hào): 907551
注冊(cè): 2009-11-19
性別: GG
專業(yè): 應(yīng)用地球物理學(xué)

引用回帖:

6樓: Originally posted by heshuangrong at 2015-09-01 17:28:16
我的是6.5,去下個(gè)高點(diǎn)的版本看看
...

好噠！

發(fā)自小木蟲(chóng)Android客戶端

回復(fù)此樓

lovelinan

7樓2015-09-01 19:57:00

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

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

[考研] 求調(diào)劑 +4	青春裁為三截 2026-03-29	4/200	2026-03-29 15:01 by 唐沐兒
[考研] 348求調(diào)劑 +5	小懶蟲(chóng)不懶了 2026-03-28	5/250	2026-03-29 10:34 by 唐沐兒
[考研] 學(xué)碩274求調(diào)劑 +9	Li李魚(yú) 2026-03-26	9/450	2026-03-28 21:42 by bymhappy
[考研] 求化學(xué)調(diào)劑 +4	wulanna 2026-03-28	4/200	2026-03-28 13:37 by 唐沐兒
[考研] 調(diào)劑 +3	好好讀書(shū)。 2026-03-28	3/150	2026-03-28 12:04 by 王保杰33
[考研] 328求調(diào)劑 +7	嗯滴的基本都 2026-03-27	7/350	2026-03-28 04:19 by fmesaito
[考研] 求調(diào)劑 +4	零八# 2026-03-27	4/200	2026-03-27 18:07 by yu221
[考研] 272求調(diào)劑 +7	腳滑的守法公民 2026-03-27	7/350	2026-03-27 17:23 by laoshidan
[考研] 322求調(diào)劑 +4	我真的很想學(xué)習(xí) 2026-03-23	4/200	2026-03-27 13:51 by 楊楊楊紫
[考研] 調(diào)劑 +3	李嘉圖·S·路 2026-03-27	3/150	2026-03-27 11:19 by wangjy2002
[考研] 一志愿鄭大085600，310分求調(diào)劑 +5	李瀟可 2026-03-26	5/250	2026-03-27 11:14 by 不吃魚(yú)的貓
[考研] 286求調(diào)劑 +4	lim0922 2026-03-26	4/200	2026-03-27 10:28 by guoweigw
[考研] 325求調(diào)劑 +5	李嘉圖·S·路 2026-03-23	5/250	2026-03-27 00:42 by wxiongid
[考研] 調(diào)劑求收留 +7	果然有我 2026-03-26	7/350	2026-03-27 00:26 by wxiongid
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3	崔wj 2026-03-26	3/150	2026-03-26 19:57 by nihaoar
[考研] 機(jī)械學(xué)碩總分317求調(diào)劑！�。�！ +4	Acaciad 2026-03-25	4/200	2026-03-25 19:59 by hanserlol
[考研] 求b區(qū)院校調(diào)劑 +4	周56 2026-03-24	5/250	2026-03-25 17:12 by yishunmin
[考研] 【2026考研調(diào)劑】制藥工程 284分求相關(guān)專業(yè)調(diào)劑名額 +4	袁奐奐 2026-03-25	8/400	2026-03-25 14:32 by lbsjt
[考研] 一志愿北化315 求調(diào)劑 +3	akrrain 2026-03-24	3/150	2026-03-24 19:35 by 了了了了。。
[考研] 277分求調(diào)劑，跨調(diào)材料 +3	考研調(diào)劑lxh 2026-03-24	3/150	2026-03-24 13:52 by JourneyLucky