《数学软件Matlab—— 矩阵运算》由会员分享，可在线阅读，更多相关《数学软件Matlab—— 矩阵运算（18页珍藏版）》请在金锄头文库上搜索。
1、数学软件数学软件 Matlab 矩阵运算1本讲主要内容本讲主要内容n Kronecker 乘积乘积n 数组运算数组运算n 函数的使用函数的使用n 变量的输入输出变量的输入输出 2Kronecker 乘乘积积l 矩阵矩阵 Kronecker 乘积乘积的定义的定义l 设设 A 是是 nm 矩阵，矩阵，B 是是 pq 矩阵，则矩阵，则 A 与与 B 的的 kronecker 乘积为：乘积为：l 是是 npmq 矩阵矩阵l 任意两个矩阵都有任意两个矩阵都有 Kronecker 乘积乘积l 通常通常l Matlab 中计算两个矩阵的中计算两个矩阵的 Kronecker 乘积乘积C=kron(A,B);3
2、本讲主要内容本讲主要内容n Kronecker 乘积乘积n 数组运算数组运算n 函数的使用函数的使用n 变量的输入输出变量的输入输出 4矩阵的数组运算矩阵的数组运算l 数组运算：数组运算：对应元素进行运算对应元素进行运算l 数组运算包括：数组运算包括：点乘、点除、点幂点乘、点除、点幂l 相应的四个数组运算符为：相应的四个数组运算符为： .* ./ . . 点与算术运算符之间不能有空格点与算术运算符之间不能有空格！参与运算的对象必须具有相同的形状参与运算的对象必须具有相同的形状！A=1,2,3; 4,5,6; B=3,2,1; 6,5,4;C=A.*B D=A./BE=A.B例例：5本讲主要内容
3、本讲主要内容n Kronecker 乘积乘积n 数组运算数组运算n 函数的使用函数的使用n 变量的输入输出变量的输入输出 6函数取值函数取值设设 x 是变量，是变量， f 是一个函数是一个函数l 当当 x = a 是标量时，是标量时，f(x) = f(a) 也是一个标量也是一个标量l 当当 x = x1, x2, , xn 是向量时，则是向量时，则 f(x) = f(x1), f(x2), , f(xn) 是一个与是一个与 x 长度相同的向量长度相同的向量l 函数作用在矩阵上的取值函数作用在矩阵上的取值l 若若 A 是矩阵，则是矩阵，则 f (A) 是一个与是一个与 A 同形状的矩阵同形状的矩
4、阵f 作用在作用在 x 的的每个分量每个分量上上！x=0:pi/4:pi; A=1,2,3; 4,5,6;y1=sin(x) y2=exp(A)y3=sqrt(A)例例：7矩阵的超越函数矩阵的超越函数l Matlab 提供的矩阵函数提供的矩阵函数：expm、sqrtm、logm详情参见联机帮助详情参见联机帮助（help expm / sqrtm / logm ）l 更一般的矩阵函数：更一般的矩阵函数： funmfunm(A, fun)fun 可以是可以是 exp, log，cos，sin，. .怎样计算怎样计算 eA ?l 矩阵函数矩阵函数l 将矩阵作为一个整体参与运算将矩阵作为一个整体参与运
5、算8矩阵与数的运算矩阵与数的运算l 加减加减：矩阵的每个元素都与数作加减运算矩阵的每个元素都与数作加减运算l 数乘：数乘：矩阵的每个元素都与数作乘法运算矩阵的每个元素都与数作乘法运算l 矩阵除以一个数：矩阵除以一个数：每个元素都除以这个数每个元素都除以这个数l 数与矩阵的点幂运算数与矩阵的点幂运算：采用数组运算：采用数组运算l 矩阵与数的运算矩阵与数的运算例例：x=1 2 3; x.2=12,22,32=1,4,92.x=21,22,23=2,4,8Matlab中所有标点符号必须在英文状态下输入中所有标点符号必须在英文状态下输入!9常用数学函数常用数学函数sin(x)、cos(x)、tan(x
6、)、cot(x)、sec(x)、csc(x)、sinh(x)、cosh(x)、tanh(x) .asin(x)、acos(x)、atan(x)、acot(x)、asec(x)、acsc(x) .exp(x) % ex 自然指数自然指数 (以以 e 为底为底)pow2(x) % 2x 以以 2 为底的指数为底的指数 log(x) % ln(x) 自然对数自然对数 (以以 e 为底为底)log2(x) % 以以 2 为底的对数为底的对数log10(x) % 以以 10 为底的对数为底的对数sqrt(x) % 平方根平方根abs(x) % 绝对值绝对值10常用数学函数常用数学函数conj(z) %
7、复数的共轭复数的共轭 real(z) % 复数的实部复数的实部 imag(z) % 复数的虚部复数的虚部angle(z) % 复数的辐角复数的辐角sign(x) % 符号函数符号函数round(x) % 取整函数取整函数 (四舍五入四舍五入)fix(x)、floor(x)、ceil(x) % 另外三个取整函数另外三个取整函数 mod(x,y) % 计算计算 x 除以除以 y 的余数（结果与的余数（结果与 y 同号）同号）rem(x,y) % 计算计算 x 除以除以 y 的余数（结果与的余数（结果与 x 同号）同号）gcd(x,y) % 最大公约数最大公约数lcm(x,y) % 最小公倍数最小公
8、倍数11常用数学函数常用数学函数linspace(a,b,n) % 生成生成 a 到到 b，n 个数的等差数列个数的等差数列logspace(a,b,n) % 生成生成 10a 到到 10b，n 个数的等比数列个数的等比数列max(x) % 求求 向量向量 x 中的最大值中的最大值min(x) % 求最小值求最小值mean(x) % 求平均值求平均值sum(x) % 求和求和sort(x) % 排序排序det(A) % 矩阵行列式矩阵行列式inv(A) % 矩阵的逆矩阵的逆eig(A) % 矩阵的特征值矩阵的特征值rank(A) % 矩阵的秩矩阵的秩若参数若参数 x 是矩阵，是矩阵，则作用在其
9、各列上则作用在其各列上更多数学函数，参见课程主页更多数学函数，参见课程主页12本讲主要内容本讲主要内容n Kronecker 乘积乘积n 数组运算数组运算n 函数的使用函数的使用n 变量的输入输出变量的输入输出 13l 输出格式输出格式Matlab 的输出格式的输出格式l Matlab 以双精度执行所有的运算，运算结果可以以双精度执行所有的运算，运算结果可以在屏幕在屏幕上输出上输出，同时，同时赋给指定变量赋给指定变量；若无指定变量，则系统会自动；若无指定变量，则系统会自动将结果赋给变量将结果赋给变量 “ans” l 在屏幕上输出运算结果时，可以在屏幕上输出运算结果时，可以通过通过 format
10、 命令指命令指定输出格式定输出格式format 只只改变改变输出格式，不会改变变量的值输出格式，不会改变变量的值！14各种各种 format格式格式解释解释例例format短格式（缺省格式），同短格式（缺省格式），同 short3.1416format short短格式（缺省格式），小数点后短格式（缺省格式），小数点后 4 位位3.1416format long长格式，双精度长格式，双精度 15位，单精度位，单精度 7位位3.14159265358979format shortE短短 E 格式（科学计数格式格式（科学计数格式）3.1416e+000format longE长长 E 格式格式3.
11、141592653589793e+000format shortG短短 G 格式格式3.1416format longG长长 G 格式格式3.14159265358979format rat分数形式分数形式355/113format compact紧凑格式紧凑格式format loose宽松格式宽松格式format +， format bank， format hex (详情查看联机帮助详情查看联机帮助)15变量的存储变量的存储l 将将变量变量保存到保存到 mat 文件中文件中save将将所有变量所有变量保存到文件保存到文件 matlab.mat 中中save fnamesave fname.
12、mat 将将所有变量所有变量保存到文件保存到文件 fname.mat 中中save fname 变量变量将将指定变量指定变量保存到文件保存到文件 fname.mat 中中l 可同时保存多个变量，各变量之间用可同时保存多个变量，各变量之间用空格空格隔开隔开例例：x=2:5; A=magic(3); a=2.4;save mydata A x;l 从从 mat 文件中读取文件中读取变量变量load fname读取读取 fname.mat 中的中的所有变量所有变量load fname 变量变量从从 fname.mat 中读取中读取指定的变量指定的变量16上机上机要求要求l将完成每题所用的命令写入一个
13、文件，取名为将完成每题所用的命令写入一个文件，取名为 m02.ml然后将该文件作为附件，通过然后将该文件作为附件，通过 foxmail 以邮件形式发给以邮件形式发给 mhjssystem.maill邮件主题为：邮件主题为：机号机号-学号学号-姓名姓名l其中机号为其中机号为 两位数两位数l三个字段之间用英文状态下的减号链接三个字段之间用英文状态下的减号链接q 上机要求上机要求171、自己练习，体会各种数学函数的使用、自己练习，体会各种数学函数的使用 (不用上交不用上交) 2、写出可以完成下列任务的、写出可以完成下列任务的Matlab命令命令 (需上交需上交) ：（1）生成列向量）生成列向量 x
14、= 10-10,10-9, ., 10-1, 1, 2, 4, 6, 8, , 80（2）数列数列 x 同（同（1），求数列），求数列 y 使得使得 yi=sin(xi)/xi（3）产生一个）产生一个由由 1 到到 10 组成组成的随机整数排列的随机整数排列 S1（参见（参见教材教材 P143）（4）产生一个）产生一个长度为长度为10的的序列序列 S2，每个元素均是，每个元素均是 1 到到 10 之间，且满足之间，且满足均均匀匀分布的分布的随机整数随机整数 （参见参见教材教材 144），问：），问：S1与与S2有什么区别有什么区别？（5）利用函数）利用函数 rand 生成一个含有满足均匀分布的
15、生成一个含有满足均匀分布的 5 阶随机整数方阵阶随机整数方阵A， 元素都是元素都是 0100 之间之间的整数的整数（6）生成一个与）生成一个与 A 同阶的满足正态分布的随机矩阵同阶的满足正态分布的随机矩阵 B（7）计算）计算 A 的转置的转置 与与 B 的下三角部分的的下三角部分的 kronecker 乘积乘积 C（8）生成由）生成由 B 与其同阶魔方矩阵点乘得到的矩阵与其同阶魔方矩阵点乘得到的矩阵 D（9）求出矩阵）求出矩阵 D 中绝对值最大的元素的绝对值中绝对值最大的元素的绝对值（10）删除）删除 D 的第三行和第二列的第三行和第二列（11）生成由）生成由 B 的第的第 2, 4, 5 行和第行和第 4, 1, 4 列组成的矩阵列组成的矩阵 E （12）教材第）教材第 53 页，第页，第 1 题（题（15）和第）和第 5 题题上机作业上机作业18
