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

1、浙江大学远程教育学院可视化计算课程作业（夏）姓名：王春强学 号：年级：2016春学习中心：桐庐学习中心1. 应用附件中的软件“画图”制作照片大小的灰度为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);%转

2、换成双精度数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(合成图)截图如下：2. 对第一题中合成后的图像，进行灰度变换，把灰度值少于55变换为0，把灰度值大于204变换为255，再把灰度值为55204之间映射到0255之间（10分）。答：程序如下：A1=im

3、read(D:COM.tif);%读取原来合成的图B1=imadjust(A1,0.22 0.8,0 1);%55/255=0.22,204/255=0.8，0 1表示映射到0-255subplot(2,1,1);imshow(A1)xlabel(原图)subplot(2,1,2);imshow(B1)xlabel(灰度值204变换为255，55204之间映射到0255间)截图如下：3. 信号y = sin(pi*50*t)+3*cos(2*pi*100*t)幅度为1的白噪声的干扰，请画出此信号，然后进行傅立叶变换，画出变换后的频域信号（10分）。答：cleart=0:0.001:0.3;%根

4、据信号函数，选取了0到0.3y = 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(频域信号)截图如下：4. 设计一个低通滤波器，从混合信号x(t)=sin(2*pi*10*t)+cos(2*pi*100*t)+0.2*randn(size(t)中获取10Hz的信号（10分）。答：程序如下：clear;

5、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 onplot(t,y,r.,MarkerSize,10)legend(Input,Output)信号截图如下：5. 编写一个MATLAB程序，读取一个音频数据文件，以不同频率播放（频率由键盘读入），并显示声音波形（10分）。答：程序如下：clc;close all;y=wavread(d:

6、再回首.wav);%再回首铃声，用MP3转换而来plot(y);p=input(请输入播放频率(输入0表示关闭功能）：)while(p=0) %功能是：输入0则结束这项程序，否则进行对应频率播放disp(当前播放频率为：);disp(p);wavplay(y,p);p=input(请输入播放频率(输入0表示关闭功能）：)enddisp(输入频率功能已关闭，可重新运行程序开始);测试截图如下：声音波形如下：6. 应用MATLAB求下列两组数据的y1=60 90 100 34 80 ;y2=60 70 90 80 85如求下列样本的样本方差和样本标准差、方差和标准差（10分）。答：程序如下：clc

7、;y1=60 90 100 34 80 ;y2=60 70 90 80 85;Dy11=var(y1)%y1样本方差Sy11=std(y1)%y1样本标准差Dy1=var(y1,1)%y1方差Sy1=std(y1,1)%y1标准差Dy22=var(y2)%y2样本方差Sy22=std(y2)%y2样本标准差Dy2=var(y2,1)%y2方差Sy2=std(y2,1)%y2标准差计算结果如下：（按样本方差、样本标准差、方差、标准差顺序）Dy11 = 689.2000Sy11 = 26.2526Dy1 = 551.3600Sy1 = 23.4811Dy22 = 145Sy22 = 12.0416Dy2 = 116Sy2 = 10.7703截图如下：7. 当两个互相垂直的简谐振动频率不同时，合成的轨迹与频率之比和两者的相位都有关系，请用MATLAB进行数值模拟及数值用动态图形仿真后的李萨如图形的动态运动过程（10分）。答 程序如下：clc;clear all;t=linspace(0,2*pi,200);m=1;for n=1:1:3for k=-100:1:100X=3*cos(n*t+pi*k/100)+8*i*sin(t);subplot(2,2,m);plot(X);xlabel(T1:T2=1:n);pause(0.05