MATLAB自相关编写及周期函数分析

MATLAB自相关编写及周期函数分析MATLAB自相关编写及周期函数分析班级:机自23姓名:于滨学号:21201010681.在Matlab环境下编写自相关分析程序,并给出周期信号、含随机噪声的正弦信号的相关分析结果。b5E2RGbCAP答:绘制纯净正弦函数自相关的程序,如下:dt=.1;采样间隔t=[0:dt:100];x=sin(t);生成正弦函数[a,b]=xcorr(x,'unbiased');对正弦函数进行无偏自相关subplot(4,1,1);p1EanqFDPwplot(t,x);绘制正弦函数图subplot(4,1,2);plot(b*dt,a);绘制正弦函数自相关图绘制加了噪声的正弦函数自相关程序如下:dt=0.1;采样间隔t=[0:dt:100];x=sin(t);生成正弦函数y=randn(1,1001);z=awgn(x,15,'measured');subplot(4,1,3);plot(t,z);绘制添加噪声的正弦函数Title添加噪声的正弦信号[a,b]=xcorr(z,'unbiased');求有噪声正弦函数的无偏自相关subplot(4,1,4);DXDiTa9E3dplot(b*dt,a);绘制自相关图形title绘制结果如下:分析:正弦函数的自相关图像具有明显的周期性,添加随机噪声后自相关图像的周期性下降,可以看出在机械设备运行过程中可以通过采集相关信号进行自相关分析,当自相关图像无任何周期性时,可以初步认为干扰主要是由于周边环境无规则干扰所致;如果其自相关图像具有明显的周期性,可以初步认为干扰是其内部故障,产生周期性的干扰。RTCrpUDGiT2。绘制三角波。t=0:0.0001:1.5;x=sawtooth(2*pi*30*t,0.5);subplot(3,1,1);plot(t,x),axis([00.4-11])fsnt=fftshift(fft(x));pfsnt=abs(fsnt);subplot(3,1,2);plot(pfsnt),axis([50001000008000]);绘制结果如下: