MATLAB答案

???????????????????????? MATLAB语言实验指导书 ????????????????????????

华东交通大学电气学院

张永贤

2006年2月

实验一 MATLAB工作环境熟悉及简单命令的执行

一、实验目的：熟悉MATLAB的工作环境，学会使用MATLAB进行一些简单的运算。 二、实验内容：MATLAB的启动和退出，熟悉MATLAB的桌面（Desktop），包括菜单

（Menu）、工具条 （Toolbar）、命令窗口(Command Window)、历史命令窗口、工作

空间(Workspace)等；完成一些基本的矩阵操作；学习使用在线帮助系统。

三、实验步骤：

1、启动MATLAB，熟悉MATLAB的桌面。

2、在命令窗口执行命令完成以下运算，观察workspace的变化，记录运算结果。 （1）（365-52?2-70）?3 （2）>>area=pi*2.5^2

（3）已知x=3，y=4，在MATLAB中求z：

x2y3 z?2?x?y?（4）将下面的矩阵赋值给变量m1,在workspace中察看m1在内存中占用的字节数。

?162313??511108?? m1=??97612????414151?执行以下命令 >>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、以下变量名是否合法？为什么？ （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 , : )

1

实验二 MATLAB语言矩阵运算

一、实验目的：掌握基本的矩阵运算及常用的函数。 二、实验内容：

?147??1??123??24?1????? c?0 d?852 a??b?????????456??135????360????2??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求下面的的方程组。

21?2??x1??4??7?9153?2??x??7????2???? (1)???2?2115??x3???1???????13213???x4??0??x?y?z?1?x?2y?z?w?8?(2)? ?2x?y?3w?3??3x?3y?5z?6w?521?2??7?9153?2?? 3、已知A????2?2115???13213??(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

2

(3) r3 = a == b* (-10) (4) r4 = ~b | c

三、思考题

y?n??10n?10?99102?2?2???2?2?,求y=？（用format long查看y的值）

10实验三 程序的编辑及调试

一、实验目的：掌握MATLAB程序编辑、运行及调试方法。

二、实验内容：

1、启动MATLAB后，点击File|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 \ while 1

n=input('Enter n,negative quits:'); if n<=0 break end a=n; while n>1

if rem(n,2)==0 n=n/2; else

n=3*n+1; end a=[a,n]; end a end

i2?10000的最小m值。 2、编程求满足?i?1m三、思考题

?x2e?sinx在[0，1]内的解，并验证，在程序中统计出对分次数。 用对分法求解方程

提示：先将原方程转化成

f(x)?2e?x?sinx?0的形式。

对分法的基本思想是：一个一元方程f(x)=0，若f(x1)*f(x2)<0，则在[x1,x2]区间内有实

数解。取该区间的中点xm=(x1+x2)/2，判定f(x1)和f(x2)二者中哪一个与f(xm)异号，若

3