确信罪犯藏身地址问题Word文件下载.docx
《确信罪犯藏身地址问题Word文件下载.docx》由会员分享,可在线阅读,更多相关《确信罪犯藏身地址问题Word文件下载.docx(28页珍藏版)》请在冰豆网上搜索。
5
6
25'
2°
15'
7
47'
46'
8
39'
9
27'
13'
10
42'
52'
11
12
40'
13
34'
(1)将角度制的坐标转换为弧度制坐标;
(2)将系列犯法的地址(13个)的球面坐标转换为平面坐标;
并画出这些点的图形(要求带适当宽度的网格,行宽*10^3,列宽h=*10^3)。
(3)求这13个点之间距离的最大值和最小值;
别离求这13个点相邻点横坐标和纵坐标之差绝对值的最大值和最小值;
4.对所得13个点的每一个点,在其周围找16个点(可疑点),其中格宽为d=500。
将这16个点与其余的12个犯法的地址进行比较,假设与某犯法的地址距离小于等于d,那么去掉该可疑点。
求出所有剩余的点(组成可疑点集合)。
(5)关于d=500:
500:
60000,重复上述进程,关于不同的d,求可疑点的个数,画出d与可疑点个数的图形,并确信d为何值使,可疑点个数最小。
(6)对上述所得的d,类似问题4,求可疑点YD。
(7)关于YD中每一个点(如第j个点),计算该点到问题2中的13个点
的距离
,成立该点的衰减函数:
,求
的最大值及YD中对应的点——罪犯所在的可能性最大的点,并将该点加在问题2的图形上。
(8)关于d=(21:
30)*500,重复上述进程,关于不同的d,求罪犯所在的可能性最大的点,并将这些点加在问题2的图形上。
四
编辑m文件disiti1
functionradian=disiti1(x)
d=size(x,1);
fori=1:
d;
a=deg2rad(x(i,1));
b=deg2rad(x(i,2)/60);
c=deg2rad(x(i,3)/3600);
radian(i,1)=a+b+c;
end
输入纬度:
weidu=[53,49,;
53,48,;
53,50,;
53,49,;
53,25,;
53,47,;
53,39,;
53,27,;
53,42,;
53,49,]
结果:
disiqi1(weidu)
输入经度为:
jingdu=[1,32,;
1,31,;
1,29,;
1,45,;
2,15,;
1,46,;
2,13,;
1,52,;
1,40,;
1,34,]
disiti1(jingdu)
以地心为原点,赤道平面为投影平面,其中a是坐标点的纬度,b是坐标点的经度。
地球半径R为6380000米。
x=Rcos(a)cos(b)
y=Rcos(a)sin(b)
新建M文件:
functionplane=disiti2(a,b)
c=size(a,1);
r=6400000;
c
d(i,1)=r*cos(a(i))*cos(b(i));
d(i,2)=r*cos(a(i))*sin(b(i));
x1=(min(d(:
1))-1250):
1250:
(max(d(:
1))+1250);
y1=(min(d(:
2))-2500):
2500:
2))+2500);
[X1,Y1]=meshgrid(x1,y1);
holdon;
plot(X1,Y1,'
b'
);
plot(X1'
Y1'
'
plot(d(:
1),d(:
2),'
o'
)
holdoff;
vpa(d)
输入:
a=[]
b=[]
运行m文件函数:
disiti2(a,b)
结果为:
ans=
[,]
[,]
4.3.1
a=
[[,]
[,]]
z=zeros(13,13);
fork=1:
13
forl=1:
z(k,l)=sqrt((a(k,1)-a(l,1))*(a(k,1)-a(l,1))+(a(k,2)-a(l,2))*(a(k,2)-a(l,2)))
d=z+eye(13,13)*10000
d=+004*
m=max(max(d))
n=min(min(d))
m=+004
n=
4.3.2
k=min(d)
h=zeros(13,2)
fori=1:
u=find(d(:
i)==k(i));
h(i,1)=a(i,1)-a(u,1);
h(i,2)=a(i,2)-a(u,2);
h;
v=abs(h)
m=max(v)
n=min(v)
v=
+003*
m=+003*
n=
t=[];
d=500;
e=[(a(i,1)+d),a(i,2);
(a(i,1)+2*d),a(i,2);
(a(i,1)-d),a(i,2);
(a(i,1)-2*d),a(i,2);
a(i,1),(a(i,2)+d);
a(i,1),(a(i,2)+2*d);
a(i,1),(a(i,2)-d);
a(i,1),(a(i,2)-2*d);
(a(i,1)+d),(a(i,2)+d);
(a(i,1)+2*d),(a(i,2)+2*d);
(a(i,1)-d),(a(i,2)-d);
(a(i,1)-2*d),(a(i,2)-2*d);
(a(i,1)+d),(a(i,2)-d)
(a(i,1)+2*d),(a(i,2)-2*d)
(a(i,1)-d),(a(i,2)+d);
(a(i,1)-2*d),(a(i,2)+2*d)];
fork=1:
16
forj=1:
ifj~=i
w(j,1)=sqrt((e(k,1)-a(j,1))*(e(k,1)-a(j,1))+(e(k,2)-a(j,2))*(e(k,2)-a(j,2)));
ifmin(w)>
d
h=e(k,:
t=[t;
h];
t
t=
+006*
u=size(t,1)
u=184
l=[]
ford=500:
60000
q=[];
u=size(t,1);
l=[l;
u]
l
d=500:
60000;
plot(d,l)
y=min(l)
y=119
x=find(l==119)
x=35
c=35*500
c=
17500
d=17500
(a(i,1)+d),(a(i,