常见数学函数.doc

2.1 基本数学函数
2.1.1 三角函数与双曲函数
函数 sin、sinh
功能 正弦函数与双曲正弦函数
格式 Y = sin(X) %计算参量X（可以是向量、矩阵，元素可以是复数）中每一个角度分量的正弦值Y，所有分量的角度单位为弧度。
Y = sinh(X) %计算参量X的双曲正弦值Y
注意：sin(pi)并不是零，而是与浮点精度有关的无穷小量eps，因为pi仅仅是精确值π浮点近似的表示值而已；对于复数Z= x+iy，函数的定义为：sin(x+iy) = sin(x)*cos(y) + i*cos(x)*sin(y)，，
例2-1
x = -pi:0.01:pi; plot(x,sin(x))
x = -5:0.01:5; plot(x,sinh(x))
图形结果为图2-1。

图2-1 正弦函数与双曲正弦函数图
函数 asin、asinh
功能 反正弦函数与反双曲正弦函数
格式 Y = asin(X) %返回参量X（可以是向量、矩阵）中每一个元素的反正弦函数值Y。若X中有的分量处于[-1,1]之间，则Y = asin(X)对应的分量处于[-π/2,π/2]之间，若X中有分量在区间[-1,1]之外，则Y= asin(X)对应的分量为复数。
Y = asinh(X) %返回参量X中每一个元素的反双曲正弦函数值Y
说明 反正弦函数与反双曲正弦函数的定义为：，
例2-2
x = -1:.01:1; plot(x,asin(x))
x = -5:.01:5; plot(x,asinh(x))
图形结果为图2-2。

图2-2 反正弦函数与反双曲正弦函数图
函数 cos、cosh
功能 余弦函数与双曲余弦函数
格式 Y = cos(X) %计算参量X（可以是向量、矩阵，元素可以是复数）中每一个角度分量的余弦值Y，所有角度分量的单位为弧度。我们要指出的是，cos(pi/2)并不是精确的零，而是与浮点精度有关的无穷小量eps，因为pi仅仅是精确值π浮点近似的表示值而已。
Y = sinh(X) %计算参量X的双曲余弦值Y
说明 若X为复数z= x+iy，则函数定义为：cos(x+iy) = cos(x)*cos(y) + i*sin(x)*sin(y)，，
例2-3
x = -pi:0.01:pi; plot(x,cos(x))
x = -5:0.01:5; plot(x,cosh(x))
图形结果为图2-3。

图2-3 余弦函数与双曲余弦函数图
函数 acos、acosh
功能 反余弦函数与反双曲余弦函数
格式 Y = acos(X) %返回参量X（可以是向量、矩阵）中每一个元素的反余弦函数值Y。若X中有的分量处于[-1,1]之间，则Y = acos(X)对应的分量处于[0,π]之间，若X中有分量在区间[-1,1]之外，则Y = acos(X)对应的分量为复数。
Y = asinh(X) %返回参量X中每一个元素的反双曲余弦函数Y
说明 反余弦函数与反双曲余弦函数定义为：，
例2-4
x = -1:.01:1; plot(x,acos(x))
x = -5:.01:5; plot(x,acosh(x))
图形结果为图2-4。

图2-4 反余弦函数与反双曲余弦函数图
函数 tan、tanh
功能 正切函数与双曲正切函数
格式 Y = tan(X) %计算参量X（可以是向量、矩阵，元素可以是复数）中每一个角度分量的正切值Y，所有角度分量的单位为弧度。我们要指出的是，tan(pi/2)并不是精确的零，而是与浮点精度有关的无穷小量eps，因为pi仅仅是精确值π浮点近似的表示值而已。
Y = tanh(X) %返回参量X中每一个元素的双曲正切函数值Y
例2-5
x = (-pi/2)+0.01:0.01:(pi/2)-0.01; % 稍微缩小定义域
plot(x,tan(x))
x = -5:0.01:5; plot(x,tanh(x))
图形结果为图2-5。

