平稳和非平稳信号_非平稳信号处理

%信号的平稳性分析
clc,clear;
fs=44100;
s=audioread('m.mp3');%通过一个样本来分析（平稳性是针对随机信号的）
s(:,2)=[];

%% 看看全局平稳性
%1均值
n=400;
N=length(s);
mean(s((1:n)+round(N/4)))
mean(s((N-n+1:N)-round(N/4)))
%2相关函数
tao=round(N/4);%自定义时间差tao
for n0=1:30
    tmp=corrcoef(s((1:n)+n0+round(N/4)),s((1:n)+n0+tao+round(N/4)));c(n0)=tmp(1,2);
end
figure,plot(c);axis([1 30 0 1]);grid on;title('tao1(全局)相关系数，可见非平稳');

%% 看看局部平稳性
t=0.01; %选取局部时间10ms
NN=fs*t;
x=s(round(N/3)+(1:NN));%取一段长度为NN的信号
%1均值
n=400;
mean(x(1:n))
mean(x(NN-n+1:NN))
%2相关函数
tao=3;%自定义时间差tao
for n0=1:30
    tmp=corrcoef(x((1:n)+n0),x((1:n)+n0+tao));c(n0)=tmp(1,2);
end
figure,plot(c);axis([1 30 0 1]);grid on;title('tao2(局部)相关系数，可见平稳');

%结论：从相关函数可以看出（期望看不出），此音乐信号具有短时平稳性