大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。
Jetbrains全家桶1年46,售后保障稳定
平稳和非平稳都是针对随机信号说的。
平稳信号是指分布参数或者分布律随时间不发生变化的信号。
非平稳信号是指分布参数或者分布律随时间发生变化的信号。
常见的平稳信号:一段短的语音信号(10~30ms)
常见的非平稳信号:一段长的语音信号,一首歌等
给出一个信号(可以随意选取一个.mp3文件,命名为m.mp3),对其进行平稳性分析:
%信号的平稳性分析
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(局部)相关系数,可见平稳');
%结论:从相关函数可以看出(期望看不出),此音乐信号具有短时平稳性
发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/219969.html原文链接:https://javaforall.cn
【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛
【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...