基于椭圆雷达定位的matlab仿真

本仿真采用椭圆定位原理，对自由空间目标已经墙体后面的目标实现定位。主要可用于穿墙雷达，探地雷达研究中的目标定位，此定位方法简单，计算量小。在少目标情况下，得到定位效果十分明显。

%自由空间的三角定位法
%设发射天线坐标T[00]   接受天线坐标R1[-20]   R2[40]
clear all

plot（00‘r.‘‘markersize‘20）;
hold on
plot（-20‘r.‘‘markersize‘20）;
hold on
plot（40‘r.‘‘markersize‘20）;

%设R1收到3e-8（s）收到回波信号   R2经过4.5e-8（s）收到信号
td1 = 3e-8;
td2 = 4.5e-8;

C = 3e8;%光速

%传播距离
X1 = C*td1;
X2 = C*td2;

%以T，R1为焦点画椭圆1
%以T，R2为焦点画椭圆2

%求解椭圆参数
%求a

a1 = X1/2;
a2 = X2/2;

%椭圆参数C为焦距/2
c1 = 1;
c2 = 2;

b1 = sqrt（a1^2-c1^2）;
b2 = sqrt（a2^2-c2^2）;
hold on
ezplot（‘（x+1）^2/（4.5^2） + （y）^2/（4.387482193696061^2） = 1‘）;
hold on
ezplot（‘（x-2）^2/（6.75^2）+（y）^2/（6.446898479113814^2） = 1‘）;
hold on

line（[-66][00]）;
line（[00][-66]）;

%求解目标点
[xy]=solve（‘（x+1）^2/（4.5^2） + （y）^2/（4.387482193696061^2） = 1（x-2）^2/（6.75^2）+（y）^2/（6.446898479113814^2） = 1‘）;
x = double（x）;
y = double（y）;

hold on 
plot（-3.6250000000000003.563657219455566 + 0.000000000000000i‘r.‘‘markersize‘20）;
hold on
plot（-3.625000000000000-3.563657219455566 + 0.000000000000000i‘r.‘‘markersize‘20）;

%连接路径
A1 = [-3.625000000000000 0];
B1 = [3.563657219455566 + 0.000000000000000i 0];
hold on 
plot（A1B1‘k‘）;

A2 = [-3.625000000000000 4];
B2 = [3.563657219455566 + 0.000000000000000i 0];
hold on
plot（A2B2‘k‘）;

A3 = [-3.625000000000000 -2];
B3 = [3.563657219455566 + 0.000000000000000i 0];
hold on 
plot（A3B3‘k‘）;

title（‘自由空间三角定位‘）;

%///////////////////////////////////////////////////////////////////////////////
%加入墙体后的三角定位法
figure
plot（00‘r.‘‘markersize‘20）;
hold on
plot（-20‘r.‘‘markersize‘20）;
hold on
plot（40‘r.‘‘markersize‘20）;

%设置墙高d = 1;
hold on
plot（01‘b.‘‘markersize‘20）;
hold on
plot（-21‘b.‘‘markersize‘20）;
hold on
plot（41‘b.‘‘markersize‘20）;


d = 1;%墙体高度
%墙体厚度和墙体的相对介电常数较小。研究电磁波和墙体的相互作用时，忽略了电磁波的折射
%设相对介电常数Er = 4;
Er = 4; %相对介电常数
Vd = C/（sqrt（Er））;  %电磁波在墙体中的传播速度

td = d/Vd;%电磁波在墙体中的传播时间
%设R1收到3e-8（s）收到回波信号   R2经过4.5e-8（s）收到信号
td1 = 3e-8;
td2 = 4.5e-8;

%在自由空间内的真实时延为
td11 = td1 - 2 * td;
td22 = td2 - 2 * td;

%建立椭圆方程，利用椭圆性质求椭圆方程参数a
aq1 = C*（td11）/2;
aq2 = C*（td22）/2;

%焦点距离24.得出椭圆参数c
cq1 = 1;
cq2 = 2;

bq1 = sqrt（aq1^2-cq1^2）;
bq2 = sqrt（aq2^2-cq2^2）;

hold on
ezplot（‘（x+1）^2/（2.5^2） + （y-1）^2/（2.291287847477920^2） = 1‘）;
hold on
ezplot（‘（x-2）^2/（4.75^2）+（y-1）^2/（4.308421984903522^2） = 1‘）;

line（[-66][00]）;
line（[00][-66]）;
line（[-66][11]）;
line（[-6-6][01]）;
line（[66][01]）;

%求解目标点
[xqyq]=solve（‘（x+1）^2/（2.5^2） + （y-1）^2/（2.291287847477920^2） = 1（x-2）^2/（4.75^2）+（y-1）^2/（4.308421984903522^2） = 1‘）;
xq = double（xq）;
yq = double（yq）;

plot（-2.2019230769230773.009109248003752 + 0.000000000000000i‘b.‘‘markersize‘20）;
title（‘加入墙体后的三角定位仿真‘）;

%连接路径
A4 = [-2.201923076923077 -2];
B4 = [3.009109248003752 + 0.000000000000000i 1];
hold on 
plot（A4B4‘k‘）;

A5 = [-2 -2];
B5 = [1 0];
hold on 
plot（A5B5‘k‘）;

A6 = [-2.201923076923077 0];
B6 = [3.009109248003752 + 0.000000000000000i 1];
hold on 
plot（A6B6‘k‘）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       3493  2018-02-28 19:02  sanjiaodingwei.m

----------- ---------  ---------- -----  ----

                 3493                    1