预测表面混乱力量Range-Doppler空间

频率= 10 e9;%赫兹alt = 3 e3;% mhelperSurfaceDopplerLimitsChart(频率、alt)

dopres = 78;%赫兹速度= 120;% m / s潜水= 10;%度

R = 5 e3:40:10e3;

放牧= grazingang (alt, R,“平”);反射= surfaceReflectivityLand;nrc =反射(放牧、频率);

[rcs,夹住]= clutterSurfaceRangeDopplerRCS (nrc, R,频率,dopres, alt、速度、潜水);helperPlotSurfaceRCS (R,夹住,rcs)

beamw = 10;%度部= 30;%度helperBeamDopplerLimitsChart(频率、alt、速度、beamw dep)

阿兹= 0;%度(R1、R2 d1, d2) = helperConicalBeamLimits(频率、alt、速度、潜水、beamw az, dep)

helperPlotSurfaceRCS (R,夹住,rcs, R1, R2, d1, d2)

阿兹= -60;%度(R1、R2 d1, d2) = helperConicalBeamLimits(频率、alt、速度、潜水、beamw az, dep);helperPlotSurfaceRCS (R,夹住,rcs, R1, R2, d1, d2)

bandw = 0(1,元素个数(R));为gIdx = 1:元素个数(R) I1 =找到(rcs (gIdx:), 1,“第一”);I2 =找到(rcs (gIdx:), 1,“最后一次”);bandw (gIdx) = (I2 - I1 + 1) * dopres;结束

pf = polyfit (R, bandw 2);图绘制(bandw e3 / 1, R / 1 e3)在情节(polyval e3 (pf, R) / 1, R / 1 e3);持有从包含(的带宽(赫兹))ylabel (的距离(公里))传说(带宽估计的,“二次健康”,“位置”,“西北”甘氨胆酸)组(,“XLimitMethod”,“紧”)

helperDopplerWrappingExample

numdop = 128;

脉冲重复频率= dopres * numdop

[rcs,夹住]= clutterSurfaceRangeDopplerRCS (nrc, R,频率,dopres, alt,速度,潜水,“NumDopplerBins”,numdop);

helperPlotSurfaceRCS (R,夹住,rcs)

频率= 10 e9;%赫兹rngres = 80;% mdopres = 31;%赫兹startrange = 10 e3;% mnumrange = 400;R = startrange + (0: numrange-1) * rngres;

alt = 4 e3;% m速度= 240;% m / s潜水= 0;%度

tgtrng = 14 e3;% mtgtgndrng =√tgtrng ^ 2-alt ^ 2);tgtpos = [0 tgtgndrng 0];tgtrcs = 6;% dBsm

rdrpos = [0 0 alt];rdrvel =速度* [cosd(潜水)0信德(潜水)];λ= freq2wavelen(频率);tgtdop = 2 /λ*点(rdrvel tgtpos-rdrpos) / tgtrng

放牧= grazingang (alt, R,“平”);nrc =反射(放牧、频率);[rcs,夹住]= clutterSurfaceRangeDopplerRCS (nrc, R,频率,dopres, alt、速度、潜水);

rcs = rcs / 2;

rcs0Hz = rcs(:,夹住= = 0);

[~,tgtrngbin] = min (abs (tgtrng-R));surfrcs = pow2db (rcs0Hz (tgtrngbin));tgtscr = tgtrcs - surfrcs

rngres = 30;% mR = startrange + (0: numrange-1) * rngres;放牧= grazingang (alt, R,“平”);nrc =反射(放牧、频率);[rcs,夹住]= clutterSurfaceRangeDopplerRCS (nrc, R,频率,dopres, alt、速度、潜水);rcs = rcs / 2;rcs0Hz_highRes = rcs(:,夹住= = 0);[~,tgtrngbin] = min (abs (tgtrng-R));surfrcs = pow2db (rcs0Hz_highRes (tgtrngbin));tgtscr = tgtrcs - surfrcs

helperPlotSurfTgtRCS (R, rcs0Hz rcs0Hz_highRes、tgtrng tgtrcs)

频率= 2.4 e9;%赫兹rngres = 80;% m脉冲重复频率= 6 e3;%赫兹numdop = 128;

dopres =脉冲重复频率/ numdop;

alt = 4 e3;% m速度= 100;% m / s潜水= 10;%度

startrange = alt;numrange = 120;R = startrange + (0: numrange-1) * rngres;

