实验解答20110530精要.docVIP

实验一 MATLAB工作环境熟悉及简单命令的执行 一、esktop），包括菜单（Menu）、工具条 （Toolbar）、命令窗口(Command Window)、历史命令窗口、工作空间(Workspace)等；完成一些基本的矩阵操作；学习使用在线帮助系统。 三、实验步骤： 1、workspace的变化，记录运算结果。 （1）（365-52(2-70）(3 (365-52*2-70)/3 ans = 63.6667 （2）area=pi*2.5^2 area = 19.6350 （3）已知x=3，y=4，在MATLAB中求z： x=3 y=4 z = x ^2 * y ^3 / (x - y) ^2 z = 576 （4）将下面的矩阵赋值给变量m1,在workspace中察看m1在内存中占用的字节数。 m1= 执行以下命令 m1 =[16 2 3 13 ; 5 11 10 8 ; 9 7 6 12 ; 4 14 15 1 ] m1( 2 , 3 ) ans = 10 m1( 11 ) ans = 6 m1( : , 3 ) ans =3 10 6 15 m1( 2 : 3 , 1 : 3 ) ans =5 11 10 9 7 6 m1( 1 ,4 ) + m1( 2 ,3 ) + m1( 3 ,2 ) + m1( 4 ,1) ans = 34 （5）执行命令help abs 查看函数abs的用法及用途，计算abs( 3 + 4i ) （6）执行命令 x=0:0.1:6*pi; y=5*sin(x); plot(x,y) （6）运行MATLAB的演示程序，demo，以便对MATLAB有一个总体了解。 四、思考题 1、x2 合法 （2）3col 不合法，首字符是数字 （3）_row 不合法，首字符必须是字母 （4）for 不合法，for为MATLAB的保留关键字 2、求以下变量的值，并在MATLAB中验证。 （1）a = 1 : 2 : 5 ; a = 1 3 5 （2）b = [ a a a ] ; b = 1 1 1 3 3 3 5 5 5 （3）c = a + b ( 2 , : ) c =4 6 8 实验二 MATLAB语言矩阵运算 一、 1、下列运算是否合法，为什么？如合法，结果是多少？ (1) result1 = a (2) result2 = a * b 不合法，两矩阵的内阶不相等。 (3) result3 = a + b (4) result4 = b * d (5) result5 = [b ; c ] * d (6) result6 = a . * b (7) result7 = a . / b (8) result8 = a . * c 不合法，两矩阵的维数不相等，点乘相求两矩阵的维数相同。 (9) result9 = a . \ b (10) result10 = a . ^2 (11) result11 = a ^2 不合法，a矩阵不是方阵，方阵才能取平方。 (12) result12 = 2 . ^ a 以上运算的关键不是结果是多少，而是结果怎么得来的，比如result5 = [b ; c ] * d，要学会将其分解后去理解，可在命令窗口输入以下语句： c t=[b;c] r5=t * d 2、用MATLAB求下面的的方程组。 (1) A=[7,2,1,-2;9,15,3,-2;-2,-2,11,5;1,3,2,13] b=[4;7;-1;0] x=A\b x = 0.4979 0.1445 0.0629 -0.0813 (2) A=[1,1,1,0;1,2,1,-1;2,-1,0,-3;3,3,5,-6] b=[1;8;3;5] xyzw=inv(A)*b xyzw= 1.0000 5.0000 -5.0000 -2.0000 3、已知 A=[7,2,1,-2;9,15,3,-2;-2,-2,11,5;1,3,2,13] (1)求矩阵A的秩(rank) rank(A) (2)求矩阵A的行列式(determinant) det(A) (3)求矩阵A的逆(inverse) inv(A) (4)求矩阵A的特征值及特征向量(eigenvalue and eigenvector) [v,d]=eig(A)