現(xiàn)代信號處理實驗一平穩(wěn)信號產(chǎn)生及時頻分析實驗結(jié)果-

1、 ylabel('y(t' title('orignal signal' grid; subplot(3,1,3 plot(lags/Fs,c,'r' xlabel('t' ylabel('Rxy(t' title('correlation' grid; orignal signal 1 x(t 0 -1 0 0.2 0.4 0.6 0.8 1 1.2 t orignal signal 1.4 1.6 1.8 2 0.5 y(t 0 -0.5 0 0.2 0.4 0.6 0.8 1 1.2 t co

2、rrelation 1.4 1.6 1.8 2 0.5 Rxy(t 0 -0.5 -0.4 -0.3 -0.2 -0.1 0 t 0.1 0.2 0.3 0.4 圖1.5(a 信號x(t的波形圖 （b）信號y(t的波形圖 （c）信號x(t和信號y(t的互相關(guān)波形圖 3、 協(xié)方差實現(xiàn) 在 MATLAB 中計算平穩(wěn)信號序列協(xié)方差的函數(shù)為 xcov，該函數(shù)的調(diào)用格式為： v=xcov(x,y v=xcov(x v=xcov(x,option c,lags=xcov(x,y,maxlags c,lags=xcov(x,maxlags c,lags=xcov(x,y,maxlags,options 協(xié)方

3、差函數(shù) xcov 的參數(shù)定義與相關(guān)函數(shù) xcorr 的一致。 例 9產(chǎn)生兩個正弦信號加白噪聲的樣本數(shù)據(jù)，并繪出其互相關(guān)函數(shù)與互協(xié)方差的估計值曲線圖。 解：程序為： % exa9.m % to test xcorr function; clear N=256; f=0.1; 16 a1=5; a2=3; Mlag=N/4; %產(chǎn)生兩個正弦加白噪聲的數(shù)據(jù)； x=a1*sin(2*pi*f*(0:N-1+2*randn(1,N; y=a2*sin(2*pi*f*(0:N-1+randn(1,N; subplot(2,2,1; plot(x(1:N/2; grid on; title('原始信

4、號1' subplot(2,2,2; plot(y(1:N/2; grid on; title('原始信號2' %求這兩個數(shù)據(jù)向量的互相關(guān)函數(shù)； corrxy=xcorr(x,y,Mlag,'unbiased' subplot(2,2,3; plot(-Mlag:1:Mlag,corrxy; grid on; title('互相關(guān)函數(shù)' %求這兩個數(shù)據(jù)向量的互協(xié)方差函數(shù)； covxy=xcov(x,y,Mlag,'biased' subplot(2,2,4; plot(-Mlag:1:Mlag,covxy; title(&

5、#39;互協(xié)方差函數(shù)' 運行結(jié)果為： 原始信號1 10 5 0 -5 -10 0 10 原始信號2 5 0 50 100 150 -5 0 50 100 150 互相關(guān)函數(shù) 10 5 0 -5 -10 -100 10 5 0 -5 -10 -100 互協(xié)方差函數(shù) -50 0 50 100 -50 0 50 100 圖1.6 (a原始信號1x(n波形圖 （b）原始信號2y(n波形圖 (c信號1與信號2的互相關(guān)函數(shù)波形圖 （d）信號1與信號2的互協(xié)方差函數(shù)波形圖 17 5、相關(guān)函數(shù)的應(yīng)用 試編程序完成。 在某個音樂廳內(nèi)，原始音頻信號的回音由于墻壁和天花板等的反射而產(chǎn)生，聽眾所感受到的音頻信 號是 x(n 和它的回音的合成。令 y (n = x(n + ax(n - k 其中， x(n =