(d1, d2) = helperSurfaceDopplerLimits(频率、alt、速度、潜水、R(结束));d2 - d1 <脉冲重复频率

放牧= grazingang (alt, R,“平”);nrc =反射(放牧、频率);

(反射)反射。斑纹=“威布尔”;反射。SpeckleShape = 0.8;spck = 0 (numrange numdop);为印第安纳州= 1:numdop [~, spck(:,印第安纳州)]=反射(放牧、频率);结束

[rcs,夹住]= clutterSurfaceRangeDopplerRCS (nrc, R,频率,dopres, alt,速度,潜水,“NumDopplerBins”,numdop);rcs = rcs。* spck;rcs = rcs / 2;

azLook = 90;%度,相对速度矢量elLook = -30;%度

λ= freq2wavelen(频率);[残雪,cy] = helperRangeDopplerCellCenters (R,夹住,λ,alt、速度、潜水);

cy helperPlotCellCentroids (cx)

sensorFrame = rotz (azLook) * roty (-elLook);%的视觉重心在场景框架losx =残雪;losy = cy;losz = alt;%转换为传感器坐标系[losx, losy losz] =交易(sensorFrame (1,1) * losx + sensorFrame (2, 1) * losy + sensorFrame (3,1) * losz,…sensorFrame (1、2) * losx + sensorFrame (2, 2) * losy + sensorFrame (3 2) * losz,…sensorFrame (1、3) * losx + sensorFrame (2、3) * losy + sensorFrame (3,3) * losz);%得到阿兹/ el角度重心(阿兹,el, ~) = cart2sph (losx、losy losz);阿兹= rad2deg (az);el = rad2deg (el);

数组=分阶段。(14(14),λ/ 2);

antresp = phased.ArrayResponse (“SensorArray”、数组);nonnan = ~ isnan (az) & ~ isnan (el);antGain = 0 (numrange numdop);antGain (nonnan) = antresp(频率(阿兹(nonnan) el (nonnan)]。');antGain = mag2db (abs (antGain));

τ= range2time (rngres);Ts = 290;noisePwr = 10 * log10 (physconst (玻耳兹曼的)* Ts /τ);

txPwr = 10 e3;

clutPwr =无穷(numrange numdop);为rIdx = 1: numrange为dIdx = 1: numdop如果rcs (rIdx dIdx) > 0 clutPwr (rIdx dIdx) = radareqsnr(λ,R (rIdx) txPwr,τ,RCS的rcs (rIdx dIdx),“获得”antGain (rIdx dIdx),“t”Ts) + noisePwr;结束结束结束

显示亮度图像(计划/ 1 e3, R / 1 e3, clutPwr)组(gca),“ydir”,“正常”)包含(“多普勒(千赫)”)ylabel (的距离(公里))标题(“表面反射(瓦分贝)的力量”)colorbar

helperPlotCellCentroids (cx, cy clutPwr)

clutMag = db2mag (clutPwr);

PH = helperFormSlowTimeSignal (clutMag,脉冲重复频率,R);

numchan = getNumElements(数组);noisePwrSum = noisePwr + 10 * log10 (numchan);倪= randn(大小(PH));nq = randn(大小(PH));酸碱= + db2mag (noisePwrSum) *我* nq (ni + 1) /√(2);

RDM = fftshift (fft (PH值,[],2),2);

显示亮度图像(计划/ 1 e3, R / 1 e3, 20 * log10 (abs (RDM)))集(gca),“ydir”,“正常”)包含(“多普勒(千赫)”)ylabel (的距离(公里))标题(“预测RDM(瓦分贝))colorbar爬([noisePwrSum-20 + 20 * log10 (numdop) -10])

% helperCompareToScenario(数组、numdop脉冲重复频率,反射,rngres,频率,azLook, elLook, txPwr, noisePwr, alt,速度,潜水,R,夹住)

