代码与象征性的初始条件求解一阶的颂歌
7视图(30天)
显示旧的评论
我很新的MATLAB,编写MATLAB代码解决这个问题# 20 1.1节的Kreyszig (2011)。
问题# 20如下:
Exponentialdecay。Subsonicflight。
亚音速飞机的引擎的效率取决于空气压力和通常是最大的附近35000英尺。找到空气压力
在这个高度。
物理信息。
变化的速度
压力成正比。在18000年英国《金融时报》是其价值的一半
在海平面。
![](http://www.tianjin-qmedu.com/ch/matlabcentral/answers/uploaded_files/1207263/image.png)
![](http://www.tianjin-qmedu.com/ch/matlabcentral/answers/uploaded_files/1207268/image.png)
![](http://www.tianjin-qmedu.com/ch/matlabcentral/answers/uploaded_files/1207273/image.png)
我写了下面的代码来解决这个问题,但是得到一个错误的答案。
清晰的
信谊y (t) k
ode = diff (y) = = k * y
气孔导度= y (18000) = = y (0) / 2
ySol (t) = dsolve(颂歌,气孔导度)
我有一个简单的解决方案。
如何让一个非平凡的解决这个问题吗?
如果我解决问题用铅笔和橡皮擦,简单的解决方案
。
![](http://www.tianjin-qmedu.com/ch/matlabcentral/answers/uploaded_files/1207288/image.png)
0评论
接受的答案
Torsten
2022年11月25日
编辑:Torsten
2022年11月25日
信谊y (h) k
ode = diff (y) = = k * y;
sol_ode_general = dsolve (ode)
var = symvar (sol_ode_general)
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])
Pressure_at_36000_feet =潜艇(sol_ode_specific一样,var (2), 36000 * 0.3048)在m %压力Pa,距离
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()和这一步后通过一个初始条件解决()为每一个适当的变量。