Matlab基础（三）——数值计算

已于 2024-04-10 22:19:16 修改
阅读量783 收藏 29
点赞数 15
分类专栏： MATLAB学习笔记 文章标签： matlab 开发语言
于 2024-04-10 20:56:14 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/m0_74970386/article/details/137220373
版权

Matlab基础（三）——数值计算

一、矩阵运算

（1）矩阵的加减法

矩阵与标量的加减等于矩阵中所有元素与此标量的加减.

c =
     1     2     3
     4     5     6
>> c+2
ans =
     3     4     5
     6     7     8

（2）矩阵的乘法

①正常乘（*）

C = A*B 是 A 和 B 的矩阵乘积。如果 A 是 m×p 矩阵，B 是 p×n 矩阵，则 C 是通过以下公式定义的 m×n 矩阵：
在这里插入图片描述

②数乘（*）

用这个常数去乘矩阵的每一个元素

c =
     1     2     3
     4     5     6
>> 3*c
ans =
     3     6     9
    12    15    18

③点乘（.*）

要求两个矩阵有相同的维数，对应元素相乘

d =
     3     6     9
    12    15    18
>> c
c =
     1     2     3
     4     5     6
>> x=d.*c
x =
     3    12    27
    48    75   108

④标量积（dot）

返回 A 和 B 的标量点积（内积）
语法：

C = dot(A,B)
C = dot(A,B,dim)

dim：按指定维数dim求两矩阵的标量积
沿其运算的维度，指定为正整数标量。如果未指定值，则默认值是大小不等于 1 的第一个数组维度。

考虑两个二维输入数组：A 和 B：

dot(A,B,1) 将 A 和 B 的列视为向量，并返回对应列的点积。

dot(A,B,2) 将 A 和 B 的行视为向量，并返回对应行的点积。
在这里插入图片描述

如果 dim 大于 ndims(A)，则 dot 返回 conj(A).*B。

⑤向量积（cross）

返回 A 和 B 的叉积
cross(A,B)
cross(A,B,dim)

⑥混合积

dot(A,cross(B,C))

⑦卷积（conv）

conv(A,B)

⑧解卷积

[q,r]=deconv(B,A)

⑨张量积（kron）

kron(A,B)
 

（3）矩阵的除法

通过矩阵的除法可以求解线性方程组的解

①左除（ \ ）

当A、B是线性方程组AX=B中的系数矩阵和常数矩阵时，A \ B是方程组AX=B的解

②右除

当A、B是线性方程组XA=B中的系数矩阵和常数矩阵时，B \ A是方程组XA=B*的解

③元素左除

A .\ B 是指B的个元素除以A的对应元素

④元素右除

A ./ B 是指A的个元素除以B的对应元素

（4）矩阵的乘方

  • 矩阵的乘法:
    调用格式：A^n,n为整数，A是方阵，计算矩阵A的n次方
 若n>0，则计算A自乘n次
 若n=0，则计算结果为与A同阶的单位阵
 若n<0，则计算A的逆矩阵自乘n次

  若n为非整数，A为方阵
A = p   λ t p − 1 A=p\ \lambda^{t}p^{-1} A=p λtp−1
 其中 lambda为A的特征对角矩阵，p为A的特征向量矩阵

  • 矩阵元素的乘方
    调用格式：A.^b,b可为非整数，A为方阵，计算A中每一个元素的b次方组成的矩阵。
  • 数的矩阵元素的乘方
    调用格式：b.^A, b为非整数，A为方阵，计算常数b的A的每一个元素次方排成的矩阵
  • 数的矩阵乘方:
    调用格式： b.^A ,b为非整数，A为方阵
    A = p   b λ p − 1 A=p\ b^{\lambda}p^{-1} A=p bλp−1
     其中 lambda为A的特征对角矩阵，p为A的特征向量矩阵

（5）矩阵的操作与变换

在这里插入图片描述
inv(A):计算非奇异矩阵的逆矩阵
det(A):计算方阵的行列式
rank(A):计算方阵的秩
trace(A):计算方阵的迹
eig(A):计算n阶方阵A的特征值
circshift(A,[m,n]):对A矩阵进行m行，n列的圆周移位