函数helperSurfaceDopplerLimitsChart(频率、alt)%创建一个图来可视化平面的多普勒极限。社民党= linspace (90150、80);maxRange = linspace (3 e3, 10 e3, 80);潜水= 0;d1_0 = 0(元素个数(maxRange),元素个数(spd));d2_0 = d1_0;为spdIdx = 1:元素个数(spd)为rIdx = 1:元素个数(maxRange) [d1_0 (rIdx spdIdx) d2_0 (rIdx spdIdx)] = helperSurfaceDopplerLimits(频率、alt、社民党(spdIdx),潜水,maxRange (rIdx));结束结束潜水= 30;d1_30 = 0(元素个数(maxRange),元素个数(spd));d2_30 = d1_30;为spdIdx = 1:元素个数(spd)为rIdx = 1:元素个数(maxRange) [d1_30 (rIdx spdIdx) d2_30 (rIdx spdIdx)] = helperSurfaceDopplerLimits(频率、alt、社民党(spdIdx),潜水,maxRange (rIdx));结束结束critR = alt /信德(潜水);h =图;次要情节(1、2、1)轮廓(社民党,maxRange / 1 e3, d2_0/1e3 (2 4 5 6 7 8 9),“黑”,“showtext”,“上”)标题({“最大的多普勒频率(赫兹)”;“潜水= 0度”})包含(“雷达速度(米/秒)”)ylabel (的最大距离(公里)次要情节(1、2、2)轮廓(社民党,maxRange / 1 e3, d2_30/1e3 (2 4 5 6 7 8 9),“黑”,“showtext”,“上”)标题({“最大的多普勒频率(赫兹)”;“潜水= 30度”})包含(“雷达速度(米/秒)”)ylabel (的最大距离(公里))线(xlim [critR critR] / 1 e3,“颜色”,“红色”,“线型”,“——”)h。位置= h。位置+ [0 0 560 0];集(gca),“XLimitMethod”,“紧”)结束