图2-5 正切函数与双曲正切函数图
函数 atan、atanh
功能 反正切函数与反双曲正切函数
格式 Y = atan(X) %返回参量X（可以是向量、矩阵）中每一个元素的反正切函数值Y。若X中有的分量为实数，则Y = atan(X)对应的分量处于[-π/2,π/2]之间。
Y = atanh(X) %返回参量X中每一个元素的反双曲正切函数值Y。
说明 反正切函数与反双曲正切函数定义为：，
例2-6
x = -20:0.01:20; plot(x,atan(x))
x = -0.99:0.01:0.99; plot(x,atanh(x))
图形结果为图2-6。

图2-6 反正切函数与反双曲正切函数图
函数 cot、coth
功能 余切函数与双曲余切函数
格式 Y = cot(X) %计算参量X（可以是向量、矩阵，元素可以是复数）中每一个角度分量的余切值Y，所有角度分量的单位为弧度。
Y = coth(X) %返回参量X中每一个元素的双曲余切函数值Y
例2-7
x1 = -pi+0.01:0.01:-0.01; % 去掉奇点x = 0
x2 = 0.01:0.01:pi-0.01; % 做法同上
plot(x1,cot(x1),x2,cot(x2))
plot(x1,coth(x1),x2,coth(x2))
图形结果为图2-7。

图2-7 余切函数与双曲余切函数图
函数 acot、acoth
功能 反余切函数与反双曲余切函数
格式 Y = acot(X) %返回参量X（可以是向量、矩阵）中每一个元素的反余切函数Y
Y = acoth(X) %返回参量X中每一个元素的反双曲余切函数值Y
例2-8
x1 = -2*pi:pi/30:-0.1; x2 = 0.1:pi/30:2*pi; % 去掉奇异点x = 0
plot(x1,acot(x1),x2,acot(x2))
x1 = -30:0.1:-1.1; x2 = 1.1:0.1:30;
plot(x1,acoth(x1),x2,acoth(x2))
图形结果为图2-8。

图2-8 反余切函数与反双曲余切函数图
函数 sec、sech
功能 正割函数与双曲正割函数
格式 Y = sec(X) %计算参量X（可以是向量、矩阵，元素可以是复数）中每一个角度分量的正割函数值Y，所有角度分量的单位为弧度。我们要指出的是，sec(pi/2)并不是无穷大，而是与浮点精度有关的无穷小量eps的倒数，因为pi仅仅是精确值π浮点近似的表示值而已。
Y = sech(X) %返回参量X中每一个元素的双曲正割函数值Y
例2-9
x1 = -pi/2+0.01:0.01:pi/2-0.01; % 去掉奇异点x = pi/2
x2 = pi/2+0.01:0.01:(3*pi/2)-0.01;
plot(x1,sec(x1),x2,sec(x2))
x = -2*pi:0.01:2*pi;
plot(x,sech(x))
图形结果为图2-9。

图2-9 正割函数与双曲正割函数图
函数 asec、asech
功能 反正割函数与反双曲正割函数
格式 Y = asec(X) %返回参量X（可以是向量、矩阵）中每一个元素的反正割函数值Y
Y = asech(X) %返回参量X中每一个元素的反双曲正割函数值Y
例2-10
x1 = 5:0.01:-1; x2 = 1:0.01:5;
plot(x1,asec(x1),x2,asec(x2))
x = 0.01:0.001:1; plot(x,asech(x))
图形结果为图2-10。

图2-10 反正割函数与反双曲正割函数图
函数 csc、csch
功能 余割函数与双曲余割函数
格式 Y = csc(X) %计算参量X（可以是向量、矩阵，元素可以是复数）中每一个角度分量的余割函数值Y，所有角度分量的单位为弧度。
Y = csch(X) %返回参量X中每一个元素的双曲余割函数值Y
例2-11
x1 = -pi+0.01:0.01:-0.01; x2 = 0.01:0.01:pi-0.01; % 去掉奇异点x=0
plot(x1,csc(x1),x2,csc(x2))
plot(x1,csch(x1),x2,csch(x2))
图形结果为图2-11。

图2-11 余割函数与双曲余割函数图
函数 acsc、acsch
功能 反余割函数与反双曲余割函数。
格式 Y = asec(X) %返回参量X（可以是向量、矩阵）中每一个元素的反余割函数值Y
Y = asech(X) %返回参量X中每一个元素的反双曲余割函数值Y
例2-12
x1 = -10:0.01:-101; x2 = 1.01:0.01:10; % 去掉奇异点x = 1
plot(x1,acsc(x1),x2,acsc(x2))
x1 = -20:0.01:-1; x2 = 1:0.01:20;
plot(x1,acsch(x1),x2,acsch(x2))
图形结果为图2-12。

