matlab的計算過程模板

1、題目：用MATLAB計算樣品C17H15O2N3S2的焓H、熵S、吉布斯自由能G學(xué)生姓名: 何 建 云 學(xué) 號: 院 系: 化工學(xué)院 專 業(yè): 制藥工程 年 級: 2010級 用MATLAB計算樣品C17H15O2N3S2的焓H、熵S、吉布斯自由能G(化工學(xué)院 制藥工程 何建云 )經(jīng)過前期的熱分析計算已經(jīng)得出樣品C17H15O2N3S2的比熱Cp與溫度T的關(guān)系式如下：Cp=2.9423-1.531810-2T-+2.986810-5T2。現(xiàn)在需要計算樣品的焓、熵、吉布斯自由能，已知它們與比熱的關(guān)系式： a(0), a(1), a(2)是從圖上讀出。Cp=a(0)+a(1)*T+a(2)*T2+

2、 a(3)*T3= a(0)T+ a(1)*T2/2+ a(2)*T2/3+ a(3)*T3/4= a(0)log(T)+ a(1)*T+ a(2)*T2/2+ a(3)*T3/3有基本積分公式：將Cp=2.9423-1.531810-2T+2.986810-5T2帶入上述計算式子可得到：=2.9423.*(t-298.15)-0.*(t.2-298.15.2)+0.*(t.3-298.15.3)=2.9423.*(log(t)-log(298.15)-0.*(t-298.15)+0.*(t.2-298.15.2) 就代表自然對數(shù)ln 代表t的平方用MATLAB軟件計算樣品的H,S,G；其中T

3、以5K為間隔，從283K計算到353K。樣品的相對分子質(zhì)量為357。MATLAB計算過程如下： 下列紅色部分代表需要輸入或者重要的部分。藍色是說明文字。To get started, select MATLAB Help or Demos from the Help menu.The element type name must be terminated by the matching end-tag .Could not parse the file: d:matlab7toolboxccslinkccslinkinfo.xml t=283:5:353 起點溫度：間隔溫度：終點溫度。t =

4、 283 288 293 298 303 308 313 318 323 328 333 338 343 348 353 h=2.9423.*(t-298.15)-0.*(t.2-298.15.2)+0.*(t.3-298.15.3)h = Columns 1 through 10 -15.3578 -10.3396 -5.2744 -0.1545 5.0274 10.2789 15.6075 21.0205 26.5254 32.1298 Columns 11 through 15 37.8411 43.6667 49.6141 55.6908 61.9043 a=h.*0.357 注：0.

5、357代表的是相對分子質(zhì)量縮小1000倍；另外a就是最后所需要的焓值。a = Columns 1 through 10 -5.4827 -3.6913 -1.8830 -0.0552 1.7948 3.6696 5.5719 7.5043 9.4696 11.4703 Columns 11 through 15 13.5093 15.5890 17.7122 19.8816 22.0998 s=2.9423.*(log(t)-log(298.15)-0.*(t-298.15)+0.*(t.2-298.15.2)s = Columns 1 through 9 -0.0455 -0.0303 -0

6、.0153 -0.0004 0.0143 0.0289 0.0435 0.0580 0.0725 Columns 10 through 15 0.0870 0.1015 0.1161 0.1307 0.1454 0.1602 b=s.*0.357 注：0.357代表的是相對分子質(zhì)量縮小1000倍；另外b就是最后所需要的熵值。b = Columns 1 through 9 -0.0162 -0.0108 -0.0055 -0.0002 0.0051 0.0103 0.0155 0.0207 0.0259 Columns 10 through 15 0.0311 0.0362 0.0414 0.0467 0.0519 0.0572 c=a-t.*b 注：0.357代表的是相對分子質(zhì)量M縮小1000倍；另外c就是最后所需要的吉布斯自由能值。c = Columns 1 through 9 -0.8865 -0.5755 -0.2828 -0.0080 0.2493 0.4895 0.7128 0.9194 1.1097 Columns 10 through 15 1.2838 1.4417 1.5837 1.7097 1.8