heshuangrong

新蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 784.2
散金: 25
帖子: 96
在線: 93.5小時(shí)
蟲號: 1498593
注冊: 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論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有130人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請問還有沒有用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à)值貼推薦，對您同樣有幫助:

向大神求教一個(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é)級數(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í)
蟲號: 907551
注冊: 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í)
蟲號: 1498593
注冊: 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í)
蟲號: 907551
注冊: 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í)
蟲號: 907551
注冊: 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í)
蟲號: 1498593
注冊: 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í)
蟲號: 907551
注冊: 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ā)表

[考研] 本人考085602 化學(xué)工程專碩 +21	不知道叫什么！ 2026-03-15	23/1150	2026-03-22 00:12 by BruceLiu320
[考研] 313求調(diào)劑 +4	肆叁貳壹22 2026-03-19	4/200	2026-03-21 17:33 by ColorlessPI
[考研] 求調(diào)劑 +3	.m.. 2026-03-21	4/200	2026-03-21 16:25 by barlinike
[考研] 332求調(diào)劑 +3	鳳凰院丁真 2026-03-20	3/150	2026-03-21 10:27 by luoyongfeng
[考研] 306求調(diào)劑 +4	chuanzhu川燭 2026-03-18	4/200	2026-03-21 08:25 by laoshidan
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6	77chaselx 2026-03-20	6/300	2026-03-21 07:24 by JourneyLucky
[考研] 265求調(diào)劑 +9	梁梁校校 2026-03-17	9/450	2026-03-21 02:17 by JourneyLucky
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +5	守候夕陽CF 2026-03-18	5/250	2026-03-21 01:06 by JourneyLucky
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +9	Ncdx123456 2026-03-19	9/450	2026-03-20 23:41 by lovewei0727
[考研] 294求調(diào)劑材料與化工專碩 +15	陌の森林 2026-03-18	15/750	2026-03-20 23:28 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
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9	Doleres 2026-03-19	9/450	2026-03-20 22:36 by JourneyLucky
[考研] 324求調(diào)劑 +5	lucky呀呀呀鴨 2026-03-20	5/250	2026-03-20 22:30 by 促天成
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3	@taotao 2026-03-19	3/150	2026-03-20 21:34 by JourneyLucky
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3	葵梓衛(wèi)隊(duì) 2026-03-20	3/150	2026-03-20 19:28 by zhukairuo
[考研] 工科材料085601 279求調(diào)劑 +7	困于星晨 2026-03-17	9/450	2026-03-20 17:38 by 無懈可擊111
[考研] 0703化學(xué)調(diào)劑 +4	18889395102 2026-03-18	4/200	2026-03-19 16:13 by 30660438
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4	我不是一根蔥 2026-03-18	4/200	2026-03-19 09:11 by 浮云166
[考研] 收復(fù)試調(diào)劑生 +4	雨后秋荷 2026-03-18	4/200	2026-03-18 14:16 by elevennnne
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6	太想進(jìn)步了0608 2026-03-16	6/300	2026-03-16 16:13 by kykm678