matlab2019中solve求解析式,solve求解方程组，方程中含有赋过值的字母

这是我的程序

load('CO.mat')

load('CO2.mat')

load('CH4.mat')

load('H2.mat')

load('N2.mat')

load('O2.mat')

load('H2O.mat')

l=[0.5 0.6 0.7 0.8 0.9];

for i=1:5

disp('>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>');

n=l(i)

L_O2=(0.5*CO+0.5*H2+2*CH4-O2)*0.01

K1=0.333

K2=550

syms b1 b2 b3 b4 a b c d e f;

b1=(CO+CO2+CH4)*0.01

b2=(H2+2*CH4+H2O)*0.01

b3=(0.5*CO+CO2+O2+0.5*H2O)*0.01+n*L_O2

b4=0.01*N2+3.76*n*L_O2

a='VCO2+VCO+VCH4=b1';

b='VH2+VH2O+2*VCH4=str2num(b2)';

c='VCO2+0.5*VCO+0.5*VCH4=str2num(b3)';

d='VN2=str2num(b4)';

e='(VCO2*VH2)/(VCO*VH2O)=0.333';

f='(VH2^2)/((VCO2+VH2+VCO+VH2O+VCH4+VN2)*VCH4)=550';

[z1,z2,z3,z4,z5,z6]=solve(a,b,c,d,e,f,'VCH4','VCO','VCO2','VH2','VH2O','VN2');

vpa([z1,z2,z3,z4,z5,z6],3)

end

2017-11-16 07:40 上传

点击文件名下载附件

175 Bytes, 下载次数: 0

2017-11-16 07:40 上传

点击文件名下载附件

177 Bytes, 下载次数: 0

2017-11-16 07:40 上传

点击文件名下载附件

174 Bytes, 下载次数: 0

2017-11-16 07:40 上传

点击文件名下载附件

177 Bytes, 下载次数: 0

2017-11-16 07:40 上传

点击文件名下载附件

173 Bytes, 下载次数: 0

2017-11-16 07:40 上传

点击文件名下载附件

178 Bytes, 下载次数: 0

2017-11-16 07:40 上传

点击文件名下载附件

177 Bytes, 下载次数: 0