计算机在化工中的应用实验报告

《计算机在化工中的应用实验报告》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、《计算机在化学化工中的应用》实验报告学院：化学与化工学院班级：12级硕勋励志班姓名：徐凯杰学号120702028实验一传热实验中多变量的曲线的拟合一、实验目的1)熟悉VB编程平台2)掌握多变S:曲线拟合的算法3)编拟合所给的传热实验模型的VB程序4)通过实验数据求出模型数据、并掌握解线性方程组的克拉默法则二、运行环境1)MicrosoftWindowsXP2)VB6.0三、实验原理略四、vb代码PrivateSubCommandl_Click()DimmAsInteger,m=inputbox(“实验次数”)m=7DimxlO,x20,yO

2、Dimi，j，kAsIntegerDima(lTo10,1To10)，y(lTo10)，y1(1To10)，a0,al，a2Dims,SI,S2,S3,b(lTo10,1To10),xxDimxl(lTo10),x2(lTo10),YY,sd,open"dem.dat,'forinputas#lTori=1tom’input#l，xx，YY'xl(i)=xx•x2(iO=xxA21y(i)=YY’nexti•close#17组努塞尔准数、雷诺数及普兰德准数，数据最大时应采用直接从文件读取方法xl0=Array(0,100,200,300,5

3、00,100,700,800)1注意下标的起点处理(加0)x20=Array(0,2,4,1，0.3,5,3,4)'注意下标的起点处理(加0)y0=Array(0,1.127,2.416,2.205,2.312,1.484,6.038,7.325)•注意下标的起点处理(加0)Fori=1Tomxl(i)=Log(xl0(i))x2(i)=Log(x20(i))y(i)=Log(yO(i))Nexti1求解法A程系数矩阵a(1，1)=ma(l,2)=0Fori=1Toma(l，2)=a(l，2)+xl(i)Nextia(2,1)=a(l,2)

4、a(l，3)=0Fori=1Toma(l，3)=a(l，3)+x2(i)Nextia(3,1)=a(l,3)a(2,2)=0Fori=1Toma(2,2)=a(2,2)+xl(i)*xl(i)Nextia(3,3)=0Fori=1Toma(3,3)=a(3,3)+x2(i)*x2(i)Nextia(2,3)=0Fori=1Toma(2,3)=a(2,3)+xl(i)*x2(i)Nextia(3,2)=a(2,3)*求解法方程常数向量yl⑴=0Fori=1Tomyl(l):yl⑴+y(i)Nextiyi⑵Fori=1Tomyl(2)=yl(2

5、)+xl(i)*y(i)Nextiyl⑶=0Fori=1Tomyl(3)=yl(3)+x2(i)*y(i)Nexti利用克拉默法则解法方程/线性lh常组)s=a(l,1)*a(2,2)*a(3,3)+a(l,2)*a(2,3)*a(3,l)+a(l,3)*a(2,1)*a(3,2)s=s-a(l,1)*a(2,3)*a(3,2)-a(l,2)*a(2,1)*a(3,3)-a(l,3)*a(2,2)*a(3,1)Forj=1To3b(j,l)=a(j,1)a(j,l)=yl①NextjSI=a(l,1)*a(2,2)*a(3,3)+a(l,

6、2)*a(2,3)*a(3,1)+a(l,3)*a(2,1)*a(3,2)51=SI-a(l,1)*a(2,3)*a(3,2)•a(l,2)*a(2,1)*a(3,3)-a(l,3)*a(2,2)*a(3,1)Forj=1To3a(j，1)=b(j,1)NextjForj=1To3b(j,2)=a(j,2)a(j，2)=yl①Nextj51=a(l,1)*a(2,2)*a(3,3)+a(l,2)*a(2,3)*a(3,1)+a(l,3)*a(2,1)*a(3,2)52=S2-a(l,1)*a(2,3)*a(3,2)-a(l,2)*a(2,1

7、)*a(3,3)-a(l,3)*a(2,2)*a(3,1)Forj=1To3a(j,2)=b(j,2)NextjForj=1To3b(j,3)=a(j,3)a(j，3)=yl⑴Nextj53=a(l,1)*a(2,2)*a(3,3)+a(l,2)*a(2,3)*a(3,1)+a(l,3)*a(2,1)*a(3,2)S3=S3-a(l,1)*a(2,3)*a(3,2)-a(l,2)*a(2,1)*a(3,3)-a(l,3)*a(2,2)*a(3,1)aO=Sl/sal=S2/sa2=S3/sText1.Text=Int(lOOO*Exp(aO

8、)+0.5)/1000•四舍五入保留三位Text2.Text=Int(1000*al+0.5)/1000Text3.Text=Int(1000*a24-0.5)/1000sd=O