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

浙江大学远程教育学院 《可视化计算》课程作业2016年（夏） 姓名： 学 号： 年级： 学习中心： ————————————————————————————— 注意：所有图像的标题必须呈现足够你本人信息 （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 （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(变换后信号) ： 答案： 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) 图像： 答案： 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 图像： 答案： 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) 像： 答案： clear; I=imread(C:\Users\tu\Desktop\×÷òμ\2.jpg); J=imread(C:\Users\tu\Desktop\×÷òμ\1.jpg); [m1,n1]=size(I); [m2,n2]=size(J); I1=(I); J1=(J); subplot(1,2,1); imshow(I1); subplot(1,2,2); imshow(J1); 像： ?二张照片分别进行傅立叶变换： FFT_I = fft2(double(I1)); FFT_J = fft2(double(J1)); subplot(1,2,1); i