二、多项式运算

（1）多项式的系数与根

① poly2sym

由多项式的系数得到多项式的表达式

>> p=[1 0 -5 7 -14 47]
p =
     1     0    -5     7   -14    47
>> poly2sym(p)
ans =
x^5 - 5*x^3 + 7*x^2 - 14*x + 47

②poly（r）

由此 MATLAB 函数（其中 r 是向量）返回多项式的系数，其中多项式的根是 r 的元素

③fzero(‘fun’,x0):

计算由函数fun=0决定的方程的根。
x0为常数时，计算x0附近的根；
如果x0为区间，计算区间内的根，若此区间无根，则输出错误信息。

求余弦函数在 1 和 2 之间的零点。

fun = @cos; % function
x0 = [1 2]; % initial interval
x = fzero(fun,x0)
x = 1.5708

请注意，cos(1) 和 cos(2) 的符号不同。

④roots( p )

计算多项式p(X)=0的解，结果可以为复数

 
 

（2）多项式的值

①polyval(p,x)

计算多项式 p 在 x 的每个点处的值。

②[y,delta]=polyval(p,x,s)

同①，y是多项式的值，delta是根据polyfit命令给出的矩阵s的误差估计向量

③polyvalm(p,A)

以矩阵方式返回多项式 p 的计算值，是计算：
P ( A ) = a 0 A n + a 1 A n − 1 + a 2 A n − 2 + … … + a n I P(A)=a_0A^n+a_1A^{n-1}+a_2A^{n-2}+……+a_nI P(A)=a0​An+a1​An−1+a2​An−2+……+an​I

 

（3）多项式的拟合与插值

①ployfit

此 MATLAB 函数 返回次数为 n 的多项式 p(x) 的系数，该阶数是 y 中数据的最佳拟合（基于最小二乘指标）。
p 中的系数按降幂排列，p 的长度为 n+1
语法：p = polyfit(x,y,n)

x = 1:50; 
y = -0.3*x + 2*randn(1,50); 
p = polyfit(x,y,1);% 进行线性回归
poly2sym(p)

输出：

ans =
 
8881560245387251/18014398509481984 - (5883261650978903*x)/18014398509481984

②interp1(x,y,xi,‘method’)

数据插值是通过已知的有限个数据点，推导出数据点之间未知位置上的数值。
vq = interp1(x,v,xq,method) 指定备选插值方法：'linear'、'nearest'、'next'、'previous'、'pchip'、'cubic'、'v5cubic'、'makima' 或 'spline'。默认方法为 ‘linear’。

示例:

x = 0:pi/4:2*pi; 
v = sin(x);
xq = 0:pi/16:2*pi;
vq1 = interp1(x,v,xq);
plot(x,v,'o',xq,vq1,':.');

输出
在这里插入图片描述
在这里插入图片描述

vq = interp1(x,v,xq,method,extrapolation) 用于指定外插策略，来计算落在 x 域范围外的点。如果希望使用 method 算法进行外插，可将 extrapolation 设置为 'extrap'。您也可以指定一个标量值，这种情况下，interp1 将为所有落在 x 域范围外的点返回该标量值。

 
 

（4)多项式的运算

①conv(A,B)

计算多项式A与B的乘积，仍以向量形式输出

②[q,r]=deconv(A,B)

计算多项式A与B的商，q为多项式的系数向量，r为A、B相除后的残值多项式系数向量

③polyder(A)

计算多项式A的导数
polyder(A,B)：计算多项式A与B乘积的导数
[q,r]=polyder(A,B)：计算多项式A与B的商的导数，其中，q为导数分子多项式的系数向量，r为导数分母多项式的系数向量

④ployint(A,k)

计算多项式A的积分，参数k为积分常数，省略时积分常数为零。

a =
     1     2     1

>> b=poly2sym(a)
b =
x^2 + 2*x + 1
 
>> c=poly2sym(polyint(a))
c =
x^3/3 + x^2 + x

⑤residue

[r,p,k] = residue(b,a)

