matlab学习笔记第四章

最新推荐文章于 2023-10-29 21:16:11 发布

speop

最新推荐文章于 2023-10-29 21:16:11 发布

阅读量2.3k

点赞数 2

文章标签： matlab 开发语言

本文链接：https://blog.csdn.net/m0_52024881/article/details/121305507

版权

本文详细介绍了MATLAB中的数值计算方法，包括矩阵运算、线性方程组求解、向量除法、矩阵求逆以及条件数的概念。此外，还讲解了矩阵的乘方运算、特征值计算、微分方程求解和不同插值方法，如线性、最近点和三次样条插值。同时，提到了多项式拟合的polyfit和polyval函数，以及如何进行图形可视化。

摘要由CSDN通过智能技术生成

第四讲数值计算

一.矩阵

在这里插入图片描述

example:

./每个元素相除

在这里插入图片描述

求解线性方程组

在这里插入图片描述

X=A\B
X=inv(A)*B

在这里插入图片描述

两个向量相除
右除就是模长相除

左除
在这里插入图片描述

求逆运算A没有逆，就添加零元素

help
mldivide - 对线性方程组 Ax = B 求解 x

此 MATLAB 函数 对线性方程组 A*x = B 求解。矩阵 A 和 B 必须具有相同的行数。如果 A 未正确缩放或接近奇异值，MATLAB
将会显示警告信息，但还是会执行计算。

x = A\B
x = mldivide(A,B)

help /
mrdivide - 对线性方程组 xA = B 求解 x

此 MATLAB 函数 对线性方程组 x*A = B 求解 x。矩阵 A 和 B 必须具有相同的列数。如果 A 未正确缩放或接近奇异值，MATLAB
将会显示警告信息，但还是会执行计算。

x = B/A
x = mrdivide(B,A)

x = B/A 对线性方程组 x*A = B 求解 x。矩阵 A 和 B 必须具有相同的列数。如果 A 未正确缩放或接近奇异值，MATLAB® 将会显示警告信息，但还是会执行计算。

如果 A 是标量，那么 B/A 等于 B./A。

如果 A 是 n×n 方阵，B 是 n 列矩阵，那么 x = B/A 是方程 x*A = B 的解（如果存在解的话）。

如果 A 是矩形 m×n 矩阵，且 m ~= n，B 是 n 列矩阵，那么 x=B/A 返回方程组 x*A = B 的最小二乘解。

x = mrdivide(B,A) 是执行 x=B/A 这一操作的替代方法，但很少使用。它可以启用类的运算符重载。

在这里插入图片描述

矩阵的乘方运算

在这里插入图片描述
运算规则

矩阵在工作精度下奇异的意思就是矩阵的条件数乘上工作精度达到Θ (1)或者更大，也就是说条件数超过10^16 9

在这里插入图片描述

矩阵的特征值

eig
eig(A)
在这里插入图片描述

6特征值多项式

在这里插入图片描述

7矩阵的秩

在这里插入图片描述

8求矩阵元素的个数

在这里插入图片描述

方阵的行列式

在这里插入图片描述

方阵的迹（对角线元素之和）

在这里插入图片描述

summary

在这里插入图片描述

常用的符号运算

极限

在这里插入图片描述

显示声明自变量

在这里插入图片描述

极限

在这里插入图片描述

在这里插入图片描述
左趋近
右趋近

第四个参数注意right left

微分diff

在这里插入图片描述
diff第二个参数是求几次微分

diff 求差分和求微分
第三个参数为1时,则变为行差分运算
第三个参数为2时,则变为列差分运算

diff(A,2,2)第2个参数为2表示为二阶差分变换.即在diff(A,1,2)的基础上再进行一次列差分变换.

积分int

函数int

int(f,v,a,b)，其中，a、b分别表示定积分的下限和上限。该函数求被积函数f在区间[a,b]上的定积分。a 和 b 可以是两个具体的数，也可以是一个符号表达式，还可以是无穷（inf(极限))。当函数f关于变量x在闭区间[a,b]上可积时，函数返回一个定积分结果。当a,b中有一个是inf 时，函数返回一个广义积分；当a、b 中有一个有符号表达式时，函数返回一个符号函数。
中文名符号函数的定积分和不定积分
在这里插入图片描述

对什么变量求积分

在这里插入图片描述

微分方程求解

在这里插入图片描述
dsolve求解微分方程

默认变量是t

在这里插入图片描述

研究一下常微分方程和偏微分方程的求解

pde ode

插值

在这里插入图片描述

linear线性插值
默认插值方式，把插值点靠近的两个数据点用直线连接，然后在直线上个选取对应插值点的数据

nearest最近点插值

cubic 3次多项式插值求出 pchip
spline3次样条插值构造光滑

在这里插入图片描述

在什么位置
x，y
插什么值x1
采用什么方法插值
线性方法，最近方法，三次样条，三次插值
Y1=interp1(X,Y,X1,method)
在这里插入图片描述
默认的插值方法是linear线性插值法

在这里插入图片描述

依赖于函数的表达式哪个是最好的
找最近的点模型和函数

插值和拟合之间的区别:
polyfit三次方拟合

在这里插入图片描述

p=polyfit(x,y,n) 多项式曲线拟合函数
[p,s]= polyfit(x,y,n)
说明：x,y为数据点，n为多项式阶数，返回p为幂次从高到低的多项式系数向量p。x必须是单调的。矩阵s用于生成预测值的误差估计。

多项式曲线求值函数：polyval( )
调用格式： y=polyval(p,x)
[y,DELTA]=polyval(p,x,s)
说明：y=polyval(p,x)为返回对应自变量x在给定系数P的多项式的值。
[y,DELTA]=polyval(p,x,s) 使用polyfit函数的选项输出s得出误差估计Y DELTA。它假设polyfit函数数据输入的误差是独立正态的，并且方差为常数。则Y DELTA将至少包含50%的预测值。
polyfit多项式曲线拟合
polyval