浙江大学远程教育学院可视化计算离线作业.doc

浙江大学远程教育学院 《可视化计算》课程作业2015年（冬） 姓名： 学 号： 年级： 学习中心： ————————————————————————————— 注意：所有图像的标题必须呈现足够你本人信息 （10分）求解下列线性方程组的解： 源程序： clear; 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]; C=A\B C = 3.5056 -0.8979 -0.2745 0.1438 0.0137 图形： （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(2,2,1); plot(y(1:100)); xlabel(变换前的信号); x=fft(y,256); subplot(2,2,2); x=real(x); plot(x(1:256)); xlabel(变换后频域信号); （10分）在空间有一个带正电的点电荷，请画出此点电荷的空间电位分布与电场强度的空间分布图。 源程序： 程序1： clear; x=-2:0.1:2; y=-2:0.1:2; [x,y]=meshgrid(x,y); z=1./sqrt(x.^2+y.^2+0.01); mesh(x,y,z); xlabel(电荷空间分布); 程序2： clear;clf; q=2e-6; k=9e9; x=-6:0.6:6; y=x; [X,Y]=meshgrid(x,y); r = sqrt(X.^2+Y.^2); V=q*k./r; [Ex,Ey]=gradient(-V); AE=sqrt(Ex.^2+Ey.^2); Ex=Ex./AE; Ey=Ey./AE; contour(X,Y,V,r-); hold on; quiver(X,Y,Ex,Ey,.5,b) 图形： （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); plot(t,x,b-) hold on plot(t,y,r.,MarkerSize,10) legend(Input,Output) xlabel(第四题); 图形： 可视化作业第四题- （20分）设计一个程序，应用函数subplot(1,2,1)、subplot(1,2,2)分别显示您本人的二张照片，然后对二张照片分别进行傅立叶变换，并分别画出变换后的频域信号。再把2个频域信号相加，经傅立叶逆变换后，显示时域信号的图像。 源程序： clear; ben=imread(ben.jpg); ben2=imread(ben2.jpg); figure(1); subplot(1,2,1) imshow(ben); title(照片1-); subplot(1,2,2) imshow(ben2); title(照片2-); F1=fft2(ben); FS1=fftshift(F1); F2=fft2(ben2); FS2=fftshift(F2); figure(2); subplot(1,2,1) imshow(FS1,[]);title(傅立叶变换之后1-); subplot(1,2,2) imshow(FS2,[]);title(傅立叶变换之后2-); k=min(min(FS1),min(FS2)):max(max(FS1),max(FS2)); s1=zeros(1,length(k));s2=s1s1(find((k=min(FS1))(k=max(k1))==1))=FS1; s2(find((k=min(FS2))(k=max(k2))==1))=FS2; y=s1+s2; z=ifft(y); imshow(z);title(逆变换后图像); 图形： 可视化作业第五题- （30分）小论文 　　 根据工作中的实际需要，请设计一个实际工程