MATLAB的这个程序每一段分别是什么意思？

这是一个循梁搜环计算。

s1=s(Start:Start+N-1)-s(Start-2*N:Start-2*N+N-1)

s2=s(Start+N:Start+2*N-1)-s(Start-2*N:Start-2*N+N-1)

s3=s(Start+2*N:Start+3*N-1)-s(Start-2*N:Start-2*N+N-1)

分别在s中计算出s1，s2，s3的数值

deltaS=zeros(1,3*N)

设定一个1*3N（矩阵一列，长度指盯为3*N）的矩阵deltaS

deltaS(1:N)=s1

deltaS(N+1:2*N)=s2

deltaS(2*N+1:3*N)=s3

将所有前面计算橡逗历出的数据放到矩阵deltaS

% 将以下程序copy到matlab编辑器中羡裤运行，或衡颂直接在工作区运行即可

fs = 20e3 % 采样频率

fn = 3e3 % 固有频率

y0 = 5 % 位移常数

g = 0.1% 阻尼系咐派郑数

T = 0.01 % 重复周期

N = 4096 % 采样点数

NT = round(fs*T) % 单周期采样点数

t = 0:1/fs:(N-1)/fs % 采样时刻

t0 = 0:1/fs:(NT-1)/fs % 单周期采样时刻

K = ceil(N/NT)+1 % 重复次数

y = []

for i = 1:K

y = [y,y0*exp(-g*2*pi*fn*t0).*sin(2*pi*fn*sqrt(1-g^2)*t0)]

end

y = y(1:N)

Yf = fft(y) % 频谱

figure(1)

plot(t,y)

axis([0,inf,-4,5])

title('轴承故障仿真信号时域波形图')

xlabel('Time(s)')

ylabel('Amplitude')

figure(2)

f = 0:fs/N:fs-fs/N

plot(f/1e3,abs(Yf))

xlabel('Frequency(KHz)')

ylabel('\itY\rm(\itf\rm)')

title('轴承故障仿真信号幅度谱图')

1、出图后输搏敬入第1行代码

2、在图片中点选要测距的两点

3、输入第2行代码（输出的基链慎距离单位是km）以上就是Matlab在唤稿两条线上取点测距的方法。

---