函数helperPlotSurfaceRCS (R,夹住,rcs, R1, R2, d1, d2)%的表面情节RCS range-Doppler空间。可选地,放大%指定时间间隔。h =图;如果输入参数个数> 3次要情节(1、2、1)结束显示亮度图像(计划/ 1 e3, R / 1 e3, 10 * log10 (rcs))集(gca,“ydir”,“正常”)包含(“多普勒(千赫)”)ylabel (的距离(公里))标题(“表面RCS (dBsm)”)colorbar如果输入参数个数> 3补丁(d1 d2 d2 d1 / 1 e3, (R1 R1, R2 R2) / 1 e3,”,“FaceColor”,“没有”,“EdgeColor”,“红色”次要情节(1、2、2)RIdx1 =找到(R < = R1, 1,“最后一次”);RIdx2 =找到(R > = R2, 1,“第一”);DIdx1 =找到(夹住< = d1, 1,“最后一次”);DIdx2 =找到(计划> = d2, 1,“第一”);显示亮度图像(夹住(DIdx1: DIdx2) / 1 e3, e3 (RIdx1: RIdx2) / 1, 10 * log10 (rcs (RIdx1: RIdx2 DIdx1: DIdx2)))集(gca,“ydir”,“正常”)包含(“多普勒(千赫)”)ylabel (的距离(公里))标题({“表面RCS (dBsm)”;“放大”})colorbar h。位置= h。位置+ [0 0 560 0];结束结束

函数helperBeamDopplerLimitsChart(频率、alt、社民党、bw dep)%创建一个图来可视化多普勒波束的极限。阿兹= linspace (-90、90、1 e3);潜水= 0;d1_0 = 0(1,元素个数(az));d2_0 = d1_0;为印第安纳州= 1:元素个数(az) [~, ~, d1_0(印第安纳州),d2_0(印第安纳州)]= helperConicalBeamLimits(频率、alt、社民党、潜水bw,阿兹(印第安纳州),部);结束潜水= 30;d1_30 = 0(1,元素个数(az));d2_30 = d1_30;为印第安纳州= 1:元素个数(az) [~, ~, d1_30(印第安纳州),d2_30(印第安纳州)]= helperConicalBeamLimits(频率、alt、社民党、潜水bw,阿兹(印第安纳州),部);结束h =图;次要情节(1、2、1)情节(az, d1_0/1e3 az, d2_0/1e3)标题({“梁多普勒极限”;“潜水= 0度”})包含(“方位相对侧向(度)”)ylabel (“多普勒极限(千赫)”)传说(分钟多普勒的,“马克斯·多普勒”网格)在集(gca),“XLimitMethod”,“紧”次要情节(1、2、2)情节(az, d1_30/1e3 az, d2_30/1e3)标题({“梁多普勒极限”;“潜水= 30度”})包含(“方位相对侧向(度)”)ylabel (“多普勒极限(千赫)”网格)在传奇(分钟多普勒的,“马克斯·多普勒”甘氨胆酸)组(,“XLimitMethod”,“紧”)h。位置= h。位置+ [0 0 560 0];结束

函数(R1、R2 d1, d2) = helperConicalBeamLimits(频率、alt、社民党、潜水bw, az,放牧)%返回锥形光束的距离和多普勒极限足迹。如果吃草< = bw / 2 R1 =正;R2 =正;elseif吃草< = bw / 2 R1 = alt /信德(放牧+ bw / 2);R2 =正;elseif吃草< = 90 - bw / 2 R1 = alt /信德(放牧+ bw / 2);R2 = alt /信德(graze-bw / 2);其他的R1 = alt;R2 = alt /信德(graze-bw / 2);结束λ= freq2wavelen(频率);azVel = az + 90;%利用方位相对速度矢量斜视= acosd (cosd(潜水)* cosd (azVel) * cosd(放牧)+信德(潜水)*信德(放牧));如果斜视< = bw t1 / 2 =斜视+ bw / 2;t2 = 0;其他的t1 =斜视+ bw / 2;t2 =斜视- bw / 2;结束d1 = 2 /λ*社民党* cosd (t1);d2 = 2 /λ*社民党* cosd (t2);结束

函数helperDopplerWrappingExample%情节展示多普勒包装alt = 1 e3;频率= 2 e9;脉冲重复频率= 1.4 e3;ndop = 256;dopres =脉冲重复频率/ ndop;潜水= 70;速度= 120;R = 2 e3:40:12e3;nrc = 1(大小(R));[RCS1, dopbins1] = clutterSurfaceRangeDopplerRCS (nrc, R,频率,dopres, alt、速度、潜水); [RCS2,dopbins2] = clutterSurfaceRangeDopplerRCS(nrcs,R,freq,dopres,alt,speed,dive,“NumDopplerBins”,ndop);h =图;次要情节(1、2、1)显示亮度图像(dopbins1 R / 1 e3 10 * log10 (RCS1))集(gca,“ydir”,“正常”)标题(打开多普勒频谱的)包含(“多普勒(Hz)”)ylabel (的距离(公里))ch = colorbar;ch.Title.String =“RCS (dBsm)”;次要情节(1、2、2)显示亮度图像(dopbins2 R / 1 e3 10 * log10 (RCS2))集(gca,“ydir”,“正常”)标题(包装的多普勒频谱的)包含(“多普勒(Hz)”)ylabel (的距离(公里))ch = colorbar;ch.Title.String =“RCS (dBsm)”;h。位置= h。位置+ [0 0 560 0];结束

函数helperPlotSurfTgtRCS (R, rcsLow rcsHigh、tgtrng tgtrcs)%情节RCS对范围和注释和RCS目标范围。情节(R / 1 e3, 10 * log10 (rcsLow))在情节(R / 1 e3, 10 * log10 (rcsHigh))情节(tgtrcs tgtrng / 1 e3,*黑的)举行从包含(的距离(公里))ylabel (“RCS (dBsm)”)传说(“表面RCS(低分辨率)”,的表面RCS(高清),“目标”,“位置”,“最佳”)标题(表面和目标RCS的)结束

函数(d1, d2) = helperSurfaceDopplerLimits(频率、alt、社民党、潜水maxRange)%返回一个平面上的多普勒极限。%俯角最大范围如果maxRange < alt maxRangeDep = 90;其他的maxRangeDep = 90 - acosd (alt / maxRange);结束%规范化最大关闭率如果潜水< maxRangeDep maxClosing = cosd (maxRangeDep-dive);其他的maxClosing = 1;结束%的最低收盘汇率如果那里潜水< maxRangeDep minClosing = -cosd (maxRangeDep +潜水);其他的minClosing = 1;结束%的因素多普勒速度和转换λ= freq2wavelen(频率);d1 = 2 /λ*社民党* minClosing;d2 = 2 /λ*社民党* maxClosing;结束

函数(x, y) = helperRangeDopplerCellCenters (R,夹住,λ,alt、速度、潜水)%的x / y坐标返回的一组range-Doppler细胞中心%在表面上R = R (:);夹住=计划(:)。”;r =√r ^ 2-alt ^ 2);ctheta =(λ* r *计划/(2 *速度)- alt *信德(潜水)。/ (r * cosd(潜水));:ctheta (r = = 0) = 0;ctheta (ctheta > 1) =南;ctheta (ctheta < 1) =南;θ= acosd (ctheta);x = r。* cosd(θ);y = r。*信德(θ);结束

函数helperPlotCellCentroids (cy残雪,cdata)%的情节单元重心在笛卡尔空间。选择指定的颜色数据。如果输入参数个数< 3情节(cx (:) / 1 e3, cy (:) / 1 e3,“布莱克”,“MarkerSize”,1)其他的散射(cx (:) / 1 e3, cy (:) / 1 e3, 10, cdata (:),“填充”)colorbar结束标题(“Range-Doppler单元形心”)包含(“X(公里))ylabel (“Y”(公里))轴平等的集(gca),“XLimitMethod”,“紧”甘氨胆酸)组(,“YLimitMethod”,“紧”)结束

