技术标签: 二阶有限差分三对角matlab
摘要:本课题主要研究运用有限差分法来数值求解线性及非线性的二阶常微分方程的边值问题,以差商公式替换微分方程中的各阶导数,得到差分方程。经过分析差分方程的本质,联系解三对角线性方程组的Crout分解法和解非线性方程组的Newton法,求解差分方程,从而给出解二阶常微分方程边值问题的有限差分算法。最后,我们利用MATLAB数学应用软件上机编写程序,并用他们求解一些具体的二阶常微分方程的边值问题。我们将实验的数值结果同理论的精确解进行对比和进行误差分析,由此说明了这些方法的有效可行。同时,我们对同一道题目选择不同的实验步长,比较在不同情况下的误差大小,进一步明确本课题所研究运用的算法的优缺点。
关键词 常微分方程;边值问题;有限差分;数值解
目录
摘要
Abstract
1 绪论-1
1.1 研究背景及意义-1
1.2 研究的现状-1
1.3 本文的主要工作-2
2 有限差分法的原理-3
2.1 差分概念-3
2.2 构造差分的方法-4
2.3 建立差分格式的步骤-4
3 求解线性及非线性方程组的核心算法-5
3.1 解三对角线性方程组的Crout分解法-5
3.2 解非线性方程组的Newton方法-6
4 解两类二阶常微分方程的边值问题的方法-9
4.1 线性问题的有限差分法-9
4.2 非线性问题的有限差分法-11
5 数值实验及应用-15
5.1 线性问题的数值求解及应用-15
5.1.1 验证型数值实验-15
5.1.2 应用型实例解析-17
5.2 非线性问题的数值实验及应用-19
5.2.1 验证型数值实验-19
5.2.2 应用型实例解析-22
结论-25
致谢-26
参考文献-27
附录-28
最近实施的一个项目,除开linux操作系统不用自己安装,所有操作都需要自己亲力亲为,对于linux命令不是很熟练的我,通过各种查资料终于把最开始的磁盘挂载和扩充搞定,留下笔记,供大家参考。 逻辑磁盘挂载方法一(PVS分盘) [[email protected] ~]# pvs PV VG Fmt...
logback.xml配置文件修改<springProperty scope="context" name="log.consolePattern" source="hcf.log.console-pattern" defaultValue="%red(%d{yyyy-MM-dd HH:mm:ss.SSS})|%highlight(%level)|%green(%thread)|%X{clientIp}|%X{traceId}|%X{rpcId}|%boldMagenta(%c.%M[%L])|
DP–m处理器问题- m processors(FZU - 1442)DescriptionThere are n data packets and m processors in a network communication system. In this problem, you need to distribute the n data packets to the m processors so that the workload of each processor is as balanc
头条号的微头条你们都发些什么,该如何写?心灵相约头条号。呼唤着,努力着,拼搏着,它属于努力地人们,是传播正能量的家园。今天起要根扎头条号,写出优秀地原创作品。头条号的朋友们努力吧!用心血和汗水浇灌头条号之花,让它生根开花结满丰收的果实。三农领域的微头条怎么写?有大神能够指导下吗?微头条,尤其是三农领域的微头条近来受到越来越多三农创作者的喜爱。但是在实际的操作中,很多三农创作者对于微头条的认知还停留...
JEECG 3.7.7 闪电版本发布,提供5套主流UI代码生成器模板!此版本为性能和表单UI深化加强版本,简称闪电版本 (闪电般的速度,主流的Bootstrap表单风格)。平台性能访问速度提升至少3倍,表单提供Bootstrap风格,新增加5套Bootstrap风格代码生成器模板,极大提升美感和代码可控性,降低了表单的开发成本,新鲜出炉,哇咔咔!!
51CTO下载中心为大家收集了一个批IOS免豆资料,做了一个下载地址目录,方便大家选择性下载,希望大家喜欢。完整目录下载地址:http://down.51cto.com/data/452500资料列表:iPhone 开发笔记实录(代码讲解)http://down.51cto.com/data/164304Phone Clube 第六期:高效手机程序开发经验分
上几次涂鸦写了几篇博文,没想到引起大家很大的反响。在这里,首先谢谢各位的厚爱。也接到许多来自学生朋友和刚出校门1-2年的朋友的邮件和QQ交流。有人说,我写的太意识流,看了很多当年流行的技术,都被看晕了,不知道如何去选择自己适合的技术。有人说,说的有些随意,不够严谨,有些地方有偏见,没有真正去查数据论证。确实是,因为只是自己偶尔心得体会,就在网页上直接敲的文字 ,只是希望给大
人生苦短,编程苦手,不妨学起Python,感受一飞冲天的快乐。不要害怕学习的过程枯燥无味,这里有程序员jackzhenguo打造的一份中文Python“糖果包”:147个代码小样,60秒一口,营养又好玩,从Python基础到机器学习尽皆囊括。入门简单如十进制转二进制,尽显Python简洁之美:In [1]: bin(10)Out[1]: 0b1010冬天到了,就算没有点亮手绘技能,也能用简单几行...
JNI Android NKD初学者遇到问题不要气馁!
【ElasticSearch】精确匹配text字段 用match加.keyword 或 term1.错误示范由于记忆混淆,记成了使用match_phrase对text字段精确匹配。#测试match_phraseGET /test/external/_search{ "query":{ "bool": { "must": [ { ...
前端开发VS Code下载地址https://code.visualstudio.com/初始设置1、中文界面配置首先安装中文插件:Chinese (Simplified) Language Pack for Visual Studio Code右下角弹出是否重启vs,点击“yes”重启后如果界面没有变化,则 点击 左边栏Manage -> Command Paletet…【...
题目描述Description有 N 堆纸牌,编号分别为 1,2,…, N。每堆上有若干张,但纸牌总数必为 N 的倍数。可以在任一堆上取若于张纸牌,然后移动。 移牌规则为:在编号为 1 堆上取的纸牌,只能移到编号为 2 的堆上;在编号为 N 的堆上取的纸牌,只能移到编号为 N-1 的堆上;其他堆上取的纸牌,可以移到相邻左边或右边的堆上。 现在要求找出一种移动方法,用最少...