原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
地球物理信息处理基础实验2
地球物理信息处理基础实验二
指导老师：邢磊
姓名：申鹏
学号：12040032032
班级：2012级地球信息科学与技术
地球物理信息处理基础实验二
实验题目：
已知两个序列：X=（1,2,3,4,5,6,7），Y=（7,6,5,4），做褶积和相关，并画出图像。理解褶积与相关的区别和联系。
实验内容：
褶积
褶积是两个序列先倒序，再相乘相加，即：
1 2 3 4 5 6 7
4 5 6 7
4 5 6 7
4 5 6 7
4 5 6 7
4 5 6 7
4 5 6 7
4 5 6 7
4 5 6 7
4 5 6 7
4 5 6 7
根据此原理编程序：
N=7;
M=4;
k=0;
x=zeros(1,7);
y=zeros(1,4);
zj=zeros(1,10);
ZJ=zeros(1,10);
for n=1:N
x(n)=n;
end
for m=1:M
y(m)=8-m;
end
for i=0:M+N-2
for j=0:N-1
if (i-j)=0(i-j)=M-1
ZJ(i+1)=ZJ(i+1)+x(j+1)*y(i-j+1)
end
end
end
plot(ZJ)
我们得到了这样的图像：
横坐标所对应纵坐标的值正好和我们用褶积原理计算出来的值
利用matlab自带的褶积函数：
zj=conv(x,y)
plot(zj)
得到:
我们发现，对于褶积来说，利用自带的conv函数和编的程序所得到的图像是一样的，因为这二者的原理是一样的：对于第二个序列先反转，再相乘相加。
做DFT分析其频谱：
for k=1:10
s=0;
for i=1:10
s=s+ZJ(i)*exp(-j*2*pi*(k-1)*(i-1)/N);
end
q(k)=s;
end
plot(abs(q))
b=angle(a)
plot(b)
这样，我们就得到了上述序列做DFT的振幅谱和相位谱：
振幅谱：
相位谱：
相关
相关是两个序列直接相乘相加，不用倒序。即：
1 2 3 4 5 6 7
7 6 5 4
7 6 5 4
7 6 5 4
7 6 5 4
7 6 5 4
7 6 5 4
7 6 5 4
据此编写程序：
xg=zeros(1,7);
XG=zeros(1,7);
for i=0:6
for j=0:3
if i+j=6
XG(i+1)=XG(i+1)+x(i+j+1)*y(j+1)
end
end
end
plot(XG)
我们得到相关的图像：
利用matlab自带的xcorr函数（自相关函数）来实现两个序列的相关：
xg=xcorr(x,y)
plot(xg)
我们可以发现，对于相关来说，自己编写的程序得到的图像和自带函数xcorr函数得到的图像是不同的。首先最明显的就是图形的不同；其次就是自己编的相关，横坐标是7而xcorr函数的横坐标是14；再看图像，自己编的程序我们得到7个值而xcorr函数我们得到13个值 。这是因为xcorr函数是将短的序列先用0补齐长度，再移位相乘相加，而我们这里的相关是没有补零，没有移位，直接相乘相加，这就造成了相关图像的差异。
xcorr函数的工作原理：
1 2 3 4 5 6 7
7 6 5 4 0 0 0
7
文档评论（0）