数字信号处理&高阶(胡广书)
文章平均质量分 79
沈子恒
计算机视觉 自然语言处理 推荐系统 架构设计
展开
-
我所理解的快速傅里叶变换(FFT)
1.历史放在最前头首先FFT是离散傅立叶变换(DFT)的快速算法,那么说到FFT,我们自然要先讲清楚傅立叶变换。先来看看傅立叶变换是从哪里来的?傅立叶是一位法国数学家和物理学家的名字,英语原名是Jean Baptiste Joseph Fourier(1768-1830), Fourier对热传递很感兴趣,于1807年在法国科学学会上发表了一篇论文,运用正弦曲线来描述温原创 2016-10-22 18:01:58 · 175278 阅读 · 13 评论 -
深入探讨傅立叶变换、拉普拉斯变换、Z变换的联系与应用
1.前言一个信号,通常用一个时间的函数来表示,这样简单直观,因为它的函数图像可以看做信号的波形,比如声波和水波等等。很多时候,对信号的处理是很特殊的,比如说线性电路会将输入的正弦信号处理后,输出仍然是正弦信号,只是幅度和相位有一个变化(实际上从数学上看是因为指数函数是线性微分方程的特征函数,就好像矩阵的特征向量一样,而这个复幅度对应特征值)。因此,如果我们将信号全部分解成正弦信号的线性原创 2016-11-27 11:57:57 · 24474 阅读 · 5 评论 -
数字滤波器回顾与解析
1.综述数字滤波器是一个离散的系统。它可以对输入的离散信号进行一系列运算处理,从输入的信号中获得所需要的信息。数字滤波器的系统函数通常表示为数字滤波器分为有限冲激响应数字滤波器,即FIR数字滤波器和无限冲激响应,即IIR数字滤波器。从公式的角度来看,FIR数字滤波器的a(i)始终为零;IIR数字滤波器a(i)至少有一个非零。2.IIR滤波器IIR原创 2016-11-27 18:38:57 · 7102 阅读 · 0 评论 -
FIR数字滤波器设计_窗函数法
1.前言在数字信号处理中 ,由于信号中经常混有各种复杂成分,所以很多信号分析都是基于滤波器而进行的, 因此数字滤波器占有极其重要的地位 。数字滤波器是具有一定传输选择特性的数字信号处理装置,其输入与输出均为数字信号,实质上是一个由有限精度算法实现的线性时不变离散系统。它的基本工作原理是利用离散系统特性对系统输入信号进行加工和变换,改变输入序列的频谱或信号波形,让有用频率的信号分量通过,原创 2016-11-27 20:21:48 · 70671 阅读 · 0 评论 -
FIR数字滤波器设计——频率抽样法
1.频率抽样法设计线性相位FIR滤波器的思想频率抽样法是从频域出发,在频域直接设计,把给定的理想频率响应加以等间隔抽样,并以此作为实际FIR滤波器的频率响应。设所需滤波器的频率响应为。现要求设计一个M阶的FIR滤波器h[k],使得在M+1个抽样点上,FIR滤波器的频率响应与所需的频率响应相等,即由设计的要求给定,h[k]通过设计来确定。如果M+1个方程是线性无关原创 2016-11-27 21:07:26 · 35579 阅读 · 2 评论 -
FIR数字滤波器设计窗函数法MATLAB仿真(Hamming)
1 matlab源码f1=100;f2=200;%待滤波正弦信号频率fs=2000;%采样频率m=(0.3*f1)/(fs/2);%定义过度带宽M=round(8/m);%定义窗函数的长度N=M-1;%定义滤波器的阶数b=fir1(N,0.5*f2/(fs/2));%使用fir1函数设计滤波器%输入的参数分别是滤波器的阶数和截止频率figure(1)[h,f]=freqz(b,原创 2016-11-27 21:50:52 · 36752 阅读 · 2 评论 -
FIR数字滤波器设计频率抽样法MATLAB仿真
1.MATLAB源码M=63; Wp=0.5*pi;%所需频率采样点个数及通带截止频率m=0:(M+1)/2; Wm=2*pi*m./(M+1);%通频带上的采样点及阻带截止频率mtr=floor(Wp*(M+1)/(2*pi))+2;%向负方向入floor(3.5)=3;floor(-3.2)=-4Ad=[Wm<=Wp];Ad(mtr)=0.38;Hd=Ad.*exp(-1原创 2016-11-27 22:06:20 · 13993 阅读 · 2 评论 -
多采样率信号处理 ——信号的抽取与插值
1.引言一般多速率变换传输系统结构如上图。这种结构在传输时进行抽取,可以减少传输的数据量,接收端进行恢复时进行插值,将信号恢复成原始抽样率的信号。但是如果是在数字信号处理中只需要改变数字采样率,一般是先进行插值在进行抽取,以免造成频率混叠。进行抽取时,原始信号的频谱被周期拓展,如果原始信号最高频率大于抽取之后采样率的一半就会发生混叠(如果是复信号,则信号最高频率大于抽取之后原创 2016-11-28 10:19:52 · 45697 阅读 · 7 评论 -
信号分解:双正交、完备性、对偶向量
1.信号分解及完备性设是X由一组向量所张成,即:如果线性独立,我们则称它们为空间中的一组基”。那么信号x可以离散表示如下:若是一组两两互相正交的向量,展式称为x的正交展开。分解系数是在各个基向量上的投影。设向量和向量满足如下双正交关系:那么,我们对原始信号就行投影变换(內积):看看,我们把最关心的分解系数给弄出来了原创 2017-03-07 00:19:23 · 7551 阅读 · 1 评论 -
能量谱与功率谱
1.能量谱与功率谱能量谱与功率谱分别是针对能量有限的信号和功率有限的信号。在进行信号的谱分析的时候一定更要看准这是一个能量信号还是一个功率信号,应用不同的谱进行分析会使问题的解决思路更加明确。1.1 能量谱对于周期功率信号我们使用傅里叶级数的方法进行分析。指数形式的傅里叶级数系数的求解用下面的公式。对于非周期的能量信号,我们主要使用傅里叶分析方法进行分析。原创 2016-12-08 23:22:28 · 21085 阅读 · 3 评论 -
IIR+双向滤波实现系统零相位相移_MATLAB仿真
1.双向滤波实现零相移的思想Matlab软件有一个m文件filtfilt.m,可以实现零相位数字滤波。它先将输入序列按顺序滤波(forward filter),然后将所得结果逆转后反向通过滤波器(reverse filter),再将所得结果逆转后输出(reverse output),即得精确零相位失真的输出序列,为方便起见将这种滤波方法取名为FRR滤波。FRR滤波的原理可作如下推证:原创 2016-11-29 22:06:53 · 13199 阅读 · 8 评论 -
我所理解的离散傅里叶变换_DFT
1.闲话放在前面扯什么是频域?从我们出生,我们看到的世界都以时间贯穿,股票的走势、人的身高、汽车的轨迹都会随着时间发生改变。这种以时间作为参照来观察动态世界的方法我们称其为时域分析。而我们也想当然的认为,世间万物都在随着时间不停的改变,并且永远不会静止下来。但如果我告诉你,用另一种方法来观察世界的话,你会发现世界是永恒不变的,你会不会觉得我疯了?我没有疯,这个静止的世界就叫做频域。原创 2016-10-19 20:59:05 · 43601 阅读 · 2 评论 -
离散余弦变换_原理及应用
1.预备知识1.1可分离变换二维傅立叶变换可用通用的关系式来表示:式中:x, u=0, 1, 2, …, M-1;y, v=0, 1, 2, …, N-1;g(x,y,u,v)和h(x,y,u,v)分别称为正向变换核和反向变换核。 如果满足 :则称正、反变换核是可分离的。进一步,如果g1和g2,h1和h2在函数形式上一样原创 2016-10-29 17:54:41 · 75407 阅读 · 5 评论 -
Matlab生成各种伪随机信号_idinput指令
1.为何不是冲激信号?只要是在信号系统这条道上混过的童鞋都知道,冲激信号是测试系统的最理想信号。因为时域的卷积相当于频域的乘积,而冲击信号的频域拥有最广泛的频谱,可以“暴露”系统最完整的信息,那么问题来了?这么好的东西,你咋就不用呢?傻~还真不是这样,能让我们计算机产生一个理想的冲激信号那是多难的事啊!!!想象狄拉克,就被难为人了......那咱办呢?这时候一大群数学家来了....原创 2016-10-29 23:32:30 · 48189 阅读 · 10 评论 -
离散时间系统的相位响应
1.前言离散事件系统经过时频变换方法都会落到频域,表征频域的特性两个指标即可实现:幅频特性,相频特性。这样我们不仅要思考一个问题,到底是系统的幅频失真造成的后果严重,还是系统的相位失真造成的后果严重?下面我们可以用超声波束合成的一个例子来解释一下:实验中,我们发射一束超声(假设为3MHz)扫描人体某个部位,那么反射回来的超声波必定发生各种复杂的相位移动,并携带了人体原创 2016-11-12 00:21:23 · 10266 阅读 · 0 评论 -
Windows下FFTW的安装_史上最全
努力写好每一篇博客,做CSDN平台上的一股清流1.FFTW简介FFTW是用来计算一维或者多维的离散傅里叶变换,输入可以为实数序列也可以为复数序列的C语言的子函数库,FFTW是免费软件,是作为fft函数库的各种应用的上佳选择。这是MIT两位老爷爷开发的,据说是史上最快FFT变换。2.一步步安装1.安装包下载这一步很简单,只要搜索fftw,很容易就可以定原创 2016-11-13 11:05:59 · 12569 阅读 · 2 评论 -
FFTW在VS中集成与测试详解
1.头文件、库文件、动态链接库的处理1.1 头文件fftw3.h文件是在解压文件包时候得到的,我们需要把它配置到VS安装目录的#include内,如下图:2.库文件库文件就是我们辛辛苦苦弄出来的那三个家伙,包括:libfftw3-3.lib libfftw3f-3.lib libfftw3l-3.lib,它们不同之处在于出炉数据的精度。我们要把他们装在#lib文件下原创 2016-11-13 11:24:09 · 3579 阅读 · 0 评论 -
Savitzky-Golay 滤波器详解及C/matlab语言程序设计
1.Savitzky-Golay 滤波器Savitzky-Golay滤波器(通常简称为S-G滤波器)最初由Savitzky和Golay于1964年提出,发表于Analytical Chemistry 杂志。之后被广泛地运用于数据流平滑除噪,是一种在时域内基于局域多项式最小二乘法拟合的滤波方法。这种滤波器最大的特点在于在滤除噪声的同时可以确保信号的形状、宽度不变。使用平滑滤波器对原创 2016-11-29 13:37:59 · 53679 阅读 · 11 评论 -
IIR+全通滤波器实现相位平衡_matlab仿真
1.全通系统实现相位均衡本文从逆向角度出发,基于搜索最优化策略对全通系统相位平衡问题进行了初步研究。主要思想:在频域中对它在单位圆上抽样N个点,按照所要求的频率特性,找出一个全通滤波器,使其满足一定的相位特性。最后利用最优化(例如遗传算法)策略进行迭代,找到满足要求的最优解。适用于低通系统相位补偿的二阶全通滤波系统如下:2.设计满足特定相位条件的全通滤波器2原创 2016-11-30 17:02:56 · 16274 阅读 · 23 评论 -
IIR+全通滤波器级联实现系统零相位相移_matlab仿真
1.前言前面详细的介绍了如何通过优化的思想逆向设计符合要求的全通相位平衡系统!实际上,线性相位的要求要比零相位相移的要求苛刻的多。晚上和好友解释了一下如何利用优化思想实现线性相位,好友感觉很难实现零相移。为此,一步一步又做了一次仿真实验,并进行记录。说白了,y=0就是y=kx的一个特例!明白了这一点,还是建议大家看上一排帖子,平衡后达到线性相位更重要!2.一步一步实现零原创 2016-12-01 00:05:32 · 11291 阅读 · 4 评论 -
信号分解:标架、对偶标架、紧标架
1.前言信号分解或信号变换的基本思路是将信号x(t)和一组函数(或向量)做內积,从而得到一组分解系数an。分解(或变换)的目的是研究原始信号中有哪些有哟用的信息,并讨论如何抽取这些有用的信息。我能能够理解,正交基具有很多优点(信息不冗余,对偶基是本身),实际应用中也是最广泛的,可惜的是,在实际工作中,发现并得到一组好的正交基往往是不容易的。正式正交基,或者更广泛地说,分解对研究原创 2017-03-07 09:16:56 · 2466 阅读 · 0 评论