[未答复] 怎么将循环解方程的结果都放到一个矩阵里面??

[复制链接]
east1111 发表于 2022-6-19 18:37:02
以下是我的程序,参数b(1)b(2)b(3)b(4)一共18组,最后可以得到18组的解,但是现在只能显示出最后一组的解,想问问大家怎么把18组的结果都显示出来呢?

syms x1 x2 x3 x4
filename1='m2.xls';
M=xlsread(filename1,'sheet2');
for i=1:18
b(1)=M(i,1);b(2)=M(i,2);b(3)=M(i,3);b(4)=M(i,10);b(5)=M(i,11);
eq1=(3.55e-10)-x1.*x2.*exp(2.*b(2).*b(5)+2.*b(1).*b(4)).*exp(-0.5547.*((b(4)+b(5)+4.*x4+x2+x3)./2).^0.5./(1+1.2.*((b(4)+b(5)+4.*x4+x2+x3)./2).^0.5))./(b(4).*exp(-0.5547.*((b(4)+b(5)+4.*x4+x2+x3)./2).^0.5./(1+1.2.*((b(4)+b(5)+4.*x4+x2+x3)./2).^0.5)+2.*b(3).*b(5)+2.*b(1).*x1));
eq2=(4.60e-07)-x2.*b(5).*exp(-0.5547.*((b(4)+b(5)+4.*x4+x2+x3)./2).^0.5./(1+1.2.*((b(4)+b(5)+4.*x4+x2+x3)./2).^0.5)).*exp(-0.5547.*((b(4)+b(5)+4.*x4+x2+x3)./2).^0.5./(1+1.2.*((b(4)+b(5)+4.*x4+x2+x3)./2).^0.5))./(0.004941889.*1);
eq3=(4.66e-11)-x2.*x4.*exp(-0.5547.*((b(4)+b(5)+4.*x4+x2+x3)./2).^0.5./(1+1.2.*((b(4)+b(5)+4.*x4+x2+x3)./2).^0.5)).*exp(-0.5547.*((b(4)+b(5)+4.*x4+x2+x3)./2).^0.5./(1+1.2.*((b(4)+b(5)+4.*x4+x2+x3)./2).^0.5))./(b(5).*exp(-0.5547.*((b(4)+b(5)+4.*x4+x2+x3)./2).^0.5./(1+1.2.*((b(4)+b(5)+4.*x4+x2+x3)./2).^0.5)));
eq4=(1.0577e-14)-x2.*x3.*exp(-0.5547.*((b(4)+b(5)+4.*x4+x2+x3)./2).^0.5./(1+1.2.*((b(4)+b(5)+4.*x4+x2+x3)./2).^0.5)).*exp(-0.5547.*((b(4)+b(5)+4.*x4+x2+x3)./2).^0.5./(1+1.2.*((b(4)+b(5)+4.*x4+x2+x3)./2).^0.5));
end
S=vpasolve(eq1,eq2,eq3,eq4,x1,x2,x3,x4);
x1=double(S.x1)
x2=double(S.x2)
x3=double(S.x3)
x4=double(S.x4)

您需要登录后才可以回帖 登录 | 注册

本版积分规则

相关帖子
相关文章
热门教程
站长推荐
快速回复 返回顶部 返回列表