利用互相关测试系统单位冲击响应.docVIP

随机信号分析实验 -------利用互相关测试系统的单位冲击响应 实验小组成员： 小组成员班级： 目录 一、 实验目的 3 二、 实验原理 3 三、 实验任务与要求 4 四、 实验步骤 5 ? 高斯白噪声的导入 5 ? 通过系统 6 ? 通过被测系统后的信号与理想高斯白噪声进行互相关 8 ? 通过低通滤波器得输出信号 8 五、 计算x(t)、noise(t)、y(t)信号的均值（数学期望）、方差、相关函数、概率密度、频谱及功率谱密度等 9 1. noise(t)（白噪声） 9 2. x(t) 11 3. y(t) 15 实验六 利用互相关测试系统单位冲击响应（实验报告） 实验目的 实验目的 ⑴ 研究分析电子系统受随机信号激励后的响应及测量方法。 ⑵ 了解随机信号的特性，包括均值（数学期望）、方差、相关函数、概率密度、频谱及功率谱密度等。 实验原理 利用互相关算法可以求取线性时不变系统的冲击响应。通过被测系统后的理想高斯白噪声信号与理想高斯白噪声进行互相关运算，产生相应的输出通过一个低通滤波器，获得线性系统单位冲激响应h(t)。其原理框图如图1所示： 图1 利用互相关测量线性系统单位冲击响应 实验任务与要求 所有电路均用matlab或c/c++编程实现。 输入信号为：高斯白噪声 图 高斯白噪声 被测系统： ① 被测系统是一个低通滤波器。低通滤波器的通带为0KHz-1KHz、通带衰减小于1db、阻带衰减大于35db。 ② 被测系统是一个微分器。低通滤波器的通带为0KHz-500Hz、通带衰减小于1db、阻带衰减大于35db。 如果被测系统是低通滤波器，它的冲击响应h(τ)为撒函数。实验结果y(t)应该为图12所示： 图12 的输出 同理，如果被测系统是微分器，z输出是它的冲击响应h(τ)。 实验步骤 高斯白噪声的导入 %*************导入高斯白噪声***************% noise=randn(10000,1); t=0:1/10000:1; t=t(1:10000); figure; plot(t,noise,-r); xlabel(t);ylabel(幅度(v)); title(白噪声); Noise=fft(noise);%得noise(t) cmo=abs(Noise); f2=(0:length(Noise)-1)*10000/length(Noise); figure; plot(f2(1:length(f2)/2),cmo(1:length(f2)/2),-b); title(白噪声的频谱); %**************系统为低通滤波器***************% Fs=10000;T=1/Fs; %下面设计低通滤波器并滤波 wp=2*600/Fs; %边界频率关于pi归一化 ws=2*1000/Fs; Rp=1;As=40; [N,wc]=buttord(wp,ws,Rp,As);%计算数字滤波器阶数N和3dB截止频率wc [B,A]=butter(N,wc); %设计数字滤波器 [H,w]=freqz(B,A); figure; plot(w/pi,abs(H)); title(系统的频谱); xlabel(\omega/\pi); h=ifft(H);%求得h(t) x2=conv(h,noise);%卷积得x(t) X2=fft(x2); cmo=abs(X2); f2=(0:length(X2)-1)*10000/length(X2); figure; plot(f2(1:length(f2)/2),cmo(1:length(f2)/2),-b); title(x(t)的频谱); %**************系统为微分器***************% Fs=20000;fp=1000;rp=0.1;fs=5000;rs=30; wp=2*fp/Fs;ws=2*fs/Fs; [N,wpo]=ellipord(wp,ws,rp,rs); [Bz,Az]=ellip(N,rp,rs,wpo,high); wk=0:pi/512:pi; [Hz,w]=freqz(Bz,Az,wk); Hx=angle(Hz); figure; plot(wk,abs(Hz));%grid on; xlabel(\omega^pi);ylabel(幅度); title(系统的频谱特性); h=ifft(H);%求得h(t) x2=conv(h,noise);%卷积得x(t) X