西安建筑科技大学考试试卷A.DOC

PAGE PAGE 1 西安建筑科技大学考试试卷 考试科目： 数学实验 考试时间： 2012年11月 专业班级： 给水1101 学号 110320120学生姓名： 周博 题号 一 二 三 总分 分数 一、综合实验题（每小题4分，共60分） 选取一幅个人的生活照片（若分辨率比较大，可裁小，但必须包含个人完整的正面面部），将照片文件复制到work文件夹或MATLAB有效路径下，并重命名为不含空格的本人的拼音名字，使用MATLAB完成下列实验： (1). 写出显示该图像的命令，并将所显示的图像复制到本试卷上。提示：参考命令格式imshow(*.jpg)。 imshow(heliheng.jpg) (2). 上述彩色图像可转化为灰度图像，将灰度图像对应的矩阵赋值为某个变量(不含空格的本人的拼音名字), 不需要输出此变量。提示：参考命令格式 double(rgb2gray(imread(*.jpg))) heliheng= double(rgb2gray(imread(heiheng.jpg))) (3). 求该矩阵的行数、列数、最大值、最小值、平均值的命令。 [heliheng]=double(rgb2gray(imread(heliheng.jpg))); [i,j]=size(heliheng) max(max(heliheng)) min(min(heliheng)) mean(heliheng(:)) (4). 求该矩阵元素值大于等于200的元素个数。 sum(sum(heliheng=200)) (5). 求该矩阵的秩，并输出该矩阵转置后的图像。提示：imagesc(*)或imshow(uint8(*))。 rank(heliheng) (6). 用255减去该矩阵，并显示新矩阵对应的图像。 (7). 该矩阵乘以它的转置得到另一矩阵，求所得矩阵的行列式和所有的特征值。 [a]=ctranspose(heliheng); [b]=a*heliheng; det(b) [m n]=eig(b) (8). 对该矩阵的第一列和最后一列绘制两条曲线，横坐标为从1开始的整数，x轴上标注“row”，y轴上标注“intensity”，标题为自己的名字，对第一条、第二条曲线分别标注“First Column”和“Last Column”。需要写出命令和输出图形。 [c]=heliheng(:,1); [d]=heliheng(:,335); [e]=[1:1:542]; plot(e,d,r); hold plod(e,c, b) (9). 提取该矩阵的前3列，它们可看做3维空间中的点，绘出这些离散点，不需要直线段相连。 (10). 对该矩阵加均值为0、标准差为10的高斯噪声，显示加噪后的图像。提示：10*randn(2,3)为2x3的均值为0、标准差为10的高斯噪声矩阵 (11). 对该矩阵做空间曲线或曲面。提示：使用mesh(X,Y,Z)或surf(X,Y,Z), Z对应该矩阵，X，Y需要根据行数和列数来确定。 (12). 使用rand命令随机生成和该矩阵同样大小的随机矩阵，若随机矩阵的第i行第j列元素大于等于0.8，则原图像矩阵的第i行第j列元素赋值为255。写出命令，并显示新图像。在使用rand前，要求先执行rand(state,？)，其中“？”为本人学号。 rand(state,110330206) a=rand(size(heliheng)) heliheng(a=0.8)=225 imshow((heliheng)) (13). 将原矩阵大致等分为4块，使用subplot显示各块图像。 (14). 对原矩阵选取奇数行、奇数列，显示子矩阵的对应的图像。 (15). 对上述子矩阵进行2维网格插值，得到和原图像大小相同的矩阵，并显示此矩阵对应的图像。 二、解答题 （前两小题各5分，第3小题20分，共30分） (1). 使用数值积分求, 其中，l=dblquad(@(x,y)sin(x.^2+y.^2+eps)./(x.^2+y.^2+eps).*(x.^2+3*y.^2=4), -2,2,-sqrt(4/3),sqrt(4/3)) l = 4.9082 (2). Henon吸引子是混沌和分形的著名例子，迭代模型为: 取初值, 进行N次迭代，对于, 在处画一点(注意不要连线)可得所谓Henon引力线图，其中N为本人的学号。试写出程序，画出图形。x(1)=0;y(1)=0; k=110330206; for n=2:k x(n)=1+0.3*y(n-1)-1.4*(x(n-1))^2; y(n)=x(n-1); end plot(x(1001:end),y(1001:e