1、一元线性拟合
求HNO 3的正常沸点温度T b 及摩尔汽化热。

程序如下：
>> t=[0 20 40 50 70 80 90 100]; >> t=t+273.15;
>> p=[1919.52 6385.07 17728.9 27726.4 62251.1 89311 124902.1 170890.6] p =
1.0e+005 *
0.0192 0.0639 0.1773 0.2773 0.6225 0.8931 1.2490 1.7089 >> subplot 121 >> plot(t,p,'o',t,p) >> t1=1./t;p2=log(p); >> pp=polyfit(t1,p2,1) pp =
1.0e+003 *
-4.5691 0.0243 >> subplot 122
>> plot(t1,p2,'o',t1,p2)
>> gtext('p/pa'),gtext('T/K'),GTEXT('lnP/Pa'),gtext('T^-^1/K')
由克拉贝龙-克劳修斯方程式，~
ln v
H P C RT
∆=-+ 作1
ln ~P T
-得一直线：3
1
ln 4.5691024.30P T -=-⨯+
斜率为：~
3
4.56910v
H R
∆-⨯=-
所以摩尔汽化热为：~
314.569108.31437.99()v H kJ mol -∆=⨯⨯=⋅ 并根据拟合方程，求得一大气压时 1
32.8010T
--=⨯
则正常沸点为：357b T K =
2、多元线性拟合：
某气体混合物由四种气体组成，在常压或低压下其粘度η与各组分摩尔分数x 1,x 2,x 3,x 4之间有如下线性关系：011223344b b x b x b x b x η=++++
试根据下表所列实验数据用最小二乘法确定上式中的各个系数，并计算其复相关系数。

Matlab 程序如下：
>> a=[1.0 0.402 0.153 0.058 0.387;1.0 0.503 0.301 0.183 0.013; 1.0 0.306 0.109 0.224 0.361; 1.0 0.296 0.365 0.009 0.330; 1.0 0.309 0.405 0.109 0.177; 1.0 0.055 0.153 0.506 0.289] a =
1.0000 0.4020 0.1530 0.0580 0.3870 1.0000 0.5030 0.3010 0.1830 0.0130 1.0000 0.3060 0.1090 0.2240 0.3610 1.0000 0.2960 0.3650 0.0090 0.3300 1.0000 0.3090 0.4050 0.1090 0.1770 1.0000 0.0550 0.1530 0.5060 0.2890
>> y=[0.00625 0.00826 0.01182 0.01944 0.02372 0.03243]' y =
0.0063 0.0083 0.0118 0.0194 0.0237 0.0324
>> b=a.'*a
b =
6.0000 1.8710 1.4860 1.0890 1.5570
1.8710 0.6944 0.4879 0.2481 0.4408
1.4860 0.4879 0.4466 0.2132 0.3388
1.0890 0.2481 0.2132 0.3550 0.2742
1.5570 0.4408 0.3388 0.2742 0.5040 >> c=a.'*y
c =
0.1019
0.0252
0.0264
0.0237
0.0268
>> coeff=b\c
coeff =
0.2267
-0.2627
-0.1620
-0.1889
-0.2058
>> regress(c,b,0.05)
ans =
0.2267
-0.2627
-0.1620
-0.1889
-0.2058。