參與評論

評論

登錄后參與評論

全部評論(1條)

• 

  FyfyY5H1h1 2016-12-25 00:00:00

  下面matlab程序，按照你的要求編寫的。fs=100;N=200;lag=100;%randn('state'，0); %設(shè)置產(chǎn)生隨機數(shù)的初始狀態(tài)（噪聲）n=0:N-1;t=n/fs;x1=sin(2*pi*10*t)+2*sin(2*pi*15*t) y=fft(x1，N); %計算頻譜mag=abs(y);f=n*fs/N;figure(1)subplot(2，2，1)，plot(t，x1);title('原始信號')，xlabel('時間/s');subplot(2，2，2)，plot(f(1:N/2)，mag(1:N/2)*2/N);title('頻譜圖')，xlabel('頻率/Hz');[c，lags]=xcorr(x1，lag，'unbiased'); %自相關(guān)subplot(2，2，3)，plot(lags/fs，c);title('自相關(guān)信號')，xlabel('時間/s');m=length(c)-1; z=fft(c，m); mag_z=abs(z);ff=(0:m-1)*fs/m;subplot(2，2，4)，plot(ff(1:m/2)，mag_z(1:m/2)*2/m);title('自相關(guān)頻譜')，xlabel('頻率/Hz');figure(2)[c1，lags1]=xcorr(x1，randn(1，length(t))，lag，'unbiased'); m1=length(c1)-1; z1=fft(c1，m1); mag_z1=abs(z1);ff1=(0:m1-1)*fs/m1;subplot(2，1，1)，plot(lags1/fs，c1);title('互相關(guān)信號')，xlabel('時間/s');subplot(2，1，2)，plot(ff1(1:m1/2)，mag_z1(1:m1/2)*2/m1);title('互相關(guān)頻譜')，xlabel('頻率/Hz');

  贊(15)

  回復(fù)(0)

  評論

  評論

  登錄后參與評論