难度 3/5
检查一维衰减方程的物理残差
方程 u'(x)+u(x)=0 满足 u(0)=1。取候选族 uₐ(x)=exp(-ax),推导内部残差和边界残差,并求使两者同时为零的 a。说明自动微分在这一步计算什么。
查看提示
先对 exp(-ax) 关于 x 求导;边界残差直接代入 x=0,内部残差保留为 x 的函数。
展开分步解答
uₐ'(x)=-a exp(-ax),所以内部残差 r(x)=(1-a)exp(-ax)。边界残差 b=uₐ(0)-1=0。由于指数函数不为零,r 对所有 x 为零当且仅当 a=1。自动微分计算的是候选函数对输入 x 的导数,并把该导数代入方程残差。
结果核验:令 a=1 后,r(x)=0 且 b=0;再令 a=0 可得 r(x)=1,说明残差公式能识别不满足微分方程的常数候选。