使用 sympy 求解差分方程

2744557306

实现了解递归方程并绘制其图形的过程。让我解释一下：
2 K* N" ^4 L' {
2 G1 h8 ?/ f4 p5 Y4 R5 {/ K1.定义方程：
   x = sp.symbols('x')
$ [6 o1 j; x. n' u  g1 w* h   y = sp.Function('y')
   f = y(x+2) - y(x+1) - y(x)
   con = {
       y(1):1,
       y(2):1,
' U# m- P' S8 X5 t, C   }
2.使用 SymPy 定义了一个递归方程 f，该方程描述了函数 y 在不同点之间的关系。
3.设置了两个初始条件 y(1)=1 和 y(2)=1。
4.解递归方程：
9 y0 e, [1 x8 e1 J8 o) P
   solve = sp.rsolve(f, y(x), con)
5.使用 rsolve 函数解析地求解了递归方程，得到了其解析解。
6.画图：
   x1 = np.linspace(1,10,10)
, ?$ a# c6 o. d) o, d, o- L# m   y1 = []
   for each in x1:
" [# f5 S7 e$ l  L* s: q       y1.append((solve.subs(x,each).evalf()))
; e* n. w: E& A) Q8 a2 G1 O, D
   import matplotlib.pyplot as plt
5 \9 b- s' k) I$ o) u. B
   plt.plot(x1,y1)
( U/ U0 n0 \, `5 U$ ]2 s   plt.scatter(x1,y1)


" i* _( {0 o7 t4 r7 i7.生成了一组横坐标 x1，并通过代入解析解中的解来计算相应的纵坐标 y1。
) Z3 x- r+ {1 z' w+ E8.使用 Matplotlib 绘制了解析解的图形，并用散点表示离散点。
9 G; [/ J. e0 {. t1 l6 w
: Q' s$ n; i' I; X& h# L这样，代码就完成了对递归方程的解析求解，并将结果可视化的过程。