南华大学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 3441 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乘积，计算A&B，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 p/3A p/6 p/2(1) 求矩阵B1,B1中每一元素为对应矩阵中每一元素的正弦函数(2) 求矩阵B2, B2中每一元素为对应矩阵中每一元素的余弦函数(3) 求 B12+B22(4) 求矩阵的特征值与特征矢量：称特征矢量为，而特征值矩阵为(5) 求Msin(L)M-1 (6) 使用funm命令求矩阵A的正弦函数（结果应该与（5）同）(7) 求cosA(8) 证明 sin2A+cos2A=I I是全1阵6 某专业有三名研究生,本学期选修了四门课程,若这些研究生的姓名,学号,性别,出生年月,课程名称,考试成绩可任意假定, (1)分别用结构型变量和细胞型变量表示以上信息; (2)举例说明查阅以上任何一条信息的方法; (3)求每一个研究生的平均成绩.7. 已知多项式，试求：（1）p(x)的根；（2）由其根生成一个多项式q(x)并与p(x)比较；（3）计算p(1.5)，p(-2)，p(5)的值。8 已知矩阵A=1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4,求 (1)A的特征多项式 (2)特征多项式中未知数为20 时的值 (3)特征多项式的根 (4)特征多项式的导数9. 在实验中测得如下10组数据: X 1 2 3 4 5 6 10 12 15 16 Y 7 15 19 30 38 37 23 69 39 28（1）求最多能拟合多项式的次数是多少？并求出各项系数。（2）将数据点和拟合曲线在同一图中绘出。五、实验过程及结果（含源代码）1、在MATLAB中新建.m文件，输入实验代码如下并保存：clc,clearA=11 12 13 14;21 22 23 24;31 32 33 34;41 42 43 44;A1=A(:,1)A2=A(2,:)A3=A(:,2:3)A4=A(2:3,2:3)A5=A(:,1:2:3)A6=A(2:3)A7=A(:)A8=A(:,:)A9=ones(2,2)A10=eye(2)A11=A,ones(2,2);eye(2)A12=diag(A)A13=diag(A,1)A14=diag(A,-1)A15=diag(A,2)点击“run”运行得到结果如下：A1 = 11 21 31 41A2 = 21 22 23 24A3 = 12 13 22 23 32 33 42 43A4 = 22 23 32 33A5 = 11 13 21 23 31 33 41 43A6 = 21 31A7 = 11 21 31 41 12 22 32 42 13 23 33 43 14 24 34 44A8 = 11 12 13 14 21 22 23 24 31 32 33 34 41 42 43 44A9 = 1 1 1 1A10 = 1 0 0 1A11 = 11 12 13 14 1 1 21 22 23 24 1 1 31 32 33 34 1 0 41 42 43 44 0 1A12 = 11 22 33 44A13 = 12 23 34A14 = 21 32 43A15 = 13 242、 在命令窗口输入代码如下：A=2 3 -5 4;-1 2 7 -3;5 -4 2 1;7 6 -2 -3;inv(A)ans = 0.0198 0.0010 0.1069 0.0610 0.1199 0.1116 -0.0636 0.0271 0.0709 0.1616 0.0667 -0.0448 0.2388 0.1178 0.0777 -0.1069det(A)ans = -1918trace(A)ans = 3E=eig(A)E = -5.6545 4.6888 1.9828 + 8.2711i 1.9828 - 8.2711iV,D=eig(A)V = 0.3850 0.4655 -0.0889 + 0.4794i -0.0889 - 0.4794i -0.1196 0.4083 -0.1447 - 0.5624i -0.1447 + 0.5624i -0.1972 0.4872 0.5903 0.5903 -0.8936 0.6157 -0.1444 + 0.2359i -0.1444 - 0.2359iD = -5.6545 0 0 0 0 4.6888 0 0 0 0 1.9828 + 8.2711i 0 0 0 0 1.9828 - 8.2711i3、 在命令窗口键入实验代码如下： A=magic(3),B=7*rand(3)A = 8 1 6 3 5 7 4 9 2B = 6.6509 3.4019 3.1953 1.6180 6.2391 0.1295 4.2479 5.3347 5.7499 A*Bans = 80.3126 65.4622 60.1908 57.7778 78.7438 50.4824 49.6611 80.4287 25.4465 A&Bans = 1 1 1 1 1 1 1 1 1 A|Bans = 1 1 1 1 1 1 1 1 1 Aans = 0 0 0 0 0 0 0 0 0 A=Bans = 0 0 0 0 0 0 0 0 0 ABans = 1 0 1 1 0 1 0 1 0 A=magic(3)A = 8 1 6 3 5 7 4 9 2 B=rand(3)B = 0.9501 0.4860 0.4565 0.2311 0.8913 0.0185 0.6068 0.7621 0.8214 C=zeros(1,2)C = 0 0 D=ones(3,4)D = 1 1 1 1 1 1 1 1 1 1 1 1 E=eye(3)E = 1 0 0 0 1 0 0 0 14、在命令窗口键入实验代码及运行结果如下： H=hilb(4)H = 1.0000 0.5000 0.3333 0.2500 0.5000 0.3333 0.2500 0.2000 0.3333 0.2500 0.2000 0.1667 0.2500 0.2000 0.1667 0.1429 Hans = 1.0000 0.5000 0.3333 0.2500 0.5000 0.3333 0.2500 0.2000 0.3333 0.2500 0.2000 0.1667 0.2500 0.2000 0.1667 0.1429 rot90(H)ans = 0.2500 0.2000 0.1667 0.1429 0.3333 0.2500 0.2000 0.1667 0.5000 0.3333 0.2500 0.2000 1.0000 0.5000 0.3333 0.2500 fliplr(H)ans = 0.2500 0.3333 0.5000 1.0000 0.2000 0.2500 0.3333 0.5000 0.1667 0.2000 0.2500 0.3333 0.1429 0.1667 0.2000 0.2500 flipud(H)ans = 0.2500 0.2000 0.1667 0.1429 0.3333 0.2500 0.2000 0.1667 0.5000 0.3333 0.2500 0.2000 1.0000 0.5000 0.3333 0.25005、在命令窗口键入实验代码及运行结果如下： A=0 pi/3;pi/6 pi/2A = 0 1.0472 0.5236 1.5708 B1=sin(A)B1 = 0 0.8660 0.5000 1.0000 B2=cos(A)B2 = 1.0000 0.5000 0.8660 0.0000 B12+B22ans = 1.8660 1.3660 1.3660 1.8660 M,L=eig(A)M = -0.9628 -0.4896 0.2703 -0.8719L = -0.2940 0 0 1.8648 M*sin(L)*inv(M)ans = -0.1200 0.6048 0.3024 0.7873 funm(A,sin)ans = -0.1200 0.6048 0.3024 0.7873 funm(A,cos)ans = 0.7873 -0.6048 -0.3024 -0.1200 (funm(A,sin)2+(funm(A,cos)2ans = 1.0000 -0.0000 -0.0000 1.00006、 实验代码及运行结果如下：(1) 分别用细胞型变量和结构型变量表示题示信息如下：% 细胞型变量A1=李明,20130001,男,19890214,MATLAB,99,ENGLISH,88,Matrix,96,Math,94;A2=张远,20130006,男,19900425,MATLAB,98,ENGLISH,80,Matrix,90,Math,96;A3=王云,20130009,女,19901001,MATLAB,97,ENGLISH,89,Matrix,88,Math,90;% 结构型变量B1=struct(Name,李明,StuNumber,20130001,Sex,男,Birth,19890214,. MATLAB,99,ENGLISH,88,Matrix,96,Math,94);B2=struct(Name,张远,StuNumber, 20130006,Sex,男,Birth, 19900425,. MATLAB,98,ENGLISH,80,Matrix,90,Math,96);B3=struct(Name,王云,StuNumber, 20130009,Sex,女,Birth, 19901001,. MATLAB,97,ENGLISH,90,Matrix,88,Math,90);C=B1 B2 B3;(2) 查询信息，如查询李明的学号，实验代码如下：name=李明information=StuNumber;for i=1:3 if strcmp(C(i).Name,name)=1 disp(information); getfield(C(i),information) break; endend运行结果如下：name =李明StuNumberans = 20130001(3) 求每位研究生的平均成绩：for i=1:3Average(i)=(C(i).MATLAB+C(i).ENGLISH+C(i).Matrix+C(i).Math)/4; endAverage结果如下：Average = 94.2500 91.0000 91.25007、在命令窗口键入实验代码及运行结果如下： p=3 -2 7 9; x=roots(p)x = 0.7513 + 1.7392i 0.7513 - 1.7392i -0.8358 q=poly(x)q =1.0000 -0.6667 2.3333 3.0000 b=1.5 -2 5; polyval(p,b)ans = 25.1250 -37.0000 369.00008、在命令窗口键入实验代码及运行结果如下： A=1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4; p=poly(A); poly2sym(p) a