二阶线性微分方程的幂级数解法

内容来源

常微分方程(第四版) (王高雄,周之铭,朱思铭,王寿松) 高等教育出版社


考虑二阶齐次线性微分方程

d 2 y d x 2 + p ( x ) d y d x + q ( x ) y = 0 \frac{\mathrm{d}^2y}{\mathrm{d}x^2}+ p(x)\frac{\mathrm{d}y}{\mathrm{d}x}+q(x)y=0 dx2d2y+p(x)dxdy+q(x)y=0

满足初值条件 y ( x 0 ) = y 0 , y ′ ( x 0 ) = y 0 ′ y(x_0)=y_0,y'(x_0)=y'_0 y(x0)=y0,y(x0)=y0 的情况

不失一般性,设 x 0 = 0 x_0=0 x0=0,否则设 t = x − x 0 t=x-x_0 t=xx0,经此变换后,方程形式不变。

定理

若上式中的系数 p ( x ) p(x) p(x) q ( x ) q(x) q(x) 都能展开成 x x x 的幂级数,且收敛区间为 ∣ x ∣ < R |x|<R x<R,则方程有形如

y = ∑ n = 0 ∞ a n x n y=\sum^\infty_{n=0}a_nx^n y=n=0anxn

的特解,也以 ∣ x ∣ < R |x|<R x<R 为收敛区间。

例1

y ′ ′ − x y = 0 y''-xy=0 y′′xy=0

设通解为

y = a 0 + a 1 x + ⋯ + a n x n + ⋯ y=a_0+a_1x+\cdots+a_nx^n+\cdots y=a0+a1x++anxn+

求导两次得

y ′ ′ = 2 ⋅ 1 a 2 + 3 ⋅ 2 a 3 x + ⋯ + n ( n − 1 ) a n x n − 2 + ⋯ y''=2\cdot1a_2+3\cdot2a_3x+\cdots+n(n-1)a_nx^{n-2}+\cdots y′′=21a2+32a3x++n(n1)anxn2+

y , y ′ ′ y,y'' y,y′′ 代入方程,比较同次幂的系数得

a 2 = 0 , n ( n − 1 ) a n − a n − 3 = 0 a_2=0,n(n-1)a_n-a_{n-3}=0 a2=0,n(n1)anan3=0

所以

a 3 k = a 0 ∏ i = 1 k 3 i ( 3 i − 1 ) a 3 k + 1 = a 1 ∏ i = 1 k 3 i ( 3 i + 1 ) a 3 k + 2 = 0 \begin{align*} &a_{3k}=\frac{a_0}{\prod^k_{i=1}3i(3i-1)}\\ &a_{3k+1}=\frac{a_1}{\prod^k_{i=1}3i(3i+1)}\\ &a_{3k+2}=0 \end{align*} a3k=i=1k3i(3i1)a0a3k+1=i=1k3i(3i+1)a1a3k+2=0

其中 a 0 , a 1 a_0,a_1 a0,a1 是任意的,因而

y = a 0 [ 1 + ∑ k = 1 ∞ x 3 k ∏ i = 1 k 3 i ( 3 i − 1 ) ] + a 1 [ x + ∑ k = 1 ∞ x 3 k + 1 ∏ i = 1 k 3 i ( 3 i + 1 ) ] y=a_0 \left[1+\sum^\infty_{k=1}\frac{x^{3k}}{\prod^k_{i=1}3i(3i-1)}\right] +a_1 \left[x+\sum^\infty_{k=1}\frac{x^{3k+1}}{\prod^k_{i=1}3i(3i+1)}\right] y=a0[1+k=1i=1k3i(3i1)x3k]+a1[x+k=1i=1k3i(3i+1)x3k+1]

例2

y ′ ′ − 2 x y ′ − 4 y = 0 y''-2xy'-4y=0 y′′2xy4y=0

初值条件 y ( 0 ) = 0 , y ′ ( 0 ) = 1 y(0)=0,y'(0)=1 y(0)=0,y(0)=1

还是设通解形式为

y = a 0 + a 1 x + ⋯ + a n x n + ⋯ y=a_0+a_1x+\cdots+a_nx^n+\cdots y=a0+a1x++anxn+

由初值条件可得

a 0 = 0 , a 1 = 1 a_0=0,a_1=1 a0=0,a1=1

所以

