當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

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)有54人回復(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的回帖

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

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

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

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

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

[考研] 328求調(diào)劑，英語六級(jí)551，有科研經(jīng)歷 +5	生物工程調(diào)劑 2026-03-17	9/450	2026-03-21 23:32 by zhujy1982
[考研] 326求調(diào)劑 +5	諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15	8/400	2026-03-21 19:33 by ColorlessPI
[考研] 326求調(diào)劑 +4	mlpqaz03 2026-03-15	4/200	2026-03-21 19:10 by ColorlessPI
[考研] 工科0856求調(diào)劑 +3	沐析汀汀 2026-03-21	3/150	2026-03-21 18:30 by 學(xué)員8dgXkO
[考研] 求調(diào)劑 +3	.m.. 2026-03-21	4/200	2026-03-21 16:25 by barlinike
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6	77chaselx 2026-03-20	6/300	2026-03-21 07:24 by JourneyLucky
[考研] 機(jī)械專碩299求調(diào)劑至材料 +3	kkcoco25 2026-03-16	4/200	2026-03-21 03:52 by JourneyLucky
[考研] 307求調(diào)劑 +3	wyyyqx 2026-03-17	3/150	2026-03-21 03:20 by JourneyLucky
[考研] 303求調(diào)劑 +5	睿08 2026-03-17	7/350	2026-03-21 03:11 by JourneyLucky
[考研] 265求調(diào)劑 +9	梁梁校校 2026-03-17	9/450	2026-03-21 02:17 by JourneyLucky
[考研] 324分 085600材料化工求調(diào)劑 +4	llllkkkhh 2026-03-18	4/200	2026-03-21 01:24 by JourneyLucky
[考研] 321求調(diào)劑 +9	何潤(rùn)采123 2026-03-18	11/550	2026-03-20 23:19 by JourneyLucky
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3	我想要讀研究生 2026-03-20	3/150	2026-03-20 20:13 by JourneyLucky
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生（接收計(jì)算機(jī)專業(yè)調(diào)劑） +3	123阿標(biāo) 2026-03-17	3/150	2026-03-20 15:58 by 飛行琦
[考研] 085601材料工程專碩求調(diào)劑 +10	慕寒mio 2026-03-16	10/500	2026-03-19 15:26 by 丁丁*
[考研] 材料專碩306英一數(shù)二 +10	z1z2z3879 2026-03-16	13/650	2026-03-18 14:20 by 007_lilei
[考研] 318求調(diào)劑 +3	Yanyali 2026-03-15	3/150	2026-03-16 16:41 by houyaoxu
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6	太想進(jìn)步了0608 2026-03-16	6/300	2026-03-16 16:13 by kykm678
[考研] 327求調(diào)劑 +6	拾光任染 2026-03-15	11/550	2026-03-15 22:47 by 拾光任染
[考研] 0856專碩279求調(diào)劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015