在这里插入图片描述
定义较为复杂，可在Matlab中查看示例（help residue）

  • 例如：使用 residue 求以下多项式之比 F(s) 的部分分式展开式
    在这里插入图片描述
		b = [-4 8];
		a = [1 6 8];
		[r,p,k] = residue(b,a)

输出：

r =
   -12
     8 
p =
    -4
    -2
k =
     []

在这里插入图片描述

也可以使用 residue 将部分分式展开转换回多项式系数。

[b,a] = residue(r,p,k)

此处部分分式展开的计算可以与“信号与系统”中求解单位冲激响应联系起来。

 
 
 

三、数理统计

1.数据的最值与排序

(1)sort(x)

B = sort(A) 按升序对 A 的元素进行排序。

  • 如果 A 是向量，则 sort(A) 对向量元素进行排序。
  • 如果 A 是矩阵，则 sort(A) 会将 A 的列视为向量并对每列进行排序。
  • 如果 A 是多维数组，则 sort(A) 沿大小不等于 1 的第一个数组维度计算，并将这些元素视为向量。

B = sort(A,dim) 返回 A 沿维度 dim 的排序元素。
如果 A 是一个矩阵，则 sort(A,2) 对每行中的元素进行排序。

[B,I] = sort(___) 还会为上述任意语法返回一个索引向量的集合。
例如，如果 A 是一个向量，则 B = A(I)。

[C,I]=sort(A,'descend')—— 降序排列
 

(2)sortrows(x,col)

B = sortrows(A) 基于第一列中的元素按升序对矩阵行进行排序。当第一列包含重复的元素时，sortrows 会根据下一列中的值进行排序，并对后续的相等值重复此行为。

B = sortrows(A,column) 基于向量 column 中指定的列对 A 进行排序。例如，sortrows(A,4) 会基于第四列中的元素按升序对 A 的行进行排序。sortrows(A,[4 6]) 首先基于第四列中的元素，然后基于第六列中的元素，对 A 的行进行排序。
 

(3)max()

M = max(A) 返回数组的最大元素，如果元素是复数，则求最大模。

  • 如果 A 是向量，则 max(A) 返回 A 的最大值。

  • 如果 A 为矩阵，则 max(A) 是包含 A 的每一列的最大值的行向量。

  • 如果 A 是多维数组，则 max(A) 沿大小大于 1 的 A 的第一个维度计算，并将这些元素视为向量。此维度中 M 的大小变为 1，而所有其他维度的大小仍与在 A 中相同。如果 A 是第一个维度长度为零的空数组，则 M 是与 A 大小相同的空数组。

  • 如果 A 是表或时间表，则 max(A) 返回包含每个变量的最大值的一个单行表。 (自 R2023a 起)

M = max(A,[],"all") 查找 A 的所有元素的最大值。

M = max(A,[],dim) 返回维度 dim 上的最大元素。
例如，如果 A 为矩阵，则 max(A,[],2) 返回包含每一行的最大值的列向量。

[M,ind]=max(A)求矩阵A中每一列中最大值组成的行向量M，并给出最大值的行下标向量。

max(A,B)求矩阵A、B中对应元素最大值组成的矩阵
 

2.求和、乘积与差分

(1)sum()

S = sum(A) 返回沿大小大于 1 的第一个数组维度计算的元素之和。

  • 如果 A 是向量，则 sum(A) 返回元素之和。

  • 如果 A 是矩阵，则 sum(A) 将返回包含每列总和的行向量。
    S = sum(A,"all") 返回 A 的所有元素的总和。

  • 如果 A 是多维数组，则 sum(A) 沿大小大于 1 的第一个数组维度计算，并将这些元素视为向量。此维度中 S 的大小变为 1，而所有其他维度的大小仍与在 A 中相同。

  • 如果 A 是表或时间表，则 sum(A) 返回一个包含每个变量总和的单行表。 (自 R2023a 起)

(2)prod()

