资源描述：
《matlab实验报告(期末作业)》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、阿布力克木.奥布力计算机在生命科学当中的应用学号：1040407114计算机在生命科学中的应用《MATLAB》实验报告课题名称：计算机在生命科学中的应用专业：信息管理与信息系统班级：10404071学号：1040407114学生姓名：阿卜力克木.奥不力学校：江苏科技大学2011年11月18日阿布力克木.奥布力计算机在生命科学当中的应用学号：1040407114试验名称MATLAB试验（一）试验目的1.矩阵的加(+)2.矩阵的减(-)3.矩阵的数乘(*)参加试验时间2011-10-11试验地点江苏科技大学西校区图书馆四楼试验代码和试验结果
2、1.打开MATLAB程序在代码区里输入如图所示的代码，并分析试验结果注意：A=[12;34];B=[5,6;7,8];C=A+B,D=A-B则回车后，得到的结果就只有C和D的值了。由此可见其中的两个分号(一个紧接着A=[12;34]的后面，另一个紧接着B=[5,6;7,8]的后面)的作用2．两个数的相乘运算在MATLAB程序在代码区里输入如图所示的代码并仔细观察试验结果阿布力克木.奥布力计算机在生命科学当中的应用学号：1040407114试验代码和试验结果矩阵的相乘(*)矩阵的数乘试验结果：如图所示阿布力克木.奥布力计算机在生命科学当中
3、的应用学号：1040407114试验名称MATLAB试验（二）试验目的1.矩阵乘方(^)2.除法运算3.矩阵转置参加试验时间2011-10-25试验地点江苏科技大学西校区图书馆四楼试验代码和试验结果1.MATLAB提供了两种除法运算：左除（）和右除（/）。当矩阵A可逆时，X=AB是方程A*X=B的解(这里当然还要求A的行数=B的行数)，而X=C/A是方程X*A=C的解(这里当然还要求A的列数=C的列数)。2.在命令框（commandwindow）中输入如图所示的代码并查看结果3.矩阵乘方(^)注意：上面的[1,2;2,1]^(-2)
4、要求[1,2;2,1]是可逆的。阿布力克木.奥布力计算机在生命科学当中的应用学号：1040407114试验代码和试验结果4.MATLAB程序中乘方(^)示这样实现的注意：上面的[1,2;2,1]^(-2)要求[1,2;2,1]是可逆的。5.矩阵转置(’)注意：若A为复数矩阵，则A’表示A的共轭转置。对于复数矩阵A，若仅希望转置，则用如下命令：A.’。见图阿布力克木.奥布力计算机在生命科学当中的应用学号：1040407114试验名称MATLAB试验（三）试验目的1.求线性方程的唯一解或特解2.求线性方程的通解参加试验时间2011-11-2
5、试验地点江苏科技大学西校区图书馆四楼试验代码和试验结果（1）求方程组的一个特解。解：用MATLAB把该方程组的增广矩阵化为行最简形(见图，从中可以看出该方程组有无数多解，而且X=[1.2500–0.250000]’就是该方程组的一个特解。用矩阵的初等变换解线性方程组（2）求方程组的解.[分析]记该方程组的增广矩阵列为a_1，a_2，a_3，a_4，a_5，b，并且令系数矩阵的行列式为D，依次用b替换D所得到的行列式分别记为D_1，D_2，D_3，D_4，D_5。根据克拉默法则，当D¹0时，该方程组有唯一解：x_1=D_1/D，x_2=D
6、_2/D，x_3=D_3/D，x_4=D_4/D，x_5=D_5/D。见图阿布力克木.奥布力计算机在生命科学当中的应用学号：1040407114用克拉默法则解线性方程组（3）求线性方程的通解在MATLAB中，函数null用来求解矩阵A的零空间的一组基，即齐次线性方程组A*X=0的解空间的一组基(基础解系)，由此可得齐次线性方程组的通解。命令格式：null(A)或null(A,’r’)其中null(A)的返回值是一个矩阵，其列向量组为A的零空间的一组标准正交基，而null(A,’r’)的列向量组为A的零空间的一组基(一般不是单位向量组，也
7、未必是两两正交的)。（4）求方程组的通解。解：先用函数null求系数矩阵A=的零空间的一组基阿布力克木.奥布力计算机在生命科学当中的应用学号：1040407114试验名称MATLAB试验（四）试验目的1认识MATLAB很强的绘图功能1.MATLAB平面图的绘制2.MATLAB三维图的绘制参加试验时间2011-11-16试验地点江苏科技大学西校区图书馆四楼试验代码和试验结果（1）画出椭圆在区域[-3,3]´[-4,4]内的图形。分析：在MATLAB的命令窗口输入如下命令：ezplot('x^2/4+y^2/5-1',[-3,3,-4,4]
8、)运行后得：用ezplot命令绘制的x2/4+y2/5=1的图象(2)画出曲线在区间[0,p]内的图形。解：在MATLAB的命令窗口输入如下命令：ezplot('sin(3*t)*cos(t)','sin(
