matlab实验指导答案详解(非常详细正确)汇编

合集下载
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

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

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

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

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

三、实验步骤:

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

2、在命令窗口执行命令完成以下运算,观察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 :

()

2

3

2y x y x z -= >>x=3 >>y=4

>>z = x ^2 * y ^3 / (x - y) ^2 z = 576

(4)将下面的矩阵赋值给变量m1,在workspace 中察看m1在内存中占用的字节数。

m1=⎥⎥⎥⎥

⎤⎢⎢⎢

⎢⎣⎡11514412679810115133216 执行以下命令

>>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、以下变量名是否合法?为什么? (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 语言矩阵运算

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

⎥⎦⎤⎢⎣⎡=654321a ⎥⎦⎤⎢⎣⎡-=531142b ⎥⎥

⎥⎦

⎢⎢⎢⎣⎡-=201c ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=063258741d

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)⎥⎥⎥⎥

⎤⎢⎢⎢⎢⎣⎡-=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⨯⎥⎥⎥⎥⎦⎤⎢⎢⎢

⎢⎣⎡----01741323151122231592127

4321x x x x 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)⎪⎪⎩⎪⎪⎨⎧=-++=--=-++=++5

6533332821

w z y x w y x w z y x z y x 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、已知⎥⎥⎥⎥

⎤⎢⎢⎢

⎢⎣⎡----=1323151122231592127

A 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)

相关文档
最新文档