Using for loop to change different input while solving silmultaneous equations

2020 2 27
0 个回答
5 次查看（30 天）

0 个投票

I want to solve the simultaneous equation with one different input value. Help me to solve this!
Below is the command that I used:
clear
clc
%Energy balance in hot stream
%IN-OUT-ACCUMULATION=0
%m_HW_IN*c*T_HW_IN - m_HW_OUT*c*T_HW_OUT - h*A*(T_HW_OUT-T_CW_OUT)=0
%m_HW_IN*c*(T_HW_IN-T_HW_OUT) - h*A*(T_HW_OUT-T_CW_OUT)=0
%Energy balance in cold stream
%IN-OUT+ACCUMULATION=0
%m_CW_IN*c*T_CW_IN - m_CW_OUT*c*T_CW_OUT + h*A*(T_CW_OUT-T_CW_OUT)=0
%m_CW_IN*c*(T_CW_IN-T_CW_OUT) + h*A*(T_HW_OUT-T_CW_OUT)=0
%Constant
m_HW_IN=132.93;%unit=g/s
c=4.186;%unit=kJ/kg*°C
T_HW_IN=60;%unit=°C
T_CW_IN=33;%unit=°C
h=9000;%unit=W/m^2*°C
A=16*2*pi*0.0079375*4.06;%unit=m^2
syms T_CW_OUT T_HW_OUT
T_HW_OUT=zeros(5,1);
T_CW_OUT=zeros(5,1);
sol=zeros(5,1);
for m_CW_IN=[83.08;99.70;116.32;132.93;149.55];
eqn1=m_CW_IN.*c.*(T_CW_IN-T_CW_OUT)+h.*A.*(T_HW_OUT-T_CW_OUT)==0;
eqn2=m_HW_IN.*c.*(T_HW_IN-T_HW_OUT)-h.*A.*(T_HW_OUT-T_CW_OUT)==0;
sol = solve([eqn1,eqn2],[T_CW_OUT,T_HW_OUT]);
T_HW_OUT= sol.T_HW_OUT;
T_CW_OUT= sol.T_CW_OUT;
T_HW_OUT=vpa(T_HW_OUT,4)
T_CW_OUT=vpa(T_CW_OUT,4)
end

请先登录，再回答此问题。

回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Symbolic Math Toolbox 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by