tshuru1.asv

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

clear all;clc
% 输入
% 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;                                      % 粒度区间离散化

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;                                                            % 请输入停留时间
%n=input('请输入n(取1-50):n=');
for n=1:10;
TJ=importdata('1tiaojian.xls');
T=TJ(n,2);Nk=TJ(n,3);ak0=TJ(n,4);Sg=TJ(n,5);
A0=1.645*Nk/ak0;
sim('testmultitank1.mdl');                                          % 运行模型
% 模型计算后的粒度转化及与测试粒度的对比
n1=[A;N(end,:)];                                                    % 模型运行后粒度转化质量百分比
C=LEFT(:,1);
n2=interp1(n1(1,:),n1(2,:),C)                                       % 插值恢复到初始的粒度分布区间
for j=1:length(C)
    v2(j)=n2(j)*pi/6*C(j)^3*2.42e-12;                               % 粒度转化为质量百分比
   
end
v2=v2';
v0(1)=v2(1);                                                        % 首末计算为nan特殊处理
v0(2)=v2(2);
for i=3:length(C)
    v0(i)=(v0(i-1)+v2(i))                                               % 转化为累计百分比
end
v0=v0/(sum(v2(2:end-1)))*100 ;                                     % 累计大于100，重新计算百分比
v0(1)=0;v0(end)=100;
v0=v0';
v3=v2-LEFT2(:,2);   
v4=abs(v3)./v2*100

plot(C,v2,'r',C,LEFT2(:,2),'b')
legend('计算','初始')

% xlabel('粒度');ylabel('质量百分数');title('种分槽的粒度分布');