浙江大学远程教育2016年(夏)可视化计算离线作业答案.docVIP

浙江大学远程教育学院 《可视化计算》课程作业（夏） 姓名： 王春强 学 号： 年级： 2016春 学习中心： 桐庐学习中心 ————————————————————————————— 应用附件中的软件“画图”制作照片大小的灰度为255级的图像，图像的内容为您的名字，再读取您本人的照片进行合成，请用MATLAB编写程序，显示原先的两幅图像，最后显示合成后的图像（10分）。 答： 程序如下： A=imread(D:\name.tif);%灰度为255级的名字图片不太理解，直接用画图做了个 B=imread(D:\myphoto.jpg); SA=size(A); m=SA(1),n=SA(2); B1=imresize(B,[m n]);%将数组B变为和A一样大小 A1=double(A);%转换成双精度数 C=double(B1); D=A1.*C/128; COM=uint8(D); subplot(2,2,1);imshow(A) xlabel(姓名) subplot(2,2,2);imshow(B) xlabel(我的照片) subplot(2,2,3);imshow(COM) imwrite(COM,D:\COM.tif,tif);%为了完成第二题，导出合成照到D盘，以tif格式保存 xlabel(合成图) 截图如下： 对第一题中合成后的图像，进行灰度变换，把灰度值少于55变换为0，把灰度值大于204变换为255，再把灰度值为55－204之间映射到0－255之间（10分）。 答： 程序如下： A1=imread(D:\COM.tif);%读取原来合成的图 B1=imadjust(A1,[0.22 0.8],[0 1]);%55/255=0.22,204/255=0.8，[0 1]表示映射到0-255 subplot(2,1,1);imshow(A1) xlabel(原图) subplot(2,1,2);imshow(B1) xlabel(灰度值55变换为0，204变换为255，55－204之间映射到0－255间) 截图如下： 信号y = sin(pi*50*t)+3*cos(2*pi*100*t)幅度为1的白噪声的干扰，请画出此信号，然后进行傅立叶变换，画出变换后的频域信号（10分）。 答： clear t=0:0.001:0.3;%根据信号函数，选取了0到0.3 y = sin(pi*50*t)+3*cos(2*pi*100*t); subplot(2,1,1);%设定成2行1列 plot(y(1:100));%绘制前100个点的图像 xlabel(原信号图); bianhuan=fft(y,500); Y=real(bianhuan(1:500)); subplot(2,1,2); plot(Y(1:500)); xlabel(频域信号) 截图如下： 设计一个低通滤波器，从混合信号x(t)=sin(2*pi*10*t)+cos(2*pi*100*t)+0.2*randn(size(t))中获取10Hz的信号（10分）。 答： 程序如下： clear; ws=1000; t=0:1/ws:0.4; x=3*sin(2*pi*10*t)+2*cos(2*pi*100*t)+0.8*randn(size(t)); wn=ws/2; [B,A]=butter(10,30/wn); y=filter(B,A,x); plot(t,x,b-) title(提取的10HZ信号); hold on plot(t,y,r.,MarkerSize,10) legend(Input,Output) 信号截图如下： 编写一个MATLAB程序，读取一个音频数据文件，以不同频率播放（频率由键盘读入），并显示声音波形（10分）。 答： 程序如下： clc; close all; y=wavread(d:\再回首.wav);%再回首铃声，用MP3转换而来 plot(y); p=input(请输入播放频率(输入0表示关闭功能）：) while(p~=0) %功能是：输入0则结束这项程序，否则进行对应频率播放 disp(当前播放频率为：); disp(p); wavplay(y,p); p=input(请输入播放频率(输入0表示关闭功能）：) end disp(输入频率功能已关闭，可重新运行程序开始); 测试截图如下： 声音波形如下： 应用MATLAB求下列两组数据的y1=[60 90 100 34 80 ];y2=[60 70 90 80 85]如求下列样