数字信号处理
文章平均质量分 66
thefutureisour
ZTE基带部码农,平时工作很忙,不再更新opencv相关内容,也不解答大家问题了。
展开
-
1 生成基本序列
数字信号处理的对象,是对模拟信号采样、量化后形成的数字信号,在一般的书中,列举的信号包括,单位冲击序列、单位阶跃序列、矩形序列、正弦序列、指(复)数序列。下面演示如何通过matlab生成这些序列。首先要考虑的是,在课本的学习中,我们是通过“函数”的角度认识他们的,即给出它们的表达式;而在计算机中,能够处理的是他们的值。(这里并不考虑matlab的符号代数系统)。所以,要想完整的表示序列,需要两原创 2014-01-29 15:15:17 · 3693 阅读 · 0 评论 -
2.序列的运算
序列的运算主要包括基本的算术运算,以及对于翻转、平移、卷积、自相关、互相关。先说代数运算。对两个序列做加减乘除用笔算都是一件很容易的事情,写程序算也比较容易,关键是能够将序列“对齐”:假如一个序列的下标是-3~4;另一个是1~5。那么就得把它们长度补的一样长,并且覆盖二者的范围,然后再进行运算。下面是对应的matlab程序:function [ xe1,xe2,n ] = dataAlig原创 2014-01-30 15:37:46 · 4400 阅读 · 1 评论 -
4 利用计算系统的输出
这里主要介绍滤波函数filter。比如如果系统函数是y(n)+0.62y(n-1)+0.13y(n-2)=x(n-2)输入是x(n)=u(n)-u(n-15)那么求输出对应的程序是:n = -5:25;N = 15;step = [n>=0];stepN = [(n-N)>=0];x = step-stepN;b = [0,0,1];a = [1,0.62,0.13原创 2014-02-01 20:36:37 · 1308 阅读 · 0 评论 -
3 利用传输函数函数分析系统性能
分为两大部分,1是计算系统的幅度响应和相位响应,2是计算系统的零点、极点。举个例子,如果系统的差分方程为:y(n)-0.53y(n-1)+0.73y(n-2) = 0.14[x(n)-x(n-2)]matlab提供了计算频率响应的函数freqz,使用时只用输入传输函数的系数即可:a = [1,-0.53,0.73];b = [1,0,-1];[h,w] = freqz(b,a)原创 2014-01-31 10:56:27 · 2714 阅读 · 0 评论 -
5 DFT与功率谱估计
首先说明频率盒问题。假设数据有N点,采样频率为fs,那么DFT的分辨率为fs/N,在没有噪声的情况下;如果信号恰好位于fs/N*k的频率上,那么就会在这个频点处为最大值,其他频点处为0。例如,采样频率=8000,信号采样点为64个。那么各个频点处的值为8000/64=125的整数倍。如果信号频率恰好为1500(12倍),那么在该点处幅度为最大值,其他点出为0,如图所示:但是,如原创 2014-03-02 23:01:32 · 5693 阅读 · 1 评论 -
6 DFT估计功率谱的增益
前一小节介绍了基本的使用DFT估计功率谱的知识,这一小节对于一些具体问题进行探讨。如何用DFT获得增益?总的来说,有两种方式:相关增益和累积增益。相关增益的获得,是通过增加信号的采样点数获得的。因为对于DFT,可以将每个X(k)看做是一个窄带滤波器,增加DFT点数,使得滤波器的通带变小,这样就能滤去更多的噪声,使得原来淹没在噪声中的信号被检测出来。例如,对于一个单音信号,进行64点DFT时原创 2014-03-03 23:46:49 · 3707 阅读 · 1 评论 -
7 数字滤波器的结构
有人可能会很纠结,滤波不就是输入信号和一组抽头系数做卷积么,滤波器的结构是来干什么的?总体上说,是为了减少抽头系数和时延器。比如,由于某些特殊要求,可以使得滤波器的抽头的系数可以有某种对称性质,(比如线性相位的FIR滤波器的系数就是这样),为了节省乘法单元,就可以使用一些特殊的结构来巧妙地实现原来较多的运算。再比如,对于IIR滤波器,可以让输入信号的时延与反馈信号的时延公用时延器。先说FI原创 2014-03-30 23:02:51 · 3396 阅读 · 0 评论