数学学院《数字号处理》课程上机考试题.docVIP

2014数学学院《数字信号处理》课程上机考试题 说明 1，开卷考试。 2，试题中的文件均在D:\A下。所有的文件均要保存在D：\A下，保存在其余位置得0分 3，文件名统一为“班级代码”+ “座号”+“姓名”。例如：“110766张三”，表示1107班66号张三。 4，交卷方法：关闭所有打开的窗口＝＞双击桌面的“交卷.exe”＝＞在打开的对话框中选择试卷类型、选择班级、输入学号＝＞交卷。 班级： 姓名： 座号： 把lena图像转化为灰度图像。（5分） clear;clc; I=imread(F:\大三二资料\数字图像\1107-08\lena.jpg); I1=rgb2gray(I); subplot(121),imshow(I);title(原图); subplot(122),imshow(I1);title(灰度图像); 把1得到的灰度图像的大小扩大到。将图像顺时针旋转45度，显示旋转后的图像。 对扩大后的图像添加椒盐噪声图像进行傅里叶变换与余弦变换并显示它们的频谱。（10分） clear;clc; I=imread(F:\大三二资料\数字图像\1107-08\lena.jpg); I1=rgb2gray(I); I2=imresize(I1,[512,512]);%将图像扩大到512*512 I3=fftshift(fft2(I2));%计算傅里叶变换并移位 I4=dct2(I2);%使用dct2对图像进行DCT变换 subplot(231),imshow(I2);title(扩大到512*512的图像); subplot(232),imshow(I3);title(傅立叶变换的图像); subplot(233),imshow(log(abs(I3)),[]);title(傅立叶变换频谱); subplot(234),imshow(I4);title(DCT变换); subplot(235),imshow(log(abs(I4)),[]);title(DCT变换频谱); 编程实现对2得到的图像进行二次小波分解并进行重构。（10分） clear;clc; I=imread(F:\大三二资料\数字图像\1107-08\lena.jpg); I1=rgb2gray(I); I2=imresize(I1,[512,512]);%将图像扩大到512*512 [X,map]=gray2ind(I2);%将灰度图像转化为索引图像 nbcol=size(map,1); cod_X=wcodemat(X,nbcol);%返回数据矩阵I2的编码矩阵 cod_X, nbcol 伪编码的最大值，即编码范围为 0～ nbcol figure,image(cod_X); colormap(pink(nbcol)); axis(square); [ca1,ch1,cv1,cd1]=dwt2(X,db1);%对图像进行分解 cod_ca1=wcodemat(ca1,nbcol); cod_ch1=wcodemat(ch1,nbcol); cod_cv1=wcodemat(cv1,nbcol); cod_cd1=wcodemat(cd1,nbcol); figure, image([cod_ca1,cod_ch1;cod_cv1,cod_cd1]); colormap(pink(nbcol)); axis(square); [ca2,ch2,cv2,cd2]=dwt2(ca1,db1);%对图像进行二次分解 cod_ca2=wcodemat(ca2,nbcol); cod_ch2=wcodemat(ch2,nbcol); cod_cv2=wcodemat(cv2,nbcol); cod_cd2=wcodemat(cd2,nbcol); figure, image([cod_ca2,cod_ch2;cod_cv2,cod_cd2]); colormap(pink(nbcol)); axis(square); a0=idwt2(ca1,ch1,cv1,cd1,db1,size(X));%用idwt2重构图像 a0=wcodemat(a0,nbcol); figure, image(a0); colormap(pink(nbcol)); axis(square); 6、绘制1得到的灰度图像中直线路径的强度值; 绘制football.jpg图像中直线路径的强度值. 其中直线路径为{A(12,23);B(4,6);C(45,34);D(200,64)、 AB;BC;CD}（10分） clear;clc; I=imread(F:\大三二资料\数字图像\1107-08\len