数值分析期末复习

第一章 科学计算

误差

解题步骤

  1. 先求绝对误差:
    ∣ x − x ∗ ∣ |x - x^*| xx
  2. 求相对误差限:
    ∣ x − x ∗ ∣ x ∗ \frac{|x\,\,-\,\,x^*|}{x^*} xxx
  3. 求有效数字
    ∣ x − x ∗ ∣ 需要小于它自身的半个单位 |x-x^*|\text{需要小于它自身的半个单位} xx需要小于它自身的半个单位,然后算小数点后一共有多少数字

举个例子:
相减得出结果为0.0000345则小于0.0005,则有效数字为4

例题1:
ea1dba3fae17eef446cea04935ac79c.jpg

第二章 线性代数直接法

高斯消去法

高斯顺序消去法解题步骤

(假设是一个三行三列的矩阵):

  1. 先用第一行消去2,3行
  2. 再用第二行消去第三行

例题1:
987c473889804cafbae9f9f1718c044.jpg
例题2:
cc0434961ba0e72df953fa4ee6a96c4.jpg

高斯列主元消去法解题步骤

  1. 比较哪一行的绝对值最大,然后交换
  2. 用第一行消去第2、3行
  3. 再次比较哪一行绝对值最大,交换
  4. 重复步骤

例题1:
52e744960581f51f8b2ce30cd385017.jpg
例题2:
2579dce62a578c83fc05cc25928d28a.jpg
a661f62d416d4d48b9a94a9df97078d.jpg

LU分解

LU分解又称为:杜利特尔 (Doolittle)分解法,直接三角分解法

解题步骤

  1. 将A矩阵分解成L、U矩阵
    1. L矩阵:下三角矩阵,对角线全为1,其他元素为x
    2. U矩阵:上三角矩阵,第一行元素和A矩阵相同,其他元素为x
  2. 从A中矩阵逆向推导,L、U剩下的元素逐一相乘得出结果
    1. 按照顺序一行一行的元素去算

例题1:
2e32d42f8f43c0ae85d2720b61fbba7.jpg

3915868f5273582e7f11b0177bac7c8.jpg
31b5216620ff5a2f33ec9ef7ec72958.jpg

追赶法

追赶法又称为:克劳特分解

解题步骤

  1. 将A矩阵分解为L、U矩阵
  2. L矩阵的特点:下三角矩阵,对角线为未知数 α \alpha α,其他元素对A照抄
  3. U矩阵的特点:上三角矩阵,对角线为1,对角线上面的元素为 β \beta β
  4. α , β \alpha,\beta α,β全部算出来
  5. L y = b Ly=b Ly=b -> U x = y Ux=y Ux=y

例题:
899fe9ec4c43c1e320cf522730c939e.jpg
212cd76b186e8a5d88ad977e8fb26bf.jpg
84492d380115456b36a99289cdba506.jpg

第三章 线性代数方程组的迭代法

范数和条件数

  1. 1范数(列范数):每一列元素的绝对值之和的最大值 ∣ ∣ A ∣ ∣ 1 ||A||_1 ∣∣A1
  2. 无穷范数(行范数):每一行元素的绝对值之和的最大值
  3. 2范数:
    1. 向量:向量元素平方的和的平方根
    2. 矩阵(又称为谱范数):null
  4. 无穷范数条件数:
    c o n d ∞ ( A ) = ∣ ∣ A ∣ ∣ ∞ ∣ ∣ A − 1 ∣ ∣ ∞ cond_{\infty}\left( A \right) \,\,=\,\,||A||_{\infty}||A^{-1}||_{\infty} cond(A)=∣∣A∣∣A1

例题1:
f690ea10aedc5224444911545e27b35.jpg
例题2:
32ddcfdd090dda264110dcbf51eabc5.jpg

A − 1 A^{-1} A1的方法

  1. 初等变换法
    image.png

雅可比迭代法

解题步骤

  1. 整体思路:将 Ax=b ->x=Bx+g 的形式
  2. 先将第一行转换为 x 1 = . . . x_1=... x1=...
  3. 第二行 x 3 = . . . . x_{3}= .... x3=....
  4. 以此类推
  5. 画出表格
    image.png

