无法在Jupyter和联机编译器中执行Matlab代码_Matlab

无法在Jupyter和联机编译器中执行Matlab代码

matlab

无法在Jupyter和联机编译器中执行Matlab代码,matlab,Matlab,我需要看到以下matlab代码的代码输出 %Defining the arrays i=1200; XbFe=zeros(1,i); XbCr=zeros(1,i); XbO2=zeros(1,i); XbFeO=zeros(1,i); XbCr2O3=zeros(1,i); CbFe=zeros(1,i); CbCr=zeros(1,i); CbO2=zeros(1,i); CbFeO=zeros(1,i); CbCr2O3=zeros(1,i); keq_FeO=zeros(1,i);

我需要看到以下matlab代码的代码输出


%Defining the arrays

i=1200;
XbFe=zeros(1,i);
XbCr=zeros(1,i);
XbO2=zeros(1,i);
XbFeO=zeros(1,i);
XbCr2O3=zeros(1,i);
CbFe=zeros(1,i);
CbCr=zeros(1,i);
CbO2=zeros(1,i);
CbFeO=zeros(1,i);
CbCr2O3=zeros(1,i);
keq_FeO=zeros(1,i);
keq_Cr2O3=zeros(1,i);
val=zeros(1,i);
XiFe=zeros(1,i);
XiCr=zeros(1,i);
CiFe=zeros(1,i);
CiCr=zeros(1,i);
JFe=zeros(1,i);
JCr=zeros(1,i);
MolesTransFe=zeros(1,i);
MolesTransCr=zeros(1,i);
MolesTransfFeO=zeros(1,i);
MolesTransfCr2O3=zeros(1,i);
MolesTransO2=zeros(1,i);
MFe=zeros(1,i);
MCr=zeros(1,i);
MO2InMetal=zeros(1,i);
Moxy=zeros(1,i);
MCr2O3=zeros(1,i);
MFeO=zeros(1,i);
MassOfFe=zeros(1,i);
MassOfCr=zeros(1,i);
MassofO2=zeros(1,i);
MassOfFeO=zeros(1,i);
MassOfCr2O3=zeros(1,i);
TotalMolesInMetal=zeros(1,i);
TotalMolesInSlag=zeros(1,i);
TotalMassOfMetal=zeros(1,i);
TotalMassOfSlag=zeros(1,i);
VM=zeros(1,i);
VS=zeros(1,i);
MolarVM=zeros(1,i);
MolarVS=zeros(1,i);    
PerFeInMetal=zeros(1,i);
PerCrInMetal=zeros(1,i);
PerO2InMetal=zeros(1,i);
PerFeOInSlag=zeros(1,i);
PerCr2O3InSlag=zeros(1,i);
time=zeros(1,i);
time(1)=0;

%User Input

TotalMassOfMetal(1) = input("Input Total Mass of Metal:" );
TotalMassOfSlag(1) = input("Input Total Mass of Slag: ");
DensityOfMetal(1) = input("Input Density of Metal: ");
DensityOfSlag(1) = input("Input Density of Slag: ");
PerFeInMetal(1) = input("Input Percentage of Fe in Metal: ");
PerCrInMetal(1) = input("Input Percentage of Cr in Metal: ");
PerO2InMetal(1)= input("Input Percentage of O2 in Metal: ");
PerFeOInSlag(1) = input("Input Percentage of FeO in Slag: ");
PerCr2O3InSlag(1) = input("Input Percentage of Cr2O3 in Slag: ");
O2Top(1) = input("Input Oxygen added from Top: ");
O2Bottom(1) = input("Input Oxygen added from bottom: ");
dt = input("Input time step of each calculation: ");
Temp=input("Initial Temperature:" );
MassTransCoeffFe=input("Mass transfer Coefficients of Fe: ");
MassTransCoeffCr=input("Mass transfer Coefficients of Cr: ");
A=input("Area Assumed: ");

%Standard Known Values/ Given Inputs

AtomicMassofFe= 56;
AtomicMassofCr= 52;
AtomicMassofO2= 16;
AtomicMassofFeO= 72;
AtomicMassofCr2O3= 152;

