使用持久性频谱查找干扰
可视化嵌入在宽带信号中的干扰窄带信号。
以1 kHz的形式产生一个chirp,持续500秒。在测量过程中,CHIRP的频率从180 Hz增加到220 Hz。
FS = 1000;t =(0:1/fs:500)';x = chirp(t,180,t(end),220) + 0.15*randn(size(t));
该信号还包含210 Hz正弦曲线。正弦曲线的振幅为0.05,仅在总信号持续时间的1/6中存在。
idx =落地(长度(x)/6);x(1:idx)= x(1:idx) + 0.05*cos(2*pi*t(1:idx)*210);
将信号保存为MATLAB®时间表。
s =时间表(秒(t),x);
打开信号分析仪并将时间表从工作区浏览器显示。点击时频按钮添加频谱图。在频谱图标签,下时间分辨率, 选择指定并输入1秒钟的时间分辨率。设置频率限制至100 Hz和290 Hz。两个信号组件都是可见的。
回到展示标签。点击时间按钮删除时间视图,然后单击光谱按钮以添加电源谱视图。频率范围继续从100 Hz到290 Hz。弱的正弦曲线被鸣叫掩盖。
单击箭头光谱按钮更改光谱
查看持续频谱
看法。在持续频谱标签,下时间分辨率, 选择指定并输入1秒钟的时间分辨率。指定相邻段之间的零重叠。设置功率限制到–50 dB和0 dB,密度极限到0.1和4.现在两个信号组件都清晰可见。
在展示标签,下分享, 点击生成脚本并选择持久性频谱脚本
。该脚本出现在MATLAB编辑器中。
%计算持续频谱由MATLAB(R)9.7和信号处理工具箱8.2生成的%。%生成:26-DEC-2018 16:07:45% 参数timelimits =秒([0 500]);%秒频率limits = [100 290];%Hz计时度= 1;%秒重叠= 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,...“持久性”,,,,...“频率限制”,频率限制,...“计时解决”,计时解决,...“超贴”,夸张);
也可以看看
应用
功能
相关示例
- 在相关信号之间找到延迟
- 通过不同的窗口泄漏来解决音调
- 使用复杂的包膜调节和解调
- 使用重新签名的光谱图查找和跟踪脊
- 从音乐信号中提取声音
- 重新采样并过滤不均匀采样的信号
- 使用您自己的功能拒绝饱和信号
- 计算振动信号的包络光谱
- 从鲸鱼的摘录区域提取区域