（标准作业指导）MATLAB实验指导书加程序+上机实例.doc

(((((((((((((((((((((((( MATLAB语言实验指导书 (((((((((((((((((((((((( 华东交通大学电气学院 张永贤 2006年2月 实验一 MATLAB工作环境熟悉及简单命令的执行 一、esktop），包括菜单（Menu）、工具条 （Toolbar）、命令窗口(Command Window)、历史命令窗口、工作空间(Workspace)等；完成一些基本的矩阵操作；学习使用在线帮助系统。 三、实验步骤： 1、workspace的变化，记录运算结果。 （1）（365-52(2-70）(3 （2）area=pi*2.5^2 （3）已知x=3，y=4，在MATLAB中求z： （4）将下面的矩阵赋值给变量m1,在workspace中察看m1在内存中占用的字节数。 m1= 执行以下命令 m1( 2 , 3 ) m1( 11 ) m1( : , 3 ) m1( 2 : 3 , 1 : 3 ) m1( 1 ,4 ) + m1( 2 ,3 ) + m1( 3 ,2 ) + m1( 4 ,1) （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 2、求以下变量的值，并在MATLAB中验证。 （1）a = 1 : 2 : 5 ; （2）b = [ a a a ] ; （3）c = a + b ( 2 , : ) 实验二 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 (12) result11 = 2 . ^ a 2、用MATLAB求下面的的方程组。 (1) (2) 3、已知 (1)求矩阵A的秩(rank) (2)求矩阵A的行列式(determinant) (3)求矩阵A的逆(inverse) (4)求矩阵A的特征值及特征向量(eigenvalue and eigenvector) 4、关系运算与逻辑运算 已知a=20,b=-2,c=0,d=1 (1) r1 = a b (2) r2 = a b c d (3) r3 = a == b* (-10) (4) r4 = ~b | c 三、思考题 ,求y=？（用format long查看y的值） 实验三 程序的编辑及调试 一、ile|New|M-File，启动MATLAB的程序编辑及调试器（Editor/Debugger），编辑以下程序，点击File|Save保存程序，注意文件名最好用英文字符。点击Debug|Run运行程序，在命令窗口查看运行结果，程序如有错误则改正。 注：数论中一个有趣的题目：任意一个正整数，若为偶数，则用2除之，若为奇数，则与3相乘再加上1。重复此过程，最终得到的结果为1。如： 2(1 3(10(5(16(8(4(2(1 6(3(10(5(16(8(4(2(1 运行下面的程序，按程序提示输入n=1,2,3,5,7等数来验证这一结论。 %classic 3n+1 problem from number theory. while 1 n=input(Enter n,negative quits:); if n=0 break end a=n; while n1 if rem(n,2)==0 n=n/2; else n=3*n+1; end a=[a,n]; end a end 2、编程求满足的最小m值。 三、思考题 用对分法求解方程在[0，1]内的解，并验证，在程序中统计出对分次数。 提示：先将原方程转化成的形式。 对分法的基本思想是：一个一元方程f(x)=0，若f(x1)*f(x2)0，x1,x2]区间内有实数解。取该区间的