有哪些信誉好的足球投注网站- 1、本文档共15页，可阅读全部内容。
- 2、有哪些信誉好的足球投注网站（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
化学反应工程_四釜串连数据处理_matlab程序
四釜串联停留时间的分布的测定数据处理：
化工03-1 方培林2003031116
数据处理：
将实验数据导入并作图： data1
x=data1(:,1); y1=data1(:,2); y2=data1(:,3); y3=data1(:,4); y4=data1(:,5);
plot(x,y1,x,y2,x,y3,x,y4)
以第四个反应釜为例：
将实验数据V(t)-t曲线拟合：
（1）拟合曲线：
用高斯三阶拟合求曲线表达式：
y=0.534*exp(-((x-352.8)/222.3)^2)+0.4597*exp(-((x-218.4)/127.8)^2)+ 0.483*exp(-((x-710.6)/999.4 )^2)
积分求面积：
A1=int(0.534*exp(-((x-352.8)/222.3)^2)+0.4597*exp(-((x-218.4)/127.8)^2)+ 0.483*exp(-((x-710.6)/999.4 )^2),0, 846.72)
得：A1= 669.46898835250308069810477106541
减去空白面积：A=A1-846.72*0.45=288.4449 由实验数据分析终点时间为：846.72 在第385组数据
ET=(y4(1:385,:)-0.45)/A;
pt=sum(x(1:385,:).*ET)/sum(ET) %平均停留时间
sigmaf=sum(x(1:385,:).^2.*ET)/sum(ET)-pt^2
sig=sigmaf/pt^2 n=1/sig
A=288.4449 pt =344.33458 sigmaf =23259.5671 sig=.19617 n=5.0975
面积：A=288.4449 平均停留时间：tm= 344.33458 方差：
无因次方差： 小于1符合实际流型：相当于釜数N=5.0975
检验归一性
jifen=int((0.534*exp(-((x-352.8)/222.3)^2)+0.4597*exp(-((x-218.4)/127.8)^2)+ 0.483*exp(-((x-710.6)/999.4 )^2)-.45)/288.4449,0,846.72)
jifen=1.0000003063063450964052571949284
计算Peclet准数:取其中一个釜进行分析为开放体系：
解方程：per=solve(8*x^2+2*x-0 .19617=0) 得两组解：
per =[ -.32536529140547271941324376838209]（舍去） [ .75365291405472719413243768382092e-1]
Peclet准数为：0. 75365
作E(t)-t关系图：
fplot((0.534*exp(-((x-352.8)/222.3)^2)+0.4597*exp(-((x-218.4)/127.8)^2)+ 0.483*exp(-((x-710.6)/999.4 )^2)-.45)/288.4449,[0 846.72])
作F(t)-t关系图：
syms x
FT=int((0.534*exp(-((x-352.8)/222.3)^2)+0.4597*exp(-((x-218.4)/127.8)^2)+ 0.483*exp(-((x-710.6)/999.4 )^2)-.45)/288.4449)
function y=myfun(x)
y=217535077020598272/1057161738408951875*pi^(1/2)*erf(10/2223*x-392/247)+538299452482387968/5285808692044759375*pi^(1/2)*erf(5/639*x-364/213)+126367970891923456/151023105486993125*pi^(1/2)*erf(5/4997*x-187/263)-13194139533312/8457293907271615*x%将积分结果，定义为函数myfun
fplot(myfun,[0 846.72]) 第四个反应釜F(t)-t曲线：
将第三个实验数据V(t)-t曲线拟合：
用高斯四阶拟合求曲线表达式：
f=0.7926*exp(-((x-212.2)/126.5)^2)+0.7935*exp(-((x-304.4)/256)^2+0.5084*exp(-((x-127)/72)^2)+ 0.5974*exp(-((x-867.4)/ 825.4)^2)
积分求面积：
A1=int(0.7926*exp(-((
文档评论（0）