可视化计算离线作业资料.doc

精品文档 精品文档 PAGE 精品文档 浙江大学远程教育学院 《可视化计算》课程作业 2016年（夏） 姓名： 学 号： 年级： 学习中心： ————————————————————————————— 注意：所有图像的标题必须呈现足够你本人信息 （10分）求解下列线性方程组的解： u1.5v 2x 9y 7z3 3.6v 0.5x 4y 4z 4 7u 10v 3x 22y 33z 20 3u 7v 8.5x 21y 6z 5 3u 8v 90y 20z 16 答案： A=[11.5297;03.60.5-4-4;710-32233;378.5216;38090-20]B=[3-420516] X=inv(A)*B X=3.5056-0.8979-0.27450.14380.0137 2.（10分）信号y=5*sin(pi*20*t)+3*cos(2*pi*50*t) 幅度为1的白噪声的干 扰，请画出此信号，然后进行傅立叶变换，画出变换后的频域信号。 节欖讪匦執紀镳鳧区謳簖较猻鎔鈦。 答案： clear; 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) xlabel( '变换前信号') Y=fft(y,512) subplot(1,2,2) Y=real(Y) plot(Y(1:512)) 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) 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.5,1,3,5,8,12]) hold on quiver(x,y,dx,dy)秃卢簽骅叢阀幣鋸鋨癱钻營霽树铿。 图像： （10分）仿照课本第11章的太阳｜地球｜月亮｜卫星，绕转动画实例，呈现地球绕太阳运转的动画。 答案： figure( 'name' ,'μ??òè?ì???×a?ˉ'); 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'  ,[-2020],  'ylim'  ,[-2020]); p1=plot(r1*cos(sita1),r1*sin(sita1),  'color'  ,'b'  ,'marker'  ,'.'  ,'markers ize'  ,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=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/