2016年浙江大学远程教育学院可视化计算课程作业附的答案.docVIP

浙江大学远程教育学院 《可视化计算》课程作业2016年（夏） 姓名： 学 号： 年级： 2017春 学习中心： ————————————————————————————— 注意：所有图像的标题必须呈现足够你本人信息 （10分）求解下列线性方程组的解： clc; A=[ 1,1.5,2,9,7; 0,3.6,0.5,-4,-4; 7,10,-3,22,33; 3,7,8.5,21,6; 3,8,0,90,-20 ] B=[3;-4;20;5;16] X=B/A （10分）信号y = 5*sin(pi*20*t)+3*cos(2*pi*50*t)幅度为1的白噪声的干扰，请画出此信号，然后进行傅立叶变换，画出变换后的频域信号。 clear; title(袁洋 717129012012) 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(717129012012) xlabel(变换后的频域信号) （10分）在空间有一个带正电的点电荷，请画出此点电荷的空间电位分布与电场强度的空间分布图。 clc; 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(袁洋 717129012012) 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(717129012012) （10分）仿照课本第11章的太阳｜地球｜月亮｜卫星，绕转动画实例，呈现地球绕太阳运转的动画。 clc; figure(name,地球绕太阳转动); title(袁洋 717129012012) 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 （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-) title(袁洋 717129012012 提取的10HZ信号); hold on plot(t,y,r.,MarkerSize,10) legend(Input,Output) （20分）设计一个程序，应用函数subplot(1,2,1)、subplot(1,2,2)分别显示您本人的二张照片，