可视化计算离线作业2015冬讲述.doc

浙江大学远程教育学院 《可视化计算》课程作业2015年（冬） 姓名： 夏林 学 号： 715099012001 年级： 2015年秋 学习中心： 洞头学习中心 ————————————————————————————— 注意：所有图像的标题必须呈现足够你本人信息 （10分）求解下列线性方程组的解： （10分）信号y = 5*sin(pi*20*t)+3*cos(2*pi*50*t)幅度为1的白噪声的干扰，请画出此信号，然后进行傅立叶变换，画出变换后的频域信号。 代码：clear; title(夏林715099012001) t=0:0.001:0.6; y=5*sin(pi*20*t)+3*cos(2*pi*50*t); y=y+randn(1,length(t)); subplot(1,2,1) plot(t,y) title(夏林) xlabel(变换前信号) Y=fft(y,512); subplot(1,2,2) Y=real(Y); plot(Y(1:512)) title(715099012001) xlabel(变换后的频域信号) 图像 （10分）在空间有一个带正电的点电荷，请画出此点电荷的空间电位分布与电场强度的空间分布图。 clear; subplot(1,2,1) [x,y]=meshgrid(-4:0.2:4,-4:0.2:4); z=1./sqrt(x.^2+y.^2+0.01); mesh(x,y,z) title(夏林) subplot(1,2,2) [x,y]=meshgrid(-2:0.1:2,-2:0.1:2); z=1./sqrt((x).^2+y.^2+0.01); [dx,dy]=gradient(z); contour(x,y,z,[-12,-8,-5,-3,-1,-0.5,-0.1,0.1,0.5,1,3,5,8,12]); hold on quiver(x,y,dx,dy) title(715099012001) （10分）仿照课本第11章的太阳｜地球｜月亮｜卫星，绕转动画实例，呈现地球绕太阳运转的动画。 figure(name,地球绕太阳转动); title(夏林715099012001) 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 End 图在下面-----截的图 不会动 在matlab中是逆时针转 （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=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-) hold on plot(t,y,r,MarkerSize,10) legend(input,Output) xlabel(715099012001) ylabel(夏林) 图： （20分）设计一个程序，应用函数subplot(1,2,1)、subplot(1,2,2)分别显示您本人的二张照片，然后对二张照片分别进行傅立叶变换，并分别画出变换后的频域信号。再把2个频域信号相加，经傅立叶逆变换后，显示时域信号的图像。