《合成孔径雷达成像算法与实现》Figure3.7

代码复现如下:

Matlab 复制代码
clc
clear all
close all

%参数设置
TBP = 100;                                              %时间带宽积
T = 10e-6;                                              %脉冲持续时间

%参数计算
B = TBP/T;                                              %信号带宽
K = B/T;                                                %信号调频频率
alpha_os = 8;                                           %过采样率
F = alpha_os*B;                                         %采样率
N = 2*ceil(F*T/2);                                      %采样点数
dt = T/N;                                               %采样时间间隔
df = F/N;                                               %采样频率间隔

%变量设置
t = -T/2:dt:T/2-dt;                                     %时间变量
f = -F/2:df:F/2-df;                                     %频率变量
t_out = linspace(2*t(1),2*t(end),2*length(t)-1);        %循环卷积后的信号长度

%信号表达
st = exp(1j*pi*K*t.^2);                                 %chirp信号复数表达式
ht = conj(fliplr(st));                                  %时域匹配滤波器表达式
%在信号中加入高斯噪声的两种方式
st_noise = awgn(st,-2.5,'measured');
%st_noise = st+0.75*randn(1,N);
s_out_noise = conv(st_noise,ht);
s_out_noise_nor = s_out_noise/max(s_out_noise);
s_out_noise_log = 20*log10(abs(s_out_noise_nor)/max(abs(s_out_noise_nor))+eps);

%绘图
figure
subplot(221),plot(t*1e+6,real(st_noise))
subplot(222),plot(t_out*1e+6,s_out_noise_log),axis([-1 1,-30 5])
subplot(223),plot(t_out*1e+6,real(s_out_noise_nor)),axis([-3 3,-0.5 1.5])
subplot(224),plot(t_out*1e+6,angle(s_out_noise_nor)),axis([-1 1,-5 5])

1. st_noise = awgn(st,-2.5,'measured');

y = awgn(x,snr)将白高斯噪声添加到向量信号x中。标量snr指定了每一个采样点信号与噪声的比率,单位为dB。如果x是复数的,awgn将会添加复数噪声。这个语法假设x的能量是0dBW。 y = awgn(x,snr,'measured')和y = awgn(x,snr)是相同的,除了agwn在添加噪声之前测量了x的能量。

相关推荐
დ旧言~23 分钟前
【高阶数据结构】图论
算法·深度优先·广度优先·宽度优先·推荐算法
张彦峰ZYF28 分钟前
投资策略规划最优决策分析
分布式·算法·金融
The_Ticker44 分钟前
CFD平台如何接入实时行情源
java·大数据·数据库·人工智能·算法·区块链·软件工程
爪哇学长1 小时前
双指针算法详解:原理、应用场景及代码示例
java·数据结构·算法
Dola_Pan1 小时前
C语言:数组转换指针的时机
c语言·开发语言·算法
繁依Fanyi2 小时前
简易安卓句分器实现
java·服务器·开发语言·算法·eclipse
烦躁的大鼻嘎2 小时前
模拟算法实例讲解:从理论到实践的编程之旅
数据结构·c++·算法·leetcode
xrgs_shz2 小时前
MATLAB读入不同类型图像并显示图像和相关信息
图像处理·计算机视觉·matlab
C++忠实粉丝2 小时前
计算机网络socket编程(4)_TCP socket API 详解
网络·数据结构·c++·网络协议·tcp/ip·计算机网络·算法
用户37791362947553 小时前
【循环神经网络】只会Python,也能让AI写出周杰伦风格的歌词
人工智能·算法