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

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)有239人回復(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é)

引用回帖:

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

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

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é)

【答案】應(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ā)表

[基金申請] 山東省面上項(xiàng)目限額評審 +4	石瑞0426 2026-03-19	4/200	2026-03-22 08:50 by Wei_ren
[考研] 一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑 +4	vv迷 2026-03-21	5/250	2026-03-22 08:38 by ChemPharm
[考研] 0703化學(xué)調(diào)劑，六級已過，有科研經(jīng)歷 +14	曦熙兮 2026-03-15	14/700	2026-03-21 19:12 by ColorlessPI
[考研] 求調(diào)劑 +4	要好好無聊 2026-03-21	4/200	2026-03-21 18:57 by 學(xué)員8dgXkO
[考研] 0703化學(xué)297求調(diào)劑 +3	Daisy☆ 2026-03-20	3/150	2026-03-21 17:45 by ColorlessPI
[考研] 材料與化工（0856）304求 B區(qū) 調(diào)劑 +3	邱gl 2026-03-21	3/150	2026-03-21 13:47 by lature00
[考研] 材料學(xué)學(xué)碩080502 337求調(diào)劑-一志愿華中科技大學(xué) +4	順順順mr 2026-03-18	5/250	2026-03-21 10:22 by luoyongfeng
[考研] 316求調(diào)劑 +6	梁茜雯 2026-03-19	6/300	2026-03-21 06:32 by Ecowxq666！
[考研] 求調(diào)劑 +3	Ma_xt 2026-03-17	3/150	2026-03-21 02:05 by JourneyLucky
[考研] 296求調(diào)劑 +6	www_q 2026-03-18	10/500	2026-03-20 23:56 by JourneyLucky
[考研] 321求調(diào)劑 +9	何潤采123 2026-03-18	11/550	2026-03-20 23:19 by JourneyLucky
[考研] 一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +4	晨昏線與星海 2026-03-19	4/200	2026-03-20 22:15 by JourneyLucky
[考研] 260求調(diào)劑 +3	朱芷琳 2026-03-20	3/150	2026-03-20 20:35 by 學(xué)員8dgXkO
[考研] 319求調(diào)劑 +3	小力氣珂珂 2026-03-20	3/150	2026-03-20 19:47 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
[考研] 一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +5	1孫悟空 2026-03-17	6/300	2026-03-19 23:46 by zcl123
[考研] 085601材料工程專碩求調(diào)劑 +10	慕寒mio 2026-03-16	10/500	2026-03-19 15:26 by 丁丁*
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +10	Liwangman 2026-03-15	10/500	2026-03-19 10:25 by 無際的草原
[考研] 344求調(diào)劑 +6	knight344 2026-03-16	7/350	2026-03-18 20:13 by walc