计算器解题步骤

  1. 先将A、B、C、D、E、F设置为0 ( A 代表 x 1 , B 代表 x 2 , C 代表 x 3 A代表x_1,B代表x_2,C代表x_3 A代表x1,B代表x2,C代表x3)
    1. 0 sto A
    2. 0 sto B
    3. 0 sto C
    4. 0 sto D
    5. 0 sto E
    6. 0 sto F
  2. 将每一行公式输入到计算器中,使用 : : :进行分割
    1. D = …:E = …:F = …:A=D:B=E:C=F
      1. 这里是因为一开始不迭代,所以要设置DEF

高斯迭代法

解题步骤

  1. 与雅可比迭代类似
  2. 但是每次都会迭代前面那个值

计算器解题步骤

  1. 先将A、B、C设置为0 ( A 代表 x 1 , B 代表 x 2 , C 代表 x 3 A代表x_1,B代表x_2,C代表x_3 A代表x1,B代表x2,C代表x3)
    1. 0 sto A
    2. 0 sto B
    3. 0 sto C
  2. 将每一行公式输入到计算器中,使用 : : :进行分割
    1. A = …:B = …:C = …

第四章 多项式插值和样条插值

拉格朗日插值

一共 2 个部分:

  1. 插值多项式
  2. 插值余项

插值多项式

  • l n ( x ) = [ ∏ i = 0 , i ≠ j n x − x i x j − x i ] y i l_n(x) =[ \prod_{i=0,i\ne j}^{n}\frac{x-x_i}{x_j-xi}] y_i ln(x)=[i=0,i=jnxjxixxi]yi
  • L n ( x ) = ∑ j = 0 n L j ( x ) y j L_n(x)=\sum_{j=0}^{n}L_j(x)y_j Ln(x)=j=0nLj(x)yj

线性 n=1,以此类推后面就是 2 次、3 次

插值余项

  • ∣ R n ( x ) ∣ = M n + 1 ( n + 1 ) ! ∣ W n + 1 ( x ) ∣ |R_n(x)|=\frac{M_{n+1}}{(n+1)!} |W_{n+1}(x)| Rn(x)=(n+1)!Mn+1Wn+1(x)
  • M n + 1 = max ⁡ a ≤ x ≤ b ∣ f n + 1 ( x ) ∣ M_{n+1} = \max_{a\le x\le b}|f^{n+1}(x)| Mn+1=maxaxbfn+1(x)
  • W n + 1 ( x ) = ( x − x 0 ) ( x − x 1 ) . . . ( x − x n ) W_{n+1}(x) = (x-x_0)(x-x_1)...(x-x_n) Wn+1(x)=(xx0)(xx1)...(xxn)

牛顿插值

插值多项式

  • f ( x 0 , . . . , x n ) = f n − f n − 1 x n − x 0 f(x_0,...,x_n) = \frac{f_n - f_{n-1}}{x_n-x_0} f(x0,...,xn)=xnx0fnfn1

解题步骤

  1. 列差商表
xf(x)一阶差商
x 0 x_0 x0 f 0 f_0 f0
x 1 x_1 x1 f 1 f_1 f1 f ( x 0 , x 1 ) f(x_0,x_1) f(x0,x1)

以此类推,有 n 个 x 的值就有多少次 n-1 阶差商

  1. 最后的结果公式
    N n ( x ) = f 0 + f ( x 0 , x 1 ) ( x − x 0 ) + . . . + f ( x 0 , . . . , x n + 1 ) ( x − x 0 ) ( x − x 1 ) . . . ( x − x n ) N_{n}(x)=f_0 + f(x_0,x_1)(x-x_0) +...+f(x_0,...,x_{n+1})(x-x_0)(x-x_1)...(x-x_n) Nn(x)=f0+f(x0,x1)(xx0)+...+f(x0,...,xn+1)(xx0)(xx1)...(xxn)

牛顿插值余项

需要补充

第九章 常微分方程初边值问题数值解

龙格-库塔公式

基本概念

一般问题会有 y ′ , h , f ( x ) = y y', h , f(x) = y y,h,f(x)=y等参数
将其转换为
image.png

注意h的值,一般是在 0 ≤ x ≤ 1 0 \le x \le 1 0x1之间,逐渐相加之后递增到1结束计算

四阶四段龙格库塔公式如下:
image.png

解题步骤

  1. x 0 , y 0 , h x_0,y_0,h x0,y0,h写在旁边
  2. 将将题目中给出的已知信息代入 k 1 , k 2 , k 3 , k 4 k_1,k_2,k_3,k_4 k1,k2,k3,k4
  3. 更新 y n y_n yn的值
  4. 重复过程

