Control 版 (精华区)

发信人: baibai (哑琴弦), 信区: Control
标  题: [合集]如何在Matlab里建一个m函数计算e^At
发信站: 哈工大紫丁香 (2003年11月25日06:52:57 星期二), 站内信件


────────────────────────────────────────
 cattiger (画虎类猫)                  于 发信站: 紫 丁 香 (Wed Apr 26 20:03:29 2000) WWW-POST
【 在 wentao (涛) 的大作中提到: 】
:  本人在做毕业设计时遇到一个难题.要求在Matlab里建一个M函数计算
: e^At.不知那位朋友肯帮忙,不胜感激!

::y=exp(At);就可以了，若想了解更多的东西，可利用命令help象

?help exp

 EXP    Exponential.
    EXP(X) is the exponential of the elements of X, e to the X.
    For complex Z=X+i*Y, EXP(Z) = EXP(X)*(COS(Y)+i*SIN(Y)).
 
    See also LOG, LOG10, EXPM, EXPINT.

 Overloaded methods
    help sym/exp.m
    help demtseries/exp.m

────────────────────────────────────────
 matrix (火鸟)                        于 发信站: 紫 丁 香 (Wed Apr 26 20:16:33 2000), 转信
【 在 cattiger (画虎类猫) 的大作中提到: 】
【 在 wentao (涛) 的大作中提到: 】
:  本人在做毕业设计时遇到一个难题.要求在Matlab里建一个M函数计算
: e^At.不知那位朋友肯帮忙,不胜感激!
::y=exp(At);就可以了，若想了解更多的东西，可利用命令help象
   ~~~~~~~~~
    商量着说，好像不能这么算把？此式相当于对（At）的每一元素求
    指数运算，至少盗版MATLAB5.2是这样。
?help exp
 EXP    Exponential.
    EXP(X) is the exponential of the elements of X, e to the X.
    For complex Z=X+i*Y, EXP(Z) = EXP(X)*(COS(Y)+i*SIN(Y)).
    See also LOG, LOG10, EXPM, EXPINT.
 Overloaded methods
    help sym/exp.m
    help demtseries/exp.m

────────────────────────────────────────
 protel (tango)                       于 发信站: 紫 丁 香 (Thu Apr 27 19:25:44 2000), 转信
    MATLAB里有现成的矩阵指数函数，说明如下：
    expm       矩阵指数函数
    expm1      expm函数的.M文件实现
    expm2      Taylor级数法求矩阵指数
    expm3      特征值特征向量法求矩阵指数
    如果按矩阵指数函数的定义，应采用expm2。我
把它的程序拷贝如下：
E = zeros(size(A));
F = eye(size(A));
k = 1;
while norm(E+F-E,1) > 0
   E = E + F;
   F = A*F/k;
   k = k+1;
end
程序中norm（x，1）计算矩阵（或向量）x的一范数，norm（x，1）=
max（sum（abs（x）））。需要注意的是：
    设B=A*t，t取具体的数值时，才能计算expm（B）=expm（A*t）。
可以动脑筋编一个验证程序，一定会很有意思。这个快乐留给你了。
【 在 wentao (涛) 的大作中提到: 】
:  本人在做毕业设计时遇到一个难题.要求在Matlab里建一个M函数计算
: e^At.不知那位朋友肯帮忙,不胜感激!

────────────────────────────────────────