实验二十四 MATLAB语言基础实验报告.doc

实验二十四 MATLAB语言基础 年级专业:10级计算机科学与技术2班 姓名: 完成时间:2011.11.23 (一)、实验目的 了解Matlab语言及其在理工科课程中的应用。 （二）、实验仪器 装有Matlab软件的计算机一台。 （三）、实验内容 1、矩阵运算。 (1）矩阵的输入。输入矩阵最简单的方法是把矩阵的元素直接排列在括号中。每行内的元素间用空格或逗号公开，行与行的内容间用分号隔开。 例如输入:A=[1 1 1;2 2 2;3 3 3] 截图如下： (2）矩阵转置。矩阵的转置用符号“ ’ ”来表示和实现。 例如输入:A=[1 1 1;2 2 2;3 3 3] 截图如下： (3）矩阵加减。矩阵的加减运算使用的是“+”和“-”运算符。而矩阵必须具有相同阶数才进行加、减运算。例如A是3X3矩阵，D是3X1矩阵，就不能进行A+D运算。 例如输入A=[1 3 5;7 9 11;13 15 17] B=[2 4 6;8 10 12;14 16 18]则截图如下： 加法： 减法： (4）矩阵乘法用*表示，当两矩阵中前一矩阵的列数和后一矩阵和行数相同，可以进行乘法运算，这与数学上的形式是一致的，两个相同维数向量的内积也可用这种乘法来实现。 例如输入A=[1 2 3]’,B=[4 5 6]’。截图如下： (5) 矩阵的行列式。N阶矩阵A的行列式由det(A)给出，如 A=[1 2 3;4 5 6;7 8 9]截图如下： (6）矩阵求逆，非奇异矩阵A的逆矩阵由inv(A)给出。例如A=[1 2 3;4 5 6;7 8 0]结果为： (7)特征性。如果A是n阶矩阵，满足式Ax＝λxλ 称为A的特征值，x为A的特征向量。例如A=[0 1;-1 0]截图如下： 2、绘图 （1）二维数据作图 如果y是一个向量，那么plot(y)绘制一个y元素和y元素序号之间关系的线性坐标图，例如要画出y元素的序号1，2，3，4，5，6，7，对应的y元素分别为0，0.48，0.84，1，0.91，0.6，0.14的图形，则输入命令 y[0，0.48，0.84，1.0，0.91，0.6，0.14] plot(y) 坐标x、y的刻度是自动给出的。 MATLAB显示图形是以给出图形窗口来显示的。若想对图形加上标注，则将当前窗口切换成MATLAB命令窗口并使之成为流动窗口，再从命令窗口中输入相应的图形加注命令。 例如：输入下列命令 title(‘My first plot’); xlabel(‘fortnights’); ylabel(‘furions’); grid 这里，第一条语句是为了给出图的题头，第二条、第三条语句分别为了给出x轴与y轴标注，第四条是为了图形加上网格。最后结果截图如下： （2)二维函数作图 如果x和y是同样长度的向量，plot(x，y)命令可画出x元素对于y元素的y-x坐标图。x＝0:0.05:4*pi; y＝sin(x); plot(x,y) x所取得区间0至4π，步长为0.05。第二、三条语句绘制正弦函数图。 (3）三维函数作图 为显示两变量函数Z＝f(x,y)的图形,第一步是产生特定期的X和Y矩阵,它们在整个函数定义域中分别由重复的行和列构成。在这个基础上,函数便可以直接计算和绘图。 考察sin(r)/r函数,它产生像”阔边帽”形状的图形,对此图形用以下方法: x＝-8:0.5:8; y＝x’ X＝ones(size(y))*x; Y＝y*ones(size(x)); R＝sqrt(X.^2+Y.^2)+eps; Z＝sin(R)./R; mesh(Z)x，y）的距离；之后形成sin(R)./R函数并用mesh命令给出结果。截图如下： 3、数值积分 函数语句 [Q，CNT]＝QUAD(’FUNFCN’，A，B，TOL，TRACE） 参数说明： FUNFCN包含函数名的字符串，给出被积函数f A 输入参数，函数f的积分上限 B 输入参数，函数f的积分下限 TOL 输入参数，积分精度，为小的正数 TRACE 输入参数，如果trace 不为零，则作出积分函数f的积分图；否则不作图 Q 输出参数， 函数的积分值 CNT 参数输出 ，函数f 的估计值 例 计算定积分。利用下列命令 I＝quad(‘log(1+x)’,0,1) 4、roots(a),其中a是多项式的系数组成的向量。 例：求三次方程式的根。截图如下：