可视化计算离线作业2015冬解析.pdf

浙江大学远程教育学院 《可视化计算》课程作业 2015 年（冬） 姓名： 夏林 学 号： 715099012001 年级： 2015 年秋 学习中心： 洞头学习中心 ————————————————————————————— 注意：所有图像的标题必须呈现足够你本人信息 1. （10 分）求解下列线性方程组的解： u 1.5v 2x 9y 7z 3 3.6v 0.5x 4 y 4z 4 7u 10v 3x 22 y 33z 20 3u 7v 8.5x 21y 6z 5 3u 8v 90y 20z 16 2. （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(' 变换后的频域信号 ') 图 像 3. （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') 4. （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','.','markers ize',30); while 1 set(p1,'xdata',r1*cos(sita1),