B = prod(A) 返回 A 的数组元素的乘积。

  • 如果 A 是向量，则 prod(A) 返回元素的乘积。

  • 如果 A 为非空矩阵，则 prod(A) 将 A 的各列视为向量，并返回一个包含每列乘积的行向量。
    B = prod(A,"all") 返回 A 的所有元素的乘积。

  • 如果 A 为 0×0 空矩阵，prod(A) 返回 1。

  • 如果 A 为多维数组，则 prod(A) 沿第一个非单一维度运算并返回乘积数组。此维度中 B 的大小将减少至 1，而所有其他维度的大小仍与在 A 中相同。

  • 如果 A 是表或时间表，则 prod(A) 返回每个变量乘积的单行表。 (自 R2023a 起)

(3)diff()

Y = diff(X) 计算沿大小不等于 1 的第一个数组维度的 X 相邻元素之间的差分：

  • 如果 X 是长度为 m 的向量，则 Y = diff(X) 返回长度为 m-1 的向量。Y 的元素是 X 相邻元素之间的差分。
    Y = [X(2)-X(1) X(3)-X(2) ... X(m)-X(m-1)]

  • 如果 X 是不为空的非向量 p×m 矩阵，则 Y = diff(X) 返回大小为 (p-1)×m 的矩阵，其元素是 X 的行之间的差分。
    Y = [X(2,:)-X(1,:); X(3,:)-X(2,:); ... X(p,:)-X(p-1,:)]

  • 如果 X 是 0×0 的空矩阵，则 Y = diff(X) 返回 0×0 的空矩阵。

3.平均值、中值与标准差

  • mean(x) —— 计算向量x元素的算数平均数
  • median(x) —— 计算向量x元素的中位数
  • std(x) —— 计算向量x元素的标准差

(各类用法与sum等十分相似)

4.协方差与相关系数

  • cov(x) —— 计算向量 x 的协方差
  • cov(x,y) —— 计算两个向量x,y的协方差
  • cov(A) —— 计算矩阵A的协方差矩阵
  • correoef(A) —— 计算矩阵A的相关系数矩阵
  • corrcoef (x,y) —— 计算两个向量x,y的相关系数

5.数据处理

  • ceil —— 向上取整
  • floor —— 向下取整
  • fix —— 向零取整
  • round —— 四舍五入
  • mod(x,y) —— 取余数
  • rem(x,y)—— 取余数
    (当x和y的正负号一样的时候，两个函数结果是等同的；当x和y的符号不同时，rem函数结果的符号和x的一样，而mod和y一样。)
  • abs(x) —— 对数据取模
  • angle(x) —— 对数据取相位
  • phase(x) —— 对数据取相位并解卷绕
  • real(x) —— 对数据取实部
  • imag(x) —— 对数据取虚部
