查看: 79|回复: 0|关注: 0
打印 上一主题 下一主题

[未答复] 请问一组序列的自相关需不需要先对树立归一化再算

[复制链接]

新手

10 麦片

财富积分


050


85

主题

325

帖子

0

最佳答案
跳转到指定楼层
1#
发表于 2019-6-23 23:27:27 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式
请问一组序列的自相关需不需要先对树立归一化再算?谢谢
%测试测震数据在没有地震情况下的自相关
clf;
clc;
clear all;
pid=fopen('newtest00ud.txt','wt+');
fid=fopen('test00ud.txt','r');
data_dfa=textscan(fid,'%f\r','HeaderLines',1);
fclose(fid);
row=length(data_dfa{1});
for i=1:row
fprintf(pid,'%f\n',data_dfa{1,1}(i));
end
x=load('newtest00ud.txt');
N=length(x);
%xmax=max(x);
%xmin=min(x);
%xnew=(x-xmin)/(xmax-xmin);  %需要归一化吗?
Fs=100;%数据长度及采样频率
n=0:N-1;t=n/Fs;%时间序列
Lag=100;%延迟样点数
%[c,lags]=xcorr(xnew,Lag,'unbiased');%对原始信号进行无偏相关估计
[c,lags]=xcorr(x,Lag,'unbiased');%对原始信号进行无偏相关估计
subplot(2,1,1),plot(t,x);%绘制原始信号x
xlabel('时间/s');ylabel('x(t)');title('原始信号');grid on;
subplot(2,1,2);plot(lags/Fs,c);%绘制x信号自相关,lags/Fs为时间序列
xlabel('时间/s');ylabel('R(t)');title('原始信号的自相关');grid on;



序列
-28793.000000
-28833.000000
-28792.000000
-28812.000000
-28838.000000
-28824.000000
-28813.000000
-28840.000000
-28813.000000
-28796.000000
-28857.000000
-28800.000000
-28746.000000
-28810.000000
-28817.000000
-28809.000000
-28850.000000
-28815.000000
-28781.000000
-28798.000000
-28813.000000
-28830.000000
-28838.000000
-28826.000000
-28860.000000
-28815.000000
-28802.000000
-28852.000000
-28762.000000
-28901.000000
-28970.000000
-28842.000000
-28936.000000
-28967.000000
-28922.000000
-28911.000000
-28929.000000
-28970.000000
-28985.000000
-29000.000000
-28973.000000
-28951.000000
-28924.000000
-28910.000000
-28927.000000
-28934.000000
-28989.000000
-28979.000000
-28889.000000

.........


计算结果为啥范围不是-1到1呢?

您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /3 下一条

快速回复 返回顶部 返回列表