![](https://img-blog.csdnimg.cn/20201014180756928.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
数字信号处理
文章平均质量分 69
dsp的一些理解
季马宝宝
侠之小者,为情为私
展开
-
FIR滤波器设计(Matlab代码)
之前已经提过了,IIR级数低,元器件少,成本低,但起不到线性相位的作用。而FIR滤波器由于没有反馈,比较不稳定,消耗资源比较多,但是却可以产生线性相位。线性相位有什么用,然后我简单说一下FIR滤波器原理(图说搬来的),原理很简单,我们设计一个低通滤波器,就是用一个窗函数去截频谱,频域上就是信号和一个抽样函数卷积。抽样函数是无限长的,但是离y轴越远衰减越厉害,我们计算机只能先对它进行截断。截断完傅里叶回去,窗函数就不是理想窗函数了。但是我们不能直接使用截断的抽样函数进行卷积,我们对窗函数还是有要求原创 2021-12-21 22:16:32 · 6074 阅读 · 0 评论 -
相移、群延迟、线性相位
相移很简单,就是经过系统后,信号中特定的频率分量,相位会落后多少。线性相位很简单,就是θ=ω∗k\theta=\omega *kθ=ω∗k。信号中相位滞后(只考虑因果信号,不可能先于输入)和频率成正比。群延迟delay=dθdωdelay=\frac{d\theta}{d\omega}delay=dωdθ,显然线性相位群延迟就等于k,是一个常数,这个常数代表,系统中任何频率分量的延迟时间都是常数。你去查线性相位很多都会告诉到你这里,但是这是为什么呢?其实很简单,因为这段频率相角转动θ\thetaθ原创 2021-12-21 22:15:02 · 2795 阅读 · 0 评论 -
FFT频谱分析(matlab代码)
实验要求在Matlab软件平台,利用函数wavrecord函数录制一段语音,采样频率fs=10000Hz,长度为5s。用sound函数播放声音。用wavread对语音信号进行采样,记住采样频率和采样点数。 或者在Python环境下,利用Pyaudio语音处理包实现录音。语音信号的频谱分析首先画出语音信号的时域波形;然后对语音号进行快速傅里叶变换,得到语音信号的频谱特性 。加入单频噪声原始的语音信号加上一个高频余弦噪声,频率为4kHz。画出加噪后的语音信号时域和频谱图,与原始信号对比,原创 2021-12-21 20:08:04 · 12599 阅读 · 3 评论 -
IIR滤波器设计代码(巴特沃斯+脉冲响应不变法/双线性变换法) Matlab代码
clc;clear;%% 参数确定Wp=0.1*pi; %通带截至频率Ws=0.4*pi; %阻带截止频率ap=1; %通带最大衰减as=25; %阻带最大衰减fs=1; %采样频率%% 脉冲响应不变法+3dB通带截止频率wp=Wp*fs;ws=Ws*fs;N=buttord(wp,ws,ap,as,'s');wc=wp/(10^(0.1*ap)-1)^(1/2/N);[b,a]=butter(N,wc,'s');[c,d]=impinvar(b,a,fs);w=linspa原创 2021-12-21 20:03:57 · 2161 阅读 · 2 评论 -
频域抽样与频谱泄露
上次的文章谈论了DFT是怎么来的,为什么我们需要DFT。但是其中还有一部分细说——频率抽样。我们的频谱在计算机内也是以离散的形式存储的,那么我们是如何对频谱进行离散化的呢?时域上因为我们采集到的模拟信号,可以进行量化,但是我们并不会得到频谱图的连续信号啊,也就是说整个采样不是直接在频谱上量化的,而是在时域上完成的。频域的采样(离散化)其实就是时域的周期化,我们在求DFT的时候不是进行了周期延拓嘛,这就相当于完成了频域的采样,而时域周期越长(通过对信号补0实现),那频谱分辨率(频谱分辨率就是频谱抽样的间隔原创 2021-12-02 23:30:12 · 486 阅读 · 0 评论 -
傅立叶级数?变换?FFT?
引言本人专业是自动化,离不开信号与系统这门课,提到信号与系统,自然是离不开大名鼎鼎的傅立叶,傅立叶级数、傅立叶变换都是常用的计算工具。同时,我也搞搞算法竞赛,众所周知,FFT(快速傅立叶变换)是一种计算多项式乘法的算法。虽然我都接触过这两个东西,但是快速傅立叶变换和傅立叶变换,有什么关系吗?当然有,在此之前,你得先知道有个东西叫离散傅立叶变换(DFT),FFT的最早提出其实是为了快速计算出DFT。也就是说,FFT是一种算法,它的命名是因为它最早被用于求解离散傅立叶变换,而同时它也可以解决一些其他问题原创 2021-09-13 22:16:50 · 955 阅读 · 0 评论 -
傅立叶变换再理解——离散傅里叶变换DFT到底是用来干什么的?它有什么局限性?
FT、FS,DTFT、DFS第一行,自然界信号其实都是模拟的,都是连续的非周期信号,它的频谱也应该是像右侧一样的连续非周期频谱。但是我们是数字信号处理,什么是数字信号,就是采样量化后的信号,也就是说我们实际采集到的信号是像第三行一样的。而第三行的频谱是连续周期的,很显然我们永远都可以找到无线高频的信号恰好在那几个点被采样到,那假如说我们知道目标信号频率不会最多不会超过某个值,那频谱高于这个频率的值全都可以当作0,奈奎斯特采样定理告诉我们只要采样频率大于两倍信号最高频率,就可以保证频谱只有-fs(fs是.原创 2021-12-02 22:31:22 · 3891 阅读 · 0 评论 -
数字信号处理理解
傅里叶变换啥的写在单独的博客里了。泰勒级数展开存在收敛域的限制,收敛半径是到最近极点的距离。洛朗级数可以延拓到整个定义域。洛朗级数的计算方法有两种。一种是计算计算围线积分。极点把定义域分成了一圈又一圈,分别给每个圈算积分,其中包含了很多复变知识,比较难算。第二种是展开,0-1可以映射到1-无穷,取倒数用泰勒公式展开就行了。洛朗级数就是Z变换,形式是一样的。Z变换是对于序列而言的,序列的值对应的就是幂级数的系数。虽然形式一样,但代表的含义不一样,如果你能把含义联系起来记得告诉我。留数是复变原创 2021-10-12 19:27:49 · 1484 阅读 · 0 评论