Chapter4-数字相关和数字卷积概述.ppt

4.7 相关技术的应用…… （2）信号为白噪声的自相关函数 设有一功率谱为 ，则自相关函数为： * “白”的含义指两点之间没有任何相关性，带限白噪声，功率谱为矩形波。sinc函数性状，就是m=0有最大值，m足够大时趋近0。 一个观测信号x(n)实际上是周期信号s(n)和随机信号w(n)的叠加加过程，x(n)= s(n)+ w(n)。 如果信号和噪声互不相关，则自相关函数： rxx(m)= rss(m)+ rww(m) m足够大， rxx(m)不为0，可以判定周期信号存在！ */58 4.7 相关技术的应用…… 【例4-13】设周期信号s(n)=0.8sin(πn/5)，噪声w(n)为随机产生的白噪声，观测信号x(n)= s(n)+ w(n)。 * %产生s信号 clear m=1:300; for n=1:300 s(n)=0.8*sin(pi*n/5); end %产生随机白噪声 w=randn(1,300); %生成观测信号 x=s+w; %产生自相关函数 rww=xcorr(w); V=circlel(s); rss=s*V;%或者rss=xcorr(s); for k=301:599 rss(k)=rss(k-300+1); end rxx=xcorr(x); %生成原始信号图 figure,subplot(3,1,1),plot(m,w) subplot(3,1,2),plot(m,s) subplot(3,1,3),plot(m,x) %生成自相关函数图 figure,subplot(3,1,1),plot(-299:299,rww) subplot(3,1,2),plot(-299:299,rss) subplot(3,1,3),plot(-299:299,rxx) %扩大观测信号的自相关函数 figure,plot(120:299,rxx(120:299)) */58 4.7 相关技术的应用…… 结果： * */58 4.7 相关技术的应用…… 【例4-14】设x(n)= e-0.05ncos(πn/6)，y(n)=1.2x(n-n0),它们的波形如下图所示，试估计延迟n0。 * %生成信号x和y m=0:23; for n=1:24 x(n)=exp(-0.05*n)*cos(pi*n/6); end for n=5:24 y(n)=1.2*x(n-4); end y(1:4)=0; figure,subplot(2,1,1),stem(m,x),title(信号x) subplot(2,1,2),stem(m,y),title(信号y) %求信号x和y的互相关函数rxy和ryx rxy=xcorr(x,y); ryx=xcorr(y,x); figure,subplot(2,1,1),stem(-23:23,rxy),title(信号x和y的互相关函数rxy) subplot(2,1,2),stem(-23:23,ryx),title(信号x和y的互相关函数ryx) */58 4.7 相关技术的应用…… 结果： * */58 本章小结 * 1、掌握：线性相关、线性卷积、相干函数和相干系数； 2、熟悉：循环相关和循环卷积； 3、了解：相关技术的应用。 */58 本章习题 * 用直接法和表格法做如下例题： 1、设x(n)和y(n)是有限长的序列，x(n)=[2,0.1,-2,0.1], y(n)=[0.1,2,0.1,-2], 求两序列的线性相关、循环相关、线性卷积和循环卷积函数。 */58 * 下集预告 第五章 维纳滤波 */58 * 实验二详解…… % 选择信号类型并设定参数，产生信号x(n) clear; clc; disp(请选择信号); disp(1 ---- 伪随机序列randn()); disp(2 ---- 实际测量的心电信号); disp(3 ---- 实际测量的脑电信号); b = input(信号：); switch b % 输入序号，产生相应信号，L=128，N=8 case 1 L = input(每段数据长度 L \n); N = input(数据共多少段 N \n); x = randn(1, L*N); case 2 load ecgdata; display([数据总长度,num2str(length(ecgdata)),点]); L = input(每段数据长度 L \n); N = in