MATLAB实验六。2..docxVIP

一：实验题目6：（提取淹没在噪声中的周期信号）提取淹没在噪声中的周期信号要求自相关（以判断周期）和互相关（以恢复信号自身）。（a）生成信号x[n] =sin(0.1*n*pi),0n499.增加一些均匀分布的随机噪声（其噪声幅度为2，平均值为0）以得到噪声信号s[n]。作出每个信号随参数n变化的图形。是否能够从所画的x[n]中判断出任何的周期性？如果可以，求周期N。能否由所画得s[n]判断出任何周期性？(b)求x[n]的周期自相关rpx[n]并作图。能否由所作的rpx[n]判断出任何周期性？如果可以，求周期N。这个N是否与x[n]的周期相同？(c)利用上面所求的N值(如果没有，则由x[n]判断)生成一个500样本的冲击串, 0≤n≤499。求出s[n]和i[n]的周期互相关。二：实验目的：（1）熟悉一些常用函数的生成方法。（2）练习用matlab绘制图形。（3）用自相关和互相关提取淹没在噪声中的周期信号。三：实验原理：周期信号和随机信号都是功率信号（功率信号是信号功率为有限正值而信号能量无限大）。自相关函数反映了一个信号与延迟k后的同一信号的相关程度。互相关函数反映了一个信号和延迟k后的另一信号的相关程度。一个功率信号的功率谱密度与信号的自相关函数是一对傅立叶变换关系。由周期函数的巴赛法尔定理，功率谱密度函数在离散的采样点上与信号的频谱有一定的关系。如果能够在功率谱密度函数上离散采样，构造出频谱，然后再对其进行傅立叶逆变换，即可得到时域的函数曲线。同理，两个相同周期的功率信号的互功率谱与他们的周期互相关是一对傅立叶变换。四：实验步骤简述：（1）生成信号x[n] =sin(0.1*n*pi),0n499。和增加一些均匀分布的随机噪声（其噪声幅度为2，平均值为0）以得到噪声信号s[n]。（2）求x[n]的周期自相关（自相关函数xcorr()）rpx[n]并作图。并读出周期。（3）生成一个500样本的冲击串, 0≤n≤499。计算是s(n)与i(n)的互相关xcorr(s,i)。五：程序框图：六：源程序：close all;clear all;clc;%生成信号x[n]n=0:499;x=sin(0.1*n*pi);subplot(6,1,1);plot(n,x,k);title(正弦信号x[n]);%生成随机噪声r[n]r=2*rand(1,500)-1;subplot(6,1,2);plot(n,r,k);title(随机噪声r[n]);%生成噪声信号s[n]s=x+r;subplot(6,1,3);plot(n,s,k);title(噪声信号s[n]);%求x[n]的自相关函数rpx[n]rx=xcorr(x,x);subplot(6,1,4);plot(rpx);title(x[n]的自相关函数rpx[n]);%生成500样本冲击串i=ones(1,500); %产生500样本的冲击串Ifor m=1:500if mod(m,20)==0 i(1,m)=1;else i(1,m)=0;endendsubplot(6,1,5);stem(i,.);axis([0,500,-0.2,1.2]);title(冲击串i[n]);%求s和i的互相关rsi[n]rsi=xcov(s,i);subplot(6,1,6);plot(rsi);title(s[n]和i[n]的互相关函数rsi[n]);七：程序结果及图表：八：实验总结(a)由图像可知，x[n]具有周期性，且周期N=20；不能从s[n]中判断出周期性。(b) x[n]的自相关函数具有周期性，周期N=20，与x[n]周期相同。(c)s[n]和i[n]的互相关函数具有周期性，且周期N=20，与x[n]的周期相同。可以估计，冲击信号可以从噪声信号中提取出原始信号。1、此实验考查周期函数和随机噪声的产生，常见二维平面图形的绘制方法，以及自相关函数和互相关函数的用法。需要一定Matlab基础知识，以及它的操作，画出图形。而且要对自相关函数和互相关函数了解，通过学习，进一步复习了周期信号的周期性的判断，知道了函数自相关和互相关的内涵，以及它们所实现的算法实质。2、作为初次使用Matlab实现信号处理的题目，在实验过程中，我还是遇到了许多问题，例如对某些Matlab函数的功能不够清晰，对数字信号处理的一些相关知识也有些遗忘，所以在复习数字信号处理后，我积极的向同学和老师请教，在他们的帮助下，我顺利完成了题目要求。一：实验题目7：该题目的目的是说明一个PN扩频信号在抑制正弦干扰中的有效性。现考虑下图所示的二进制通信系统，对信号发生器的输出乘上一个二进制（1）PN序列。同一个二进制PN序列用来与解调器输入相乘，因此消除了这个PN序列在期望信号上的