matlab实验报告数字信号的产生和频谱分析.doc

实验一数字信号的产生和频谱分析 实验目的 掌握采样定理内容，熟悉连续信号经理想采样前后的频谱变化关系，加深时域 采样定理的理解。 熟悉FFT算法原理和FFT子程序的应用。 学习用FFT对连续信号和吋域离散信号进行谱分析的方法，了解可能出现的分 析误差及其原因，以便在实际中正确应用FFT。掌握利用FFT进行信号谱分析 的原理。 理解信号通过高斯0噪声信道的原理。 实验要求 按照采样定理生成CW信号和LFM信号； 画出信号时域波形图和频谱图； 生成高斯分布的白噪声； 生成一定信噪比的带噪信号，并对其进行谱分析。 基本原理 所谓抽样。就是对时间连续的信号隔一定的时间间隔T抽取一个瞬时幅 度值在一个频带限制在(0，t)内的时间连续信号f (t)，如果以fs的抽 样频率对它进行抽样，那么根据这些抽样值就能完全恢复原信号。根据抽 样定理只有当发生fs^2fm,抽样后的信号才能含原连续的全部信息。 A 对于一个连续信号Aw进行理想采样的过程可用下式表示： Xa(t) = Xa(t)p(t) A 其屮是对的理想采样，为周期冲激脉冲，即 ”?= 一oo A A 的傅里叶变换x (人仰为 A 上式表明的周期延拓，其延拓周期为采样角频率 。只有满足采样定理时，才不会发生频率混叠失真。 将②代入①进行傅里叶变换可得 之(》= 卯-nT)]e-?t =-00 oo ZZ1II II ZZ1 II II xa(nT)e~jQnT 式中的久就是采样后得到的序列x ，即^) = xa(nT) x(n)的傅里叶变换为 x(n)的傅里叶变换为 OO X(ej(0)= X*細 /2=-oo z\ 比较两式可得 Xa(jco) = X(eJ ) |^Qr 这说明两者之间只在频率度量上差一个常数因子T 为了在计算机上观察分析各种序列的频域特征，通常对％^；在 [0, 2 JT]上进行M点采样来观察分析。对长度为N的有现场序列x(n), 有 N—\ X(e*) = Zx(m)e’kn =0 -K中 M ，k=o, 1,…，M-1 通常M应取得大一些，以便观察谱的细节变化。取模\(^)可绘出幅 频特性曲线。 一个时域离散线性非移变系统的输入/输出关系为 y(n) = x(n)* h(n) = V x(m)h(n - m) /=-oo 如果h(n)和x(n)的长度分别为N和M，则y (n)的长度为L=M+N_1。这样 上式的卷积运算就是序列移位、相乘和累加的过程。 上述卷积过程也可以在频域实现 Y(ejM) = X(ej0))H(ejM) 上式右边的相乘是在各频点｛o k｝上的频谱值相乘 在本实验中选取抽样频率为10倍的信号频率。本试验中所用到的函数 分别为 CW 信号表达式：s(t)=cos (2 Ji ft+l) |t|T/2 LFM 信号表达式:s(t)=cos[2H (ft+kt2/2 )+ (1)]; |t|T/2 ;k=B/T;B 为 带宽。 我们用十倍于信号频率的采样频率对信号进行采样，再利用FFT变换， 使其变换为离散的功率谱函数。随机产生高斯白噪声，叠加在原噪声上， 再进行频谱分析，两者可进行比较。 实验流程 按题目要求初始化参数f=100Hz，Fs=1000Hz. 产生CW信号并画出信号的时域图和频谱图。 3. 图1-1生成CW信号 4.产生LFMs (t) =cos(2*pi*(f*t+k*t*t/2)+pi/3)信号并画出信号的吋域图和频谱 图O 5. 图1-2生成LFM信号 调用matlab的rand函数生成高斯分布的白噪声。 通过给定的信噪比后改变CW信号的幅值，SNR=a,由于SNR=101g (Pa/o2 ) 高斯白噪声的功率为o2=l,因此计算出Pa=10 (SNR/10),就为CW信号的 幅值。 将产生的高斯白噪声叠加在CW信号上，再画出加入噪声后的CW信号的时 域图与频谱图。 图1-2画带噪信号频谱图 图1-3程序流程图 五.实验结果及分析（在matlab下仿真的波形） 15 cS 0 5 1 ? ■ -O 0.05 0.1 0.15 0.2 时间 60 40 20 屮心频率f=100F 300 线性调频信竹时域波形 0.05 0.1 0.15 0.2 时问 1 5 o 5 -O. 30 lb -300 -200 ?100 0 100 频率/Hz LFM信号频谱图 200 ?1000 0 频率/Hz 1000 2 O -2000 2000 图1-1 CW侣号和LFM fe号波形及频谱图 如上图可以看出在采样频率为1000时能够准确的岡出正弦信号的波形图， 从频谱图中可以看出CW信号的中心频率为100Hz，LFM波形在时域上形成一 个有规则的包络，从频谱可以看出屮心频率在600Hz,与计算的理论值： cos(2*pi*(f*t+k*t*