奈奎斯特采样率与压缩感知学习报告

数字信号处理第一次大作业奈奎斯特采样定理与信号稀疏采样学习报告专业：信息对抗技术学生姓名：石星宇02123010指导教师：吕雁目录奈奎斯特采样定理与信号稀疏采样学习报告 (1)一、奈奎斯特采样定理 (1)1、奈奎斯特采样定理说明 (1)2、信号的采样与恢复 (1)3、相关代码 (3)4、关于奈奎斯特采样定理的一些问题 (5)二、信号稀疏采样 (5)1、为什么要提出信号的稀疏采样 (5)2、压缩感知概述 (6)3、压缩感知基本概念 (6)4、压缩感知仿真 (7)5、压缩感知仿真程序 (8)三、总结 (9)四、参考资料 (10)奈奎斯特采样定理与信号稀疏采样学习报告一、奈奎斯特采样定理1、奈奎斯特采样定理说明采样过程所应遵循的规律，称为取样（采样）定理、抽样定理。

采样定理说明采样频率与信号频率之间的关系，是连续信号离散化的基本依据。

在进行模拟/数字信号的转换过程中，当采样频率s f 大于等于信号中最高频率c f 的2倍时，采样之后的数字信号完整地保留了原始信号中的信息，可由采样得到的数字信号恢复原来的模拟信号。

一般实际应用中保证采样频率为信号最高频率的5～10倍。

采样定理又称奈奎斯特采样定理。

将c s f f 2=称为奈奎斯特频率。

2、信号的采样与恢复结合实例，说明奈奎斯特采样定理与内插恢复的应用。

假设有模拟信号()t f t f t x a 212cos 2cos ππ+=，其中Hz f Hz f 50,2021==。

该信号波形及频谱如下图所示：对信号()t f t f t x a 212cos 2cos ππ+=以采样频率为Hz f f s 10022==进行采样，得到如下所示的离散时间信号，即序列()s s nT f nT f nT x 212cos 2cos ππ+=，其中s s f T /1=。

该序列的频谱如下：由此可见，采样过程对原始信号的频谱有一定的影响。

但是随着采样频率的逐渐增加，会使得采样信号的频谱与原始信号的频谱逐渐接近。

现在利用内插公式对采样得到的离散时间信号进行恢复。

定义内插函数为()TtT t t ππsin=Φ则()()()()()()()[]()TnT t T nT t nT x nT t nT x t nT x t x n n a //sin --=-Φ=Φ*=∑∑∞-∞=∞-∞=ππ根据上式，便可由采样得到的序列()nT x 完整的恢复出原始信号()t x a 。

下面给出利用MATLAB 计算的结果：从上图可以看出，利用内插公式，完整地将原始信号恢复了出来。

3、相关代码 close all clear allclcdf=0.5;%频率分辨率tp=1/df;%保证df所需的信号持续时间t=linspace(0,tp,1024);%连续时间变量f1=20;f2=50;%信号频率fc=max(f1,f2);%信号最高频率fs=2*fc;%采样率ts=1/fs;%采样间隔N=2^ceil(log2(fs/df));n=1:N;xa=cos(2*pi*f1*t)+cos(2*pi*f2*t);%模拟信号xn=cos(2*pi*f1*n*ts)+cos(2*pi*f2*n*ts);%模拟信号plot(t,xa,'r')hold onstem(n*ts,xn,'b')legend('模拟信号','采样信号')title('模拟信号和采样信号')xlabel('t'),ylabel('x(t)')axis([0 tp/2 min(xa)-0.2 max(xa)+0.2])figure(2)subplot(211)fftxa=fft(xa);fa=0.5*(-length(t)/2:length(t)/2-1*fs/length(t));plot(fa,fftshift(abs(fftxa)),'r')%模拟信号频谱title('模拟信号频谱')xlabel('f/Hz')subplot(212)fftxn=fft(xn);fn=0.5*(-length(n)/2:length(n)/2-1*fs/length(n));plot(fn,fftshift(abs(fftxn)))%采样信号频谱title('采样信号频谱')xlabel('f/Hz')xaa=zeros(1,length(t));for tt=1:length(t)%计算采样内插值xaa(tt)=0;for n=1:Nxaa(tt)=xn(n)*(sin(pi*(t(tt)-n*ts)/ts)/(pi*(t(tt)-n*ts)/ts))+xaa(tt);endendfigure(3)plot(t,xaa,'b')title('采样内插恢复信号')xlabel('t/s'),ylabel('x(t)')axis([0 tp/2 min(xaa)-0.2 max(xaa)+0.2])4、关于奈奎斯特采样定理的一些问题假设模拟信号为()Hz f t f t x a 202cos 11==，π，用奈奎斯特频率对其采样，发现采样点处的取值均为零（如下图），因此用这些采样点是无法恢复原始信号的。

这也就是为什么实际中采用的采样频率要大于奈奎斯特频率的原因。

此外，实际中我们处理的信号不可能是简单的正弦信号，因此遇到采样点均为零的情况几乎不可能，上述只是一个特例。

二、信号稀疏采样1、为什么要提出信号的稀疏采样 首先考虑奈奎斯特采样定理的几点缺陷：（1）采样率不得低于信号最高频率的两倍，这使得硬件系统面临很大的采样速率压力；（2）在压缩编码过程中，为了降低存储、处理和传输的成本，大量变换计算得到的小系数被丢弃，造成了数据计算和内存资源的浪费。

