tshuru.m

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

% 输入
% 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=');
TJ=importdata('jztiaojian.xls');
T=TJ(n,2);Nk=TJ(n,3);ak=TJ(n,4);Sg=TJ(n,5);
for i=1:1:30
    A(i)=0.439*2^(1/3*i);
end

LEFT=importdata('jzlidu.xls');                                    % 读取初始晶种的粒度分布
B=importdata('Bshuzu.xls'); 

yi=interp1(LEFT(:,1),LEFT(:,n+1),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=2e7;                                                         % 输入成核速率,Rn=2e5

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

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


%%%                                                             % 画出粒度分布图
% P=[y11(170,2) y12(170,2) y13(170,2) y14(170,2) y15(170,2) y16(170,2) y17(170,2) y18(170,2) y19(170,2) y110(170,2) y111(170,2) y112(170,2) y113(170,2) y114(170,2) y115(170,2) y116(170,2) y117(170,2) y118(170,2) y119(170,2) y120(170,2)];
% P1=[y11(170,:) y12(170,:) y13(170,:) y14(170,:) y15(170,:) y16(170,:) y17(170,:) y18(170,:) y19(170,:) y110(170,:) y111(170,:) y112(170,:) y113(170,:) y114(170,:) y115(170,:) y116(170,:) y117(170,:) y118(170,:) y119(170,:) y120(170,:)];
% plot(A,P);
% surf([A;P1]);
% hold on
% plot(A,yi)
% xlabel('粒度');
% ylabel('粒子数');
% title('种分槽的粒度分布');