代码与象征性的初始条件求解一阶的颂歌

7视图(30天)
我很新的MATLAB,编写MATLAB代码解决这个问题# 20 1.1节的Kreyszig (2011)。
问题# 20如下:
Exponentialdecay。Subsonicflight。 亚音速飞机的引擎的效率取决于空气压力和通常是最大的附近35000英尺。找到空气压力 在这个高度。 物理信息。 变化的速度 压力成正比。在18000年英国《金融时报》是其价值的一半 在海平面。
我写了下面的代码来解决这个问题,但是得到一个错误的答案。
清晰的
信谊y (t) k
ode = diff (y) = = k * y
ode (t) =
气孔导度= y (18000) = = y (0) / 2
气孔导度=
ySol (t) = dsolve(颂歌,气孔导度)
ySol (t) =
0
我有一个简单的解决方案。
如何让一个非平凡的解决这个问题吗?
如果我解决问题用铅笔和橡皮擦,简单的解决方案

接受的答案

Torsten
Torsten 2022年11月25日
编辑:Torsten 2022年11月25日
信谊y (h) k
ode = diff (y) = = k * y;
sol_ode_general = dsolve (ode)
sol_ode_general =
var = symvar (sol_ode_general)
var =
eqn1 =潜艇(sol_ode_general一样,var (2), 0) = = 1.013 e5;在m %压力Pa,距离
eqn2 =潜艇(sol_ode_general一样,var (2), 18000 * 0.3048) = = 1.013 e5/2;在m %压力Pa,距离
(溶胶。c₁sol.k] =解决([eqn1 eqn2], [var (1), var (3)));
sol_ode_specific =潜艇(sol_ode_general, var (1) var (3)], [sol.C_1 sol.k])
sol_ode_specific =
Pressure_at_36000_feet =潜艇(sol_ode_specific一样,var (2), 36000 * 0.3048)在m %压力Pa,距离
ans =
25325年
3评论
沃尔特·罗伯森
沃尔特·罗伯森 2022年11月25日
dsolve()仅仅是为了解决一个初始条件每级的导数。它构造的一般公式解决方案(t) = = f (t) + C,然后花时间的第一个初始条件,如y的0 (0),subsitutes,在得到解决方案(0)= = f (0) + C,替代品的已知值初始情况下,得到y0 = f (0) + C和C,解决和back-substitutes得到解决(t) = f (t) + specific_constant。当第二个初始条件y(18000),在溶液中没有剩余常数求解(t) = f (t) + specific_constant——在同一水平,至少不是一个常数,如果y0 发生了 涉及多个未知变量,第二个初始条件不知道哪一个是旨在解决。
dsolve()不是为了专门寻找周期解。万博 尤文图斯如果你有未知数的方程,其价值是由多个初始条件,然后求出初始条件使用dsolve()和这一步后通过一个初始条件解决()为每一个适当的变量。

登录置评。

更多的答案(0)

社区寻宝

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

开始狩猎!