数学实验-实验指导书2[1].5.docVIP

实验一 MATLAB初步 实验目的 掌握矩阵运算和函数作图，掌握函数文件和命令文件的编写。 实验内容 矩阵运算、函数作图、编写函数文件、编写命令文件。 实验准备 一、向量的初始化 v=[a b c] 各元素之间用逗号或者空格分开。 v=first:increment:last first是向量的初始值，increment是步长。v从first出发，每隔increment步长取一个元素，直至不超过last的最大值。若省略increment，则步长默认为1。 v=linspace(first，last，n) 从first到last等间距的取n个元素。 x=[1 2 3 4 5 8 7 18] y=1:7 z=3:2:9 v=[y z] u=linspace(2,9,11) 二、矩阵的初始化 逗号或者空格分隔同一行的元素，分号表示一行的结束。 m=[1 2 3 4;5 6 7 8;9 10 11 12] 几个常用矩阵函数 ones(m,n) 产生一个m行n列的元素全为1的矩阵 zeros(m,n) 产生一个m行n列的零矩阵 eye(m,n) 产生一个m行n列的单位矩阵 三、子矩阵的提取 A(v1,v2) v1表示要提取的行号构成的向量 v2表示要提取的列号构成的向量 v1若为:，则表示提取所有行 v2若为:，则表示提取所有列 关键字end表示最后一行或者最后一列 a=[1 2 3;4 5 6;7 8 9] a1=a(2,:) a2=a(:,2) a3=a(1:2,2:3) a4=a(2:-1:1,:) a5=a(:,3:-1:2) a6=[a a2] a7=[a;a1] 四、矩阵的代数运算 A+B 矩阵的代数加法 A-B 矩阵的代数减法 A*B 矩阵的代数乘法 A’ 矩阵的转置 Ax=B的解 x=A\B（左除） xA=B的解 x=B/A（右除） A^x（矩阵乘方）比如A^2 A^3 说明：A必须为方阵 五、矩阵的点运算（MATLAB一大特色） 两个矩阵之间的点运算是它们对应元素的直接运算。 A.*B A.^x（A不一定是方阵） A.^A A./B 说明：函数也是以点运算方式进行。 六、MATLAB函数 三角函数、指数函数、整值函数、绝对值函数、矩阵操作函数、矩阵函数 七、选择结构 if (条件) 语句组一; else 语句组二; end 如果条件为真，就执行语句组一，否则执行语句组二。 八、循环结构 for i=v 循环体语句; end while (条件式) 循环体语句; end 1、for循环允许一组命令以固定的和预定的次数重复。在for和end语句之间的循环体语句按向量v中的每一个元素执行一次。在每一次循环后，i被指定为向量的下一个元素。 例：求 for n=1:10 x(n)=sin(n*pi/10); end x 2、while循环以不定的次数求一组语句的值。只要条件式为真，就执行while和end语句之间的循环体语句。 例：设银行年利率为11.25%。将10000元钱存入银行，问多长时间会连本带利翻一番？ money=10000; years=0; while money20000 years=years+1; money=money*(1+11.25/100); end 例：求100以内的勾股数（编写命令文件） G=[ ]; for a=1:98 for b=(a+1):99 for c=(b+1):100 if a^2+b^2==c^2 G(size(G,1)+1,:)=[a,b,c]; end end end end G 九、二维作图 1、plot(x,y) 直角坐标，x，y是向量,分别表示点集的横坐标和纵坐标。 例：函数作图 输入以下命令： x=-2:0.05:2; y=sin(x)./x; plot(x,y) 2、polar(t,r) 极坐标，t，r是向量,分别表示点集的角坐标和极半径坐标。 例：作图 输入以下命令： t=0:0.01:6*pi; r=5*sin(4*t./3); polar(t,r) 三维作图 1、plot3(x,y,z) 空间曲线，x，y，z是向量 例：函数作图 输入以下命令： t=0:0.02:10*pi; plot3(sin(t),cos(t),t) 2、mesh(X,Y,Z) 空间曲面，X，Y，Z是矩阵 例：函数作图 输入以下命令： x=-2:0.02:2; y=x; [X,Y]=meshgrid(x,y); R=sqrt(X.^2+Y.^2)+eps; Z=sin(R)./R; mesh(X,Y,Z) 符号函数作图 ezplot(‘f(x)’,[a,b]) 在[