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

浙江大学远程教育学院 《可视化计算》课程作业（夏） 姓名： 学 号： 年级： 15年春 学习中心： ————————————————————————————— 应用附件中的软件“画图”制作照片大小的灰度为255级的图像，图像的内容为您的名字，再读取您本人的照片进行合成，请用MATLAB编写程序，显示原先的两幅图像，最后显示合成后的图像（10分）。 答： I1=imread(C:\Documents and Settings\Administrator\My Documents\My Pictures \姓名.bmp); I2=imread( C:\Documents and Settings\Administrator\My Documents\My Pictures \照片.JPG); I3= imread( C:\Documents and Settings\Administrator\My Documents\My Pictures \合成.JPG); subplot(1,3,1);imshow(I1);title(姓名); subplot(1,3,2);imshow(I2);title(照片); subplot(1,3,3);imshow(I3);title(合成); 对第一题中合成后的图像，进行灰度变换，把灰度值少于55变换为0，把灰度值大于204变换为255，再把灰度值为55－204之间映射到0－255之间（10分）。 答： I3= imread(C:\Documents and Settings\Administrator\My Documents\My Pictures\合成.bmp); I4=imadjust(I3,[55/255 204/255],[0,1]); subplot(1,2,1); imshow(I3); title(合成); subplot(1,2,2); imshow(I4); title(灰度变换后) 信号y = sin(pi*50*t)+3*cos(2*pi*100*t)幅度为1的白噪声的干扰，请画出此信号，然后进行傅立叶变换，画出变换后的频域信号（10分）。 答： clear t=0:0.001:0.6; y=sin(2*pi*25*t)+3*cos(2*pi*100*t); y=y+randn(1,length(t)); subplot(2,2,1); plot(y(1:100)); xlabel(变换前的信号); Y=fft(y,256); subplot(2,2,2); Y=real(Y); plot(Y(1:256)); xlabel(变换后频域信号) Pyy=Y.*conj(Y)/256; f=1000*(0:255)/256; subplot(2,2,3); plot(Pyy(1:256)); 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=sin(2*pi*10*t)+cos(2*pi*100*t)+0.2*randn(size(t)); wn=ws/2; [B,A]=butter(10,30/wn); y=filter(B,A,x); plot(t,x,b-) hold on plot(t,y,r.,MarkerSize,10) legend(Input,Output) 编写一个MATLAB程序，读取一个音频数据文件，以不同频率播放（频率由键盘读入），并显示声音波形（10分）。 答： clc close all; y=wavread(C:\Documents and Settings\Administrator\My Documents\My Pictures\Windows XP 登录音.wav); plot(y); wavplay(y,44100); wavplay(y,22050); wavplay(y,11025); 应用MATLAB求下列两组数据的y1=[60 90 100 34 80 ];y2=[60 70 90 80 85]如求下列样本的样本方差和样本标准差、方差和标准差（10分）。 答： y1=[60 90 100 34 80 ]; y2=[60 70 90 80 85]; DY_1=var(y1, 1) %y1方差 sigma_1=std(y1, 1) %y1标准差 DY1_1=var(y1) %y1样本方差 sigma1