2016年（夏）浙江大学远程教育学院《可视化计算》课程作业附答案

《2016年（夏）浙江大学远程教育学院《可视化计算》课程作业附答案》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、浙江大学远程教育学院《可视化计算》课程作业2016年（夏）姓名：学号：年级：2017春学习中心：—————————————————————————————注意：所有图像的标题必须呈现足够你本人信息1.（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'/A1.（10分）信号y=5*sin(pi*20*t)+3*cos(2*pi*50*t)幅度为1的白噪声的干扰，请画出此信号，然后进行傅立叶变换，画出变换

2、后的频域信号。clear;title('袁洋2')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('2')xlabel('变换后的频域信号')1.（10分）在空间有一个带正电的点电荷，请画出此点电荷的空间电位分布与电场强度的空间分布图。clc;subplot(1,2,1)

3、[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('袁洋2')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]);holdonquiver(x,y,dx,dy)title('2')1.（10分）仿照课本第11章

4、的太阳｜地球｜月亮｜卫星，绕转动画实例，呈现地球绕太阳运转的动画。clc;figure('name','地球绕太阳转动');title('袁洋2')s1=[0:.01:2*pi];holdon;axisequal;axisoffr1=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','m

5、arker','.','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','.','markersize',30);while1set(p1,'xdata',r1*cos(sita1),'ydata',r1*s

6、in(sita1));sita1=sita1+w1*pausetime;pause(pausetime);drawnowend1.（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);p

7、lot(t,x,'b-')title('袁洋2提取的10HZ信号');holdonplot(t,y,'r.','MarkerSize',10)legend('Input','Output')1.（20分）设计一个程序，应用函数subplot(1,2,1)、subplot(1,2,2)分别显示您本人的二张照片，然后对二张照片分别进行傅立叶变换，并分别画出变换后的频域信号。再把2个频域信号相加，经傅立叶逆变换后，显示时域信号的图像。clc;A=imread('D:Me正常血压范围.jpg');subplot(1,2,1)ims