函数PH = helperFormSlowTimeSignal (specmag,脉冲重复频率,R)%形式的时域信号幅度谱[numr, numd] =大小(specmag);%平移谱级specmag = ifftshift (specmag 2);%初始化时域信号PH值= 0(大小(specmag));%样品时间tslow = 0:1 /脉冲重复频率(numd-1): /脉冲重复频率;tfast = 2 * (R (:) - R (1)) / 299792458;t = tslow + tfast;%的随机相位偏移量为每一个垃圾箱φ=兰德(numr numd) 1/2;%多普勒轴df =脉冲重复频率/ numd;夹住= 0:df (prf-df);在每本%随机多普勒偏移量fc =计划+(兰德(numr numd) 1/2) * df;为印第安纳州= 1:numd%的整个信号在当前的多普勒频率y = specmag(:,印第安纳州)。* exp(1 * 2 *π* (fc(:,印第安纳州)。* t +φ(:,印第安纳州)));%积累酸碱= + y;结束结束

函数helperCompareToScenario(数组、numdop脉冲反射,rngres, fc, azLook, elLook, txPwr, noisePwrDb, alt,速度,潜水,R,夹住)%模拟杂波返回radarScenario和阴谋的结果%使场景场景= radarScenario (“UpdateRate”0,“StopTime”(numdop-1/2) /脉冲重复频率);landSurface(场景中,“RadarReflectivity”,反射);%使雷达c = physconst (“光速”);rdr = radarTransceiver (“MountingAngles”(azLook -elLook 0),“NumRepetitions”,numdop);rdr.TransmitAntenna。传感器=数组;rdr.ReceiveAntenna。传感器=数组;rdr.TransmitAntenna。OperatingFrequency = fc;rdr.ReceiveAntenna。OperatingFrequency = fc;rdr.Waveform。脉冲重复频率=脉冲重复频率;sampleRate = c / (2 * rngres);sampleRate =脉冲重复频率*轮(sampleRate /脉冲重复频率);%调整与脉冲重复频率匹配约束rdr.Receiver。SampleRate = SampleRate;rdr.Waveform。SampleRate = SampleRate;rdr.Waveform。脉冲宽度= 2 * rngres / c;rdr.Transmitter。获得= 0;rdr.Transmitter。LossFactor = 0;rdr.Transmitter。PeakPower = txPwr * sqrt (2); rdr.Receiver.Gain = 0; rdr.Receiver.LossFactor = 0; rdr.Receiver.NoiseMethod =“噪音力量”;rdr.Receiver。NoisePower = 10 ^ (noisePwrDb / 10);%增加雷达平台平台(场景中,“传感器”rdr,“轨迹”kinematicTrajectory (“位置”[0 0 alt),“速度”、速度* [cosd(潜水)0信德(潜水))));%配置杂乱的一代clut = clutterGenerator (rdr的场景,“UseBeam”假的,“决议”rngres / 2,“RangeLimit”R(结束));ringClutterRegion (clut 0 sqrt (R(结束)^ 2-alt ^ 2), 180年,azLook);%收集返回信号和总和梁形式iqsig =接收(场景);PH =排列(总和(iqsig {1}, 2), [1 3 2]);%形式RDMRDM = fftshift (fft (PH值,[],2),2);%截断范围uar = c /(2 *脉冲重复频率);拉尔= 0:rngres: (uar-rngres);RDM = RDM(拉尔> = R(1) &拉尔< = R(结束):);%绘制结果图显示亮度图像(计划/ 1 e3, R / 1 e3, 20 * log10 (abs (RDM)))集(gca),“ydir”,“正常”)包含(“多普勒(千赫)”)ylabel (的距离(公里))标题(“模拟RDM(瓦分贝))colorbar numchan = getNumElements(数组);noisePwrSum = noisePwrDb + 10 * log10 (numchan);爬([noisePwrSum-20 + 20 * -10 (log10 (numdop)))结束