日志
发现NI的COT算法存在不足
热度 1
一、原始信号包含一个基频从5Hz到30Hz线性变化的信号，以及一个2倍频信号及一个65倍频信号。重采样为128倍基频，所以如果低通滤波处理不好就会存在频率混叠。
clc;
clear;
Ls=1024*64*4; %样本数
fs=4096*8; %采样频率
t0 =0:1/fs:(Ls-1)/fs;
t1=t0(length(t0));
f0=5; %起始频率
f1=30; %终止频率
x_1x = 1*chirp(t0,f0,t1,f1,[],280); %一倍频信号
x_65x = 0.65*chirp(t0,65*f0,t1,65*f1,[],115); %65倍频信号
x_2x = 0.75*chirp(t0,2*f0,t1,2*f1,[],20); %2倍频信号
x_noise=0.15*rand(1,length(x_1x)); %噪声
x=x_1x+x_2x+x_65x;%+x_noise; %合成信号
clear;
Ls=1024*64*4; %样本数
fs=4096*8; %采样频率
t0 =0:1/fs:(Ls-1)/fs;
t1=t0(length(t0));
f0=5; %起始频率
f1=30; %终止频率
x_1x = 1*chirp(t0,f0,t1,f1,[],280); %一倍频信号
x_65x = 0.65*chirp(t0,65*f0,t1,65*f1,[],115); %65倍频信号
x_2x = 0.75*chirp(t0,2*f0,t1,2*f1,[],20); %2倍频信号
x_noise=0.15*rand(1,length(x_1x)); %噪声
x=x_1x+x_2x+x_65x;%+x_noise; %合成信号
1、原始频谱
2、NI的COT重采样信号阶比谱
3、我的COT重采样阶比谱
二、如果将上述的65倍频信号改为100倍频，结果NI的COT重采样信号阶比谱中混叠频率就没有了
全部作者的其他最新日志
- • Circular Plot
- • pyFFTW
- • 轴心轨迹信号处理
- • CNN在轴承故障诊断上的应用
- • 阶比包络分析
- • 关于谱峭度选择包络解调分析的中心频率和带宽
发表评论 评论 (7 个评论)
- 回复 ckakaqiang
- 为什么在阶次谱中没有65阶次的信号呢？
- 回复 ckakaqiang
- 那为什么要重采样为128倍基频呢？要是显示65倍频信号的话，提高重采样就行了吧？但你是为了说明阶比谱中混叠频率现象，那么是不是重采样越高越好呢？就是等角域采样时插的值越密越好呢？
- 回复 impulse
- 阶比太高也没意义，主要看你振动分析需求，太高的话对信号分析时间也不利，比如很多国外商业软件对于转速超出一定范围以后（比如18000rpm），不允许你选择256插值，需要选择128、64等等
- 回复 ckakaqiang
- 那么如果我的重采样很高的话，是不是能从一定程度上减小混叠频率现象呢？
- 回复 ckakaqiang
- 嗯嗯，也对，谢谢指导！
