自动化18级MATLAB大作业.docVIP

常州大学大作业 科目 MATLAB程序设计 成绩 卷 共 12 页，第 PAGE 1 页 PAGE 1 PAGE 1 系（部） 自动化 教研室 2019 年 10 月 28 日 ……………………………………………………………………………………………………………………………… 装 订 线 …………………………………………………………………………………………………………………………tc ……………………………………………………………………………………………………………………………… 装 ……………………………………………………………………………………………………………………………… 装 订 线 …………………………………………………………………………………………………………………………tc ……………………………………………………………………………………………………………………………… 装 订 线 ………………………………………………………………………………………………………………………… 班级 自动化18 学号 姓名 。 （注意事项，作业中有些代码中字符是全角字符，这在MATLAB中是非法的，为了识别是否为半角字符，可将有关代码粘贴至脚本程序中观察其顔色； 或者你可手工以半角字符英文方式重新输入有关代码。 蓝色部份是参考性知识介绍或方法应用背景） 给定长度为9的行向量number,其各元素顺次存放你自己的学号各位(例如若你学号为164630528,则number＝[1,8,4,6,3,0,5,2,8]) 正定矩阵 二阶正定矩阵的特征值和特征向量 1.考虑矩形平板在其所处平面内两个相互垂直方向上的拉伸引起的均匀变形产生的坐标变换： 或者写成 这里和分别为变形前和变形后的坐标。经此变形（线性变换），原来的单位圆 被变换为椭圆 若将坐标轴逆时针旋转过角度，令变形前后的坐标分别为和，，则可推知变形前后的坐标关系为 可以推出， 其中，是的逆矩阵。 假定将角度换为，编写符号计算程序并运行之，计算出相应的矩阵的特征值和特征向量。 （显然，虽然从另一方向观察变形情况，原来的圆仍旧是变换为长半轴和短半轴分别为的椭圆，这正是的两个特征值，也是的两个特征值； 观测方向虽然不同，观测到的变形却应是相同的；和之间有内在联系，它们不过是同一变换在不同正交坐标系下的表示形式而已） 数组的生成和操作 2. number=[用空格分隔的你的9个学号数字]; A=[number(1),number(3)]; B=[number(8), number(9)]; C=[B; A]; D=[A, B]; 问： 程序运行后 C, D各为何？ 3.. rand(state,65868); dx=rand(1,9); dy=90*sin(dx); dz=sum(dy); number=[用空格隔开的你的9个学号数字]; v=number+dz; 问： 程序运行后, dz, v 各为何？ 4.clear; clc; number=[用空格隔开的你的9个学号数字]; V=mean(number); A = magic(3)+V; [vals, rows]=min(A); [val, col]=min(vals); row=rows(col) 问： 观察运行结果，并说明程序运行后 row, , col,val 为何？ 5. clear; clc; number=[用空格隔开的你的9个学号]; a=number(3:6)+1; b= -exp(1./a).’; c=a.^3; 程序运行之后，a, b,c 各为何 6. 已知程序： clc; clear; number=[用空格隔开的你的9个学号]; odds=1:2:length(number); b(odds)=number(1:2:end); 问程序运行之后，b为何？ 7. 用eig函数计算矩阵　A=number(end-2:end).*number(4:6); 　的特征值和它们各自所对应的特征向量; 用两个返回参数：[v d]=eig(A) A=ones(3,3)+eye(3); B=A(A1); 运行后，B为何？ 8.clear; clc; number=[用空格隔开的你的