可视化计算离线作业_716060014001卢子强.doc

浙江大学远程教育学院 《可视化计算》课程作业2016年（夏） 姓名： 卢子强 号： 716060014001 年级： 2016春 学习中心： 深圳 716060014001) s1=[0:.01:2*pi]; hold on;axis equal; axis off r1=10; w1=1; t=0; pausetime=.002; sita1=0;sita2=0; set(gcf,doublebuffer,on) plot(-20,18,color,r,marker,.,markersize,40); text(-17,18,太阳); p1=plot(-20,16,color,b,marker,.,markersize,20); text(-17,16,地球); plot(0,0,color,r,marker,.,markersize,60); plot(r1*cos(s1),r1*sin(s1)); set(gca,xlim,[-20 20],ylim,[-20 20]); p1=plot(r1*cos(sita1),r1*sin(sita1),color,b,marker,.,markersize,30); while 1 set(p1,xdata,r1*cos(sita1),ydata,r1*sin(sita1)); sita1=sita1+w1*pausetime; pause(pausetime); drawnow （10分）设计一个低通滤波器，从混合信号： 　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) xlabel(第5题); （20分）设计一个程序，应用函数subplot(1,2,1)、subplot(1,2,2)分别显示您本人的二张照片，然后对二张照片分别进行傅立叶变换，并分别画出变换后的频域信号。再把2个频域信号相加，经傅立叶逆变换后，显示时域信号的图像。 源程序： clear; ben=imread(ben.jpg); ben2=imread(ben2.jpg); figure(1); subplot(1,2,1) imshow(ben); title(照片1-); subplot(1,2,2) imshow(ben2); title(照片2-); F1=fft2(ben); FS1=fftshift(F1); F2=fft2(ben2); FS2=fftshift(F2); figure(2); subplot(1,2,1) imshow(FS1,[]);title(傅立叶变换之后1-); subplot(1,2,2) imshow(FS2,[]);title(傅立叶变换之后2-); k=min(min(FS1),min(FS2)):max(max(FS1),max(FS2)); s1=zeros(1,length(k));s2=s1s1(find((k=min(FS1))(k=max(k1))==1))=FS1; s2(find((k=min(FS2))(k=max(k2))==1))=FS2; y=s1+s2; z=ifft(y); imshow(z);title(逆变换后图像); 图形： （30分）小论文 　　 根据工作中的实际需要，请设计一个实际工程问题的可视化。可以选择以下之一：（1）工程动画的可视化；（2）大数据处理中的可视化；（3）算法与模型计算的可视化；（4）实际生产流程的可视化；（5）或其它有创新意义的可视化科学计算。要求： （1）题目有实际意义。 （2）有分析、算法描述 （3）程序源代码设计。 （4）问题结果有可视化显示。 （4）题目的问题有一定的新意。 小论文的字数不能少于2500字，格式由下列各部分组成： 中文题目 摘要： 中文关键词： 英文题目： 英文摘要： 英文关键词： 引言 算法基础 程序代码 结果分析（必须有可视化图） 结论 参考文献 基于MATLAB关于交混回响 的仿真研究 摘要： MATLAB是进行科学研究和数值分析的首选软件，本次对于交混回响的研究正是借助MATLAB的强大功能，通过编写程序实现了模