MATLAB的答案

通过建设难以建模的温度变化

1视图(30天)
洛根帕里什
洛根帕里什 2022年3月6日
回答: 大卫。马塞罗说道 2022年3月6日
我试图通过一个简单的两层小屋模型温度变化但我保持与preallocating遇到问题我的温度数组或循环模型和欧拉方法的温度。模型我试图做的小屋,有一个壁炉,输出一定量的能量在一楼。有三个通量方程:通量从1楼到外面,从1楼到2楼,从二楼到外面。这些通量方程使用牛顿冷却定律创建两个一阶常微分方程1和2层的温度。
通量方程是:
温度常微分方程是: 和
在我的代码我包括所有常数所示的方程。
%确定常数
k1 = 0.2;% W / (m ^ 2 * K)
k2 = 3;%的范围从1到10 W / (m ^ 2 * K)
k3 = 0.5;% W / (m ^ 2 * K)
秦= 2500;% 2500 Max, W / (m ^ 2 * k)
是= 0.5;% m ^ 2
A1 = 30;% 1楼的面积,m ^ 2
A2 = 30;% 2楼面积,m ^ 2
rho_air = 1.2;% ^ 3公斤/米
以下简称= 1005;% 1005 J /(公斤* K)
V1 = 90;%的低地板上
V2 = 38.97;%的楼上
h = 1;% stepsize
t = linspace (0, 2 * 86400);
N = t / h;
%初始化数组
T1 = 0 (1, N);
T2 = 0 (1, N);
%集成电路
T1 (1) = 5;% C
T2 (1) = 7;% C
T_out = -10 * sin((2 *π。* t) / 86400);% C
%的功能
Q1 = k1 * (T1 - T_out);
Q2 = k2 * (T1, T2 + 5);
第三季度= k3 * (T2 - T_out);
dt₁= (Ain *秦)/(以下简称* rho_air * V1) - (A1 * Q1) /(以下简称* rho_air * V1) + (A1 * Q2) /(以下简称* rho_air * V1);
dT2 = (Ain *秦)/(以下简称* rho_air * V1) + (A1 * Q2) /(以下简称* rho_air * V1) + (A2 * Q3) /(以下简称* rho_air * V2);
为我= 1:n - 1
t (i + 1) = t(我)+ h;
T1 (i + 1) = T1 (i) + h * dt₁(T1 (i), t(我));
T2 (i + 1) = T2 (i) + h * dT2 (T2 (i), t(我));
结束
情节(t, T1,的r -)
情节(t, T2,“g -”)
任何帮助将不胜感激得到这件事。

在回答这个问题。

答案(1)

大卫。马塞罗说道
大卫。马塞罗说道 2022年3月6日
下面的代码是
%确定常数
k1 = 0.2;% W / (m ^ 2 * K)
k2 = 3;%的范围从1到10 W / (m ^ 2 * K)
k3 = 0.5;% W / (m ^ 2 * K)
秦= 2500;% 2500 Max, W / (m ^ 2 * k)
是= 0.5;% m ^ 2
A1 = 30;% 1楼的面积,m ^ 2
A2 = 30;% 2楼面积,m ^ 2
rho_air = 1.2;% ^ 3公斤/米
以下简称= 1005;% 1005 J /(公斤* K)
V1 = 90;%的低地板上
V2 = 38.97;%的楼上
h = 1;% stepsize
t = 0: h: 2 * 86400;
N =长度(t);
%初始化数组
T1 = 0 (1, N);
T2 = 0 (1, N);
%集成电路
T1 (1) = 5;% C
T2 (1) = 7;% C
T_out = -10 * sin((2 *π。* t) / 86400);% C
为我= 1:n - 1
Q1 = k1 * (T1 (i) - T_out(我));
Q2 = k2 * (T1 (i) - T2(我)+ 5);
第三季度= k3 * (T2 (i) - T_out(我));
dt₁= (Ain *秦)/(以下简称* rho_air * V1) - (A1 * Q1) /(以下简称* rho_air * V1) - (A1 * Q2) /(以下简称* rho_air * V1);
dT2 = (Ain *秦)/(以下简称* rho_air * V1) + (A1 * Q2) /(以下简称* rho_air * V1) + (A2 * Q3) /(以下简称* rho_air * V2);
T1 (i + 1) = T1(我)+ h * dt₁;
T2 (i + 1) = T2(我)+ h * dT2;
结束
图(1)
情节(t, T1,的r -t, T2,“g -”)
传奇(“T1”,《终结者2》,“位置”,“最佳”)

类别

找到更多的在不连续在帮助中心和文件交换

标签

s manbetx 845


释放

R2021b

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!