%Volume of Metal and Slag

VM(1)=TotalMassOfMetal(1)/DensityOfMetal;
VS(1)=TotalMassOfSlag(1)/DensityOfSlag;

%Initial Moles of Metal and Slag in kmoles

MFe(1)=TotalMassOfMetal(1)*(PerFeInMetal(1)/100)*(1/AtomicMassofFe);
MCr(1)=TotalMassOfMetal(1)*(PerCrInMetal(1)/100)*(1/AtomicMassofCr);
MFeO(1)=TotalMassOfSlag(1)*(PerFeOInSlag(1)/100)*(1/AtomicMassofFeO);
MCr2O3(1)=TotalMassOfSlag(1)*(PerCr2O3InSlag(1)/100)*(1/AtomicMassofCr2O3);
MO2InMetal(1)=TotalMassOfMetal(1)*(PerO2InMetal(1)/100)*(1/AtomicMassofO2);

%Moles of Oxygen Added in kmoles

MO2Added=((O2Bottom+O2Top)*2*(273/298)/(22.4*60));
Moxy(1)=MO2InMetal(1)+MO2Added;

%Total Moles in Metal and Slag in kMoles

TotalMolesInMetal(1)=MCr(1)+MFe(1)+Moxy(1);
TotalMolesInSlag(1)=MCr2O3(1)+MFeO(1);

%Molar Volume 

MolarVM(1)=TotalMolesInMetal(1)/VM(1);
MolarVS(1)=TotalMolesInSlag(1)/VS(1);

