✅作者简介:热爱科研的Matlab仿真开发者,修心和技术同步精进,matlab项目合作可私信。
🍎个人主页:Matlab科研工作室
🍊个人信条:格物致知。
⛄ 内容介绍
In the paper, a novel estimation method is proposed, concerning the wide-sense stationarity testof the signals. Background information is given about the concept of stationarity of the processes and signals. The problem of the signal stationarity estimation is a ddressed along with criticism of the availablestationarity tests. Further, a new wide-sense stationarity estimation method is described, involving estimation of the mean-, variance- and autocovariance- stationarity of a signal. Finally, a few representative signals are tested and the results clearly indicate the consistence of the proposed test method.
⛄ 部分代码
clear, clc, close all
%% generate test signal
% sampling frequency
fs = 44100;
% TS1A
t1 = 0:1/fs:1;
x1 = 1.0*sin(2*pi*1000*t1);
% TS2A
ta = 0:1/fs:0.2;
tb = 0:1/fs:0.8;
xa = 1.0*sin(2*pi*1000*ta);
xb = 2.0*sin(2*pi*1000*tb);
x2 = [xa xb];
% TS3A
tc = 0:1/fs:0.2;
td = 0:1/fs:0.8;
xc = 1.0*sin(2*pi*1000*tc);
xd = 1.0*sin(2*pi*2000*td);
x3 = [xc xd];
% TS4A
t4 = 0:1/fs:1;
x4 = chirp(t4, 1000, t4(end), 10000);
%% statiogram computation
% set the frame durations to check about stationarity
fr_dur = [5 10 20 30 40 50 100 200 300 400 500]*1e-3;
fr_len = round(fr_dur*fs);
% set the confidence level
gamma = 0.5;
% plot the signal statiogram
figure(1)
subplot(2, 2, 1)
statiogram(x1, fr_len, gamma)
xlabel('Frame duration, ms')
ylabel([])
xticklabels(num2str(fr_dur'*1000))
yticklabels({'Mean', 'Var', 'AutoCoVar'})
title('Statiogram for TS1A')
subplot(2, 2, 2)
statiogram(x2, fr_len, gamma)
xlabel('Frame duration, ms')
ylabel([])
xticklabels(num2str(fr_dur'*1000))
yticklabels({'Mean', 'Var', 'AutoCoVar'})
title('Statiogram for TS2A')
subplot(2, 2, 3)
statiogram(x3, fr_len, gamma)
xlabel('Frame duration, ms')
ylabel([])
xticklabels(num2str(fr_dur'*1000))
yticklabels({'Mean', 'Var', 'AutoCoVar'})
title('Statiogram for TS3A')
subplot(2, 2, 4)
statiogram(x4, fr_len, gamma)
xlabel('Frame duration, ms')
ylabel([])
xticklabels(num2str(fr_dur'*1000))
yticklabels({'Mean', 'Var', 'AutoCoVar'})
title('Statiogram for TS4A')
⛄ 运行结果
⛄ 参考文献
[1] H. Zhivomirov, I. Nedelchev. A Method for Signal Stationarity Estimation. Romanian Journal of Acoustics and Vibration, ISSN: 1584-7284, Vol. XVII, No. 2, pp. 149-155, 2020. (http://rjav.sra.ro/index.php/rjav/article/view/178/103).
[2] --- A research article dedicated to this estimation technique is about to come by the end of 2023. ---