matlab用自己编写的自相关法对信噪比 5dB 0dB -5dB的雷达信号进行检测

以下是使用自相关法检测信噪比为5dB、0dB、-5dB的雷达信号的Matlab代码示例：

% 生成信号
fs = 100e3; % 采样率
T = 1/fs; % 采样间隔
t = 0:T:10e-3; % 时间向量
f0 = 10e3; % 载频频率
snr = [-5, 0, 5]; % 信噪比
x1 = sin(2*pi*f0*t); % 信号
x2 = awgn(x1, snr(1), 'measured'); % 加入噪声
x3 = awgn(x1, snr(2), 'measured'); % 加入噪声
x4 = awgn(x1, snr(3), 'measured'); % 加入噪声

% 自相关检测
% 自相关函数的计算
R1 = xcorr(x1, x1, 'biased'); 
R2 = xcorr(x2, x2, 'biased'); 
R3 = xcorr(x3, x3, 'biased'); 
R4 = xcorr(x4, x4, 'biased'); 

% 画出自相关函数图像
figure;
subplot(2,2,1);
plot(R1);
xlabel('Lag');
ylabel('Amplitude');
title(['SNR = ' num2str(snr(1)) 'dB']);

subplot(2,2,2);
plot(R2);
xlabel('Lag');
ylabel('Amplitude');
title(['SNR = ' num2str(snr(2)) 'dB']);

subplot(2,2,3);
plot(R3);
xlabel('Lag');
ylabel('Amplitude');
title(['SNR = ' num2str(snr(3)) 'dB']);

subplot(2,2,4);
plot(R4);
xlabel('Lag');
ylabel('Amplitude');
title('Original Signal');

% 自相关函数的峰值位置即为信号的时延
[~,I1] = max(R1);
[~,I2] = max(R2);
[~,I3] = max(R3);
[~,I4] = max(R4);

% 画出信号和检测结果
figure;
subplot(2,2,1);
plot(t,x1);
xlabel('Time (s)');
ylabel('Amplitude');
title('Original Signal');

subplot(2,2,2);
plot(t,x2);
xlabel('Time (s)');
ylabel('Amplitude');
title(['SNR = ' num2str(snr(1)) 'dB']);

subplot(2,2,3);
plot(t,x3);
xlabel('Time (s)');
ylabel('Amplitude');
title(['SNR = ' num2str(snr(2)) 'dB']);

subplot(2,2,4);
plot(t,x4);
xlabel('Time (s)');
ylabel('Amplitude');
title(['SNR = ' num2str(snr(3)) 'dB']);

% 画出检测结果
figure;
plot(snr, [I1,I2,I3]);
xlabel('SNR (dB)');
ylabel('Delay (samples)');
title('Detection Results');

运行上述代码，即可得到自相关函数图像和检测结果图像。其中，自相关函数的峰值位置即为信号的时延，通过将不同信噪比下的检测结果绘制在一张图中，可以直观地观察信噪比对检测效果的影响。