数值分析期末复习

第一章 科学计算

误差

解题步骤

  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.hqwc.cn/news/295393.html

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

相关文章

pycharm git 版本回退

参考 https://blog.csdn.net/qq_38175912/article/details/102860195 yoyoketang 悠悠课堂

二维码初体验 com.google.zxing 实现续 - web api封装

文章目录 一、概述二、最终效果三、源码结构四、完整代码 一、概述 在 二维码初体验 com.google.zxing 实现 我们实现了二维码的生成,但是大部分情况下,二维码的相关功能是作为API接口来提供服务的。 我们下面便演示在springboot、Knife4j下封装api接口…

Prometheus+Grafana搭建Jmeter性能监控平台

📢 专注于分享软件测试干货内容,欢迎点赞 👍 收藏 ⭐留言 📝 如有错误敬请指正!📢 交流讨论:欢迎加入我们一起学习!📢 资源分享:耗时200小时精选的「软件测试…

vue3项目 - 使用 pnpm 包管理器来创建项目

创建项目 npm install -g pnpm pnpm create vue 输入项目名称、包名称、选择要安装的依赖,最后 pnpm install pnpm format #规范格式 pnpm dev #启动项目

【教程】使用ipagurd打包与混淆Cocos2d-x的Lua脚本

文章目录 摘要引言正文1. 准备工作2. 使用ipaguard处理Lua文件3. 运行ipagurd进行混淆代码加密具体步骤测试和配置阶段IPA 重签名操作步骤4. IPA重签名与发布 总结 摘要 本文将介绍如何使用ipagurd工具对Cocos2d-x中的Lua脚本进行打包与混淆,以及在iOS应用开发中的…

【华为OD机试真题2023CD卷 JAVAJS】多段线数据压缩

华为OD2023(C&D卷)机试题库全覆盖,刷题指南点这里 多段线数据压缩 知识点数组栈递归矩阵循环 时间限制:1s 空间限制:256MB 限定语言:不限 题目描述: 下图中,每个方块代表一个像素,每个像素用其行号和列号表示。 为简化处理,多段线的走向只能是水平、竖直、斜向45…

地震勘探原理---数字滤波处理

目录 一. 地震数字滤波的目标 二. 数字滤波器 2.1 数字滤波器的分类 三. 一维数字滤波器 3.1 傅里叶变换与傅里叶逆变换 3.2 滤波流程 四. 二维数字滤波 4.1为什么有二维数字滤波 4.2 f-k域滤波 4.3 τ-p域滤波 4.4 相干滤波 四. 总结 一. 地震数字滤波的目标 核心任务&am…

使用 Spring Boot + MyBatis开发需要注意的事项以及开发模版

前言: 注意,本篇不适用于有相关开发经验的开发者,作为一个在职开发者,我经常在完成从0-1的模块,也就是从数据库表开始到创建实体类,以及dao层,Service层等业务需要添加相关注解,这样…

WEB渗透—PHP反序列化(八)

Web渗透—PHP反序列化 课程学习分享(课程非本人制作,仅提供学习分享) 靶场下载地址:GitHub - mcc0624/php_ser_Class: php反序列化靶场课程,基于课程制作的靶场 课程地址:PHP反序列化漏洞学习_哔哩…

【C语言】打印内存数据

C语言,用函数封装:16进制打印unsigned char *p指向的内存,长度为int l。16个字节,换一次行。16个字节用一个字符串缓存,一次打印。 以下是一个使用函数封装的C语言代码,用于以16进制格式打印unsigned char …

Tomcat远程调试

windows环境 写一个 startup-debug.bat,指定tomcat的根目录,端口自己定义 rem *******设置Tomcat目录*******-- set CATALINE_HOMED:\asd\A8-2\tomcat d: rem 8787为可用端口,为远程调试监听端口-- cd %CATALINE_HOME%/bin set JPDA_ADDRESS8787 set J…

单片机原理及应用:流水灯的点亮

流水灯是一种简单的单片机控制电路&#xff0c;由许多LED组成&#xff0c;电路工作时LED会按顺序点亮&#xff0c;类似于流水的效果。 下面是运行在keil上的代码&#xff0c;分别使用了数组&#xff0c;移位符和库函数来表示。 //数组法 #include <reg52.h> //头文…