数学建模软件求解入门.ppt

A(:, 2) = [] % 删除第二列（：代表所有行） 例 A = 1 2 3 4 5 6 7 8 9 10 11 12 我们可以对矩阵进行各种处理： A(2,3) = 5 % 改变位於第二列，第三行的元素值 A = 1 2 3 4 5 6 5 8 9 10 11 12 B = A(2,1:3) % 取出部份矩阵B B = 5 6 5 A = [A B'] % 将B转置後以行向量并入A A = 1 2 3 4 5 5 6 5 8 6 9 10 11 12 5 第二十九页，共八十三页，2022年，8月28日 A = [A; 4 3 2 1] % 加入第四行 A = 1 3 4 5 5 5 8 6 9 11 12 5 4 3 2 1 A([1 4], :) = [] % 删除第一和第四行（：代表所有列） A = 5 5 8 6 9 11 12 5 第三十页，共八十三页，2022年，8月28日 数组操作函数 Repmat（D，1，3）命令：在水平方向“铺放”三个矩阵 D＝ 1 0 0 0 1 0 0 0 1 则Repmat（D，1，3）结果为 ans = 1 0 0 1 0 0 1 0 0 0 1 0 0 1 0 0 1 0 0 0 1 0 0 1 0 0 1 Diag(A,1)命令：对矩阵A取第一上对角元素 可运行diag(A,-1),查看运行结果 Triu(A)命令：取A的上三角阵；tril(A)命令：取A的下三角阵 第三十一页，共八十三页，2022年，8月28日 矩阵的逻辑运算 矩阵的比较运算：用＝＝表示等于关系，用>=表示大于等于关系，用~= 表示不等于关系等。 A=[1 2 3;4,5,6;7,8,0] Find(A>=5) 找出矩阵中元素大于等于5的下标 ans = 3 5 6 8 该函数先将A按列构成列向量，然后判断哪些元素 大于等于5，返回其下标。 >>[i,j]=find(A>=5);[i,j] 同时返回行和列坐标 >>all(A>=5) 判断A矩阵的某列元素是否全部大于等于5，是返回1，否则为0 >>any(A>=5) 你能猜测表示什么含义吗？ 第三十二页，共八十三页，2022年，8月28日 abs(x)：纯量的绝对值或向量的长度 angle(z)：复数z的相角(Phase angle) sqrt(x)：开平方 real(z)：复数z的实部 imag(z)：复数z的虚部 conj(z)：复数z的共轭复数 round(x)：四舍五入至最近整数 fix(x)：无论正负，舍去小数至最近整数 floor(x)：地板函数，即舍去正小数至最近整数 ceil(x)：天花板函数，即加入正小数至最近整数 rat(x)：将实数x化为分数表示 rats(x)：将实数x化为多项分数展开 sign(x)：符号函数 (Signum function)。 当x<0时，sign(x)=-1； 当x=0时，sign(x)=0; 当x>0时，sign(x)=1。 rem(x,y)：求x除以y的馀数 gcd(x,y)：整数x和y的最大公因数 lcm(x,y)：整数x和y的最小公倍数 MATLAB常用的基本数学函数 第三十三页，共八十三页，2022年，8月28日 第2章 绘图与图形处理 命令1 plot 功能 线性二维图。在线条多于一条时，若用户没有指定使用颜色，则plot循环使用由当前坐标轴颜色顺序属性（current axes ColorOrder property）定义的颜色，以区别不同的线条。在用完上述属性值后，plot又循环使用由坐标轴线型顺序属性（axes LineStyleOrder property）定义的线型，以区别不同的线条。 用法 plot(X,Y) 当X,Y均为实数向量，且为同维向量（可以不是同型向量），X=[x(i)]，Y=[y(i)]，则plot(X,Y)先描出点(x(i)，y(i))，然后用直线依次相连。 二维图形的绘制 1 散点图的绘制 画散点图一般使用命令scatter，其调用格式为： scatter(x,y) 2 二维曲线的绘制 第三十四页，共八十三页，2022年，8月28日 plot(X1,Y1,LineSpec1,X2,Y2,LineSpe