图2-12 反余割函数与反双曲余割函数图
函数 atan2
功能 四象限的反正切函数
格式 P = atan2(Y,X) %返回一与参量X和Y同型的、与X和Y元素的实数部分对应的、元素对元素的四象限的反正切函数阵列P，其中X和Y的虚数部分将忽略。阵列P中的元素分布在闭区间[-pi,pi]上。特定的象限将取决于sign(Y)与sign(X)。

例2-13
z=1+2i;
r = abs(z);
theta = atan2(imag(z),real(z))
z = r *exp(i *theta)
feather(z);hold on
t=0:0.1:2*pi;
x=1+sqrt(5)*cos(t);
y=sqrt(5)*sin(t);
plot(x,y);
axis equal; hold off
计算结果为：
theta =
1.1071
z =
1.0000 + 2.0000i
图形结果为图2-13。

图213 四象限的反正切函数图
2.1.2 其他常用函数
函数 fix
功能 朝零方向取整
格式 B = fix(A) %对A的每一个元素朝零的方向取整数部分，返回与A同维的数组。对于复数参量A，则返回一复数，其分量的实数与虚数部分分别取原复数的、朝零方向的整数部分。
例2-14
>>A = [-1.9, -0.2, 3.1415926, 5.6, 7.0, 2.4+3.6i];
>>B = fix(A)
计算结果为：
B =
Columns 1 through 4
-1.0000 0 3.0000 5.0000
Columns 5 through 6
7.0000 2.0000 + 3.0000i
函数 roud
功能 朝最近的方向取整。
格式 Y = round(X) %对X的每一个元素朝最近的方向取整数部分，返回与X同维的数组。对于复数参量X，则返回一复数，其分量的实数与虚数部分分别取原复数的、朝最近方向的整数部分。
例2-15
>>A = [-1.9, -0.2, 3.1415926, 5.6, 7.0, 2.4+3.6i];
>>Y = round(A)
计算结果为：
Y =
Columns 1 through 4
-2.0000 0 3.0000 6.0000
Columns 5 through 6
7.0000 2.0000 + 4.0000i
函数 floor
功能 朝负无穷大方向取整
格式 B = floor(A) %对A的每一个元素朝负无穷大的方向取整数部分，返回与A同维的数组。对于复数参量A，则返回一复数，其分量的实数与虚数部分分别取原复数的、朝负无穷大方向的整数部分。
例216
>>A = [-1.9, -0.2, 3.1415926, 5.6, 7.0, 2.4+3.6i];
>>F = floor(A)
计算结果为：
F =
Columns 1 through 4
-2.0000 -1.0000 3.0000 5.0000
Columns 5 through 6
7.0000 2.0000 + 3.0000i
函数 rem
功能 求作除法后的剩余数
格式 R = rem(X,Y) %返回结果X - fix(X./Y).*Y，其中X、Y应为正数。若X、Y为浮点数，由于计算机对浮点数的表示的不精确性，则结果将可能是不可意料的。fix(X./Y)为商数X./Y朝零方向取的整数部分。若X与Y为同符号的，则rem(X,Y)返回的结果与mod(X,Y)相同，不然，若X为正数，则rem(-X,Y) = mod(-X,Y) - Y。该命令返回的结果在区间[0，sign(X)*abs(Y)]，若Y中有零分量，则相应地返回NaN。
例2-17
>>X = [12 23 34 45];
>>Y = [3 7 2 6];
>>R = rem(X,Y)
计算结果为：
R =
0 2 0 3
函数 ceil
功能 朝正无穷大方向取整
格式 B = floor(A) % 对A的每一个元素朝正无穷大的方向取整数部分，返回与A同维的数组。对于复数参量A，则返回一复数，其分量的实数与虚数部分分别取原复数的、朝正无穷大方向的整数部分。
例2-18
>>A = [-1.9, -0.2, 3.1415926, 5.6, 7.0, 2.4+3.6i];
>>B = ceil(A)
计算结果为：
B =
Columns 1 through 4
-1.0000 0 4.0000 6.0000
Columns 5 through 6
7.0000 3.0000 + 4.0000i
函数 exp
功能 以e为底数的指数函数
格式 Y = exp(X) %对参量X的每一分量，求以e为底数的指数函数Y。X中的分量可以为复数。对于复数分量如，z = x +i*y，则相应地计算：e^z = e^x*(cos(y) + i*sin(y))。
例2-19
>>A = [-1.9, -0.2, 3.1415926, 5.6, 7.0, 2.4+3.6i];
>>Y = exp(A)
计算结果为：
Y =

