《MATLAB作业1》由会员分享，可在线阅读，更多相关《MATLAB作业1（8页珍藏版）》请在金锄头文库上搜索。
1、MATLAB 作业 1 1 、在你的机器上安装MATLAB 语言环境， 并键入 demo 命令， 由给出的菜单系统和对话框原型演示程序，领略MATLAB 语言在求解数学问题方面的能力与方法。解：2、启动 MATLAB 环境，并给出语句tic, A=rand(500); B=inv(A); norm(A*B-eye(500),toc，试运行该语句，观察得出的结果，并利用help 命令对你不熟悉的语句进行帮助信息查询，逐条给出上述程序段与结果的解释。解：tic 命令： 开启一个新的计时点，记录下当前的时刻，与toc 命令配合使用，用于计算两点之间的时间间隔。rand 命令：rand(N) 生成一个
2、 N维的矩阵， 矩阵中的元素都是0和1之间的服从均匀分布的随机数。inv 命令： inv(A) 求一个矩阵的转置。norm命令： norm(A, N) 求矩阵或向量A的N范数。toc 命令：读取秒表的时间，和tic 命令配合使用，用于计算两点之间的时间间隔。3、试用符号元素工具箱支持的方式表达多项式5432( )34236f xxxxxx，并令11sxs，将 f(x) 替换成 s 的函数。解：f(s)=19 - (72*s4 + 120*s3 + 136*s2 + 72*s + 16)/(s + 1)5 4、用 MATLAB 语句输入矩阵A 和B 12341423324143214132231
3、4,234123324114324132234114jjjjjjjjABjjjjjjjj前面给出的是 4 4 矩阵，如果给出A(5， 6) = 5 命令将得出什么结果？解：命令 A(5, 6) = 5, 即将原来的方阵A扩展为一个五行六列的矩阵，并且第五行第六列的元素值为 5。5、假设已知矩阵A ，试给出相应的MATLAB 命令，将其全部偶数行提取出来，赋给B 矩阵，用A =magic(8) 命令生成 A 矩阵，用上述的命令检验一下结果是不是正确。解：6、用 MATLAB 语言实现下面的分段函数,( )/,hxDyf xh DxxDhxD。解：7、用数值方法可以求出63 62630212482
4、2iiS，试不采用循环的形式求出和式的数值解。由于数值方法采用double 形式进行计算的，难以保证有效位数字，所以结果不一定精确。试采用符号运算的方法求该和式的精确值。解：8、 编写一个矩阵相加函数mat_add() ， 使其具体的调用格式为A=mat_add(A1 ,A2 ,A3 , ) ，要求该函数能接受任意多个矩阵进行加法运算。（注：varargin 变量的应用）解：9 已知 Fibonacci 数列由式12,3,4,kkkaaak可以生成， 其中初值为121aa，试编写出生成某项Fibonacci 数值的 MATLAB 函数，要求函数格式为 y=fib(k) ，给出 k 即能求出第
5、k 项ka并赋给 y 向量；编写适当语句，对输入输出变量进行检验，确保函数能正确调用；利用递归调用的方式编写此函数。（注： 递归的调用方式速度较慢，比循环语句慢很多，所以不是特别需要，解这样问题没有必要用递归调用的方式。）解：10、下面给出了一个迭代模型2 1111.40.3kkkkkxyxyx写出求解该模型的M-函数（ M-脚本文件），如果取迭代初值为000,0xy，那么请进行30000 次迭代求出一组x 和y 向量，然后在所有的kx和ky坐标处点亮一个点( 注意不要连线 ) ，最后绘制出所需的图形。（提示这样绘制出的图形又称为Henon 引力线图，它将迭代出来的随机点吸引到一起，最后得出貌似连贯的引力线图。）解：11、选择合适的步距绘制出下面的图形1sin( )t，其中( 1,1)t。( 注：合适的步距包括等距与不等距 ) 解：12、对合适的范围选取分别绘制出下列极坐标图形（注：要求把图形窗口分为4块，每块绘一个图）21.0013，cos(7/ 2)，sin( ) /，31cos (7 )解：13、请分别绘制出xy和sin()xy的三维图和等高线。解：