Payton_Von
关注
  • 15
    点赞
  • 踩
  • 29
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
专栏目录
matlab计算一组数的函数值_MATLAB函数总结——数值运算和符号运算
weixin_39731782的博客
01-12 2710
系统命令：clear 后可加变量名，也可以加‘-expect c,a’除了a，c，clear all会清除全部变量，函数，mex文件cd 切换文件夹who 查看当前变量名称； whos 查看当前变量的详细信息clc 清除命令行窗口clf 清除当前的figureclose 关闭当前的figureclose all 关闭所有的figureexit或者quit 退出关于save函数：save filen...
matlab多元函数数值积分,MATLAB基础（七）——数值计算与数据分析（3）数值积分及常微分方程的数值解...
weixin_42293387的博客
03-16 1394
七、数值计算与数据分析(三)数值积分及常微分方程的数值解1.数值积分(1)一元函数的数值积分函数1 quad、quadl、quad8功能 数值定积分，自适应Simpleson积分法。格式 q = quad(fun,a,b) %近似地从a到b计算函数fun的数值积分，误差为10-6。若给fun输入向量x，应返回向量y，即fun是一单值函数。q = quad(fun,a,b,tol) ...
Matlab基础学习笔记（三）—— 数值计算
菜鸡成长之路
09-29 284
听万门的一个Matlab实践课做的笔记 知识结构 笔记链接（点我） 知识点比较零散记得比较概括，用幕布看可能更清晰一点 顺便放个邀请注册链接，可以领一个月高级版 注册幕布 目录笔记链接（点我）矩阵计算【线型代数】常用符号计算【高数】插值拟合 矩阵计算【线型代数】 乘法 * —— 矩阵相乘 .* —— 同型矩阵对应元素相乘 除法 / —— 右除 B/A = B*inv(A) \ —— 左除 ./ —— 元素右除 .\ —— 元素左除 乘方 ^ —— 矩阵的乘方.
MATLAB基础篇——基本语法
热门推荐
yun_gao_的博客
11-02 2万+
MATLAB基础篇——基本语法一、数据类型与变量数据类型变量二、矩阵字符串三、运算四、MATLAB常用函数五、矩阵分析与处理六、程序设计 MATLAB ————Matrix laboratory 一般操作： 1.操作界面：主窗口，命令窗口，工作空间窗口，当前目录窗口和搜索路径 2.帮助系统： help命令 lookfor命令 help 函数名 help（‘函数名’） lookfor 函数名 %lookfor命令更加宽松，检索所有m文件返回包含函数名的全部函数 3.注释：%这是一条注释 4.语句后加；
专题一MATLAB基础知识——1.2数值数据
Ricardo2的博客
04-06 1142
数值数据类型的分类 整型、浮点型、复数型 一、整形（8种） 1、无符号整数：无符号8位整数、无符号16位整数、无符号32位整数、无符号64位整数 有符号整数：有符号8位整数、有符号16位整数、有符号32位整数、有符号64位整数 2、无符号8位整数数据范围：00000000 ~ 11111111（0 ~ 2^7-1） 有符号8位整数数据范围：10000000 ~ 01111111（-2^7 ~ 2^7-1），0表示正数，1表示负数 3、unit8函数将数值数据转化为无符号8位整数，int8函数将数值数据转化为
数值计算方法——matlab实现
m0_63630166的博客
08-03 3509
MATLAB数值计算
matlab应用基础——工程应用
m0_55196097的博客
04-03 1664
1. fft（）函数 快速傅里叶变换 官方参考文档：快速傅里叶变换 - MATLAB fft - MathWorks 中国 语法： （持续更新。。。）
matlab数值过小为0,MATLAB数值计算——0
weixin_28889765的博客
03-18 2128
MATLAB数值计算MATLAB中文论坛基础板块常见问题归纳(出处: MATLAB中文论坛)登录http:www.mathworks.com/moler 获取NCM文件，注册账户，下载Toolbox格式，点击安装——在MATLAB中输入ncmgui出现图案即为success！(pathtool添加默认路径)开始学习：format 表示格式 help format 可以查询format帮助forma...
MATLAB基础篇——绘图
yun_gao_的博客
11-10 1203
MATLAB基础篇——绘图二维图形三维图形隐函数绘图 二维图形 将平面坐标上的数据点连接起来的平面图形 可以采用不同的坐标系——直角坐标，极坐标等 1.绘制二维曲线的基本函数 plot plot(x,y） plot(x,y,选项） plot(x1,y1,x2,y2,…) plot(x1,y1,选项1，x2，y2，选项2,…) 选项：线型（-，：，-·，–），颜色（b，y，g，w，k，等），标记符号（*，p，d，s，+等）选项 双纵坐标plotyy(x1,y1,x2,y2) 2.绘制图形的辅助操作 1).图形
工程数值计算matlab实验报告——辛普森数值微分数值积分.docx
07-14
工程数值计算matlab实验报告——辛普森数值微分数值积分.docx工程数值计算matlab实验报告——辛普森数值微分数值积分.docx工程数值计算matlab实验报告——辛普森数值微分数值积分.docx工程数值计算matlab实验报告...
工程数值计算matlab实验报告——辛普森数值微分数值积分.pdf
07-14
工程数值计算matlab实验报告——辛普森数值微分数值积分.pdf工程数值计算matlab实验报告——辛普森数值微分数值积分.pdf工程数值计算matlab实验报告——辛普森数值微分数值积分.pdf工程数值计算matlab实验报告——...
高中MATLAB创新课程教学实践与探索——以上海中学基于MATLAB的数值计算课程为例.pdf
06-23
为推动创新人才早期培育模式的探索,上海中学尝试将MATLAB课程引入高中课程体系。从学校的实际条件及学生的需求和能力出发,明确课程定位,并将其进一步细化为教学目标。在教学实践中采用学生实时操作的模式,注重创设...
Matlab经典教程——从入门到精通
11-20
分章阐述了矩阵计算、数值计算、符号计算、数据可视化、数据分析、M文件编写、Simulink仿真、句柄图形、图形用户界面、文件读写、MATLAB编译器、应用程序接口等内容；并结合相关的函数或命令，精心编写了一些具体的...
matlab教程PPT——MATLAB 语言及其应用（共340页）
10-06
讲授MATLAB语言基础入门知识，介绍MATLAB产品的体系、MATLAB桌面工具的使用方法，重点介绍MATLAB的数据可视化、数值计算的基本步骤以及如何使用MATLAB语言编写整洁、高效、规范的程序。并涉及到一些具体的专业应用...
PLC的ST语言实现IIR butterworth低通滤波器
exat500g的博客
05-15 209
PLC的ST语言实现IIR的butter低通滤波器
20240515matlab 动态命名存储文件
HITORANGE的博客
05-15 212
用sprintf命令格式化路径. 这里面的i是变量, 是用来对文件名编号的. %s 是把后面的’D:AlexStudyNew_Experiment_Data_CorridorShort3DPointsNImages’转换为字符串, %d是把后面的i转换为字符(这里i是整形数), 因为文件名都是字符串.这里的’lineN’是我们需要存储的数据, 大部分都是矩阵或者是数组之类的. 主意这里面的格式,‘‘这个很重要.save (path,‘lineN’,‘-ascii’);%linN代表变量，可以是数组。
254 基于matlab的钢筋混凝土非线性分析
最新发布
m0_59476991的博客
05-15 253
基于matlab的钢筋混凝土非线性分析，根据梁本构关系，然后进行非线性分析，绘制弯矩-曲率曲线。可设置梁的截面尺寸、混凝土本构，钢筋截面面积等相关参数，程序已调通，可直接运行。
matlab使用2-基础绘图
Ausgelebt的博客
05-14 855
选择哪种方法取决于具体需求，例如想要创建的子图的数量、位置、样式以及是否需要与其他MATLAB功能（如工具栏、菜单等）集成等等。函数允许更精细地控制子图的位置和大小，可以通过指定位置向量和轴属性来定义轴的范围和样式。函数用于在单个图形窗口中创建多个子图。你可以指定子图的网格位置和当前激活的子图。函数提供了一种更现代的方式来创建平铺的图形布局。然后使用这些网格数据和一个函数来计算Z的值。允许定义网格布局，并在这些网格中放置轴。然后，使用不同的函数来绘制三维图形。首先创建X和Y的网格数据。
matlab数值模拟,西安交通大学——温度场数值模拟(matlab
06-01
好的，关于温度场数值模拟，我可以给你一些基本的思路和步骤，希望对你有所帮助。 首先，我们需要确定模拟的物理场、边界条件和初始条件。在这个问题中，我们需要模拟的是温度场，边界条件可能包括材料的热导率、热容、边界的温度、热通量等信息，初始条件可能包括材料的初始温度分布。 接下来，我们需要选择合适的数值方法来求解温度场。常见的数值方法包括有限差分法、有限元法、边界元法等。在这里，我们可以选择有限差分法来进行数值模拟。 有限差分法的基本思路是将物理空间离散化，将连续的物理空间划分为若干个离散的空间点，然后在这些空间点上构建差分方程，通过求解差分方程得到温度场的数值解。需要注意的是，差分方程的形式和边界条件、初始条件有关，需要根据具体情况进行选择和构建。 最后，我们可以使用matlab等数值计算软件来实现温度场的数值模拟。具体步骤包括定义物理空间和离散化网格、构建差分方程、设置边界条件和初始条件、求解差分方程得到数值解、可视化温度场等。 希望以上内容对你有所帮助，如果你有其他问题或需要更详细的信息，可以继续询问我。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值