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

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

班级：机自23
姓名：于滨
学号：2120101068

1.在Matlab环境下编写自相关分析程序,并给出周期信号、含随机噪声的正弦信号的相关分析结果。
答：绘制纯净正弦函数自相关的程序，如下：
dt=.1; 采样间隔
t=[0:dt:100];
x=sin(t); 生成正弦函数
[a,b] = xcorr(x,'unbiased'); 对正弦函数进行无偏自相关
subplot(4,1,1);
plot(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);
plot(b*dt,a); 绘制自相关图形
title
绘制结果如下：

分析：
正弦函数的自相关图像具有明显的周期性，添加随机噪声后自相关图像的周期性下降，可以看出在机械设备运行过程中可以通过采集相关信号进行自相关分析，当自相关图像无任何周期性时，可以初步认为干扰主要是由于周边环境无规则干扰所致；如果其自相关图像具有明显的周期性，可以初步认为干扰是其内部故障，产生周期性的干扰。2。绘制三角波。
t = 0:0.0001:1.5;
x = sawtooth(2*pi*30*t,0.5);
subplot(3,1,1);
plot(t,x),axis([0 0.4 -1 1])
fsnt = fftshift(fft(x));
pfsnt = abs(fsnt);
subplot(3,1,2);
plot(pfsnt),axis([5000 10000 0 8000]);
绘制结果如下：