可视化计算离线作业解决方案.doc

  1. 1、本文档共16页,可阅读全部内容。
  2. 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
浙江大学远程教育学院 《可视化计算》课程作业2016年(夏) 姓名: 学 号: 年级: 学习中心: ————————————————————————————— 注意:所有图像的标题必须呈现足够你本人信息 1.(10分)求解下列线性方程组的解: 答案: 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=inv(A)*B X=3.5056 -0.8979 -0.2745 0.1438 0.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(变换后信号) 图像: 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) 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) 图像: 4.(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,[-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 图像: 5.(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) 图像: 6.(20分)设计一个程序,应用函数subplot(1,2,1)、subplot(1,2,2)分别显示您本人的二张照片,然后对二张照片分别进行傅立叶变换,并分别画出变换后的频域信号。再把2个频域信号相加,经傅立叶逆

文档评论(0)

希望之星 + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档