heshuangrong

新蟲 (小有名氣)

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

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

程序是書上的，店Debug--Run沒法運(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)過實(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ù)此樓

» 猜你喜歡

哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有7人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有166人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請(qǐng)問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有31人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有68人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)

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

向大神求教一個(gè)用matlab擬合動(dòng)力學(xué)參數(shù)的問題，求指導(dǎo) 已經(jīng)有5人回復(fù)
給我大神求教一個(gè)用matlab擬合動(dòng)力學(xué)參數(shù)的問題，求指導(dǎo) 已經(jīng)有0人回復(fù)
Matlab擬合動(dòng)力學(xué)參數(shù)問題，求指導(dǎo)！已經(jīng)有9人回復(fù)
matlab 擬合多參數(shù)動(dòng)力學(xué)方程求程序已經(jīng)有4人回復(fù)
MATLAB擬合動(dòng)力學(xué)問題已經(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ù)問題已經(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é)方程的問題，已經(jīng)走投無路了，求大家?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é)方程問題已經(jīng)有7人回復(fù)

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

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

lvgao116

鐵桿木蟲 (著名寫手)

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

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

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

不要輸入帶有function語句的那一行。或者用%注釋。我已經(jīng)試過沒問題。

贊一下(1人)

回復(fù)此樓

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

heshuangrong

lovelinan

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

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

heshuangrong

新蟲 (小有名氣)

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

送紅花一朵

引用回帖:

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

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

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

贊一下

回復(fù)此樓

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

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

lvgao116

鐵桿木蟲 (著名寫手)

應(yīng)助: 4 (幼兒園)
金幣: 20624.4
紅花: 13
沙發(fā): 1
帖子: 1691
在線: 153小時(shí)
蟲號(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給你的提示是沒有定義該函數(shù)或變量，說明你的matlab版本中沒有這個(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

鐵桿木蟲 (著名寫手)

應(yīng)助: 4 (幼兒園)
金幣: 20624.4
紅花: 13
沙發(fā): 1
帖子: 1691
在線: 153小時(shí)
蟲號(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

新蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 784.2
散金: 25
帖子: 96
在線: 93.5小時(shí)
蟲號(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給你的提示是沒有定義該函數(shù)或變量，說明你的matlab版本中沒有這個(gè)函數(shù)。ps：我的matlab版本是R2008a，7.6.0....

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

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

贊一下

回復(fù)此樓

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

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

lvgao116

鐵桿木蟲 (著名寫手)

應(yīng)助: 4 (幼兒園)
金幣: 20624.4
紅花: 13
沙發(fā): 1
帖子: 1691
在線: 153小時(shí)
蟲號(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ā)自小木蟲Android客戶端

回復(fù)此樓

lovelinan

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

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

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

[考研] 0703化學(xué)調(diào)劑 +11	妮妮ninicgb 2026-03-15	15/750	2026-03-21 19:15 by ColorlessPI
[考研] 278求調(diào)劑 +9	煙火先于春 2026-03-17	9/450	2026-03-21 17:47 by 學(xué)員8dgXkO
[考研] 297求調(diào)劑 +11	戲精丹丹丹 2026-03-17	12/600	2026-03-21 17:47 by ColorlessPI
[考研] 317求調(diào)劑 +9	申子申申 2026-03-19	15/750	2026-03-21 17:31 by 學(xué)員8dgXkO
[考研] 336求調(diào)劑 +5	rmc8866 2026-03-21	5/250	2026-03-21 17:24 by 學(xué)員8dgXkO
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +6	自由煎餅果子 2026-03-16	6/300	2026-03-21 17:21 by 學(xué)員8dgXkO
[考研] 266求調(diào)劑 +3	哇呼哼呼哼 2026-03-20	3/150	2026-03-21 16:46 by barlinike
[考研] 0805材料320求調(diào)劑 +3	深海物語 2026-03-20	3/150	2026-03-21 15:46 by 無際的草原
[考研] 一志愿山大07化學(xué) 332分四六級(jí)已過本科山東雙非求調(diào)劑！ +3	不想理你 2026-03-16	3/150	2026-03-21 03:59 by JourneyLucky
[考研] 265求調(diào)劑 +9	梁梁校校 2026-03-17	9/450	2026-03-21 02:17 by JourneyLucky
[考研] 271材料工程求調(diào)劑 +8	.6lL 2026-03-18	8/400	2026-03-21 00:58 by JourneyLucky
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4	葵梓衛(wèi)隊(duì) 2026-03-18	6/300	2026-03-20 23:02 by JourneyLucky
[考研] 289求調(diào)劑 +6	懷瑾握瑜l 2026-03-20	6/300	2026-03-20 20:30 by 學(xué)員8dgXkO
[考研] 一志愿南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +5	@taotao 2026-03-20	5/250	2026-03-20 20:16 by JourneyLucky
[考研] 261求B區(qū)調(diào)劑，科研經(jīng)歷豐富 +3	牛奶很忙 2026-03-20	4/200	2026-03-20 19:34 by JourneyLucky
[考研] 08工學(xué)調(diào)劑 +5	用戶573181 2026-03-20	5/250	2026-03-20 15:47 by xia_2003
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +4	xbxudjdn 2026-03-18	4/200	2026-03-20 09:07 by 不168
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3	Eember. 2026-03-17	9/450	2026-03-18 10:28 by Eember.
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6	薛云鵬 2026-03-15	6/300	2026-03-17 11:05 by 學(xué)員h26Tkc
[考研] 275求調(diào)劑 +4	太陽花天天開心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