%自由空間的三角定位法
%設(shè)發(fā)射天線坐標(biāo)T[00]???接受天線坐標(biāo)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）;

%設(shè)R1收到3e-8（s）收到回波信號(hào)???R2經(jīng)過4.5e-8（s）收到信號(hào)
td1?=?3e-8;
td2?=?4.5e-8;

C?=?3e8;%光速

%傳播距離
X1?=?C*td1;
X2?=?C*td2;

%以T，R1為焦點(diǎn)畫橢圓1
%以T，R2為焦點(diǎn)畫橢圓2

%求解橢圓參數(shù)
%求a

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

%橢圓參數(shù)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]）;

%求解目標(biāo)點(diǎn)
[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）;

%設(shè)置墻高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;%墻體高度
%墻體厚度和墻體的相對(duì)介電常數(shù)較小。研究電磁波和墻體的相互作用時(shí)，忽略了電磁波的折射
%設(shè)相對(duì)介電常數(shù)Er?=?4;
Er?=?4;?%相對(duì)介電常數(shù)
Vd?=?C/（sqrt（Er））;??%電磁波在墻體中的傳播速度

td?=?d/Vd;%電磁波在墻體中的傳播時(shí)間
%設(shè)R1收到3e-8（s）收到回波信號(hào)???R2經(jīng)過4.5e-8（s）收到信號(hào)
td1?=?3e-8;
td2?=?4.5e-8;

%在自由空間內(nèi)的真實(shí)時(shí)延為
td11?=?td1?-?2?*?td;
td22?=?td2?-?2?*?td;

%建立橢圓方程，利用橢圓性質(zhì)求橢圓方程參數(shù)a
aq1?=?C*（td11）/2;
aq2?=?C*（td22）/2;

%焦點(diǎn)距離24.得出橢圓參數(shù)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]）;

%求解目標(biāo)點(diǎn)
[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‘）

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

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

-----------?---------??----------?-----??----

?????????????????3493????????????????????1