原创力文档- 1、本文档共7页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
MATLAB作业1参考解析
MATLAB作业1参考答案
1 、在你的机器上安装MATLAB 语言环境，并键入demo 命令，由给出的菜单系统和对话框原型演示程序，领略MATLAB 语言在求解数学问题方面的能力与方法。
【求解】略.
2、启动MATLAB 环境，并给出语句tic, A=rand(500); B=inv(A); norm(A*B-eye(500)),toc，试运行该语句，观察得出的结果，并利用help 命令对你不熟悉的语句进行帮助信息查询，逐条给出上述程序段与结果的解释。
【求解】在MATLAB 环境中感触如下语句，则可以看出，求解随机矩阵的逆，并求出得出的逆矩阵与原矩阵的乘积，得出和单位矩阵的差，得出范数。一般来说，这样得出的逆矩阵精度可以达到。
tic, A=rand(500); B=inv(A); norm(A*B-eye(500)), toc
ans =
1.2333e-012
Elapsed time is 1.301000 seconds.
3、试用符号元素工具箱支持的方式表达多项式，并令，将f(x) 替换成s 的函数。
【求解】可以先定义出f 函数，则由subs() 函数将x 替换成s 的函数
syms s x;f=x^5+3*x^4+4*x^3+2*x^2+3*x+6;F=subs(f,x,(s-1)/(s+1))
F =
(s-1)^5/(s+1)^5+3*(s-1)^4/(s+1)^4+4*(s-1)^3/(s+1)^3+
2*(s-1)^2/(s+1)^2+3*(s-1)/(s+1)+6
syms x s
x=(s-1)./(s+1);
y=x.^5+3*x.^4+4*x.^3+2*x.^2+3*x+6;
y
simple(y)
y=
19 - (72*s^4 + 120*s^3 + 136*s^2 + 72*s + 16)/(s + 1)^5
4、用MATLAB 语句输入矩阵A 和B
前面给出的是4 ×4 矩阵，如果给出A(5; 6) = 5 命令将得出什么结果？
【求解】用课程介绍的方法可以直接输入这两个矩阵
A=[1 2 3 4; 4 3 2 1; 2 3 4 1; 3 2 4 1]
A =
1 2 3 4
4 3 2 1
2 3 4 1
3 2 4 1
若给出A(5,6)=5 命令，虽然这时的行和列数均大于B 矩阵当前的维数，但仍然可以执行该语句，得出
A(5,6)=5
A =
1 2 3 4 0 0
4 3 2 1 0 0
2 3 4 1 0 0
3 2 4 1 0 0
0 0 0 0 0 5
复数矩阵也可以用直观的语句输入
B=[1+4i 2+3i 3+2i 4+1i; 4+1i 3+2i 2+3i 1+4i;
2+3i 3+2i 4+1i 1+4i; 3+2i 2+3i 4+1i 1+4i];
B =
1.0000 + 4.0000i 2.0000 + 3.0000i 3.0000 + 2.0000i 4.0000 + 1.0000i
4.0000 + 1.0000i 3.0000 + 2.0000i 2.0000 + 3.0000i 1.0000 + 4.0000i
2.0000 + 3.0000i 3.0000 + 2.0000i 4.0000 + 1.0000i 1.0000 + 4.0000i
3.0000 + 2.0000i 2.0000 + 3.0000i 4.0000 + 1.0000i 1.0000 + 4.0000i
5、假设已知矩阵A ，试给出相应的MATLAB 命令，将其全部偶数行提取出来，赋给B 矩阵，用A =magic(8) 命令生成A 矩阵，用上述的命令检验一下结果是不是正确。
【求解】魔方矩阵可以采用magic() 生成，子矩阵也可以提取出来
A=magic(8), B=A(2:2:end,:)
A =
64 2 3 61 60 6 7 57
9 55 54 12 13 51 50 16
17 47 46 20 21 43 42 24
40 26 27 37 36 30 31 33
32 34 35 29 28 38 39 25
41 23 22 44 45 19 18 48
49 15 14 52 53 11 10 56
8 58 59 5 4 62 63 1
B =
9 55 54 12 13 51 50 16
40 26 27 37 36 30 31 33
41 23 22 44 45 19 18 48
8 58 59 5 4 62 63 1
6、用MATLAB 语言实现下面的分段函数。
【求解】两种方法，其一，巧用比较表达式解决
y=h*(xD) + h/D*x.*(abs(x)=D) -h*(x-D);
另外一种方法，用循环语句和条件转移语句
文档评论（0）