matlab学习笔记第四章

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

speoki

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

阅读量2k

点赞数 2

文章标签： matlab 开发语言

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

版权

第四讲数值计算

一.矩阵

在这里插入图片描述

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

二维插值 interp2()

在这里插入图片描述

mesh函数三个变量

如果数据很少的话就不够平滑
在这里插入图片描述

在这里插入图片描述

为什么这里要加上20？
在这里插入图片描述
加了20的和没加的没什么差别,就是z的坐标变了而已

在这里插入图片描述
注意转置矩阵
meshgrid

三维插值

在这里插入图片描述

插值函数总结

在这里插入图片描述
sin 三角函数找傅里叶
栅格数据插值
分段多项式插值精细一些

拟合

在这里插入图片描述
m几次多项式
M次方越多计算量越多

polyfit
p有6个数 s是误差
s是结构体，里面有在这里插入图片描述
s可以用作polyval的输入来获取误差估计值

在这里插入图片描述

polyval
在这里插入图片描述

在这里插入图片描述

在这里插入图片描述

5阶拟合

在这里插入图片描述

在这里插入图片描述
图形窗口菜单可视化操作

speoki

关注

2
点赞
踩
6

收藏

觉得还不错? 一键收藏
0
评论
matlab学习笔记第四章

第四讲数值计算一.矩阵example:./每个元素相除求解线性方程组X=A\BX=inv(A)*B两个向量相除右除就是模长相除左除求逆运算A没有逆，就添加零元素help mldivide - 对线性方程组 Ax = B 求解 x此 MATLAB 函数对线性方程组 A*x = B 求解。矩阵 A 和 B 必须具有相同的行数。如果 A 未正确缩放或接近奇异值，MATLAB将会显示警告信息，但还是会执行计算。x = A\Bx = mldivide(A,
复制链接

扫一扫