复数域内的微分解法

最新推荐文章于 2024-06-27 09:07:08 发布

chenbingchenbing

最新推荐文章于 2024-06-27 09:07:08 发布

阅读量3.6k

点赞数

分类专栏：综合文章标签： c float

本文链接：https://blog.csdn.net/chenbingchenbing/article/details/6599765

版权

本文探讨了在复数域内解决微分方程的简便方法，以y'=e^x*cos x为例，通过转换到复数域简化求解过程。通过数值计算验证了解的准确性，并用Lisp程序实现计算过程，尽管存在浮点数计算精度限制导致的误差，但整体解法已非常接近实际结果。

摘要由CSDN通过智能技术生成

复数域内的微分解法
先考虑一个简单的微分形式y'=e^x*cos x ,
这种形式如果考虑分离变量法，需要采用两次局部积分法最后才能获得结果；但如果在复数域

内解答将会很简单，并且我们将验证采用数值解法的结果比较；
将上面形式变为复数域的形式有：
yy'=Re { e^x*e^ix}
Go
yy'=Re{e^{x(1+i)}}
GO
yy=Re {1/(1+i)*e^{x(1+i)}+C } (C为实数)
Go
yy=Re { (1-i)/2*e^{x(1+i)}+C }
Go
yy=Re { 1/2*e^x*(cos x+i*sin x) - 1/2*e^x*(i*cos x-*sin x) +C}
Go
yy=Re { 1/2*e^x*(cos x-sin x) + i/2*e^x*(sin x -cos x) +C}
回到实数领域内有：
y= 1/2*e^x*(cos x+sin x) +C
将y代入原微分方程可以获得解答；
现在假设初始条件y(0)=1,这个时候C为1/2

下面写程序来证明：
(setq c 0.5)

(defun pow (num count)
(if (or (> count 1) (eq count 1) )
   (* num
   (pow num