数值计算实验.docVIP

实验1 Matlab基础实验 成绩 专业班级 数学142 学号 201412010217 姓名 王倩倩 报告日期 2016.3.29 . 实验类型：●验证性实验 ○综合性实验 ○设计性实验 实验目的：熟悉MATLAB系统的启动、退出、演示系统、帮助系统、 MATLAB的运行环境，掌握Matlab矩阵、算符、表达式、数据结构、数据类型、函数与程序设计及其运行、画图等的基础知识。 实验内容： MATLAB系统的启动、退出、演示系统、帮助系统，Matlab中Notebook的应用，矩阵的建立、矩阵的拆分、特殊矩阵、矩阵算符及其表达式表示和计算，矩阵求逆，求行列式的值，矩阵特征值与特征向量计算，结构体、元胞数组，常用数学函数、画图函数应用，字符串与数值的转换，数据类型的转换。编写求三位整数中是水仙花数的程序，所谓水仙花数是指一个3位数，其各位数字的立方之和等于该数本身，例如153是一个水仙花数，13+53+33=153。编写求三位整数中是完数的程序，一个数如果恰好等于它的因子之和，这个数就称为完数。 实验报告：根据实验情况和结果撰写并递交实验报告。 实验总结(学会了......; 掌握了......; 训练了......; 发现了......; 今后学习中......有待提高。） 电子报告word文件命名规则：专业班级-学号后两位-实验X-姓名.doc, 如信息123班学号为201212030315的郭海涛同学实验1报告word文件命名则应是：信息123-15-实验1-郭海涛.doc, 其中 .doc是Word文件扩展名。特别提醒：电子报告文件命名不规范的报告将不予接收。 1、 矩阵建立 a=[1 2 3;4 5 6;7 8 9] a = 1 2 3 4 5 6 7 8 9 b=[4 6 9;4 7 12;5 8 9] b = 4 6 9 4 7 12 5 8 9 2、矩阵的拆分 （1）取出矩阵的第一列 C=a(:,1) C = 1 4 7 （2）取出矩阵后两列 D=a(:,2:3) D = 2 3 5 6 8 9 3、 矩阵运算 （1）矩阵加法 c=a+b c = 5 8 12 8 12 18 12 16 18 （2）矩阵减法 d=a-b d = -3 -4 -6 0 -2 -6 2 0 0 （3）矩阵的乘 A=a*b A = 27 44 60 66 107 150 E=eig(A): 求矩阵的A的全部特征值，构成向量E A=[1 5;6 7] E=eig(A) E = -2.2450 10.2450 或者 A=[1 5;6 7] A = 1 5 6 7 [V,D]=eig(A) V = -0.8388 -0.4757 0.5444 -0.8796 D = -2.2450 0 0 10.2450 即求矩阵A的全部特征值，构成对角阵D,并求A的特征向量构成V的列向量。 7、常用数学函数 abs(-4) ans = 4 sin(pi/4) ans = 0.7071 exp(2) ans = 7.3891 8、画图函数应用 x=0:0.2:2*pi x = Columns 1 through 7 0 0.2000 0.4000 0.6000 0.8000 1.0000 1.2000 Columns 8 through 14 1.4000 1.6000 1.8000 2.0000 2.2000 2.4000 2.6000 Columns 15 through 21 2.8000 3.0000 3.2000 3.4000 3.6000 3.8000 4.0000 Columns 22 through 28 4.2000 4.4000 4.6000 4.8000 5.0000 5.2000 5.4000 A = 27 44