可视化计算离线作业解答.doc

浙江大学远程教育学院 《可视化计算》课程作业2015年（冬） 姓名： 学 号： 年级： 学习中心： ————————————————————————————— 注意：所有图像的标题必须呈现足够你本人信息 （10分）求解下列线性方程组的解： 解：A=[1 0 7 3 3;1.5 3.6 10 7 8;2 0.5 -3 8.5 0;9 -4 22 21 90; 7 -4 33 6 -20] B=[3;-4;20;5;16] X=B/A A = 1.0000 0 7.0000 3.0000 3.0000 1.5000 3.6000 10.0000 7.0000 8.0000 2.0000 0.5000 -3.0000 8.5000 0 9.0000 -4.0000 22.0000 21.0000 90.0000 7.0000 -4.0000 33.0000 6.0000 -20.0000 B = 3 -4 20 5 16 X = 3.5056 -0.8979 -0.2745 0.1438 0.0137 （10分）信号y = 5*sin(pi*20*t)+3*cos(2*pi*50*t)幅度为1的白噪声的干扰，请画出此信号，然后进行傅立叶变换，画出变换后的频域信号。 clear t=(0:0.001:1); y = 5*sin(2*pi*10*t)+3*cos(2*pi*50*t); plot(y(1:150)) title(); 傅立叶转变后： clear t=0:0.001:0.6; y=5*sin(2*pi*10*t)+3*cos(2*pi*50*t); plot(y(1:150)); y=fft(y,512); Y=real(y); figure,plot(Y(1:512)) title( ) （10分）在空间有一个带正电的点电荷，请画出此点电荷的空间电位分布与电场强度的空间分布图。 hh=figure(numbertitle,off,name,点电荷的空间电位分布与电场强度的空间分布图2015年秋 ); r0=0.12; th=linspace(0,2*pi,13); [x,y]=pol2cart(th,r0) x=[x;0.1*x] y=[y;0.1*y] plot(x,y) grid on hold on plot(0,0,o,MarkerSize,12) k=9e9; q=1e-9; r0=0.12; u0=k*q/r0; u=linspace(1,3,7)*u0 x=linspace(-r0,r0,100) [X,Y]=meshgrid(x) r=sqrt(X.^2+Y.^2); U=k*q./r; contour(X,Y,U,u) xlabel(x,fontsize,16) ylabel(y,fontsize,16) title( ) （10分）仿照课本第11章的太阳｜地球｜月亮｜卫星，绕转动画实例，呈现地球绕太阳运转的动画。 clear; clc; close all; %相对应的圆点坐标 半径等 x0=0;y0=0;r0=80;Lmin0=0;Lmax0=0;T0=2160;w0=0*pi/T0;q0=0; x1=0;y1=0;r1=40;Lmin1=25;Lmax1=30;T1=1080;w1=pi/T1;q1=0; %初始化 hh=figure(numbertitle,off,name,地球 | 太阳转演示动画 余合哲); %设置擦除方式 sun = line(0,0,color,r,linestyle,.,erasemode,xor,markersize,r0);%太阳 earth = line(x0,y0,color,k,linestyle,.,erasemode,xor,markersize,r1);%地球 %添加标注 axis off title(太阳 | 地球 ,fontname,宋体,fontsize,9,FontWeight,demi,color,black); text(-50,50,太阳);%对太阳进行标识 line(-55,50,color,r,marker,.,markersize,80); text(-50,40,地球);%对地球进行标识 line(-55,40,color,k,marker,.,markersize,40); %绘制轨道 s1=[0:.01:2*p