nongdu.asv

氧化铝生产种分过程粒度预测

clear all;
% 输入
% J=input('请输入溶解度方程常数:J=');                             % 输入溶解度方程常数
% H=input('请输入溶解度方程常数:H=');                             % 输入溶解度方程常数
% K=input('请输入溶液分解前的指前因子:K=');                       % 输入溶液分解前的指前因子
% At=input('请输入晶种的比表面积:At=');                           % 输入晶种的比表面积
% E=input('请输入溶液分解反应活化能:E=');                         % 输入活化能
% K2=input('请输入生长速率方程常数:K2=');                         % 输入生长速率方程常数
% E2=input('请输入晶体生长活化能:E2=');                           % 输入晶体生长活化能
% CA0=input('请输入初始铝酸钠溶液浓度:CA0=');                     % 输入初始铝酸钠溶液浓度
% Nk=input('请输入苛碱浓度:Nk=');                                 % 输入苛碱浓度
% t=input('请输入模型运行时间:t=');                               % 输入模型运行时间
% T1=input('请输入槽1的温度:T1=');                                % 输入首槽温度
% u=T1;
% W=input('请输入划分的离散的粒度区间数:W=');                     % 输入划分的离散的粒度区间数
% C=(425-0.3)/W;

% A=0.6:(450-0.6)/(W-1):450;                                      % 粒度区间离散化
%n=input('请输入n(取1-50):n=');
for n=1:50;
   
TJ=importdata('jztiaojian.xls');
T=TJ(n,2);Nk=TJ(n,3);ak0=TJ(n,4);Sg=TJ(n,5);aks
for i=1:1:30
    A(i)=0.439*2^(1/3*i);
end

LEFT=importdata('LEFT1.xls');                                     % 读取1号出口的粒度分布
LEFT2=importdata('LEFT2.xls');                                    % 读取2号出口的粒度分布
B=importdata('Bshuzu.xls');                                       
% B(:,:)=1e-18;                                                    

yi=interp1(LEFT(:,1),LEFT(:,3),A');                               % 直接插值计区节点的初始的粒数值
N1=abs(yi(1));
N2=abs(yi(2));
N3=abs(yi(3));
N4=abs(yi(4));                                                  
N5=abs(yi(5));
N6=abs(yi(6));
N7=abs(yi(7));
N8=abs(yi(8));
N9=abs(yi(9));
N10=abs(yi(10));
N11=abs(yi(11));
N12=abs(yi(12));
N13=abs(yi(13));
N14=abs(yi(14));
N15=abs(yi(15));
N16=abs(yi(16));
N17=abs(yi(17));
N18=abs(yi(18));
N19=abs(yi(19));
N20=abs(yi(20));
N21=abs(yi(21));
N22=abs(yi(22));
N23=abs(yi(23));
N24=abs(yi(24));
N25=abs(yi(25));
N26=abs(yi(26));
N27=abs(yi(27));
N28=abs(yi(28));
N29=abs(yi(29));
N30=abs(yi(30));
Rn=2e5;                                                         % 输入成核速率,Rn=2e5

%B=1e-18;                                                       % 输入附聚核
S=4;                                                            % 请输入停留时间

sim('testmultitank1.mdl');                                          % 运行模型