如果温度增加17到23日号,问= 15 * 10 ^ 6,如果温度23减少untill17关掉,Q = 0时我怎么能画出温度与时间T(0) = 20度

1视图(30天)
dT / dT = (Q / 15172.5) + 29.66 + 2.966 t
4评论
山姆翟
山姆翟 2022年3月23日
你的英语没有问题。我能明白你写的什么。只是我不真正的得到你想要的,因为策划结果需要解决温度控制微分方程,和建模中的磁滞行为恒温器(Q)。然而,你没有提到任何问题ODE。这就是为什么我要求澄清。
幸运的是你了。然而,我发现了两件事:
  1. 恒温器以来表现出非线性行为(滞后),拉普拉斯和拉普拉斯逆变换不能用于非线性系统。
  2. 您的建模的恒温器(Q)似乎是不正确的,并没有行动之间的定义。请确认是否这是目的。

登录置评。

答案(1)

大卫。马塞罗说道
大卫。马塞罗说道 2022年3月23日
编辑:大卫。马塞罗说道 2022年3月23日
你可以试试这个
清楚,clc
tf = 3;%的最后一次集成
T0 = 20;%初始温度
往往= 0;%初始时间的集成
idx = 0;%迭代索引
往往<特遣部队
如果T0 < 23%制定积极的热量输入当T < 23
Tx = 23;
Q = 1.5 e + 6;
其他的%设置负热量输入当T > = 23
Tx = 17;
Q = -1.5 e + 6;
结束
idx = idx + 1;
tspan =(往往,tf);
选择= odeset (“事件”@ (t, t) heatSwitch (t, t, Tx));
索尔(idx) = ode15s (@ (t, t) tempControl (t, t, Q), tspan, T0,选择);
往往=索尔(idx)方式(结束);
T0 =索尔(idx) .y(结束);
结束
情节([索尔(:)。x]、[索尔(:).y],“k”,
(索尔(:)。x),(1,长度((sol (:)。x))) * 23日“r——”,
(索尔(:)。x),(1,长度((sol (:)。x))) * 17,“r——”)
包含(,年代)
ylabel (“温度、C”)
轴([负无穷到正无穷16 24])
函数dTdt = tempControl (t, t, Q)
dTdt = (Q / 15172.5) + 29.66 + 2.966 * T;
结束
函数(位置、isterminal方向)= heatSwitch (t, t, Tx)
位置= T-Tx;%的价值,我们希望为零
isterminal = 1;%停止集成
方向= 0;% 0可以从两个方向靠近
结束
然而,请注意,当T = 23日我开关Q - Q否则温度不会向下回到17岁。事实上,如果Q = 0, dT / dT在你espression仍然是积极的。
8的评论

登录置评。

类别

找到更多的在编程帮助中心文件交换

标签

社区寻宝

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

开始狩猎!