heshuangrong

新蟲 (小有名氣)

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

[求助] MATLAB擬合動力學問題已有1人參與

程序是書上的，店Debug--Run沒法運行，求教各位大俠

function KineticDataFit
clear all
clc

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

% 用最小二乘樣條擬合法計算微分dCA/dt--使用不經過實驗點的B樣條插值函數(shù)
knots = 3;
K = 3;    % 三次B樣條
sp = spap2(knots,K,t,CA)
sp = spap2(newknt(sp),K,t,CA);
pp = fnder(sp)       % 計算B樣條函數(shù)的導函數(shù)
dCAdt = fnval(pp,t)    % 計算t處的導函數(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)

回復此樓

» 猜你喜歡

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

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

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

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

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

lvgao116

鐵桿木蟲 (著名寫手)

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

【答案】應助回帖

引用回帖:

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

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

贊一下

回復此樓

lovelinan

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

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

查看全部 7 個回答

lvgao116

鐵桿木蟲 (著名寫手)

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

【答案】應助回帖

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

不要輸入帶有function語句的那一行�；蛘哂�%注釋。我已經試過沒問題。

贊一下(1人)

回復此樓

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

heshuangrong

lovelinan

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

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

heshuangrong

新蟲 (小有名氣)

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

送紅花一朵

引用回帖:

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

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

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

贊一下

回復此樓

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

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

lvgao116

鐵桿木蟲 (著名寫手)

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

引用回帖:

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

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

贊一下

回復此樓

lovelinan

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

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

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

[考研] 263求調劑 +4	yqdszhdap－ 2026-03-22	4/200	2026-03-22 21:20 by 1144970272
[考研] 材料與化工085600，總分304，本科有兩篇sci參與，求調劑 +4	幸運的醬醬 2026-03-22	5/250	2026-03-22 20:15 by edmund7
[考研] 306求調劑 +6	chuanzhu川燭 2026-03-18	6/300	2026-03-22 19:39 by 公瑾逍遙
[考研] 環(huán)境學碩288求調劑 +6	皮皮皮123456 2026-03-22	6/300	2026-03-22 16:52 by i_cooler
[考研] 319求調劑 +4	小力氣珂珂 2026-03-20	4/200	2026-03-22 15:53 by ColorlessPI
[考研] 求調劑院校信息 +6	CX 330 2026-03-21	6/300	2026-03-22 15:25 by 無懈可擊111
[考研] 269專碩求調劑 +6	金恩貝 2026-03-21	6/300	2026-03-22 14:31 by ColorlessPI
[考研] 一志愿北京化工大學070300 學碩336求調劑 +5	vv迷 2026-03-21	8/400	2026-03-22 14:20 by ColorlessPI
[考研] 求調劑 +5	Zhangbod 2026-03-21	7/350	2026-03-22 13:13 by Zhangbod
[考博] 招收博士1-2人 +3	QGZDSYS 2026-03-18	4/200	2026-03-22 10:25 by QGZDSYS
[考研] 化學調劑 +5	yzysaa 2026-03-21	5/250	2026-03-21 22:12 by peike
[考研] 求調劑 +3	.m.. 2026-03-21	4/200	2026-03-21 16:25 by barlinike
[考研] 295求調劑 +4	一志愿京區(qū)211 2026-03-18	6/300	2026-03-20 23:41 by JourneyLucky
[考研] 一志愿南昌大學，327分，材料與化工085600 +9	Ncdx123456 2026-03-19	9/450	2026-03-20 23:41 by lovewei0727
[考研] 350求調劑 +5	weudhdk 2026-03-19	5/250	2026-03-20 22:04 by luoyongfeng
[考研] 中南大學化學學碩337求調劑 +3	niko- 2026-03-19	6/300	2026-03-20 21:58 by luoyongfeng
[考研] 材料學碩297已過四六級求調劑推薦 +11	adaie 2026-03-19	11/550	2026-03-20 21:30 by laoshidan
[考研] 一志愿吉林大學材料學碩321求調劑 +11	Ymlll 2026-03-18	15/750	2026-03-20 19:40 by 丁丁*
[考研] 286求調劑 +6	lemonzzn 2026-03-16	10/500	2026-03-19 14:31 by lemonzzn
[考研] 【同濟軟件】軟件（085405）考研求調劑 +3	2026eternal 2026-03-18	3/150	2026-03-18 19:09 by 搏擊518