矩阵指数

计算和比较 A 的指数和 A 的矩阵指数。

A = [1 1 0; 0 0 2; 0 0 -1];
exp(A)
ans = 3×32.7183    2.7183    1.00001.0000    1.0000    7.38911.0000    1.0000    0.3679expm(A)
ans = 3×32.7183    1.7183    1.08620    1.0000    1.26420         0    0.3679

请注意，两种结果的对角线元素相等，这对于任何三角矩阵都是如此。非对角线元素（包括对角线下方的元素）则不相同。

funm

计算常规矩阵函数

语法

F = funm(A,fun) F = funm(A,fun,options) F = funm(A,fun,options,p1,p2,...) [F,exitflag] = funm(...) [F,exitflag,output] = funm(...)

说明

F = funm(A,fun) 计算在方阵参数为 A 时用户定义的函数 fun。F = fun(x,k) 必须接受向量 x 和整数 k，返回大小相同的 x 的向量 f，其中 f(i) 是在 x(i) 条件下计算的函数 fun 的第 k 个导数。fun 表示的函数必须包含具有无限收敛半径的泰勒级数，被视为特殊情况的 fun = @log 除外。

您也可以使用 funm 计算下表中列出的特殊函数在 A 处的值。

函数	计算矩阵 A 处的函数的语法
`exp`	`funm(A, @exp)`
`log`	`funm(A, @log)`
`sin`	`funm(A, @sin)`
`cos`	`funm(A, @cos)`
`sinh`	`funm(A, @sinh)`
`cosh`	`funm(A, @cosh)`

对于方阵根，请改用 sqrtm(A)。对于矩阵指数，expm(A) 和 funm(A, @exp) 哪一个更准确取决于矩阵 A。

fun 表示的函数必须包含具有无限收敛半径的泰勒级数。例外是被视为特殊情况的 @log。参数化函数解释如何在必要情况下向函数 fun 提供其他参数。

F = funm(A,fun,options) 将算法的参数设置为结构体 options 中的值。

下表列出了 options 的字段。

字段	说明	值
`options.Display`	显示级别	`'off'`（默认值）、`'on'`、`'verbose'`
`options.TolBlk`	阻止 Schur 表的容差	正标量。默认值为 `0.1`。
`options.TolTay`	计算对角线块的泰勒级数的终止容差	正标量。默认值为 `eps`。
`options.MaxTerms`	泰勒级数项的最大数目	正整数。默认值为 `250`。
`options.MaxSqrt`	计算对数时，逆缩放和二乘法中计算的最大平方根数。	正整数。默认值为 `100`。
`options.Ord`	指定 Schur 表 `T` 的排序方式。	长度为 `length(A)` 的向量。`options.Ord(i)` 是 `T(i,i)` 所放置到的块的索引。默认值为 `[]`。

F = funm(A,fun,options,p1,p2,...) 向函数传递额外的输入 p1,p2,...。

[F,exitflag] = funm(...) 返回用于描述 funm 的退出条件的 exitflag。exitflag 可以具有下列值：

0 - 算法成功。
1 - 一次或多次泰勒级数计算未收敛，在使用对数的情况下，需要的平方根太多。但是，F 的计算值可能仍然正确。

[F,exitflag,output] = funm(...) 返回包含以下字段的结构体 output：

字段	说明
`output.terms`	一个向量，其中 `output.terms(i)` 是在计算第 i 个块时所使用的泰勒级数的项数，或者在使用对数的情况下，维度大于 2 的矩阵的平方根数。
`output.ind`	重新排序的 Schur 因子 `T` 的 `(i,j)` 块为 `T(output.ind{i}, output.ind{j})` 的元胞数组。
`output.ord`	传递到 `ordschur` 时对 Schur 表排序
`output.T`	重新排序的 Schur 表

如果 Schur 表为对角线，则 output = struct('terms',ones(n,1),'ind',{1:n})。

示例

示例 1

以下命令计算 3×3 幻方矩阵的矩阵正弦值。

F=funm(magic(3), @sin)F =-0.3850    1.0191    0.01620.6179    0.2168   -0.18440.4173   -0.5856    0.8185

示例 2

以下语句

S = funm(X,@sin);
C = funm(X,@cos);

在舍入误差内生成与下面相同的结果

E = expm(i*X);
C = real(E);
S = imag(E);

