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

浙江大学远程教育学院 《可视化计算》课程作业2015年（冬） 姓名： 学 号： 年级： 学习中心： ————————————————————————————— 注意：所有图像的标题必须呈现足够你本人信息 （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.505571 -0.897945 -0.274501 0.143796 0.013740 （10分）信号y = 5*sin(pi*20*t)+3*cos(2*pi*50*t)幅度为1的白噪声的干扰，请画出此信号，然后进行傅立叶变换，画出变换后的频域信号。 cleart=0:0.001:0.6;y=5*sin(2*pi*10*t)+3*cos(2*pi*50*t);plot(y(1:100)); xlabel(变化信号) y=fft(y,512)subplot(2,2,2) y=real(y) plot(y(1:512)) 转换后 （10分）在空间有一个带正电的点电荷，请画出此点电荷的空间电位分布与电场强度的空间分布图。 （10分）仿照课本第11章的太阳｜地球｜月亮｜卫星，绕转动画实例，呈现地球绕太阳运转的动画。 clear;?clc;close?allx0=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,太阳｜地球，绕转演示动画——matlabfan);sun=line(0?,0?,color,r,linestyle,.,erasemode,xor,markersize,r0);?%太阳earth=line(x0,y0,color,k,linestyle,.,erasemode,xor,markersize,r1);?%地球axis?offtitle(太阳｜地球,fontname,宋体,fontsize,9,FontWeight,demi,Color,black);text(-20,50,?);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*pi];line(Lmax1*cos(s1),Lmin1*sin(s1),linestyle,:);?axis([-60,60,-60,60]);????t?=0;while?1if?~ishandle(hh),return,endq0=t*w0;?q1=t*w1;?t=t+1;??%设置运动规律if?t?=?4320;?t?=?0;?end??%到了一个周期就重置x0?=?Lmax0?*?cos(q1);?y1?=?Lmin0?*?sin(q1);?x1?=?x0?+?Lmax1?*?cos(q1);?y1?=?y0?+?Lmin1?*?sin(q1);?%设置地球圆心的坐标set(sun,xdata,x0,ydata,y0);????%画太阳set(earth,xdata,x1,ydata,y1);????%画地球drawnow;title(?715129012016)end （10分）设计一个低通滤波器，从混合信号： 　x(t)=sin(2*pi*10*t)　+　cos(2*pi*100*t)　+　0.2*randn(size(t)) 　中获取10Hz的信号（10分）。 clearws=1000t=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/2y=filter(B,A,x);plot(t,x,b-)hold?onplot(t,y,r,MarkerSize,10)legend(Input,Output) （20分）设计一个程序，应用函数subplot(1,2,1)、subplo