- 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)