预备实验 MATLAB使用练习.docVIP

预备实验 MATLAB使用练习 MATLAB ( MATrix LABoratory 矩阵实验室 特点： 编程效率高 计算功能强 特别是矩阵计算 使用简便 易于扩充 可建立M文件定义过程和函数以供调用，可与其它语言(C, Fortran等)混合编程。 §1 矩阵、数组与函数 矩阵的输入与运算 矩阵的直接输入 A=[1,2,3;4,5,6]( 或 A=[1 2 3;4 5 6]( 或 A=[1 2 3 4 5 6]( 显示矩阵元素： A(2,1)( ans=4 A(2,1)=7( A=1 2 3 7 5 6 A(3,4)=1( A= 1 2 3 0 7 5 6 0 0 0 0 1 函数生成矩阵 w=zeros(2,3)( u=ones(2,3)( v=eye(3,4)( 前三列为单位矩阵，后一列为零 x=rand(1,3)( 1X3 (0,1)均匀分布随机矩阵, randn(m,n) mXn标准正态分布矩阵， hilb(n), n阶Hilbert矩阵，magic(n)，n阶幻方矩阵 注：n阶Hilbert矩阵为，magic(3)= 矩阵的裁剪与拼接 A(3,:)( 显示A的第3行， A(:,2)( 显示A的第2列 A(2:5,:), A(:,1:2:5), A(2:4,1:3) C=(A,B), D=(A;B) 矩阵的运算 A+B，A-B，A, A*B，A^3, B/A(即BA-1), A\B (即A-1B), A+3(矩阵加常数，即每个元素加此常数) 数组及其运算 数组运算MATLAB中一种特定的计算，可视为两个同型矩阵间对应项的运算，运算符号为矩阵运算符号前加“.”。由于矩阵的加、减本身就是对应项间的运算，故不必再加“.”。 例如，x=t^2, y=sin(t), 0=t=pi. z=xy=t^2*sin(t)。在MATLAB中，往往用函数的一组值来表示此函数(如作图). t=0:0.1*pi:pi; x=t.^2; y=sin(t); z=x.*y (or z=t.^2.*sin(t)); 则t,x,y,z都是维数相同的向量，且有 x(k)=t(k)^2, y(k)=sin(t(k)), z(k)=x(k)*y(k)=t(k)^2*sin(t(k)). 1.3 语句、变量和表达式 1.3.1 语句形式 一般为 变量=表达式：a=[1 2 3 4]; b=a+1; c=a.*b; d=a*b; e=a*b; 1.3.2 变量 变量名由字母、数字和下划线组成，区分大小写字母，最多31个字符，第一个字符必须是字母. 不必说明类型和维数，自动分配内存空间. 几个特殊量：pi 圆周率, eps 可分辨最小正数, inf 正无穷大, NaN 不定值，i、j 复数虚单位. 以上几个量有其特殊值，但若赋予其它值后，则不再有其特殊性. 1.3.3 字符串 s1=Hello; s2=every; s3=body; s=[s1,,,s2, ,s3](，则 s= Hello,every body ss=s(1:5)(，则?ss=Hello. 1.4 函数 1.4.1 标量函数 三角函数：sin, cos, tan, cot, sec, csc, asin, acos, atan, acot, asec, acsc, sinh, cosh, tanh, asinh, acosh, atanh 其它基本函数：sqrt, exp, log, log10, abs, round(四舍五入取整), floor(向左方向取整), ceil(向右方向取整), fix(向0方向取整), sign, real(取实部), imag, angle, rats(有理逼近) rats(pi)(?ans=355/113, rats(sqrt(2))(?ans=1393/985 x=(0:0.2:1)*pi; y=sin(x)( y=0 0.5878 0.9511 0.9511 0.5878 0 round(y) 0 1 1 1 1 0; floor(y) 0 0 0 0 0 0; 另一个计算函数值的命令： feval(F,x), F:字符串，x:变量. 如上面的x：z=feval(sin,x), 同z=sin(x). 1.4.2 向量函数 max, min, sum, length, mean, median(中间值，奇数项取中间项，偶数项取中间两项的平均), prod(乘积), sort(从小到大重新排列) %以上