k 2 k_2 k2->f的 x n + h 2 x_n+\frac{h}{2} xn+2h表示 x x x,同理另外一个表示 y y y,将其代入到f(x,y)中进行化简

02a8186f69218c968ad62f27dd644e3.jpg

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

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

相关文章

亚信安慧AntDB:支撑中国广电5G业务的数据库之力

自2019年6月获得5G牌照以来,中国广电积极利用700MHz频谱资源,迅速崛起为第四大运营商,标志着其在数字通信领域取得的巨大成就。通过与中国移动紧密合作,共建共享基站已超过400万座,为实现自主运营和差异化竞争提供了坚…

直接插入排序【从0-1学数据结构】

文章目录 💗 直接插入排序Java代码C代码JavaScript代码稳定性时间复杂度空间复杂度 我们先来学习 直接插入排序, 直接排序算是所有排序中最简单的了,代码也非常好实现,尽管直接插入排序很简单,但是我们依旧不可以上来就直接写代码,一定要分析之后才开始写,这样可以提…

微软官方出品:GPT大模型编排工具,支持C#、Python等多个语言版本

随着ChatGPT的火热,基于大模型开发应用已经成为新的风口。虽然目前的大型模型已经具备相当高的智能水平,但它们仍然无法完全实现业务流程的自动化,从而达到用户的目标。 微软官方开源的Semantic Kernel的AI编排工具,就可以很好的…

设计模式(三)-结构型模式(6)-享元模式

一、为何需要享元模式(Flyweight)? 假如在网页中渲染这样的一个画面:大小不一的星星铺满了整个画布,并且都在不断的进行移动闪烁着。一批星星消失了,另一批又从另一边缘处出现。 要实现这样的渲染效果,在…

实习课知识整理2:用户登录及实现登录后用户名和头像的展示

接上一篇&#xff0c;当用户点击购买按钮后&#xff0c;还是未登录的状态&#xff0c;此时页面会跳转到登录页面&#xff0c;这时需要输入正确的用户名和密码&#xff0c;完成登录 1. 给登录按钮添加点击事件&#xff0c;并提交表单中的数据到后端 <form th:action"{/u…

Elasticsearch Reroute API 的使用

本文通过一个 Elasticsearch 集群中主分片分配不均衡的例子演示一下 Cluster reroute API 的使用。 对于 Elasticsearch 分片分配策略不了解的同学可以点一下关注&#xff0c;后面更文之后获取第一手资料。 环境信息 Windows 10 Elasticsearch 8.1 JDK17 初始集群状态 分片…

【JAVA面试题】什么是引用传递?什么是值传递?

&#x1f34e;个人博客&#xff1a;个人主页 &#x1f3c6;个人专栏&#xff1a; JAVA ⛳️ 功不唐捐&#xff0c;玉汝于成 前言 博客的正文部分可以详细介绍Java中参数传递的机制&#xff0c;强调Java是按值传递的&#xff0c;并解释了基本数据类型和对象引用在这种传…

鳄鱼目标检测数据集VOC格式100张

鳄鱼是一种生活在热带和亚热带地区的爬行动物&#xff0c;属于爬行纲鳄形目鳄鱼科。它们的体形庞大&#xff0c;有粗壮的四肢和强壮的尾巴&#xff0c;一般能长到2-6米长&#xff0c;体重可达500公斤以上。鳄鱼的皮肤粗糙&#xff0c;呈灰褐色或黑色&#xff0c;布满了坚韧的鳞…

XSKY星辰天合星海架构荣获 IT168 “2023 年度技术卓越奖”

近日&#xff0c;"2023 年度技术卓越奖"获奖名单公布&#xff0c;XSKY 星辰天合的星海架构&#xff08;XSEA&#xff0c;极速全共享架构&#xff09;获得行业 CIO/CTO 大咖、技术专家及 IT 媒体三方认可&#xff0c;成功入选&#xff01; “技术卓越奖”评选由国内著…

【Linux驱动】字符设备驱动程序框架 | LED驱动

&#x1f431;作者&#xff1a;一只大喵咪1201 &#x1f431;专栏&#xff1a;《RTOS学习》 &#x1f525;格言&#xff1a;你只管努力&#xff0c;剩下的交给时间&#xff01; 目录 &#x1f3c0;Hello驱动程序⚽驱动程序框架⚽编程 &#x1f3c0;LED驱动⚽配置GPIO⚽编程驱动…

