close all;clear all;clc
T = 1.e-6; % 信号持续时间
B = 50e6; % 信号带宽
K = B/T; % 调频率
Fc = 1.e9 ; % 信号载频
Fs = 200.e6; % 采样频率
dt = 1/Fs; % 采样间隔
N = ceil(T/dt); % 采样点数
t = ((0:N-1)-N/2)/N*T; % 时间轴
c = 300000000;
tt = linspace(-0.5,0.5,N);
st = exp(1j*pi*(2*Fc*t + K*(t).^2)); % 生成信号
%加窗
h = kaiser(length(t), 2.5);
st = st.*h';
%输出信号的实部和虚部
st_real = real(st);
st_imag = imag(st);
figure;
subplot(2,1,1),plot(st_real);
title('输入信号的实部');ylabel('幅度');
subplot(2,1,2),plot(st_imag);
title('输入信号的虚部'),ylabel('幅度');
%信号的幅度谱和相位谱
st_fft = fftshift(fft(fftshift(st)));
st_fft_range = abs(st_fft); %图像幅度谱,加log便于显示
st_fft_range = st_fft_range / max(st_fft_range);
st_fft_phase = unwrap(angle(st_fft));%*180/pi; %图像相位谱
figure;
subplot(2,1,1),plot(st_fft_range);
title('输入信号的幅度谱');ylabel('幅度');
subplot(2,1,2),plot(st_fft_phase);
title('输入信号的相位谱'),ylabel('幅度');
%脉冲压缩后波形
%脉冲压缩
ht = conj(fliplr(st)); % 匹配滤波器
st_comp = fftshift(ifft(fft(st).*fft(ht))); % 计算循环卷积
st_comp_range = abs(st_comp);
st_comp_range = st_comp_range/max(st_comp_range);
st_comp_phase = angle(st_comp)*180/pi; %图像相位谱
figure;
subplot(2,1,1),plot(st_comp_range);
title('脉冲压缩后的幅度');ylabel('幅度');
subplot(2,1,2),plot(st_fft_phase);
title('脉冲压缩后的相位');ylabel('幅度');
%目标响应的分辨率,峰值旁瓣比,积分旁瓣比
[main_peak, main_location] = max(findpeaks(st_comp_range));
IndMin = find(diff(sign(diff(st_comp_range)))>0)+1; %获得局部最小值的位置
IndMax = find(diff(sign(diff(st_comp_range)))<0)+1; %获得局部最大值的位置
figure; hold on; box on;
plot(1:length(st_comp_range),st_comp_range);
plot(IndMin,st_comp_range(IndMin),'r^')
plot(IndMax,st_comp_range(IndMax),'k*')
legend('曲线','波谷点','波峰点')
title('计算离散节点的波峰波谷信息', 'FontWeight', 'Bold');
%峰值旁瓣比
PSLR = 20*log10(st_comp_range(IndMax(main_location - 1))/st_comp_range(IndMax(main_location )))
%积分旁瓣比
first_valley = 0;
for i = 1:length(IndMin)-1;
if IndMin(i) < IndMax(main_location) && IndMin(i + 1) > IndMax(main_location);
first_valley = i;
break;
end
end
main_part_power = sum(st_comp_range(IndMin(first_valley):IndMin(first_valley + 1)));
total_power = sum(st_comp_range);
ISLR = 10*log10(( total_power - main_part_power )/ main_part_power )
%目标响应的分辨率,找3db带宽
bound_3db = 0;
for i = IndMin(first_valley):IndMax(main_location);
if st_comp_range(i) <= 1/2* main_peak && st_comp_range(i+1) >= 1/2* main_peak;
bound_3db = 2*(IndMax(main_location) - i);
break;
end
end
tao = bound_3db * T / 200;
fenbianlv = c*tao/(2)
脉冲压缩.rar_alonezet_beautifulyja_sevenci1_脉冲压缩图像_运动补偿 帧
版权申诉
54 浏览量
2022-07-15
01:48:00
上传
评论
收藏 1KB RAR 举报
weixin_42653672
- 粉丝: 93
- 资源: 1万+
最新资源
- Libraries-Comm-Controller
- 豆瓣电影爬虫 爬取top电影的评论 + 每个用户的看过的电影的评论 用于推荐系统的 协同过滤+源代码+文档说明
- 交互设计课程竞品分析内容案例设计
- c07c4b30caf2ab290c3f2eea8339b34b.mp4
- emqx服务器搭建文件
- Libraries-Comm-Controller-DOC-V2-0-1-en.pdf
- update9-20240601.5.205.slice.img.7z.003
- 9f9ae03ea06c5c991afa26c5813d8831.amr
- 计算机等级考试Python二级真题.zip
- 爱普生L301清零程序
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