南华大学MATLAB实验报告1.doc

核科学技术学院 实 验 报 告 实验项目名称 MATLAB数值计算 所属课程名称 MATLAB及应用 实 验 类 型 上机实验 实 验 日 期 指 导 教 师 班 级 学 号 姓 名 成 绩 一、实验名称 MATLAB数值计算 二、实验目的 （1）掌握MATLAB变量的使用 （2）掌握MATLAB数组的创建 （3）掌握MATLAB数组和矩阵的运算 （4）熟悉MATLAB多项式的运用? 三、实验原理 1. 矩阵分析 ???? 矩阵转置：单引号（’） ???? 矩阵的旋转：rot90(A,k)，功能是将矩阵A逆时针旋转90度的k倍，缺省值是1 ???? 矩阵的左右翻转：fliplr(A) ???? 矩阵的上下翻转：flipud(A) ???? 矩阵的逆：inv(A)，与A^(-1)等价 ???? 矩阵的行列式：det(A) ???? 矩阵的秩： rank(A) ???? 矩阵的迹：trace(A) ???? 将矩阵化为最简式：rref(A) ???? 矩阵的特征值与特征向量：(1) E = eig(A)；矩阵A的所有特征值构成向量E；(2) [V,D]=eig(A)；A的所有特征值构成对角阵D，A的特征向量构成V的列向量； ? ?? 2. 多项式 ???? 多项式的建立：若多的项的全部根构成的向量为X，则以X为根的多项式为poly(X) ???? 多项式的根：roots(p)计算以向量p为系数的多项式的根，包括重根，复根 ???? 多项式求值：polyval(p,x)，p是多项式的系数，x可以是一个数也可以是一个矩阵 多项式求拟合次数：polyfit(x,y,n)，x可以是一个数也可以是一个矩阵，y是x对应的数或矩阵 ???? 多项式的四则运算：(1)P1+P2；(2)P1-P2；(3)conv(P1,P2), (4)deconv(P1,P2) 四、实验内容 1. 已知矩阵　　　 11 12 13 14 21 22 23 24 A= 31 32 33 34 41 42 43 44 (1) 如何输出A的第1列？ (2) 如何输出A的第2行？ (3) 如何输出A的第2列和第3列？ (4) 要得到，怎么办？ (5) 如何输出A的第1列和第3列？ (6) 如何将21和31按一列输出？ (7) 如何将A所有元素按从最左列至最右列新排列为一列输出？ (8) 如何将A原阵输出？ (9) 如何输出2阶全1阵？ (10) 如何输出2阶单位阵？ (11) 试用分块矩阵的方法生成阵 11 12 13 14 1 1 21 22 23 24 1 1 31 32 33 34 1 0 41 42 43 44 0 1 (12) 求A的主对角线 (13) 求除去A阵第1列后新阵的主对角线 (14) 求除去A阵第1行后新阵的主对角线 (15) 求除去A阵从左到右前2列后新阵的主对角线 (16) 要分别单独输出A阵元素11、22、33、44，怎么办？ (17) 如何将11、31、12、32、13、33、14、34按顺序一行输出？ 2. 已知 .求（1）A的逆；（2）A的行列式；（3）A的迹； （4）A的所有特征向量和特征值。 3. A=magic(3)，B=7*rand(3)，计算数组A、B乘积，计算AB，A|B，~A，A==B，AB 如何生成一个三阶幻方、3阶随机阵、1行2列全零阵、3行4列全一阵、3阶单位阵？ 4. 生成一个4阶Hilbert矩阵H，（1）求H的转置；（2）将H旋转90度；（3）对H实行左右翻转；（4）对H实行上下翻转。 5 输入如下矩阵A 　　　　　　　　　　 0　　 ?/3　 　　　　　　　A＝ ?/6 ?/2 (1) 求矩阵B1,　B1中每一元素为对应矩阵Ａ中每一