最小二乘法简介

最小二乘法简介 1、背景描述2、最小二乘法2.1、最小二乘准则2.2、最小二乘法 3、最小二乘法与线性回归3.1、最小二乘法与线性回归3.2、最小二乘法与最大似然估计 4、正态分布&#xff08;高斯分布&#xff09; 1、背景描述 在工程应用中&#xff0c;我们通常会用一组观测数据去…

电商数仓项目----笔记六(数仓ODS层)

ODS层的设计要点如下&#xff1a; &#xff08;1&#xff09;ODS层的表结构设计依托于从业务系统同步过来的数据结构。 &#xff08;2&#xff09;ODS层要保存全部历史数据&#xff0c;故其压缩格式应选择压缩比较高的&#xff0c;此处选择gzip。 &#xff08;3&#xff09;…

C++入门-【13-C++ 多维数组】

C 多维数组 C 支持多维数组。多维数组声明的一般形式如下&#xff1a; type name[size1][size2]...[sizeN]; 例如&#xff0c;下面的声明创建了一个三维 5 . 10 . 4 整型数组&#xff1a; int threedim[5][10][4]; 二维数组 多维数组最简单的形式是二维数组。一个二维数组&am…

用23种设计模式打造一个cocos creator的游戏框架----(二十三)中介者模式

1、模式标准 模式名称&#xff1a;中介者模式 模式分类&#xff1a;行为型 模式意图&#xff1a;用一个中介对象来封装一系列的对象交互。中介者使各对象不需要显式地相互引用&#xff0c;从而使其耦合松散&#xff0c;而且可以独立地改变它们之间的交互。 结构图&#xff…

竞赛保研 基于GRU的 电影评论情感分析 - python 深度学习 情感分类

文章目录 1 前言1.1 项目介绍 2 情感分类介绍3 数据集4 实现4.1 数据预处理4.2 构建网络4.3 训练模型4.4 模型评估4.5 模型预测 5 最后 1 前言 &#x1f525; 优质竞赛项目系列&#xff0c;今天要分享的是 基于GRU的 电影评论情感分析 该项目较为新颖&#xff0c;适合作为竞…

Linux基本指令(一)

前言 基本知识 文件文件内容文件属性(对文件的操作就是对这两部分进行操作) 在Linux中以 . 开头的文件叫隐藏文件 以-开头的是普通文件 以d开头的是目录文件 几个指令 先快速认识几个指令&#xff0c;方便后续的详细介绍 whoami 查看当前使用Linux系统的用户是谁 pwd …

要参加微软官方 Copilot 智能编程训练营了

GitHub Copilot 是由 GitHub、OpenAI 和 Microsoft 联合开发的生成式 AI 模型驱动的。 GitHub Copilot 分析用户正在编辑的文件及相关文件的上下文&#xff0c;并在编写代码时提供自动补全式的建议。 刚好下周要参加微软官方组织的 GitHub Copilot 工作坊-智能编程训练营&…

【51单片机系列】C51中的中断系统扩展实验

本文是关于51单片机中断系统的扩展实验。 文章目录 一、 扩展实验一&#xff1a;使用外部中断0控制蜂鸣器&#xff0c;外部中断1控制直流电机二、扩展实验二&#xff1a;修改定时器初值&#xff0c;设定3秒钟的定时时间让LED模块闪烁三、扩展实验三&#xff1a;使用定时器1和数…

法线贴图实现地形模型皱褶、凹凸不平的纹理效果

在线工具推荐&#xff1a; 3D数字孪生场景编辑器 - GLTF/GLB材质纹理编辑器 - 3D模型在线转换 - Three.js AI自动纹理开发包 - YOLO 虚幻合成数据生成器 - 三维模型预览图生成器 - 3D模型语义搜索引擎 法线贴图在3D建模中扮演着重要的角色&#xff0c;它通过模拟表面的微…

(十七)Flask之大型项目目录结构示例【二扣蓝图】

大型项目目录结构&#xff1a; 问题引入&#xff1a; 在上篇文章讲蓝图的时候我给了一个demo项目&#xff0c;其中templates和static都各自只有一个&#xff0c;这就意味着所有app的模板和静态文件都放在了一起&#xff0c;如果项目比较大的话&#xff0c;这就非常乱&#xf…