在任一情况下，结果都满足 S*S+C*C = I，其中 I = eye(size(X))。

示例 3

要使用一个对 funm 的调用计算函数 exp(x) + cos(x) 在 A 处的值，请使用

F = funm(A,@fun_expcos)

其中 fun_expcos 是以下函数。

function f = fun_expcos(x, k)
% Return kth derivative of exp + cos at X.g = mod(ceil(k/2),2);if mod(k,2)f = exp(x) + sin(x)*(-1)^g;elsef = exp(x) + cos(x)*(-1)^g;end

$A=\bigl(\begin{smallmatrix} 2&1 & -1\\ 1& 3 &-1 \\ 1& 2& 0 \end{smallmatrix}\bigr)$ 求 $_e{At},_e{A},sin{A}$

A=[2 1 -11 3 -11 2 0];
A = sym(A);
syms t;
expm(A.*t)
expm(A)
funm(A, @exp)

和答案一样

用MATLAB计算常规矩阵函数 eAt sinA相关推荐

P2 Matlab计算基础-《Matlab/Simulink与控制系统仿真》程序指令总结
上一篇回到目录下一篇 <Matlab/Simulink与控制系统仿真>程序指令总结 Matlab_Simulink_BookExample 2. Matlab 计算基础表2.1 Ma ...
2021-05-25 传递函数阶跃响应指标的matlab计算
传递函数阶跃响应指标的matlab计算对于一个常规的传递函数,我们在进行了阶跃仿真后经常需要给出对应的响应指标, sys = tf(1,[1 2 3]) y = step(sys) figure p ...
matlab计算一组数的函数值_matlab中怎么求一组数据概率最大的数
好象直接使用就可以了,第一次出现就包含了定义顺便把些命令弄给你看下 matlab命令 matlab commands and functions list A a abs 绝对值.模.字符的ASCI ...
matlab 计算 r,MATLAB R2015b*化计算
目录第1章MATLAB R2015b概述 1.1MATLAB简介 1.1.1MATLAB的发展史 1.1.2MATLAB的优点 1.1.3MATLAB系统组成 1.1.4MATLAB应用程序 1.1 ...
matlab画复变函数,科学网—复数复变函数的Matlab计算与绘图 - 周铁戈的博文
复数复变函数的Matlab计算与绘图周铁戈复数的表示存在两种表示方法,一种是代数式,一种是指数式,在Matlab中的方式如下: >> z=1+2i #代数式,1 ...
【数字信号处理】相关函数应用 ( 使用 matlab 计算相关函数 )
文章目录一.相关函数应用场景 1.生成高斯白噪声 2.信噪比 SNR 3.根据信噪比 SNR 求信号幅度 4.产生单载波信号及最终信号 5.求自相关函数及功率 6.matlab 完整代码一.相关函 ...
Matlab计算基站覆盖面积示例
站间距用D表示:小区半径为 R: 单扇区覆盖面积一般作为一个正六边形:D和R示意如下: 计算基站覆盖面积,一般算三个正六边形:下图填充部分: 公式如下: D一般算作R的1.5倍: 用matlab计算 ...
fvdm 跟驰模型 matlab仿真_【新书推荐】【2009.01】控制系统MATLAB计算及仿真（第3版）（黄忠霖）...
本书包括MATLAB语言基础.控制系统MATLAB仿真基础.控制系统MATLAB计算与仿真3篇. 本书上篇介绍了MATLAB 7.1 系统的相关知识,包括MATLAB 7.1系统概述与活的笔记本Not ...
matlab 求复数的模长,matlab计算带有复数的函数，最后求复数函数的模，结果里面却有...
该楼层疑似违规已被系统折叠隐藏此楼查看此楼 matlab计算带有复数的函数,最后求复数函数的模,结果里面却有复数标志i,输入自变量,得到的结果确实是实数,不知道怎么回事; 代码如下: syms n2 ...
matlab计算复活节概率,复活节日期的计算方法
复活节(主复活日)是一个西方的重要节日,在每年春分月圆之后第一个星期日.基督徒认为,复活节象征着重生与希望,为纪念耶稣基督于公元30到33年之间被钉死在十字架之后第三天复活的日子. 算法复活节是西方 ...

用MATLAB计算常规矩阵函数 eAt sinA