1.0e+003 *
Columns 1 through 4
0.0001 0.0008 0.0231 0.2704
Columns 5 through 6
1.0966 -0.0099 - 0.0049i
函数 expm
功能 求矩阵的以e为底数的指数函数
格式 Y = expm(X) %计算以e为底数、x的每一个元素为指数的指数函数值。若矩阵x有小于等于零的特征值，则返回复数的结果。
说明 该函数为一内建函数，它有三种计算算法：
（1）使用文件expm1.m中的用比例法与二次幂算法得到的Pad近似值；
（2）使用Taylor级数近似展开式计算，这种计算在文件expm2.m中。但这种一般计算方法是不可取的，通常计算是缓慢且不精确的；
（3）在文件expm3.m中，先是将矩阵对角线化，再把函数计算出相应的的特征向量，最后转换过来。但当输入的矩阵没有与矩阵阶数相同的特征向量个数时，就会出现错误。
例2-20
>>A=hilb(4);
>>Y = expm(A)
计算结果为：
Y =
3.2506 1.2068 0.8355 0.6417
1.2068 1.7403 0.5417 0.4288
0.8355 0.5417 1.4100 0.3318
0.6417 0.4288 0.3318 1.2729
函数 log
功能 自然对数，即以e为底数的对数。
格式 Y = log(X) %对参量X中的每一个元素计算自然对数。其中X中的元素可以是复数与负数，但由此可能得到意想不到的结果。若z = x + i*y，则log对复数的计算如下：log (z) = log (abs (z)) + i*atan2(y,x)
例2-21 下面的语句可以得到无理数π的近似值：
>>Pi = abs(log(-1))
计算结果为：
Pi =
3.1416
函数 log10
功能 常用对数，即以10为底数的对数。
格式 Y = log10(X) %计算X中的每一个元素的常用对数，若X中出现复数，则可能得到意想不到的结果。
例2-22
>>L1 = log10(realmax) % 由此可得特殊变量realmax的近似值
>>L2 = log10(eps) % 由此可得特殊变量eps的近似值
>>M = magic(4);
>>L3 = log10(M)
计算结果为：
L1 =
308.2547
L2 =
15.6536
L3 =
1.2041 0.3010 0.4771 1.1139
0.6990 1.0414 1.0000 0.9031
0.9542 0.8451 0.7782 1.0792
0.6021 1.1461 1.1761 0
函数 sort
功能 把输入参量中的元素按从小到大的方向重新排列
格式 B = sort(A) %沿着输入参量A的不同维的方向、从小到大重新排列A中的元素。A可以是字符串的、实数的、复数的单元数组。对于A中完全相同的元素，则按它们在A中的先后位置排列在一块；若A为复数的，则按元素幅值的从小到大排列，若有幅值相同的复数元素，则再按它们在区间[-π,π]的幅角从小到大排列；若A中有元素为NaN，则将它们排到最后。若A为向量，则返回从小到大的向量，若A为二维矩阵，则按列的方向进行排列；若A为多维数组，sort(A)把沿着第一非单元集的元素象向量一样进行处理。
B = sort(A,dim) %沿着矩阵A（向量的、矩阵的或多维的）中指定维数dim方向重新排列A中的元素。
[B,INDEX] = sort(A,…) %输出参量B的结果如同上面的情形，输出INDEX是一等于size(A)的数组，它的每一列是与A中列向量的元素相对应的置换向量。若A中有重复出现的相同的值，则返回保存原来相对位置的索引。
例2-23
>>A = [-1.9, -0.2, 3.1415926, 5.6, 7.0, 2.4+3.6i];
>>[B1,INDEX] = sort(A)
>>M = magic(4);
>>B2 = sort(M)
计算结果为：
B1 =
Columns 1 through 4
-0.2000 -1.9000 3.1416 2.4000 + 3.6000i
Columns 5 through 6
5.6000 7.0000
INDEX =
2 1 3 6 4 5
B2 =
4 2 3 1
5 7 6 8
9 11 10 12
16 14 15 13
函数 abs
功能 数值的绝对值与复数的幅值
格式 Y = abs(X) %返回参量X的每一个分量的绝对值；若X为复数的，则返回每一分量的幅值：abs(X) = sqrt(real(X).^2+imag(X).^2)。
例2-24
>>A = [-1.9, -0.2, 3.1415926, 5.6, 7.0, 2.4+3.6i];
>>Y = abs(A)
计算结果为：
Y =
1.9000 0.2000 3.1416 5.6000 7.0000 4.3267
函数 conj
功能 复数的共轭值
格式 ZC = conj(Z) %返回参量Z的每一个分量的共轭复数：
conj(Z) = real(Z) - i*imag(Z)

