欧拉迭代算法#
对于标准微分方程,初值确定时,它一定有且仅有唯一解,可以计算其数值近似解
y~1=y0+y0′(x1−x0),以此类推
我们已经说明过特解,但是它有一些其他的要求。特解必须是连续的,譬如说,特解不可以是 ln∣x∣
空气净化模型#
现在有体积为 V 的空气在容器中,在0时刻,有 CCO2=ρ0 ,现在不断加入新鲜空气,注入速度是每秒 V0 体积,其 CCO2=ρ1 ,并且假定空气不被压缩,注入多少就排出多少,空气的混合会迅速完成。现在求容器中 CO2 浓度变化规律
对于从 t→t+dt 的极短时间中,注入的二氧化碳质量为 dtV0ρ1 ,排出的二氧化碳质量为 dtV0ρˉ(t,dt)
dt→0,二氧化碳变化量即为dtV0ρ1−dtV0ρ(t)=Vdρ,dρ为浓度的变化量所以浓度变化率dtdρ=VV0(ρ1−ρ(t))
可分离变量的微分方程#
微分方程大多是不可解的,可解的具有一定的特性。我们在这里讨论一阶微分方程的解法
y′=f(x,y)可分离变量的微分方程可以转化为如下形式
y′=P(x)Q(y)于是有
dxdy=P(x)Q(y)观察不难发现,若Q(y)=0有解,y1...yn 微分方程就有许多水平特解 y=y1...
若Q(y)没有零点,那么我们进行如下变换
Q(y)dy=P(x)dx此时对于两端取不定积分,就可以得到解
∫Q(y)dy=∫P(x)dxF(y)=G(x)+C之后算出常数C,最后得出 y=ϕ(x) 即可
有时也可以取定积分,对于初值条件 (x0,y0) ,若 Q(y0)=0 可以计算定积分
∫x0xQ(y)dy=∫y0yP(x)dx实际上,不存在Q(y)有零点也有非零点的情况,这会使得解不符合解的唯一性
对于通解,我们需要涵盖所有情况,因此可以出现绝对值等不连续的形式,而对于特解才有需要连续的要求
实际上,通解常常不能涵盖所有情况,因此也存在一些例外的特解,通常来自我们在求解过程中增加的一些限定,比如使某个变量不能等于0
例题:
xdy=ydxx=0,y=0,ydyln∣y∣取指数,∣y∣也就是y因此yx=0,y=0也满足这一条件同理我们也知道x=C1y是微分方程的解,=xdx=ln∣x∣+C=eC∣x∣=±eCx=Cx(C=0),也符合y=Cx因此微分方程的解是任意过原点的直线有时候,我们并不能得到可以很好地分离的微分方程,如果结构合适,我们也可以使用换元法
例题:
求dxdy=sin2(x−y+1)的通解令u则有y故dxdy=1−u′dxdu=u′cos2udutanu=x−y+1=x+1−u=sin2u=cos2u=dx=x+C