for i= 1:1200

    %Mole Fractions of elements
    
    XbFe(i)=MFe(i)/TotalMolesInMetal(i);
    XbCr(i)=MCr(i)/TotalMolesInMetal(i);
    XbO2(i)=Moxy(i)/TotalMolesInMetal(i);
    XbFeO(i)=MFeO(i)/TotalMolesInSlag(i);
    XbCr2O3(i)=MCr2O3(i)/TotalMolesInSlag(i);

    %Bulk Concentration of elements

    CbFe(i)=XbFe(i)*MolarVM(i);
    
    CbCr(i)=XbCr(i)*MolarVM(i);
    
    CbO2(i)=XbO2(i)*MolarVM(i);
    
    CbFeO(i)=XbFeO(i)*MolarVS(i);
    
    CbCr2O3(i)=XbCr2O3(i)*MolarVS(i);
  

    %Equilibriium Constant Identfication

    keq_FeO(i)=exp((-121009.9+(53.114*Temp)+ (8.314* Temp*log(0.5585/16)))/(-8.314*Temp));
    keq_Cr2O3(i)=exp((-274347+(120.55*Temp)+ (8.314* Temp*log(0.5585/16)))/(-8.314*Temp));
 
    %Interface concentration calculation of Fe and Cr
    
    val(i)=keq_Cr2O3(i)*XbO2(i);
    XiFe(i)=XbFeO(i)/(keq_FeO(i)*XbO2(i));
    XiCr(i)=(XbCr2O3(i)^(0.5))/(val(i)^(1.5));
    CiFe(i)=XiFe(i)*MolarVM(i);
    CiCr(i)=XiCr(i)*MolarVM(i);
   
    %Mass Flux Balance
    JFe(i)=-MassTransCoeffFe*(CbFe(i)-CiFe(i));
    JCr(i)=-MassTransCoeffCr*(CbCr(i)-CiCr(i));
    
    %Moles Transfered
    MolesTransFe(i)=JFe(i)*dt*A;
    MolesTransCr(i)=JCr(i)*dt*A;
    MolesTransfFeO(i)=MolesTransFe(i);
    MolesTransfCr2O3(i)=0.5*MolesTransCr(i);
    MolesTransO2(i)=MolesTransFe(i)+1.5*MolesTransCr(i);
    
    %Moles of elements after this iteration
    MFe(i+1)=MFe(i)+MolesTransFe(i);
    MCr(i+1)=MCr(i)+MolesTransCr(i);
    MO2InMetal(i+1)=MO2InMetal(i)+MolesTransO2(i);
    Moxy(i+1)=Moxy(i)+MolesTransO2(i)+MO2Added;
    MCr2O3(i+1)=MCr2O3(i)-MolesTransfCr2O3(i);
    MFeO(i+1)=MFeO(i)-MolesTransfFeO(i);

    %Mass of each elements
    MassOfFe(i+1)=MFe(i+1)*AtomicMassofFe;
    MassOfCr(i+1)=MCr(i+1)*AtomicMassofCr;
    MassofO2(i+1)=Moxy(i+1)*AtomicMassofO2;
    MassOfFeO(i+1)=MFeO(i+1)*AtomicMassofFeO;
    MassOfCr2O3(i+1)=MCr2O3(i+1)*AtomicMassofCr2O3;

    %Total Moles in metal and Slag
    TotalMolesInMetal(i+1)=MCr(i+1)+MFe(i+1)+Moxy(i+1);
    TotalMolesInSlag(i+1)=MCr2O3(i+1)+MFeO(i+1); 
    
    %Total Mass of Metal and Slag
    
    TotalMassOfMetal(i+1)=MassOfFe(i+1)+MassofO2(i+1)+MassOfCr(i+1);
    TotalMassOfSlag(i+1)=MassOfCr2O3(i+1)+MassOfFeO(i+1);
     
    %Molar Volumes
    VM(i+1)=TotalMassOfMetal(i+1)/DensityOfMetal;
    VS(i+1)=TotalMassOfSlag(i+1)/DensityOfSlag;
    MolarVM(i+1)=TotalMolesInMetal(i+1)/VM(i);
    MolarVS(i+1)=TotalMolesInSlag(i+1)/VS(i);

    %Percentage of elements 
    
    PerFeInMetal(i+1)=(MFe(i+1)*100*AtomicMassofFe)/TotalMassOfMetal(i+1);
    PerCrInMetal(i+1)=(MCr(i+1)*100*AtomicMassofCr)/TotalMassOfMetal(i+1);
    PerO2InMetal(i+1)=(Moxy(i+1)*100*AtomicMassofO2)/TotalMassOfMetal(i+1);
    PerFeOInSlag(i+1)=(MFeO(i+1)*100*AtomicMassofFeO)/TotalMassOfSlag(i+1);
    PerCr2O3InSlag(i+1)=(MCr2O3(i+1)*100*AtomicMassofCr2O3)/TotalMassOfSlag(i+1);
        
    %temperature iterations 
    Temp=Temp+0.02;
    time(i+1)=time(i)+1;
end
figure(1)
plot(PerO2InMetal,'DisplayName','PerO2InMetal');hold on;plot(PerFeOInSlag,'DisplayName','PerFeOInSlag');plot(PerFeInMetal,'DisplayName','PerFeInMetal');plot(PerCrInMetal,'DisplayName','PerCrInMetal');plot(PerCr2O3InSlag,'DisplayName','PerCr2O3InSlag');hold off;
legend({'PerO2InMetal','PerFeOInMetal','PerFeInMetal','PerCrInMetal','PerCr2O3InMetal'});
figure(2)
plot(TotalMassOfMetal,'DisplayName','TotalMassOfMetal');hold on;plot(TotalMassOfSlag,'DisplayName','TotalMassOfSlag');hold off;
legend({'TotalMassofMetal','TotalMassofSlag'});

我在conda环境中安装了matlab内核，并打开了一个Jupyter笔记本来运行代码。在my env中运行conda列表会显示安装了

matlab内核

和

matplotlib

。但是当我尝试执行代码时，没有输出。重新启动内核也不会产生任何结果。事实上，Jupyter中的单元格在[*]中读取为


因此，我切换到查看是否可以使用以下链接执行代码：

我得到以下错误：
$octave-qf--无窗口系统演示.m输入金属总质量：警告：函数。/demo.m阴影核心库函数错误：输入：读取用户输入失败！错误：demo.m中执行异常

请提供有关如何执行此代码的任何指导。
这是否回答了您的问题？