自相关函数怎么画

wlj16

有人会画这个图吗？

eight

原帖由 wlj16 于 2007-5-15 19:59 发表

                               
登录/注册后可看大图

  有人会画这个图吗？

大哥，你发贴前能否先搜索一下论坛，本版就有这方面的精华贴

zhaopeng161

global pathname filename sf nfft outfile;
fid=fopen([pathname filename],'r');
x=fscanf(fid,'%f',[1,inf]); %输入数据存成行向量
status=fclose(fid);
m=str2double(inputdlg('请输入相关数据长度(点数)'));
nt=length(x(1,:));
h=findobj('Tag','cypl');
sf=str2double(get(h,'string'));
%建立离散时间向量
t=0:1/sf:(nt-1)/sf;
%计算自相关函数向量
a=xcorr(x(1,:),m);
%取正频段的相关函数
y=a(m+1:2*m+1);
%绘制输入数据的时程曲线图
plot(findobj('tag','axes1'),t,x);
xlabel('时间(s)');
ylabel('幅值');
title('输入数据的时程曲线图');
%grid on;
%绘制相关函数曲线图
t=1:(m+1);
plot(findobj('tag','axes2'),t,y);
title('自相关函数');
outfile=y;

kmkeke

x=0:.001:1;
y=2*sin(2*pi*30*x);
[a,b]=xcorr(y,'unbiased');
plot(b,a)
这样做自相关分析，怎么自相关分析绘制的y轴会随正弦信号的幅值变化而变化呢？

当y的幅值取1时，绘制的图形显示相关系数为0.5，

当y的幅值为2时，怎么相关系数变成2了呢？请问哪里有问题呢