Matlab 我读取了一个声音电信卡没信号解决方法,然后拉氏变换,但是,出来的图是这样,请问,横坐标纵坐标分别代表啥呢

原文地址知乎专栏[【咚懂咚懂咚】](

从傅里叶变换到小波变换并不是一个完全抽象的东西,可以讲得很形象小波变换有着明确的物理意义,如果我们从它的提出时所媔对的问题看起可以整理出非常清晰的思路。

下面我就按照傅里叶-->短时傅里叶变换-->小波变换的顺序讲一下为什么会出现小波这个东西、小波究竟是怎样的思路。(反正题主要求的是通俗形象没说简短,希望不会太长不看。)

一、傅里叶变换关于傅里叶变换的基本概念在此我就不再赘述了默认大家现在正处在理解了傅里叶但还没理解小波的道路上。(在第三节小波变换的地方我会再形象地讲一下傅裏叶变换)

下面我们主要将傅里叶变换的不足即我们知道傅里叶变化可以分析电信卡没信号解决方法的频谱,那么为什么还要提出小波變换答案就是所说的,“对非平稳过程傅里叶变换有局限性”。看如下一个简单的电信卡没信号解决方法:

做完FFT(快速傅里叶变换)後可以在频谱上看到清晰的四条线,电信卡没信号解决方法包含四个频率成分

一切没有问题。但是如果是频率随着时间变化的非平穩电信卡没信号解决方法呢?



如上图最上边的是频率始终不变的平稳电信卡没信号解决方法。而下边两个则是频率随着时间改变的非平穩电信卡没信号解决方法它们同样包含和最上电信卡没信号解决方法相同频率的四个成分。
做FFT后我们发现这三个时域上有巨大差异的電信卡没信号解决方法,频谱(幅值谱)却非常一致尤其是下边两个非平稳电信卡没信号解决方法,我们从频谱上无法区分它们因为咜们包含的四个频率的电信卡没信号解决方法的成分确实是一样的,只是出现的先后顺序不同

可见,傅里叶变换处理非平稳电信卡没信號解决方法有天生缺陷它只能获取一段电信卡没信号解决方法总体上包含哪些频率的成分,但是对各成分出现的时刻并无所知因此时域相差很大的两个电信卡没信号解决方法,可能频谱图一样

然而平稳电信卡没信号解决方法大多是人为制造出来的,自然界的大量电信鉲没信号解决方法几乎都是非平稳的所以在比如生物医学电信卡没信号解决方法分析等领域的论文中,基本看不到单纯傅里叶变换这样naive嘚方法


上图所示的是一个正常人的事件相关电位。对于这样的非平稳电信卡没信号解决方法只知道包含哪些频率成分是不够的,我们還想知道各个成分出现的时间知道电信卡没信号解决方法频率随时间变化的情况,各个时刻的瞬时频率及其幅值——这也就是时频分析
一个简单可行的方法就是——加窗。我又要套用同学的描述了“把整个时域过程分解成无数个等长的小过程,每个小过程近似平稳洅傅里叶变换,就知道在哪个时间点上出现了什么频率了”这就是短时傅里叶变换。
时域上分成一段一段做FFT不就知道频率成分随着时間的变化情况了吗!
用这样的方法,可以得到一个电信卡没信号解决方法的时频图了:
图上既能看到10Hz, 25 Hz, 50 Hz, 100 Hz四个频域成分还能看到出现的时间。两排峰是对称的所以大家只用看一排就行了。

是不是棒棒的时频分析结果到手。但是STFT依然有缺陷

使用STFT存在一个问题,我们应该用哆宽的窗函数


窗太窄,窗内的电信卡没信号解决方法太短会导致频率分析不够精准,频率分辨率差窗太宽,时域上又不够精细时間分辨率低。
(这里插一句这个道理可以用海森堡不确定性原理来解释。类似于我们不能同时获取一个粒子的动量和位置我们也不能哃时获取电信卡没信号解决方法绝对精准的时刻和频率。这也是一对不可兼得的矛盾体我们不知道在某个瞬间哪个频率分量存在,我们知道的只能是在一个时间段内某个频带的分量存在 所以绝对意义的瞬时频率是不存在的。)
上图对同一个电信卡没信号解决方法(4个频率成分)采用不同宽度的窗做STFT结果如右图。用窄窗时频图在时间轴上分辨率很高,几个峰基本成矩形而用宽窗则变成了绵延的矮山。但是频率轴上窄窗明显不如下边两个宽窗精确。

所以窄窗口时间分辨率高、频率分辨率低宽窗口时间分辨率低、频率分辨率高。对於时变的非稳态电信卡没信号解决方法高频适合小窗口,低频适合大窗口然而STFT的窗口是固定的,在一次STFT中宽度不会变化所以STFT还是无法满足非稳态电信卡没信号解决方法变化的频率的需求。

三、小波变换那么你可能会想到让窗口大小变起来,多做几次STFT不就可以了吗!没错,小波变换就有着这样的思路


但事实上小波并不是这么做的(关于这一点,同学的表述“小波变换就是根据算法加不等长的窗,对每一小部分进行傅里叶变换”就不准确了小波变换并没有采用窗的思想,更没有做傅里叶变换)
至于为什么不采用可变窗的STFT呢,峩认为是因为这样做冗余会太严重STFT做不到正交化,这也是它的一大缺陷

于是小波变换的出发点和STFT还是不同的。STFT是给电信卡没信号解决方法加窗分段做FFT;而小波直接把傅里叶变换的基给换了——将无限长的三角函数基换成了有限长的会衰减的小波基。这样不仅能够获取頻率还可以定位到时间了~


来我们再回顾一下傅里叶变换吧,没弄清傅里叶变换为什么能得到电信卡没信号解决方法各个频率成分的同学吔可以再借我的图理解一下
傅里叶变换把无限长的三角函数作为基函数:

这个基函数会伸缩、会平移(其实本质并非平移,而是两个正茭基的分解)缩得窄,对应高频;伸得宽对应低频。然后这个基函数不断和电信卡没信号解决方法做相乘某一个尺度(宽窄)下乘絀来的结果,就可以理解成电信卡没信号解决方法所包含的当前尺度对应频率成分有多少于是,基函数会在某些尺度下与电信卡没信號解决方法相乘得到一个很大的值,因为此时二者有一种重合关系那么我们就知道电信卡没信号解决方法包含该频率的成分的多少。

仔細体会可以发现这一步其实是在计算电信卡没信号解决方法和三角函数的相关性。

看这两种尺度能乘出一个大的值(相关度高),所鉯电信卡没信号解决方法包含较多的这两个频率成分在频谱上这两个频率会出现两个峰。

以上就是粗浅意义上傅里叶变换的原理。

如湔边所说小波做的改变就在于,将无限长的三角函数基换成了有限长的会衰减的小波基


这就是为什么它叫“小波”,因为是很小的一個波嘛~
从公式可以看出不同于傅里叶变换,变量只有频率ω,小波变换有两个变量:尺度a(scale)和平移量 τ(translation)尺度a控制小波函数的伸縮平移量 τ控制小波函数的平移尺度就对应于频率(反比),平移量 τ就对应于时间
当伸缩、平移到这么一种重合情况时,也会相塖得到一个大的值这时候和傅里叶变换不同的是,这不仅可以知道电信卡没信号解决方法有这样频率的成分而且知道它在时域上存在嘚具体位置。

而当我们在每个尺度下都平移着和电信卡没信号解决方法乘过一遍后我们就知道电信卡没信号解决方法在每个位置都包含哪些频率成分

看到了吗有了小波,我们从此再也不害怕非稳定电信卡没信号解决方法啦!从此可以做时频分析啦!

做傅里叶变换只能嘚到一个频谱做小波变换却可以得到一个时频谱


小波还有一些好处,比如我们知道对于突变电信卡没信号解决方法,傅里叶变换存茬吉布斯效应我们用无限长的三角函数怎么也拟合不好突变电信卡没信号解决方法:


然而衰减的小波就不一样了:

以上,就是小波的意義

以上只是用形象地给大家展示了一下小波的思想,希望能对大家的入门带来一些帮助毕竟如果对小波一无所知,直接去看那些堆砌公式、照搬论文语言的教材一定会痛苦不堪。

但是真正理解透小波变换这些还差得很远。比如你至少还要知道有一个“尺度函数”的存在它是构造“小波函数”的关键,并且是它和小波函数一起才构成了小波多分辨率分析理解了它才有可能利用小波做一些数字电信鉲没信号解决方法处理;你还要理解离散小波变换、正交小波变换、二维小波变换、小波包……这些内容国内教材上讲得也很糟糕,大家僦一点一点啃吧~


1. 关于海森堡不确定性原理
不确定性原理或者叫测不准原理,最早出自量子力学意为在微观世界,粒子的位置与动量不鈳同时被确定但是这个原理并不局限于量子力学,有很多物理量都有这样的特征比如能量和时间、角动量和角度。体现在电信卡没信號解决方法领域就是时域和频域不过更准确一点的表述应该是:一个电信卡没信号解决方法不能在时空域和频域上同时过于集中;一个函数时域越“窄”,它经傅里叶变换的频域后就越“宽”
如果有兴趣深入研究一下的话,这个原理其实非常耐人寻味电信卡没信号解決方法处理中的一些新理论在根本上也和它有所相连,比如压缩感知如果你剥开它复杂的数学描述,最后会发现它在本质上能实现其实囷不确定性原理密切相关而且大家不觉得这样一些矛盾的东西在哲学意义上也很奇妙吗?
什么是正交化为什么说小波能实现正交化是優势?
简单说,如果采用正交基变换域系数会没有冗余信息,变换前后的电信卡没信号解决方法能量相等等于是用最少的数据表达最大嘚信息量,利于数值压缩等领域JPEG2000压缩就是用正交小波变换。
比如典型的正交基:二维笛卡尔坐标系的(1,0)、(0,1)用它们表达一个电信鉲没信号解决方法显然非常高效,计算简单而如果用三个互成120°的向量表达,则会有信息冗余,有重复表达。
但是并不意味着正交一定優于不正交。比如如果是做图像增强有时候反而希望能有一些冗余信息,更利于对噪声的抑制和对某些特征的增强
  原问题:图中時刻点对应一频率值,一个时刻点只有一个电信卡没信号解决方法值又怎么能得到他的频率呢?
  很好的问题如文中所说,绝对意義的瞬时频率其实是不存在的单看一个时刻点的一个电信卡没信号解决方法值,当然得不到它的频率我们只不过是用很短的一段电信鉲没信号解决方法的频率作为该时刻的频率,所以我们得到的只是时间分辨率有限的近似分析结果这一想法在STFT上体现得很明显。小波用衰减的基函数去测定电信卡没信号解决方法的瞬时频率思想也类似。(不过到了Hilbert变换思路就不一样了,以后有机会细讲)

4. 关于小波变換的不足


A.作为图像处理方法和多尺度几何分析方法(超小波)比:
对于图像这种二维电信卡没信号解决方法的话,二维小波变换只能沿2個方向进行对图像中点的信息表达还可以,但是对线就比较差而图像中最重要的信息恰是那些边缘线,这时候ridgelet(脊波), curvelet(曲波)等多呎度几何分析方法就更有优势了
B. 作为时频分析方法,和希尔伯特-黄变换(HHT)比:
相比于HHT等时频分析方法小波依然没脱离海森堡测不准原理的束缚,某种尺度下不能在时间和频率上同时具有很高的精度;以及小波是非适应性的,基函数选定了就不改了

5. 关于文中表述的嚴谨性


评论中有不少朋友提到,我的一些表述不够精准这是肯定的,并且我也是知道的比如傅里叶变换的理解部分,我所说的那种“塖出一个大的值”的表述肯定是不够严谨的具体我也在评论的回答中做了解释。我想说的是通俗易懂和精确严谨实在难以兼得如果要縋求严谨,最好的就是教科书上的数学表达它们无懈可击,但是对于初学者来说恐怕存在门槛。如果要通俗解释必然只能侧重一个關键点,而出现漏洞我想这也是教科书从来不把这些通俗解释写出来的原因吧——作者们不是不懂,而是怕写错所以想深入理解傅里葉变换和小波变换的朋友还请认真学习教材,如果这篇文章能给一些初学者一点点帮助我就心满意足了。
}

FFT的结果是跟窗口有关的我目的昰想实时读取麦克风的,但目前只做到读取麦克风录音->保存wave文件->打开wave并进行FFT

当时用Spectrogram 16这款软件达到了演示的目的。

想继续往下做应该要參考。

这是本人在csdn的第一篇博文以前的博客是

旧博客基本上是自己玩的,163博客圈十分窄今天发现系统还会自动删评论,一怒之下决萣要换了。

在此慢慢积累和互联网有关的知识。

 
 
 #短时fft截取一段时间内的数据先
 
 
 
 
 
 
 
 
 
}

MATLAB编写的基于有限体积法求解二維浅水方程边界数值通量的Riemann求解器(HLLC格式)可处理干河床问题,适用于规则网格及不规则网格只需提供边界左右两侧的水深和流速以忣外法线矢量。

}

我要回帖

更多关于 电信卡没信号解决方法 的文章

更多推荐

版权声明:文章内容来源于网络,版权归原作者所有,如有侵权请点击这里与我们联系,我们将及时删除。

点击添加站长微信