1、近代数字信号处理近代数字信号处理(Advanced Digital Signal Processing)电子电子信息工程学院信息工程学院信号与图像处理研究室简要回顾简要回顾)()(0nXtx)j()(Xtx)e(jXkxmXkx连续周期信号连续非周期信号离散非周期信号离散周期信号离散非周期频谱连续非周期频谱连续周期频谱离散周期频谱四种信号的时域与频域对应关系 什么是信号的抽样 为什么要进行抽样 如何进行信号抽样 信号抽样理论分析 抽样定理工程应用 基于信号时域分析和频域分析，以全新的方式揭示了信号时域抽样定理的本质。信号的时域抽样()t kTx kx t什么是信号时域抽样x,Fs,Bits=w
2、avread(myheart);sound(x,Fs);Fs=22,050Hz;Bits=16什么是信号时域抽样什么是信号时域抽样 Nyquist，美国物理学家，1889年出生在瑞典。他对信息论做出了重大贡献。1907年移民到美国并于1912年进入北达克塔大学学习。1917年在耶鲁大学获得物理学博士学位。19171934年在AT&T公司工作，后转入Bell电话实验室工作。1927年，Nyquist确定了对某一带宽的有限时间连续信号进行抽样，为不使原波形产生“半波损失”，抽样率至少应为信号最高频率的2倍，这就是著名的Nyquist抽样定理。Harry Nyquist(1889 1976)谁提出信
3、号抽样概念谁提出信号抽样概念V.A.Kotelnikov，科捷利尼科夫.第一个准确系统地阐述抽样定理并应用于通信工程领域的科学家，他证明了低通信号和带通信号的抽样理论。1 V.A.Kotelnikov,“On the transmission capacity of ether and wire in electro-communications,1933.Claude Elwood Shannon香农,1916-2001,出生于美国密歇根州，1940年获MIT数学博士学位和电子工程硕士学位。香农发表了两篇革命性文章，奠定了信息论的基础。l C.E.Shannon,“A mathematica
4、l theory of communication,”1948.2 C.E.Shannon,“Communication in the presence of noise”,1949.谁提出信号抽样概念谁提出信号抽样概念为什么进行信号抽样离散系统抽样重建输入x(t)xkyk输出y(t)利用数字技术处理模拟信号信号稳定性好:数据用二进制表示，受外界影响小。系统的精度高:可通过增加字长提高系统的精度。系统灵活性强:改变系统的系数使系统完成不同功能。信号可靠性高:存储无损耗，传输抗干扰。离散信号与系统的主要优点：信号处理简便:信号压缩，信号编码，信号加密等。数字信号的主要特点数字信号的主要特点u 易
5、于存储易于存储u 易于传输易于传输u 易于处理易于处理易于存储模拟存储数字存储数字存储优点：1.容量大2.稳定性好纸张相机胶卷磁带录像带DVD硬盘SD卡存储方式晶体管技术微电子技术数字信号的主要特点数字信号的主要特点容量大易于存储序号存储方式名称存储容量需要数量1模拟A4纸约2000字/页2415页2模拟胶卷36页A467卷3数字CD700MB1.4张4数字U盘8GB0.12个5数字硬盘1TB0.0009个例如：中国最长的古典小说榴花梦共483万字。存储一个汉字需要2个字节，约计943MB。700M容量CD需1.4张8G容量U盘需0.12个数字信号的主要特点数字信号的主要特点易于存储序号序号存
6、储方式存储方式名称名称存储容量存储容量需要数量需要数量1模拟胶卷36张/卷21760卷2数字CD700MB5595张3数字8G U盘8 GB478个4数字1T硬盘1 TB4个例如：中国植物图像库，783240幅高清图像。若存储一张图像需要5MB，共计约391.6万MB=3916 GB。容量大8G容量U盘需478个1T容量硬盘需4个胶卷21769卷数字信号的主要特点数字信号的主要特点1885年4月15日人类历史上第一张声音签名音盘，Bell博士(1847.3.31922.8.2)易于存储稳定性好榴莲歌 CD音质榴莲歌 磁带音质数码存储方式便于数据长期稳定保存老照片数码照片（空间离散）数字信号的主
7、要特点数字信号的主要特点模拟信号传输原理图数字信号传输原理图模拟信源A/D转换D/A转换信道模拟信宿模拟信源信道模拟信宿易于传输数字传输优点：抗干扰能力强数字信号的主要特点数字信号的主要特点模拟信号传输受到的噪声干扰数字信号传输受到的噪声干扰易于传输数字传输下，信号幅值采用二进制0、1编码。相比模拟信号而言，抗干扰能力强。抗干扰能力强数字信号的主要特点数字信号的主要特点易于传输抗干扰能力强模拟信号传输数字信号传输数字传输有效抵抗信道干扰，音质好例如：中央人民广播电台开篇录音数字信号的主要特点数字信号的主要特点易于传输传输距离远数字信号噪声消除模拟信号噪声累积数字信号的主要特点数字信号的主要特点
8、易于传输模拟电视信号数字电视信号数字信号的主要特点数字信号的主要特点模拟信号加密/解密实现非常困难数字信号传输易于加密/解密易于传输相比模拟信号而言，数字传输下，数字信号可采用加入伪随机码等多种加密方式，易于实现数字信号的保密传输。加密序列保密性好数字信号的主要特点数字信号的主要特点 信息编码加密解密信道信息解码信源信宿加密密钥解密密钥数字信号加密示意图易于传输保密性好数字信号的主要特点数字信号的主要特点语音信号加密语音伪随机码易于传输保密性好数字信号的主要特点数字信号的主要特点无法加密加密图像，窃取无法恢复易于传输保密性好数字信号的主要特点数字信号的主要特点易于处理模拟信号处理数字信号处理连
9、续系统h(t)连续信号x(t)连续信号y(t)离散系统hk离散信号xk离散信号yk计算机技术数字信号的主要特点数字信号的主要特点易于处理动画制作 帧是动画电影中的单幅图像传统动画制作：手工绘制图像帧按动画每秒播放24帧计算，1分钟的动画大概需要1440幅单帧图像！数字动画制作：计算机插值产生中间帧在拍摄动画片大闹天宫时，几十位动画工作者花了近两年的时间才完成，其中绘制了600多万张图画。提高制作效率，节约制作成本！数字信号的主要特点数字信号的主要特点wav文件(未经压缩)压缩编码易于处理压缩编码例如：音乐我的祖国参数：时长8秒，双声道，16为采样位数，44.1kHz采样频率数据量(字节/秒)=
10、(采样频率采样位数声道数)/8mp3文件(经过压缩)在基本音质保持不失真情况下，压缩率达6.5:1 207KB/1.35MB=15.3%数字信号的主要特点数字信号的主要特点易于处理识别检索采集定位特征提取特征匹配检索结果服务器可实现1秒内匹配1亿幅图像千万级规模图片特征提取可在数小时内完成！数字信号的主要特点数字信号的主要特点数字图像取证水印技术数字图像取证水印技术原始图像水印图像复制-粘贴攻击后的篡改图像定位篡改通过水印检测图像被篡改，定位篡改区域即使图像被篡改，仍可提取出可识别的版权信息嵌入水印信息“北京交大”人眼无法区别嵌入水印图像与原始图像数字信号的主要特点数字信号的主要特点a.原始图
11、像随机信息和图像随机信息和图像块平均值信息两块平均值信息两种水印信号嵌入种水印信号嵌入b.水印图像复制复制-粘贴粘贴攻击攻击利用嵌入的随利用嵌入的随机信息进行篡机信息进行篡改检测与定位改检测与定位d.定位篡改c.篡改后的图像数数字字图图像像认认证证水水印印技技术术e.恢复图像利用嵌入的平利用嵌入的平均值信息对篡均值信息对篡改的图像块进改的图像块进行恢复行恢复数字信号的主要特点数字信号的主要特点易于处理滤波去噪 含噪声图像 滤波后图像数字信号的主要特点数字信号的主要特点抽样频率对图像的影响信号抽样的基本原理信号抽样的基本原理kTttxkx)(如何选取抽样间隔如何选取抽样间隔T？信号抽样的基本原理
12、信号抽样的基本原理传统模型传统模型新模型新模型输入和输出都是连续时间信号输入是连续时间信号，输出是离散时间信号信号抽样的基本原理信号抽样的基本原理抽样定理证明模型kTkTtkTxttxtx)()()()()(sam)(*)j(21)j(samsamsamnXXn)(j 1samnXTn-j-jjsam(j)()e e(e)k TkkkXx kTx kX信号抽样的基本原理信号抽样的基本原理)()(j(1 )e(samjTnXTXn抽样定理的本质：信号时域的离散化导致其频域的周期化抽样定理的本质：信号时域的离散化导致其频域的周期化 而而 sam2 m只是上述基本结论针对只是上述基本结论针对带限信号
13、带限信号的特例。的特例。kTttxkx)(如果：如果：则有：则有：信号抽样的基本原理信号抽样的基本原理 若带限信号x(t)的最高角频率为m，则在满足一定条件下，信号x(t)可以用等间隔T的抽样值唯一表示。fsam=2fm 为最小抽样频率，称为为最小抽样频率，称为Nyquist Rate.抽样间隔抽样间隔T T 需满足：需满足：mm/1/(2)Tffsam 2fm （或或 sam 2m）信号抽样的基本原理信号抽样的基本原理单边带信号单边带信号与与窄带高频信号窄带高频信号的抽样问题的抽样问题msamsam2Bm1000k rad/s,8k rad/sB信号抽样的基本原理信号抽样的基本原理窄带高频信
14、号的抽样窄带高频信号的抽样 fsam=8kHz 抽样后的频谱。fsam=56kHz 抽样后的频谱。中心频率24kHz，带宽8kHz。解调后语音信号抽样后的语音信号(不解调)解调后语音信号信号抽样的基本原理信号抽样的基本原理fm=28 kHzfsam=8 kHz信号抽样的基本原理信号抽样的基本原理连续时间信号离散时间信号信号抽样的工程实现信号抽样的工程实现抽样保持0010010001000110011110011000011101000010数字信号量化编码()t kTx kx t抽样间隔(周期)T (s)抽样角频率 sam=2p/T (rad/s)抽样频率 fsam=1/T (Hz)通道数最高
15、抽样率表示精度信号抽样的工程实现信号抽样的工程实现许多实际工程信号不满足带限条件许多实际工程信号不满足带限条件)j(Hmm10 抗混低通抗混低通 滤波器滤波器)(tx)(1tx)(th信号抽样的工程应用信号抽样的工程应用混叠误差混叠误差 vsvs 截断误差截断误差信号抽样的工程应用信号抽样的工程应用语音信号的抗混叠滤波抽样频率fsam=44,100 Hz抽样频率fsam=5,512 Hz抽样频率fsam=5,512 Hz抽样前对信号进行了抗混叠滤波信号抽样的工程应用信号抽样的工程应用例1 已知实信号x(t)的最高频率为fm(Hz)，试计算对各信号x(2t)，x(t)*x(2t)，x(t)x(2
16、t)抽样不混叠的最小抽样频率。对信号x(2t)抽样时，最小抽样频率为 4fm(Hz)；对x(t)*x(2t)抽样时，最小抽样频率为2fm(Hz)；对x(t)x(2t)抽样时，最小抽样频率为 6fm(Hz)。解：根据信号时域与频域的对应关系及抽样定理得：x(t)xkyky(t)利用离散系统处理连续时间信号 生物医学信号处理 铁路控制信号识别5.抽样定理工程应用生物医学信号处理生物医学信号处理生物神经细胞（元）结构图 5.抽样定理工程应用生物医学信号处理生物医学信号处理AdLink PCI 9112A/D,D/A CardPersonal ComputersIn WindowOperationEn
17、vironmentsAIAODOABCBDB生物信号采集系统组成框图 5.抽样定理工程应用生物医学信号处理生物医学信号处理采集的生物信号的模式识别5.抽样定理工程应用生物医学信号处理生物医学信号处理神经元等效电路 Gion1 Gion2 Gionm Eion1 Eion2 Eionm CM Iex Ges1 Ges2 Gesn V1 V2 Vn Gcs1,1 Gcs1,2 Gcs1,p Ecs1,1 Ecs1,2 Ecs1,p Gcsn,1 Gcsn,2 Gcsn,p Ecsn,1 Ecsn,2 Ecsn,p Ionic conductances Electricalsynapses(es)C
18、hemicalsynapses(cs)+5.抽样定理工程应用 铁路控制信号识别铁路控制信号识别机车信号信号抽样的工程应用信号抽样的工程应用多制式列车控制信号的频谱铁路控制信号识别铁路控制信号识别信号抽样的工程应用信号抽样的工程应用 传统的车载信号系统，由于安全性和可靠性等技术局限，只能将车载信号作为辅助信号，司机必须瞭望地面信号机来驾驶列车。国际公认160km/h以上的高速列车运行已不能靠司机瞭望地面信号方式保证安全，而必须以车载信号作为主体信号来控制列车。信号抽样的工程应用信号抽样的工程应用铁路控制信号识别铁路控制信号识别铁路控制信号的频谱分析信号抽样的工程应用信号抽样的工程应用对连续带限信
19、号进行抽样时，只需抽样速率 fs 2fm。在工程应用中，抽样速率为何常设为 fs(35)fm？时域抽样问题的探究时域抽样问题的探究研究性课题l Nyquist抽样率是信号精确恢复的充分条件，但不是必要条件。l 信号非均匀抽样的问题。l 近年来，由D.Donoho、E.Candes及华裔科学家T.等人提出了一种新的信息获取指导理论，即，压缩感知或压缩传感压缩感知或压缩传感 (Compressive Sensing(CS)or Compressed Sampling)l CS理论主要包括信号的稀疏表示、编码测量和重构算法。相关知识拓展信号的频域抽样mnNnmNnxmXX2j2je)e(ZrNkrN
20、kn,1,1,0,令)(10rNkmNrNkWrNkxmXmkNNNkWkx10)(DFS)e(2jkxmXXNmN 信号的频域抽样即对非周期序列xk的频谱X(ej)在每个周期2p内均匀抽样N点。将将xk以以N为周期进行周期化为周期进行周期化 xk+3-4-3-2-1012345671221 x k -4-3-2-1012345671221-4-3-2-10123456721122112214kx xk-3-4-3-2-1012345671221-4-3-2-1012345672222222222223kx 当序列长度不超过当序列长度不超过N时时,周期化后的序列和原序列一个周期化后的序列和原序
21、列一个周期内的值相同。周期内的值相同。当序列长度超过当序列长度超过N时，周时，周期化后的序列会出现混叠期化后的序列会出现混叠(aliasing)。已知有限序列xk=-1,-1,4,3；k=0,1,2,3，序列xk的DTFT为X(ej)。记X(ej)在W=2p m/3;m=0,1,2的取样值为Xm，求IDFSXm。IDFTXm =xk+xk+3 =2,-1,4;k=0,1,2 X(ej)在频域的离散在频域的离散化导致对应的时域序列化导致对应的时域序列xk的周期化的周期化.信号的频域抽样X(ej)在频域的离散化导致对在频域的离散化导致对应的时域序列应的时域序列xk的周期化。的周期化。x(t)在时域
22、的离散化导致对应在时域的离散化导致对应的频谱函数的频谱函数X(j)的周期化。的周期化。和和为利用数字化方式为利用数字化方式分析和处理信号奠定了理论基础。分析和处理信号奠定了理论基础。CTFTDTFTIDTFTIDFS nTnXTX)2j(1)j(周期化)(kxtx时域抽样)e(jmXX频域抽样 nnNkxkx周期化四种信号的时域与频域对应关系tx(t)0X(j)0t0)(txkxk00kxkX(n0)0FTFSDTFTDFS抽样信号重建)/)(Sa)(rTkTtkxtxk 信号重建理想D/A输 入与输出的频谱关系)e()j(jrecTXX零阶保持零阶保持D/A零阶保持D/A模型框图 零阶保持D
23、/A输出信号的频谱 2/jze)2/(Sa)j(TTTH)(j)(j)(jzz X HX离散域进行补偿的FIR和IIR滤波器21116189161)(zzzH12125.018/9)(zzH抽样连续信号的离散处理连续信号的离散处理对如图所示系统，已知输入模拟信号对如图所示系统，已知输入模拟信号X(j)，离散系统，离散系统H(ej)，试画出，试画出X(ej)、Y(ej)和和Y(j)。1 1)(jeHp pc c p p 1)(jXsam/2sam/2对如图所示系统，已知输入模拟信号对如图所示系统，已知输入模拟信号X(j)，离散系统，离散系统H(ej)，试画出，试画出X(ej)、Y(ej)和和Yr(j)。(1)已知已知DF的截频的截频 c，等效的，等效的AF的截频的截频 c=c/T(2)要求等效要求等效AF的截频为的截频为 c，DF的截频应设计为的截频应设计为 c=T c 抽样后的谱抽样后的谱DF滤波后的谱滤波后的谱例：已知连续带通信号x(t)的频谱如下图所示,试分别画出sam1=0.5m 及sam2=0.8m时，抽样后离散序列的频谱。解：sam1=0.5m,T1=2p/sam1=4p/msam2=0.8m,T2=2p/sam2=2.5p/m
