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

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的互相關波形圖 3、 協方差實現 在 MATLAB 中計算平穩(wěn)信號序列協方差的函數為 xcov，該函數的調用格式為： 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 協方

3、差函數 xcov 的參數定義與相關函數 xcorr 的一致。 例 9產生兩個正弦信號加白噪聲的樣本數據，并繪出其互相關函數與互協方差的估計值曲線圖。 解：程序為： % exa9.m % to test xcorr function; clear N=256; f=0.1; 16 a1=5; a2=3; Mlag=N/4; %產生兩個正弦加白噪聲的數據； 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' %求這兩個數據向量的互相關函數； corrxy=xcorr(x,y,Mlag,'unbiased' subplot(2,2,3; plot(-Mlag:1:Mlag,corrxy; grid on; title('互相關函數' %求這兩個數據向量的互協方差函數； covxy=xcov(x,y,Mlag,'biased' subplot(2,2,4; plot(-Mlag:1:Mlag,covxy; title(&

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