使用持久谱发现干扰
可视化嵌入宽带信号中的干扰窄带信号。
产生一个啁啾采样在1千赫500秒。在测量过程中,啁啾的频率从180 Hz增加到220 Hz。
Fs = 1000;T = (0:1/fs:500)';x =唧唧喳喳(220 t, 180 t(结束),)+ 0.15 * randn(大小(t));
该信号还包含210赫兹的正弦波。正弦信号的幅值为0.05,仅在总信号持续时间的1/6出现。
Idx =地板(长度(x)/6);(1: idx) = x (1: idx) + 0.05 * cos(2 *π* t (1: idx) * 210);
将信号保存为MATLAB®时间表。
S =时间表(秒(t),x);
开放信号分析仪并将时间表从工作空间浏览器到一个显示器。单击时频按钮,添加光谱图视图。在光谱图选项卡,在时间分辨率中,选择指定并输入1秒的时间分辨率。设置频率的限制至100hz和290hz。两个信号组件都是可见的。
回到显示选项卡。单击时间按钮,删除时间视图,并单击光谱按钮,添加功率谱视图。频率范围继续为100hz ~ 290hz。微弱的正弦波被啁啾掩盖了。
的下面的箭头光谱按钮,以更改光谱
查看a持久性频谱
视图。在持久性频谱选项卡,在时间分辨率中,选择指定并输入1秒的时间分辨率。指定相邻段之间零重叠。设置权限到-50 dB和0 dB密度的限制到0.1和4。现在两个信号组件都清晰可见了。
在显示选项卡,在分享,点击生成脚本并选择持久性谱脚本
.该脚本出现在MATLAB编辑器中。
计算持久谱由MATLAB(R) 9.7和信号处理工具箱8.2生成。生成时间:2018年12月26日16:07:45%的参数timeLimits = seconds([0 500]);%秒frequencyLimits = [100 290];%赫兹timerresolution = 1;%秒overlapPercent = 0;% %索引到感兴趣的信号时间区域S_x_ROI = S(:,“x”);S_x_ROI = S_x_ROI(timerange(timeLimits(1),timeLimits(2),“关闭”), 1);计算光谱估计运行下面不带输出参数的函数调用来绘制结果[P,F,PWR] = pspectrum(S_x_ROI,...“坚持不懈”,...“FrequencyLimits”frequencyLimits,...“TimeResolution”timeResolution,...“OverlapPercent”, overlapPercent);
另请参阅
应用程序
功能
相关的例子
- 找出相关信号之间的延迟
- 通过改变窗户泄漏来解决音调问题
- 利用复包络的调制与解调
- 使用重新分配的光谱图寻找和跟踪脊
- 从音乐信号中提取声音
- 非均匀采样信号的重采样与滤波
- 使用自己的函数衰减饱和信号
- 计算振动信号的包络谱
- 从鲸鱼的歌声中提取感兴趣的区域