clc;clear all;close all;
n = 2000;
fs = 150; %采样频率
x = 1:1:n; % 数据点
data = sin((x/fs)*2*pi); % 频率为1Hz的正弦波
% 谱缓存 2倍采样率谱 对应谱最大周期为1/0.5Hz = 2s
spec_arr = zeros(300,n);
% 数据缓存
data_arr = zeros(301,1);
% DSP下,每次读入一个新的数据
for i = 1:1:n
data_arr(1:300) = data_arr(2:301);
data_arr(301) = data(i);
% 数据量足够做谱的话
if(i>300)
for j = 300:-1:1
spec_arr(j,i) = data_arr(301)*data_arr(j);
end
end
end
figure(1);
plot(data);
figure(2);
surf(spec_arr);
[数字信号处理]信号自相关谱算法实现
最新推荐文章于 2023-10-12 21:33:34 发布