 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 424 | 回復(fù): 1 | |||
anathma新蟲 (初入文壇)
|
[求助]
用matlab擬合積分方程- 急用！坐等！非常感謝！
|
|
方程如下， 自變量是T，因變量是Cp。Tt 和ΔH 是擬合常數(shù)。請用 β(T) 對T作圖，并根據(jù)以上方程擬合得出Tt 和ΔH的值。（Tend是數(shù)據(jù)的最后一個T值。） 以上是對于曲線只有一個峰值的情況，對于兩個峰的情況，方程有一點小改動，如下， 作圖要求跟第一種情況一樣。因為本人對于matlab知之甚少，哪怕是一點小改動對于我來說都是座高山，所以只能勞煩高手了。附件里有單峰和雙峰的數(shù)據(jù)舉例。第一列是T值，第二列是Cp值。 另外以下是我曾經(jīng)找一個朋友幫忙做的程序，當(dāng)時對方程的理解錯誤，把Cp完全積分出來作為Y值擬合的，結(jié)果在Tt轉(zhuǎn)折點附近總是積分不好，程序和擬合結(jié)果如下，（相應(yīng)的數(shù)據(jù)也在附件里面，data.txt)。希望能有參考價值。 clc clear all close all d = load('data.txt', '-ASCII'); XX = d(:,1); YY = d(:,2); %idx = (XX>280); idx = (XX>304); X = XX(idx) Y = YY(idx); R = 8.314; f = @(a,x)a(1)*(1-exp( a(2)/R * (1./x - 1/a(3))) ) ; % hyperbolic-like %f = @(a,x)a(1) * (1-(1-exp( a(2)/R * (1./x - 1/a(3))) ) ./ (1 + exp( a(2)/R * (1./x-1/a(3))) )) ; a0 = [1400 220000 304] [a1,r] = nlinfit(X, Y, f, a0); A = a1(1) B = a1(2) T = a1(3) mse = sum(r.^2)/(length(r)-length(a1)) plot(XX,YY,'k', XX, f(a1,XX), 'r') ylabel('$\Delta \hat{H}(T) \quad / \quad kCal \cdot mol^{-1} \cdot K^{-1}$','Interpreter','LaTex') xlabel('$T \quad / \quad K$','interpreter','latex') legend('data', 'fitted param') axis([XX(1), XX(end), YY(1), YY(end)]) |
回復(fù)此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : 雙峰.txt
- 附件 2 : 單峰.txt
- 附件 3 : data.txt
2014-08-21 02:26:38, 12.11 K
2014-08-21 03:16:56, 6.1 K
2014-08-21 03:17:38, 9.51 K
» 猜你喜歡
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有158人回復(fù)
-
請問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有32人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運，祝面上順利！
已經(jīng)有76人回復(fù)
-
華東師范大學(xué)芯片設(shè)計徐瓏真實水平咋樣
已經(jīng)有1人回復(fù)
-
調(diào)劑
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
中德博士后2026項目5月30日截止- 國內(nèi)博士后入站、德國兩年全職
已經(jīng)有0人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 最小二乘法的matlab 實現(xiàn)
已經(jīng)有4人回復(fù)
- 求高版本1stopt跑一個程序
已經(jīng)有7人回復(fù)
- 求數(shù)學(xué)大神用matlab幫我解決一個問題。不勝感激
已經(jīng)有6人回復(fù)
- 求助用matlab擬合動力學(xué)方程
已經(jīng)有9人回復(fù)
- 用origin8.0英文版或7.5漢化版對Langmuir進行數(shù)據(jù)的擬合，如何操作。本人就一組數(shù)據(jù)。
已經(jīng)有3人回復(fù)
- mathcad怎樣用四階龍格庫塔法求解二階微分方程？
已經(jīng)有8人回復(fù)
- 回歸分析，急用！
已經(jīng)有9人回復(fù)
- 求前輩幫忙用1stOpt你個一下數(shù)據(jù)，寫論文急用，先謝過
已經(jīng)有14人回復(fù)
- 如何用軟件繪制米氏方程曲線啊？？急求
已經(jīng)有3人回復(fù)
- MATLAB微分方程參數(shù)擬合問題，求大神
已經(jīng)有7人回復(fù)
- 求助 用matlab擬合三元線性回歸方程及分析
已經(jīng)有9人回復(fù)
- 請問matlab擬合二元方程的問題
已經(jīng)有3人回復(fù)
- 1stopt或matlab如何求解以下的非線性方程并擬合出相應(yīng)曲線？
已經(jīng)有13人回復(fù)
- matlab 擬合反應(yīng)動力學(xué)參數(shù)結(jié)果很差。大家?guī)兔匆幌?/a>
已經(jīng)有14人回復(fù)
- matlab微分方程組參數(shù)擬合，以周為單位求解，匯總后以年為單位進行數(shù)值比較
已經(jīng)有4人回復(fù)
- matlab擬合方程參數(shù)時初值的選擇
已經(jīng)有15人回復(fù)
- matlab 求指點 動力學(xué)方程擬合過程中導(dǎo)數(shù)的獲取
已經(jīng)有10人回復(fù)
- 【公式擬合問題】如何擬合下面這樣的一個公式，多謝大牛。
已經(jīng)有6人回復(fù)
- 急切求助用matlab曲線擬合色散方程
已經(jīng)有10人回復(fù)
- 【求助】使用Matlab擬合反應(yīng)動力學(xué)方程問題
已經(jīng)有7人回復(fù)
anathma
新蟲 (初入文壇)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 161.6
- 帖子: 14
- 在線: 69小時
- 蟲號: 1533615
- 注冊: 2011-12-12
- 專業(yè): 高分子物理與高分子物理化
2樓2014-08-21 22:29:33
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 anathma 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑一志愿武漢理工大學(xué)材料工程（085601） +3 | WW.' 2026-03-23 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 招08考數(shù)學(xué) +6 | laoshidan 2026-03-20 | 14/700 |
|
|
[考研] 306求調(diào)劑 +9 | chuanzhu川燭 2026-03-18 | 9/450 |
|
|
|
[考研] 070300化學(xué)求調(diào)劑 +5 | 苑豆豆 2026-03-20 | 5/250 |
|
|
|
[考研] 08工科 320總分 求調(diào)劑 +11 | 梨花珞晚風(fēng) 2026-03-17 | 11/550 |
|
|
|
[考研] 一志愿中南化學(xué)（0703）總分337求調(diào)劑 +9 | niko- 2026-03-19 | 10/500 |
|
|
|
[考研] 一志愿華中科技大學(xué)071000，求調(diào)劑 +4 | 沿岸有貝殼6 2026-03-21 | 4/200 |
|
|
|
[考研] 280求調(diào)劑 +11 | 咕嚕曉曉 2026-03-18 | 12/600 |
|
|
|
[考研] 一志愿深大，0703化學(xué)，總分302，求調(diào)劑 +4 | 七月-七七 2026-03-21 | 4/200 |
|
|
|
[考研] 材料學(xué)碩333求調(diào)劑 +3 | 北道巷 2026-03-18 | 3/150 |
|
|
|
[考研] 299求調(diào)劑 +5 | shxchem 2026-03-20 | 7/350 |
|
|
|
[考研] 316求調(diào)劑 +6 | 梁茜雯 2026-03-19 | 6/300 |
|
|
|
[考研] 二本跨考鄭大材料306英一數(shù)二 +3 | z1z2z3879 2026-03-17 | 3/150 |
|
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-18 | 6/300 |
|
|
|
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +4 | @taotao 2026-03-19 | 4/200 |
|
|
|
[考研] 353求調(diào)劑 +3 | 拉鉤不許變 2026-03-20 | 3/150 |
|
|
|
[考研] 261求B區(qū)調(diào)劑，科研經(jīng)歷豐富 +3 | 牛奶很忙 2026-03-20 | 4/200 |
|
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
|
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|
|
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|