原创力文档- 1、本文档共16页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
浙江大学远程教育学院
《可视化计算》课程作业2016年（夏）
姓名：
学 号：
年级：
学习中心：
—————————————————————————————
注意：所有图像的标题必须呈现足够你本人信息
1.（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
2.（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(变换后信号)
图像：
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)
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)
图像：
4.（10分）仿照课本第11章的太阳｜地球｜月亮｜卫星，绕转动画实例，呈现地球绕太阳运转的动画。
答案：
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
图像：
5.（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=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)
图像：
6.（20分）设计一个程序，应用函数subplot(1,2,1)、subplot(1,2,2)分别显示您本人的二张照片，然后对二张照片分别进行傅立叶变换，并分别画出变换后的频域信号。再把2个频域信号相加，经傅立叶逆
文档评论（0）