【精品】数学实验-实验指导书2[1]5.docx

《【精品】数学实验-实验指导书2[1]5.docx》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、实验_MATLAB初步实验目的掌握矩阵运算和函数作图，掌握函数文件和命令文件的编写。实验内容矩阵运算、函数作图、编写函数文件、编写命令文件。实验准备一、向量的初始化v=[abc]各元素之间用逗号或者空格分开。v=first:increment:lastfirst是向量的初始值，increment是步长。v从first出发，每隔increment步长取一个元索，直至不超过last的最人值。若省略increment,则步长默认为kv=linspace(first,last,n)从first到last等间距的取n个元素。x二[123458718]y二1:7z=3:2:9v=[yz]

2、u=linspace(2,9,11)二、矩阵的初始化逗号或者空格分隔同一行的元索，分号表示一行的结束。m=[l234;5678;9101112]几个常用矩阵函数ones(m,n)产生一个m彳亍n列的元素全为1的矩阵zeros(m,n)产生一个m行n列的零矩阵eye(m,n)产生一个m行n列的单位矩阵三、了矩阵的提取A(vl,v2)vl表示要提取的行号构成的向量v2表示要提取的列号构成的向量vl若为：，则表示提取所有行吃若为：，则表示提取所有列关键字end表示最后一行或者最后一列a=[l23;456;789]al=a(2,:)a2=a(:,2)a3=a(1:2,2:3)a4=a

3、(2：-l:1,:)a5=a(:,3：-1:2)a6=[aa2]a7=[a;al]四、矩阵的代数运算矩阵的代数减法矩阵的转置说明：A必须为方阵A+B矩阵的代数加法A-BA*B矩阵的代数乘法A'Ax二B的解x二AB（左除）xA二B的解x二B/A（右除）Vx（矩阵乘方）比如AJA,五、矩阵的点运算（MATLAB—大特色）两个矩阵之间的点运算是它们对应元素的直接运算。A.*BA.八x（A不一定是方阵）A.八AA./B说明：函数也是以点运算方式进行。六、MATLAB函数三角函数、指数函数、整值函数、绝对值函数、矩阵操作函数、矩阵函数七、选择结构if（条件）语句组一；else语句组二

4、；end否则执行语句组二。while（条件式）循环体语句；end如果条件为真，就执行语句组一,循环结构fori二v循环体语句；end1、for循环允许一组命令以固定的和预定的次数重复。在for和end语句之间的循环体语句按向量v中的每一个元索执行一次。在每一次循环后，i被指定为向量的卜'一个元素。n7tTo"10例：求£sinn=lforn=l:10x(n)=sin(n*pi/10);end2>while循环以不定的次数求一组语句的值。只要条件式为真，就执行while和end语句之间的循环体语句。例：设银行年利率为11.25%。将10000元钱存入银行，问多长时间会连本带利翻

5、一番？money=10000;years二0;whi1emoney<20000ycars=ycars+l;money二money*（1+11.25/100）;end例：求100以内的勾股数（编写命令文件）G=[];fora=l:98forb=(a+l):99forc二(b+l):100ifa^2+b^2==c^2G(size(G,1)+1,:)=[a,b,c];endendendendG九、二维作图1、plot(x,y)直角坐标，x,y是向量，分别表示点集的横坐标和纵坐标。例：函数作图沁,xw[-2,2]X输入以下命令：x二-2:0.05:2;y=sin(x)・/x;plot

6、(x,y)2、polar(t,r)极处标，t,r是向量，分别表示点集的角坐标和极半径坐标。4例:作图p=5sin(；&),T=6龙输入以下命令：t=0:0.01:6*pi;r•二5*sin(4*t./3);polar(t,r)十、三维作图1、plot3(x,y,z)空间曲线，x,y,z是向量x=sint例：函数作图＜y=cost输入以下命令：t=0:0.02:10*pi;plot3(sin(t),cos(t),t)2、mesh(X,Y,Z)空间曲面,X,Y,Z是矩阵./22例：函数作图z=，兀“―2,2],)©-2,2]输入以下命令：x二-2:0.02:2;y二x;[X,Y]

7、=meshgrid(x,y);R=sqrt(X.2+Y.2)+eps;Z=sin(R)・/R;mesh(X,Y,Z)I•一、符号函数作图在[a,b]上作f(x)的图。ezplot("f(x)‘,[a,b])例：函数作图y=cosx9xe[0,龙]输入以卜•命令：ezplot("cos(x)',[0,pi])实验作业一、设有分块矩阵4=3x3^3x2*^2x2其屮E,R,0,S分别为单位阵、随机阵、零阵和对角阵，试通过数值计算验证屮ER+RS0S2二、函数作图2、四叶玫瑰线p=sin203、x=叶形线＜)