matlab实验指导答案详解(非常详细正确)汇编
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 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)