函数 imag
功能 复数的虚数部分
格式 Y = imag(Z) %返回输入参量Z的每一个分量的虚数部分。
例2-25
>>imag(2+3i)
计算结果为：
ans =
3
函数 real
功能 复数的实数部分。
格式 Y = real(Z) %返回输入参量Z的每一个分量的实数部分。
例2-26
>>real(2+3i)
计算结果为：
ans =
2
函数 angle
功能 复数的相角
格式 P = angle(Z) %返回输入参量Z的每一复数元素的、单位为弧度的相角，其值在区间[-π，π]上。

说明 angle(z) = imag (log(z)) = atan2 (imag(z),real(z))
例2-27
>>Z =[1-i, 2+i, 3-i, 4+i;
>>1+2i,2-2i,3+2i,4-2i;
>>1-3i,2+3i,3-3i,4+3i;
>>1+4i,2-4i,3+4i,4-4i];
>>P = angle(Z)
计算结果为：
P =
-0.7854 0.4636 -0.3218 0.2450
1.1071 -0.7854 0.5880 -0.4636
-1.2490 0.9828 -0.7854 0.6435
1.3258 -1.1071 0.9273 -0.7854
函数 complex
功能 用实数与虚数部分创建复数
格式 c = complex(a,b) %用两个实数a，b创建复数c=a+bi。输出参量c与a、b同型（同为向量、矩阵、或多维阵列）。该命令比下列形式的复数输入更有用：a + i*b 或a + j*b因为i和j可能被用做其他的变量(不等于sqrt(-1))，或者a和b不是双精度的。
c = complex(a) %输入参量a作为输出复数c的实部，其虚部为0：c = a+0*i。
例2-28
>>a = uint8([1;2;3;4]);
>>b = uint8([4;3;2;1]);
>>c = complex(a,b)
计算结果为：
c =
1.0000 + 4.0000i
2.0000 + 3.0000i
3.0000 + 2.0000i
4.0000 + 1.0000i
函数 mod
功能 模数（带符号的除法余数）
用法 M = mod(X,Y) %输入参量X、Y应为整数，此时返回余数X -Y.*floor(X./Y)，若Y≠0，或者是X。若运算数x与y有相同的符号，则mod(X,Y)等于rem(X,Y)。总之，对于整数x,y，有：mod(-x,y) = rem(-x,y)+y。若输入为实数或复数，由于浮点数在计算机上的不精确表示，该操作将导致不可预测的结果。
例2-29
>>M1 = mod(13,5)
>>M2 = mod([1:5],3)
>>M3 = mod(magic(3),3)
计算结果为：
M1 =
3
M2 =
1 2 0 1 2
M3 =
2 1 0
0 2 1
1 0 2
函数 nchoosek
功能 二项式系数或所有的组合数。该命令只有对n<15时有用。
函数 C = nchoosek(n,k) %参量n,k为非负整数，返回n! / ( (n-k)! k!)，即一次从n个物体中取出k个的组合数。
C = nchoosek(v,k) %参量v为n维向量，返回一矩阵，其行向量的分量为一次性从v个物体中取k个物体的组合数。矩阵 C包含=n! / ( (n-k)! k!)行与k列。
例2-30
>>C = nchoosek(2:2:10,4)
计算结果为：
C =
2 4 6 8
2 4 6 10
2 4 8 10
2 6 8 10
4 6 8 10
函数 rand
功能 生成元素均匀分布于(0,1)上的数值与阵列
用法 Y = rand(n) %返回n*n阶的方阵Y，其元素均匀分布于区间(0,1)。若n不是一标量，在显示一出错信息。
Y = rand(m,n)、Y = rand([m n]) %返回阶数为m*n的，元素均匀分布于区间(0,1)上矩阵Y。
Y = rand(m,n,p,…)、Y = rand([m n p…]) %生成阶数m*n*p*…的，元素服从均匀分布的多维随机阵列Y。
Y = rand(size(A)) %生成一与阵列A同型的随机均匀阵列Y
rand %该命令在每次单独使用时，都返回一随机数（服从均匀分布）。
s = rand('state') %返回一有35元素的列向量s，其中包含均匀分布生成器的当前状态。该改变生成器的当前的状态，见表2-1。
表2-1
命 令
含 义
Rand(’state’,s)
设置状态为s
Rand(’state’,0)
设置生成器为初始状态
Rand(’state’,k)
设置生成器第k个状态(k为整数)
Rand(’state’,sum(100*clock))
设置生成器在每次使用时的状态都不同（因为clock每次都不同）
例：
>>R1 = rand(4,5)
>>a = 10; b = 50;
>>R2 = a + (b-a) * rand(5) % 生成元素均匀分布于(10,50)上的矩阵
计算结果可能为：
R1 =
0.6655 0.0563 0.2656 0.5371 0.6797
0.3278 0.4402 0.9293 0.5457 0.6129
0.6325 0.4412 0.9343 0.9394 0.3940
0.5395 0.6501 0.5648 0.7084 0.2206
R2 =
33.6835 19.8216 36.9436 49.6289 46.4679
18.5164 34.2597 15.3663 31.0549 49.0377
19.0026 37.1006 33.6046 39.5361 13.9336
12.4641 12.9804 35.5420 23.2916 46.8304
28.5238 48.7418 49.0843 13.0512 10.9265
函数 randn
功能 生成元素服从正态分布（N(0,1)）的数值与阵列
格式 Y = randn(n) %返回n*n阶的方阵Y，其元素服从正态分布N(0,1)。若n不是一标量，则显示一出错信息。
Y = randn(m,n)、Y = randn([m n]) %返回阶数为m*n的，元素均匀分布于区间(0,1)上矩阵Y。
Y = randn(m,n,p,…)、Y = randn([m n p…]) %生成阶数m*n*p*…的，元素服从正态分布的多维随机阵列Y。
Y = randn(size(A)) %生成一与阵列A同型的随机正态阵列Y
randn %该命令在每次单独使用时，都返回一随机数（服从正态分布）。
s = randn('state') %返回一有2元素的向量s，其中包含正态分布生成器的当前状态。该改变生成器的当前状态，见表2-2。
表2-2
命 令
含 义
randn(’state’,s)
设置状态为s
randn(’state’,0)
设置生成器为初始状态
rand(’state’,k)
设置生成器第k个状态(k为整数)
rand(’state’,sum(100*clock))
设置生成器在每次使用时的状态都不同（因为clock每次都不同）
例：
>>R1 = rand(4,5)
>>R2 = 0.6 + sqrt(0.1) * randn(5)
计算结果可能为：
R1 =
0.2778 0.2681 0.5552 0.5167 0.8821
0.2745 0.3710 0.1916 0.3385 0.5823
0.9124 0.5129 0.4164 0.2993 0.0550
0.4125 0.2697 0.1508 0.9370 0.5878
R2 =
0.4632 0.9766 0.5410 0.6360 0.6931
0.0733 0.9760 0.8295 0.9373 0.1775
0.6396 0.5881 0.4140 0.6187 0.8259
0.6910 0.7035 1.2904 0.5698 1.1134
0.2375 0.6552 0.5569 0.3368 0.3812