test.m

来自「Friedlander算法,假定服务基站的二维坐标为(0,0),基站数为4,多基」· M 代码 · 共 31 行

clear
len=100;
ff=0:1/100:1;
yy=sqrt(ff);
len=len+1;
% aaa=[1.1 1 0.2; 0 1 1.2];
aaa=[1 1 0; 0 1 1];
a(1,1)=0;
a(2,1)=0;
a(1,2:4)=aaa(1,:);
a(2,2:4)=aaa(2,:);
plot(a(1,:),a(2,:),'r.','MarkerSize',20)
for n=1:len
    bsDis=sqrt(ff(n)^2+yy(n)^2);
    rel(1,n)=sqrt((ff(n)-aaa(1,1))^2+(yy(n)-aaa(2,1))^2)-bsDis;
    rel(2,n)=sqrt((ff(n)-aaa(1,2))^2+(yy(n)-aaa(2,2))^2)-bsDis;
    rel(3,n)=sqrt((ff(n)-aaa(1,3))^2+(yy(n)-aaa(2,3))^2)-bsDis;
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %加入噪声,注释掉这一部分后为无噪的情况.
% randn( 'state',sum(100*clock) );
% voice=randn(3,len);
% voice=voice./max(max(abs(voice)));
% rel=rel+voice*0.05;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for nn=1:len
    ppp(:,nn)=getPosition2(aaa,rel(:,nn));
end
hold,
plot(ff,yy,'bo',ppp(1,:),ppp(2,:),'r.')
grid on