clc;
clear;
close all;

N=4; ? ? ? ?%參與定位的基站數?
C=3e5; ? %電磁波傳播速度300000000m/s
X=[0 5000 5000 0];
Y=[0 0 5000 5000];
x=1200; ? ?y=1600;
D(1:N)=sqrt((X(1:N)-x).^2+(Y(1:N)-y).^2);
delay(1:N)=D(1:N)./C;
T=1;
N_sample=200000;%傳輸速率200kbit/s
iTsignal=[ones(1,N_sample/2),-ones(1,N_sample/2)];%發送信號的實部
qTsignal=zeros(size(iTsignal));%發送信號的虛部
rsn=[0 3 5 10 20 30 40 50]; ?%信噪比
%rsn=40;
standard=0; ? ? ? %1為b標準，0為a標準；
for j=1:8
? ? for i=5
? ? [iRsignal1]=pathoutput(iTsignal,rsn(j),standard,delay(1));
? ? [iRsignal2]=pathoutput(iTsignal,rsn(j),standard,delay(2));
? ? [iRsignal3]=pathoutput(iTsignal,rsn(j),standard,delay(3));
? ? [iRsignal4]=pathoutput(iTsignal,rsn(j),standard,delay(4));
? ?
? ? iRsignal1(i,:)=iRsignal1;
? ? iRsignal2(i,:)=iRsignal2;
? ? iRsignal3(i,:)=iRsignal3;
? ? iRsignal4(i,:)=iRsignal4;
? ??
? ? end
? ? iRsignal1=mean(iRsignal1);
? ? iRsignal2=mean(iRsignal2);
? ? iRsignal3=mean(iRsignal3);
? ? iRsignal4=mean(iRsignal4);
? ?
? ?TDOA(1)=correlation(iRsignal2,iRsignal1,T,N_sample);
? ?TDOA(2)=correlation(iRsignal3,iRsignal1,T,N_sample);
? ?TDOA(3)=correlation(iRsignal4,iRsignal1,T,N_sample);
? ?
? ?Q=eye(N-1);
? ?for i=1:N-1
? ? Rplus(i)=TDOA(i)*C; %第一個reader為主站其余為輔站，Ri,1為第i個站到定位點的距離與第一個站到其距離的差。
? ? R(i)=(TDOA(i)*C).^2; %為了在Tdoainiial函數中使用R
? ?end ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? %Q=sigma(l).*[1 0.5 0.5;0.5 1 0.5;0.5 0.5 1]; %噪聲協方差陣隨著基站的數目變化而改變。
? ?[xu yu]=Leastsquare(X,Y,N,Q,R,Rplus,x,y);
? ? xr(j)=xu-x;
? ? yr(j)=yu-y;
? ? locationerror(j)=sqrt((xu-x).^2+(yu-y).^2);
end
figure(1);
plot(locationerror,'b*-');
xlabel('信噪比（dB）（1-8分別代表信噪比[0 3 5 10 20 30 40 50]');
ylabel('誤差');
title('信噪比對誤差的影響');
grid on
figure(2);
plot(xr,'k*-')
hold on
plot(yr,'rx--')
xlabel('信噪比（dB）（1-8分別代表信噪比[0 3 5 10 20 30 40 50]');
ylabel('x、y坐標誤差');
gtext({'xr'});
gtext({'yr'});
title('信噪比對x、y坐標誤差的影響');
grid on
?

總結

以上是生活随笔為你收集整理的基于最小二乘法的无线定位技术研究与仿真_程序的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。