非平稳信号的频谱分析方法—(短时傅立叶变换)

非平稳信号的频谱分析方法—(短时傅立叶变换)非平稳信号又称时变信号。对这一类信号,其一阶、二阶统计量和功率谱的估计显然不能简单的使用平稳信号的估计方法,必须考虑它们的时变因素。

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全家桶1年46,售后保障稳定

   非平稳信号又称时变信号。对这一类信号,其一阶、二阶统计量和功率谱的估计显然不能简单的使用平稳信号的估计方法,必须考虑它们的时变因素。

基本原理

   对非平稳信号,人们希望能有一种分析方法把时域分析和频域分析结合起来,即找到一个二维函数,它既能反映信号的频率内容,也能反映出该频率内容随时间变化的规律。研究这一问题的信号处理理论称为信号的联合时频分布。其中最重要的是以Cohen类为代表的双线性时频分布,此分布可表示为

非平稳信号的频谱分析方法---(短时傅立叶变换)   

   式中非平稳信号的频谱分析方法---(短时傅立叶变换)是一个二维的窗函数,给定不同的窗函数可以得到不同的时频分布。在上式中x(t)出现了两次,且是相乘的形式,这一特点称为双线性。

非平稳信号的频谱分析方法---(短时傅立叶变换)       

   式中w是一个一维的窗函数,则(1)式可以简化成如下的谱图

非平稳信号的频谱分析方法---(短时傅立叶变换)          

式中

非平稳信号的频谱分析方法---(短时傅立叶变换)               

称为信号x(t)的短时傅里叶变换,它反映了信号的频谱随时间和频率的分布。

Matlab仿真分析

   生成的信号两个LFM信号,一个频率随时间由小变大,另一个由大变小,中心频率都为150KHz。

Matlab代码为:

f0 = 150e3;            %定义LFM信号的中心频率
d_f = 100e3;            %定义LFM信号的频率变化大小
fs = 5*f0;             %定义采样率
points = 2^15;         %采样点数
SNR = 40;              %信噪比
 
nTs = (0:points-1)/fs;
Delta_f = (0:points-1)*fs/points;
%生成LFM信号
signal = sin(2*pi*(f0 + (d_f*fs)/(2*points)*(nTs - points/(2*fs))).*(nTs - points/(2*fs))) + ...  %频率由小变大
         sin(2*pi*(f0 + (-d_f*fs)/(2*points)*(nTs - points/(2*fs))).*(nTs - points/(2*fs)));      %频率由大变小                                                                    
%加入一定功率的白噪声
signal_noise = signal + 10^(-SNR/20)*randn(1,points);

Jetbrains全家桶1年46,售后保障稳定

N = 2^8;                                          %截取信号片段点数
window = hanning(N);                              %窗函数类型
Noverlap = N - N/2;                               %相邻两段之间重叠的点数
[S,F,T] = specgram(signal_noise,N,fs,window,Noverlap);
abs_S = abs(S);
max_abs_S = max(max(abs_S));
abs_S = abs_S/max_abs_S;
log10_abs_S = 20*log10(abs_S);
[tt,ff] = meshgrid(T,F/1000);
figure;
surf(tt,ff,log10_abs_S);
xlabel('时间/s');
ylabel('频率KHz');
zlabel('归一化功率谱P(w,t)/dB');
string = ['短时傅里叶变换,采样点数=',num2str(points),',LFM信号中心频率f0=',num2str(f0/1000),...
          'KHz,信号频率变化量d_f=',num2str(d_f/1000),'KHz,采样率fs=',num2str(fs/1000),'KHz,信噪比SNR=',...
          num2str(SNR),'dB'];
title(string);


运行一下得到如下图:

非平稳信号的频谱分析方法---(短时傅立叶变换)

调整一下得到如下的俯视图

非平稳信号的频谱分析方法---(短时傅立叶变换)

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/234168.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)
blank

相关推荐

  • 数据库恢复的三种方式

    数据库恢复的三种方式数据库恢复的三种方式

  • mysql uniqueidentifier_Uniqueidentifier数据类型[通俗易懂]

    mysql uniqueidentifier_Uniqueidentifier数据类型[通俗易懂]一.Uniqueidentifier数据类型可存储16字节的二进制值Uniqueidentifier用来存储一个全局唯一标识符,即GUID。GUID是唯一的二进制数:世界上的任何两台计算机都不会生成重复的GUID值uniqueidentifier值通常不定义为常量。您可以按下列方式指定uniqueidentifier常量:字符串格式:’6F9619FF-8B86-D011-B42D-00C0…

  • 黑马程序员——JAVA学习笔记四(继承、接口、内部类)

    黑马程序员——JAVA学习笔记四(继承、接口、内部类)

  • 图像处理中经常用的一幅图片Lena.jpg

    图像处理中经常用的一幅图片Lena.jpg在数字图像处理中,Lena(Lenna)是一张被广泛使用的标准图片,特别在图像压缩的算法研究中。(为什么用这幅图,是因为这图的各个频段的能量都很丰富:即有低频(光滑的皮肤),也有高频(帽子上的羽毛),很适合来验证各种算法)然而,这张图片背后的故事是颇有意思的,很多人都抱有学究都是呆子的看法,然而Lena对此就是一个有力的驳斥。lena(lenna)是一张于1972年11月…

  • Linux lamp_lamp搭建和配置

    Linux lamp_lamp搭建和配置LAMP架构LAMP架构是目前成熟的企业网站应用模式之一,指的是协同工作的一整套系统和相关软件,能够提供动态Web站点服务及其应用开发环境。LAMP是一个缩写词,具体包括Linux操作系统、Apache网站服务器、MySQL数据库服务器、PHP(或Perl、Python)网页编码。在构建LAMP平台时,各组件的安装顺序依次为Linux、Apache、MySQL、PHP。其中…

    2022年10月17日
  • matlab 画图的颜色

    matlab 画图的颜色plot函数代表不同颜色的标示符一共有八种:y:黄色;k:黑色;w:白色;b:蓝色;g:绿色;r:红色;c:亮青色;m:锰紫色;其他的就要自己设定了,‘color’,[000];三原色你值得拥有;…

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号