 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 422 | 回復: 1 | ||
anathma新蟲 (初入文壇)
|
[求助]
用matlab擬合積分方程- 急用！坐等！非常感謝！
|
|
方程如下， 自變量是T，因變量是Cp。Tt 和ΔH 是擬合常數(shù)。請用 β(T) 對T作圖，并根據(jù)以上方程擬合得出Tt 和ΔH的值。（Tend是數(shù)據(jù)的最后一個T值。） 以上是對于曲線只有一個峰值的情況，對于兩個峰的情況，方程有一點小改動，如下， 作圖要求跟第一種情況一樣。因為本人對于matlab知之甚少，哪怕是一點小改動對于我來說都是座高山，所以只能勞煩高手了。附件里有單峰和雙峰的數(shù)據(jù)舉例。第一列是T值，第二列是Cp值。 另外以下是我曾經找一個朋友幫忙做的程序，當時對方程的理解錯誤，把Cp完全積分出來作為Y值擬合的，結果在Tt轉折點附近總是積分不好，程序和擬合結果如下，（相應的數(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)]) |
回復此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內容負責。
本內容由用戶自主發(fā)布，如果其內容涉及到知識產權問題，其責任在于用戶本人，如對版權有異議，請聯(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
» 猜你喜歡
0702一志愿吉大B區(qū)求調劑
已經有5人回復
-
求調劑
已經有0人回復
-
物理學I論文潤色/翻譯怎么收費?
已經有239人回復
-
0702一志愿吉大B區(qū)求調劑有論文
已經有0人回復
-
請問還有沒有用Latex寫文章的小伙伴們？
已經有0人回復
-
光學工程學碩調劑信息
已經有31人回復
-
歡迎加入課題組
已經有0人回復
-
散金幣，求好運，祝面上順利！
已經有74人回復
-
華東師范大學芯片設計徐瓏真實水平咋樣
已經有1人回復
-
調劑
已經有0人回復
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經有0人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 最小二乘法的matlab 實現(xiàn)
已經有4人回復
- 求高版本1stopt跑一個程序
已經有7人回復
- 求數(shù)學大神用matlab幫我解決一個問題。不勝感激
已經有6人回復
- 求助用matlab擬合動力學方程
已經有9人回復
- 用origin8.0英文版或7.5漢化版對Langmuir進行數(shù)據(jù)的擬合，如何操作。本人就一組數(shù)據(jù)。
已經有3人回復
- mathcad怎樣用四階龍格庫塔法求解二階微分方程？
已經有8人回復
- 回歸分析，急用！
已經有9人回復
- 求前輩幫忙用1stOpt你個一下數(shù)據(jù)，寫論文急用，先謝過
已經有14人回復
- 如何用軟件繪制米氏方程曲線��？？急求
已經有3人回復
- MATLAB微分方程參數(shù)擬合問題，求大神
已經有7人回復
- 求助 用matlab擬合三元線性回歸方程及分析
已經有9人回復
- 請問matlab擬合二元方程的問題
已經有3人回復
- 1stopt或matlab如何求解以下的非線性方程并擬合出相應曲線？
已經有13人回復
- matlab 擬合反應動力學參數(shù)結果很差。大家?guī)兔匆幌?/a>
已經有14人回復
- matlab微分方程組參數(shù)擬合，以周為單位求解，匯總后以年為單位進行數(shù)值比較
已經有4人回復
- matlab擬合方程參數(shù)時初值的選擇
已經有15人回復
- matlab 求指點 動力學方程擬合過程中導數(shù)的獲取
已經有10人回復
- 【公式擬合問題】如何擬合下面這樣的一個公式，多謝大牛。
已經有6人回復
- 急切求助用matlab曲線擬合色散方程
已經有10人回復
- 【求助】使用Matlab擬合反應動力學方程問題
已經有7人回復
anathma
新蟲 (初入文壇)
- 應助: 1 (幼兒園)
- 金幣: 161.6
- 帖子: 14
- 在線: 69小時
- 蟲號: 1533615
- 注冊: 2011-12-12
- 專業(yè): 高分子物理與高分子物理化
2樓2014-08-21 22:29:33
相關版塊跳轉
我要訂閱樓主 anathma 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 276求調劑 +3 | YNRYG 2026-03-21 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 石河子大學（211、雙一流）碩博研究生長期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|
|
|
[考研] 285求調劑 +6 | ytter 2026-03-22 | 6/300 |
|
|
|
[考研] 354求調劑 +7 | Tyoumou 2026-03-18 | 10/500 |
|
|
|
[考研] 296求調劑 +4 | www_q 2026-03-20 | 4/200 |
|
|
|
[基金申請]
學校已經提交到NSFC，還能修改嗎？
 40+4
|
babangida 2026-03-19 | 9/450 |
|
|
|
[考研] 083200學碩321分一志愿暨南大學求調劑 +3 | innocenceF 2026-03-17 | 3/150 |
|
|
|
[考研] 材料專碩英一數(shù)二306 +7 | z1z2z3879 2026-03-18 | 7/350 |
|
|
|
[考研] 一志愿南京理工大學085701資源與環(huán)境302分求調劑 +4 | 葵梓衛(wèi)隊 2026-03-18 | 6/300 |
|
|
|
[考研] 288求調劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
|
[考研]
 材料與化工 322求調劑
+4
|
然11 2026-03-19 | 4/200 |
|
|
|
[考研] 329求調劑 +9 | 想上學吖吖 2026-03-19 | 9/450 |
|
|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調劑 +5 | 西南交通專材355 2026-03-19 | 5/250 |
|
|
|
[考研] 材料學求調劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
|
[考研] 261求B區(qū)調劑，科研經歷豐富 +3 | 牛奶很忙 2026-03-20 | 4/200 |
|
|
|
[考研] 一志愿南理工085701環(huán)境302求調劑院校 +3 | 葵梓衛(wèi)隊 2026-03-20 | 3/150 |
|
|
|
[考研] 085410人工智能專碩317求調劑（0854都可以） +4 | xbxudjdn 2026-03-18 | 4/200 |
|
|
|
[考研] 288求調劑，一志愿華南理工大學071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
|
[考研] 本科鄭州大學物理學院，一志愿華科070200學碩，346求調劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|