可视化计算离线作业讲解

浙江大学远程教育学院《可视化计算》课程作业2016年（夏）姓名： 贺鹏学 号： 716128012033年级：16秋专本计算机学习中心：浙大校内直属学习中心（华家池）———————————————————————————— 注意：所有图像的标题必须呈现足够你本人信息1.（10分）求解下列线性方程组的解： 1620908356215.87320332231074445.06.337925.1=-++=++++=++-+-=--+=++++z y v u z y x v u z y x v u z y x v z y x v u2.（10分）信号y = 5*sin(pi*20*t)+3*cos(2*pi*50*t)幅度为1的白噪声的干扰，请画出此信号，然后进行傅立叶变换，画出变换后的频域信号。

clear;title('贺鹏716128012033')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('716128012033')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 onquiver(x,y,dx,dy)title('716128012033')4.（10分）仿照课本第11章的太阳｜地球｜月亮｜卫星，绕转动画实例，呈现地球绕太阳运转的动画。

Clear;figure('name','地球绕太阳转动');title('贺鹏 716128012033')s1=[0:.01:2*pi];hold on;axis equal;axis offr1=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 1set(p1,'xdata',r1*cos(sita1),'ydata',r1*sin(sita1));sita1=sita1+w1*pausetime;pause(pausetime);drawnowEnd5.（10分）设计一个低通滤波器，从混合信号：x(t)=sin(2*pi*10*t) + cos(2*pi*100*t) + 0.2*randn(size(t))中获取10Hz的信号（10分）。

clear;figure('name','贺鹏716128012033');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 onplot(t,y,'r','MarkerSize',10)legend('input','Output')6.（20分）设计一个程序，应用函数subplot(1,2,1)、subplot(1,2,2)分别显示您本人的二张照片，然后对二张照片分别进行傅立叶变换，并分别画出变换后的频域信号。

再把2个频域信号相加，经傅立叶逆变换后，显示时域信号的图像。

clear;A=imread('C:\1.jpg');subplot(1,2,1)imshow(A);title('贺鹏716128012033')B=imread('C:\2.jpg');subplot(1,2,2)imshow(B);title('贺鹏716128012033')A1=fft2(A);B1=fft2(B);subplot(1,2,1);imshow(A1) title('贺鹏716128012033') subplot(1,2,2);imshow(B1) title('贺鹏716128012033')频域信号相加后经傅立叶逆变换C1=A1+B1;C=ifft2(C1);imshow(C)xlabel('716128012033');ylabel('贺鹏');7.（30分）小论文根据工作中的实际需要，请设计一个实际工程问题的可视化。

可以选择以下之一：（1）工程动画的可视化；（2）大数据处理中的可视化；（3）算法与模型计算的可视化；（4）实际生产流程的可视化；（5）或其它有创新意义的可视化科学计算。

要求：（1）题目有实际意义。

（2）有分析、算法描述（3）程序源代码设计。

（4）问题结果有可视化显示。

（4）题目的问题有一定的新意。

小论文的字数不能少于2500字，格式由下列各部分组成：Matlab在公司中如何用来计算工资摘要：随着社会的不断发展，许多大型的企业也在迅速崛起，企业员工也从以前的几十个人发展到成百上千，乃至数万人。

可想而知，单纯的靠人力去计算每个人都每月的实发工资肯定是不现实的。

最近学习了Matlab以后，发现它不仅仅能编写一些简单的程序，更可以把许多数据以图像的方式呈献给人们，这对于企业的管理以及工资的统计工作有着巨大的帮助。

因为学习时间也不是很久，所以我在根据已学的知识基础上，参考了一些其他书籍，集合书本上他人的设计经验，自己写了一个小程序，可以很快的计算出一个人每月的实际工资。

中文关键词：图像统计工资英文题目：Matlab in the company how to calculate salary英文摘要：With the continuous development of society, also is in rapid rise, many large enterprises from the dozens of personal development to enterprise employees also have had hundreds of thousands, even tens of thousands of people. Predictably, the pure by manpower to calculate each monthly net salary is not realistic. Recent learning Matlab, found that it not only can write some simple program, can put a lot more people present in the form of image data, for the enterprise management with the help of the statistical work and wages are huge. Because learning time is not long, so I am to learn knowledge, based on reference to some other books, a collection of books on the design experience of others, oneself write a small program, can quickly calculate the one real wages each month. 英文关键词：mage statistics wages1.引言Matlab是一套高性能的数值计算和可视化的工具软件。

也是一套十分强大的计算机科学计算、可视化辅助设计和可视化的教学软件。

它的应用十分的广泛，包括科学计算、信号和图像处理、通讯、控制系统计算、测试和测量、财务建模和分析，以及计算生物学等众多科学领域。

日常，我们计算一些数据量小，数据个数不多的文件，可能人工计算也能很快的得出结果，但是在大数据面前，如：每个月的工资统计。

钥匙还是使用人力计算，就需要大量的人和时间，还存在着一定的错误率，这明显是十分没有效率的；但是用Matlab设计一个简单的程序，就能解决所有的问题，更能使正确率达到100%。