solve求解方程组后结果表示为root。

Question

trmpexs 2023-5-19

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1968014-solve-root

回答： buxuber 2023-5-19

%全混流简化模型
clc
close all 
syms PA CB
global QG0 QL0 kox CA CB0 PA0 PN20 R_const T_react kla ping HA FN2 H0 S0 PM
% 基本参
%设备操作参数
                H0=12;                    %塔高 m
                S0=13.8;                    %塔截面积
                %PM=3.0;                    %进塔压力 atm 
                PM=1.5*10^5;                 %出塔压力 pa
                QL0=200/3600 ;                % 液体进料 m3/s
                QG0=8000/3600;             % 气体进料 m3/s
 %传质反应参数
                CB0=250;                  %氢化完成液氢蒽醌集合物浓度，mol/m3
                PA0=63600;                %进塔气体氧分压，pa
                PN20=234000;              %进塔气体氮分压，pa
                R_const=8.314;            % J mol-1 K-1
                T_react=325.15;          % 温度 K
                HA=10901;               % A亨利系数  [  Pa / (mol m^-3) ]  
                kox=0.002;              %  反应速率常数
                kla=100;                %  液膜体积传质系数
                ping=0.8;               % 气含率
                CA=((kla./(1-ping))*(PA./HA))./((kla./(1-ping))+kox*CB);
                FN2=(0.78*QG0*PN20)./(R_const*T_react);
                FAout=FN2*PA./(PM-PA);
                FAin=(0.21*QG0*PA0)./(R_const*T_react);
   %方程及求解
   %eq1=(0.21*QG0*PA0)./(R_const*T_react)-FN2*PA./(PMout-PA)-kox*CA*CB;
   eq1=FAin-FAout-kox*CA*CB;
   eq2=QL0*(CB0-CB)-kox*CA*CB*S0*H0;
   [solPA,solCB]=solve(eq1,eq2);