y = x + a 2 x 2 + ⋯ + a n x n + ⋯ y ′ = 1 + 2 a 2 x + ⋯ + n a n x n − 1 + ⋯ y ′ ′ = 2 a 2 + 3 ⋅ 2 a 3 x + ⋯ + n ( n − 1 ) a n x n − 2 + ⋯ \begin{align*} &y=x+a_2x^2+\cdots+a_nx^n+\cdots\\ &y'=1+2a_2x+\cdots+na_nx^{n-1}+\cdots\\ &y''=2a_2+3\cdot2a_3x+\cdots+n(n-1)a_nx^{n-2}+\cdots\\ \end{align*} y=x+a2x2++anxn+y=1+2a2x++nanxn1+y′′=2a2+32a3x++n(n1)anxn2+

代入方程,比较同次幂的系数得

a n = 2 n − 1 a n − 2 a_n=\frac{2}{n-1}a_{n-2} an=n12an2

所以

a 2 k + 1 = 1 k ! , a 2 k = 0 a_{2k+1}=\frac{1}{k!},a_{2k}=0 a2k+1=k!1,a2k=0

y = x + x 3 + x 5 2 ! + ⋯ + x 2 k + 1 k ! + ⋯ = x e x 2 y=x+x^3+\frac{x^5}{2!}+\cdots+\frac{x^{2k+1}}{k!}+\cdots=xe^{x^2} y=x+x3+2!x5++k!x2k+1+=xex2

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/pingmian/62862.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Java基础面向对象(String类)

String 特点 是内存中常量, 值在内存中一旦创建, 不可改 更改String类型引用的值本质上是将引用指向了一个新的字符串地址 String s1 "abc";String s2 s1;//引用s1的地址赋值给了s2 ​s2 "edf";//让s2指向新字符串 ​System.out.println("s1: &q…

【系统架构设计师】真题论文: 论软件质量保证及其应用(包括解题思路和素材)

更多内容请见: 备考系统架构设计师-专栏介绍和目录 文章目录 真题题目(2017年 试题4)解题思路论文素材参考软件质量保证定义和重要性软件质量保证在软件开发生命周期中的应用真题题目(2017年 试题4) 软件质量保证 (Software Quality Assurance. SQA) 是指为保证软件系统或…

LeetCode763. 划分字母区间(2024冬季每日一题 23)

给你一个字符串 s 。我们要把这个字符串划分为尽可能多的片段&#xff0c;同一字母最多出现在一个片段中。 注意&#xff0c;划分结果需要满足&#xff1a;将所有划分结果按顺序连接&#xff0c;得到的字符串仍然是 s 。 返回一个表示每个字符串片段的长度的列表。 示例 1&a…

【C++】深入优化计算题目分析与实现

博客主页&#xff1a; [小ᶻ☡꙳ᵃⁱᵍᶜ꙳] 本文专栏: C 文章目录 &#x1f4af;前言&#x1f4af;第一题&#xff1a;圆的计算我的代码实现代码分析改进建议改进代码 老师的代码实现代码分析可以改进的地方改进代码 &#x1f4af;第二题&#xff1a;对齐输出我的代码实现…

6.824/6.5840 Lab 3: Raft——Part 3B3C

芙蓉花又栖满了枝头 奈何蝶难留 漂泊如江水向东流 望断门前隔岸的杨柳 寂寞仍不休 我无言让眼泪长流 ——山外小楼夜听雨 完整代码见&#xff1a; https://github.com/SnowLegend-star/6.824 在完成Lab之前&#xff0c;务必把论文多读几遍&#xff0c;力求完全理解Leader选举、…

LeetCode - #150 逆波兰表达式求值

文章目录 前言1. 描述2. 示例3. 答案关于我们 前言 我们社区陆续会将顾毅&#xff08;Netflix 增长黑客&#xff0c;《iOS 面试之道》作者&#xff0c;ACE 职业健身教练。&#xff09;的 Swift 算法题题解整理为文字版以方便大家学习与阅读。 LeetCode 算法到目前我们已经更新…

使用go实现一个简单的rpc

什么是rpc, rpc是干什么的?几种协议的压测数据对比:tcphttp 使用tcp实现一个简单的rpc服务 什么是rpc, rpc是干什么的? rpc的作用就是实现远程的服务调用 工作流程: 客户端携带服务信息(服务名,方法名)数据 去请求服务端,服务端拿到数据,解析后执行对应的方法,将结果返回给客…

【C++】continue语句、goto语句

1、continue 语句 作用&#xff1a;在循环语句中&#xff0c;跳过本次循环中余下尚未执行的语句。继续下一次循环。 注意&#xff1a;continue只能用于循环中。 示例&#xff1a; 代码&#xff1a; //continue的用法 #include<iostream> using namespace std; int ma…

