我是数字信号处理的新手,正在阅读一篇论文。我不知道他们如何从信号中提取一些特征,如均方根 (RMS)、恒定误报率 - CFAR、平均值分散统计 - MVD、功率比 - ROP 和信号的自相关,如图所示在下图中:
谁能帮我解释一下如何,如果可能的话,你能告诉我一些代码来获取它们吗?
我试图编写一些代码来计算这些参数。你能看一下吗:
有效值:
for i= 1:floor(time(end)/tPeriod)
tBegin = (i-1)*tPeriod;
tEnd = i*tPeriod;
index = find((time>=tBegin)&(time<tEnd));
RMS(i) = sqrt(mean(data(index).^2)); %equation 1
tPlot = [tPlot tEnd];
end
CFRA:
for i= 1:floor(time(end)/tPeriod)
tBegin = (i-1)*tPeriod;
tEnd = i*tPeriod;
indexTime = find((time>=tBegin)&(time<tEnd));
y = fft(data(indexTime),M)/M; % M is the number of fft points
y = 2*abs(y);
df = fMax/M;
f = [0:M-1]*df;
indexF = find((f>=n1)&(f<n2));
CFAR(i) = sum(y(1:M/2)).^v;
tPlot = [tPlot tEnd];
end
罗普:
for i= 1:floor(time(end)/tPeriod)
tBegin = (i-1)*tPeriod;
tEnd = i*tPeriod;
indexTime = find((time>=tBegin)&(time<tEnd));
y = fft(data(indexTime),M)/M;
y = 2*abs(y);
df = fMax/M;
f = [0:M-1]*df;
indexF = find((f>=n1)&(f<n2));
ROP(i) = sum(y(indexF).^2)/sum(y.^2);
tPlot = [tPlot tEnd];
end