综合上述两点，人们便提出这样的问题：能否利用其它变换域描述信号，建立新的信号描述和处理理论框架，使得在保证信息不损失的情况下，用远低于奈奎斯特频率的频率去采样信号，同时可以完全恢复信号？答案是肯定的，这就是由 E. J. Candes 、J. Romberg 、T. Tao 和D. L. Donoho 等科学家在2004 年提出的压缩感知理论（Compressed sensing ）。

2、压缩感知概述压缩感知（Compressed sensing ），也被称为压缩采样(Compressive sampling)，稀疏采样(Sparse sampling)。

它作为一个新的采样理论，它通过开发信号的稀疏特性，在远小于Nyquist 采样率的条件下，用随机采样获取信号的离散样本，然后通过非线性重建算法完美的重建信号。

压缩感知理论的核心思想主要包括两点。

第一是信号的稀疏结构。

另外一点是不相关特性。

稀疏信号的有用信息的获取可以通过一个非自适应的采样方法将信号压缩成较小的样本数据来完成。

用数学模型对压缩感知的主要内容进行如下描述：（1）信号稀疏表示问题（稀疏变换）：对于信号N R s ∈，如何找到某个正交基ψ，使其在ψ上的表示是稀疏的（2）信号低速采样问题（非相关测量）：如何设计一个平稳的、与变换基ψ不相关的N M ⨯维的观测矩阵Φ，保证稀疏向量θ从N 维降到M 维时重要信息不遭到破坏（3）信号重构问题（重构算法）：如何设计快速重构算法从线性观测θθΘ=Φψ=y 中恢复信号。

3、压缩感知基本概念设信号R s ∈是一维实值离散信号，在正交基（稀疏基）ψ上可以稀疏表示，即()()()θψ=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=N s s s s 21 其中[]}{()N i i N ，，2,1,|||21=ψψψψ=ψ是1⨯N 的向量，稀疏系数s s T i i i ψ=ψ=,θ。

当信号在稀疏基上只有K 个非零系数时，属于严格稀疏的情况。

多数情况下信号无法满足严格稀疏的要求，但仍具有可压缩性，即信号的变换系数经排序后以指数级衰减并趋近于零时，信号是可以近似稀疏表示的。

测量矩阵有观测波形和采样方式决定。

观测波形一般包括独立同分布的高斯随机波形、伯努利随机波形和正交函数系等；采样方式包括均匀采样、随机采样和jitter 采样。

将信号s 投影到一组测量矩阵[]TM φφφ 21=Φ上，则测量值,,Tm m s y φ=即：1111||||||||⨯⨯⨯⨯⨯⨯⨯⨯Θ=ψΦ=Φ=N N M N N N N M N N M M s y θθ 也即θΘ=y重构算法是压缩感知的另一个关键因素。

目前的重构算法有贪婪算法（又称匹配追踪（Matching Pursuit ，MP ），正交匹配追踪（Orthogonal Matching Pursuit ，OMP ）、凸优化算法（最小1l 范数）和统计优化算法（Sparse Bayesian ）等。

4、压缩感知仿真与奈奎斯特采样定理仿真相同，仍然设模拟信号()t f t f t x a 212cos 2cos ππ+=，对其进行压缩感知采样并重建，并与原始信号对比如下：程序同时给出重构误差0.1032。

需要注意的是，每次运行程序所得的重构误差是不同的，是因为信号的重构过程中有随机因素在里面。

5、压缩感知仿真程序% 1-D信号压缩传感的实现(正交匹配追踪法Orthogonal Matching Pursuit)% 测量数M>=K*log(N/K),K是稀疏度,N信号长度,可以近乎完全重构% 编程人--香港大学电子工程系沙威 Email: wsha@eee.hku.hk% 编程时间：2008年11月18日% 文档下载: http://www.eee.hku.hk/~wsha/Freecode/freecode.htm% 参考文献：Joel A. Tropp and Anna C. Gilbert% Signal Recovery From Random Measurements Via Orthogonal Matching % Pursuit，IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 53, NO. 12, % DECEMBER 2007.clc;clear%% 1. 时域测试信号生成K=7; % 稀疏度(做FFT可以看出来)N=256; % 信号长度M=64; % 测量数(M>=K*log(N/K),至少40,但有出错的概率)f1=20; % 信号频率1f2=50; % 信号频率2fs=100; % 采样频率Ts=1/fs; % 采样间隔n=1:N; % 采样序列x=cos(2*pi*f1*n*Ts)+cos(2*pi*f2*n*Ts) % 完整信号,由2个信号叠加而来%% 2. 时域信号压缩传感Phi=randn(M,N); % 测量矩阵(高斯分布白噪声)64*256的扁矩阵，Phi也就是文中说的D矩阵s=Phi*x.'; % 获得线性测量，s相当于文中的y矩阵%% 3. 正交匹配追踪法重构信号(本质上是L_1范数最优化问题)%匹配追踪：找到一个其标记看上去与收集到的数据相关的小波；在数据中去除这个标记的所有印迹；不断重复直到我们能用小波标记“解释”收集到的所有数据。