最长最短单词

最长最短单词 C语言实现C实现Java实现Python实现 &#x1f490;The Begin&#x1f490;点点关注&#xff0c;收藏不迷路&#x1f490; 输入1行句子&#xff08;不多于200个单词&#xff0c;每个单词长度不超过100&#xff09;&#xff0c;只包含字母、空格和逗号。单词由至少一…

ESP32项目 --- 智能门锁(WiFi 蓝牙 OTA)

1 项目简介 1.1 项目概述 本项目是实现一款智能门锁中的智能控制部分, 可以应用在家庭, 办公室等任何使用门锁的场所. 本项目实现了以下主要功能: &#xff08;1&#xff09;通过按键配置密码 &#xff08;2&#xff09;通过按键输入密码开锁 &#xff08;3&#xff09;录…

【Qt】QTableView选中行发生变化时触发的信号

问题 QTableView选中的行发生变化时&#xff0c;使用的信号是QTableView的selectionModel()里的currentChanged信号&#xff0c;界面点击行来回切换&#xff0c;发现怎么也触发不了&#xff1f; 原因 信号槽连接放在了QTableView数据初始化前面&#xff0c;这时候QTableView…

洛谷题单-入门2-分支结构-python-下

找出出现的最早的最大值 count 0 list_number [] while True:list_number.append(list(map(int, input().split())))count 1if count 7:breaklist2_number_total []for i1,i2 in list_number:list2_number_total.append(i1i2)target max(list2_number_total)index 0 if…

Java Collections 深度探索

在 Java 编程中&#xff0c;java.util.Collections是一个非常重要的工具类&#xff0c;它提供了一系列对集合进行操作的静态方法。本文将深入探讨 Java Collections 的功能、用法、优势以及在实际编程中的应用。 一、引言 Java Collections 框架为开发者提供了一套强大而灵活…

Taro小程序开发随记

处理taro小程序显示wangeditor内的a标签跳转 Taro.options.html.transformElement (el) > {if (el.props.class h5-a) {el.__handlers.tap [() > toWebView(el.props.href)]}return el } 处理wangeditor富文本内容中图片视频到小程序中展示问题 <view class&qu…

Pytorch使用手册-What is torch.nn really?(专题九)

我们建议将本教程作为 notebook 而不是脚本运行。要下载 notebook(.ipynb)文件,请点击页面顶部的链接。 PyTorch 提供了精心设计的模块和类,如 torch.nn、torch.optim、Dataset 和 DataLoader,帮助你创建和训练神经网络。为了充分利用这些工具的强大功能并根据你的问题进…

框架模块说明 #05 权限管理_03

背景 权限设计可以分为两个主要方面&#xff1a;操作权限和数据权限。前两篇文章已经详细介绍了操作权限的设计与实现&#xff0c;以及如何将其与菜单关联起来的具体方法。本篇将聚焦于数据权限&#xff0c;为您深入讲解相关的设计与实现方式。 全局开关 Value("${syst…

Linux网络编程之---多线程实现并发服务器

下面我们来使用tcp集合多线程实现并发服务器 一.服务端 #include <stdio.h> #include <arpa/inet.h> #include <unistd.h> #include <stdlib.h> #include <string.h> #include <pthread.h>typedef struct sockinfo {char ip[16];unsigne…

二分查找!

问题描述 小明在图书馆借阅书籍&#xff0c;图书馆的书籍在系统中按序号顺次排列&#xff0c;小明在借阅后&#xff0c;需在系统中从“在馆书籍列表”中将该书删除。请帮助小明编写一个函数&#xff0c;在现有列表{1, 3, 5, 6, 7, 10, 12, 14, 26, 32, 35, 39, 42, 45, 54, 56…

按vue组件实例类型实现非侵入式国际化多语言翻译

#vue3##国际化##本地化##international# web界面国际化&#xff0c;I18N&#xff08;Internationalization&#xff0c;国际化&#xff09;&#xff0c;I11L(International&#xff0c;英特纳雄耐尔)&#xff0c;L10N&#xff08;Localization&#xff0c;本地化&#xff09;&…

Linux C/C++编程之静态库

【图书推荐】《Linux C与C一线开发实践&#xff08;第2版&#xff09;》_linux c与c一线开发实践pdf-CSDN博客《Linux C与C一线开发实践&#xff08;第2版&#xff09;&#xff08;Linux技术丛书&#xff09;》(朱文伟&#xff0c;李建英)【摘要 书评 试读】- 京东图书 (jd.com…