数值分析及其MATLAB实现第二篇MATLAB快速入门word21.docVIP

2.1 矩阵的输入与运算 MATLAB的主要数据对象是矩阵，标量、数组、行向量、列向量都是它的特例，最基本的功能是进行矩阵运算，MATLAB对于矩阵和数组的操作有一些特殊的规定。 矩阵输入有多种办法，如直接输入每个元；由语句或函数生成；在M-文件（以后介绍）中生成等。 2.1.1 矩阵的直接输入 MATLAB中不用描述矩阵的类型和维数，它们由输入的格式和内容决定。小矩阵可以用排列各个元的方法输入，同一行元用逗号或空格分开, 不同行的元用分号或回车分开。如, 在MATLAB工作窗口输入程序： 》A=[1,2,3;4,5,6] 或 》A=[1 2 3;4 5 6] 或 》A=[1 2 3 4 5 6] 都输入了一个2×3矩阵A，屏幕上显示输出变量为 A= 1 2 3 4 5 6 矩阵中的元可以用它的行、列数表示，例如, 在MATLAB工作窗口输入程序： 》a=A(2,1)%MATLAB区分大小写字母，a和A是不同的变量. 运行后输出结果如下: a= 4 或者不指定输出变量，MATLAB将回应ans（answer的缩写），如, 在MATLAB工作窗口输入程序： 》A(2,1) 运行后输出结果如下: ans= 4 A输入后一直保存在工作空间中，可随时调用，除非被清除或替代。 可以直接修改矩阵的元，如, 在MATLAB工作窗口输入程序： 》A(2,1)=7,A(3,4)=1, 运行后输出结果如下: A= 1 2 3 7 5 6 A= 1 2 3 0 7 5 6 0 0 0 0 1 原来的A没有3行4列，MATLAB自动增加行列数，对未输入的元赋值0。 2.1.2 函数生成的矩阵 MATLAB提供了一些函数来构造特殊矩阵，如, 在MATLAB工作窗口输入程序： 》w=zeros(2,3)% 2×3零矩阵. 运行后输出结果如下: w= 0 0 0 0 0 在MATLAB工作窗口输入程序： 》u=ones(3)%3×3全1矩阵，方阵只需输入行数，这几个矩阵均如此。 运行后输出结果如下: u= 1 1 1 1 1 1 1 1 1 在MATLAB工作窗口输入程序： 》v=eye(3,4)%3×4对角线为1的矩阵. 运行后输出结果如下: v= 1 0 0 0 0 1 0 0 0 0 1 0 在MATLAB工作窗口输入程序： 》x=rand(1,3)%1×3(0,1)均匀分布随机矩阵. 运行后输出结果如下: x= 0.9501 0.2311 0.6068 还有m×n标准正态分布矩阵randn(m,n)，n阶希尔伯特矩阵hilb(n)，n阶幻方矩阵magic(n)等。 2.1.3 矩阵的裁剪与拼接 从一个矩阵中取出若干行（列）构成新矩阵称为裁剪，MATLAB中“:”是非常重要的裁剪工具，如, 在MATLAB工作窗口输入程序： 》A(3,:)%A的第3行. 运行后输出结果如下: ans= 0 0 0 1 在MATLAB工作窗口输入程序： 》A(:,2)%A的第2列. 运行后输出结果如下: ans= 2 5 0 在MATLAB工作窗口输入程序： 》B=A(1:2,:)%A的第1,2行. 运行后输出结果如下: B= 1 2 3 0 7 5 6 0 在MATLAB工作窗口输入程序： 》C=B(:,2:4)%B的第2～4列. 运行后输出结果如下: C= 2 3 0 5 6 0 在MATLAB工作窗口输入程序： 》D=A(2:3,2:2:4)%A的第2,3行，第2,4列. 运行后输出结果如下: D= 5 0 0 1 请特别注意“:”的用法，想在MATLAB工作窗口输入程序： 》A(1:2:3,4:-1:2) 应该有什么结果。在MATLAB工作窗口输入程序： 》D(:,1)=[]%删除D的第1列,[]为空集符号. 运行后输出结果如下: D= 3 2 6 5 将几个矩阵接在一起称为拼接，左右拼接时行数要相同，上下拼接时列数要相同，如, 在MATLAB工作窗口输入程序： 》E=[D,zeros(2,1)] 运行后输出结果如下: E= 0 3 2 0 0 6 5 0 在MATLAB工作窗口输入程序： 》F=[A(1:2,:);eye(1,4)] 运行后输出结果如下: F= 1 2 3 0 7 5 6 0 1 1