畢業(yè)論文_基于聲音能量的無線傳感網(wǎng)絡定位算法-資料下載頁

【導讀】無線傳感器網(wǎng)絡在當今社會起著非常重要的作用。而傳感器節(jié)點的位置對整個。則所采集的數(shù)據(jù)便毫無意義。因此定位是無線傳感器網(wǎng)絡必備的功能，也是無線傳。感器網(wǎng)絡的核心技術(shù)。律和無線傳感器網(wǎng)絡的幾種算法，距離相關(guān)算法。設計主要針對三邊定位算法和極大似。法的兩種算法中極大似然估計算法更好一點的結(jié)論。

　　

【正文】 :10001 sz(i)=U1。 end for j=1:10000 sz(j+10000)=U2。 end figure(2) plot(f,10*log10(Pxx))。 window=kaiser(length(xn),)。 nfft=1024。 [Pxx,f]=periodogram(xn,window,nfft,fs)。 figure(3) plot(f,10*log10(Pxx))。 window1=boxcar(length(xn))。 window2=hamming(length(xn))。 window3=blackman(length(xn))。 [Pxx1,f]=periodogram(xn,window1,nfft,fs)。 [Pxx2,f]=periodogram(xn,window2,nfft,fs)。 [Pxx3,f]=periodogram(xn,window3,nfft,fs)。 figure(4) subplot(221) plot(f,10*log10(Pxx))。 subplot(222) plot(f,10*log10(Pxx1))。 subplot(223) 內(nèi)蒙古工業(yè)大學本科畢業(yè)設計說明書 27 plot(f,10*log10(Pxx2))。 subplot(224) plot(f,10*log10(Pxx3))。 兩定位算法程序： close all clc clear disp(39。199。235。202。228。200。235。210。209。214。170。189。218。181。227。202。253。n(n3)。39。)。 n=input(39。input n=39。)。 if n3。 disp(39。202。228。200。235。189。218。181。227。202。253。180。237。39。)。 return。 end for i=1:n x(i)=10*normrnd(0,1,1)。 % y(i)=10*normrnd(0,1,1)。 scatter(x(i),y(i),39。r39。)。 hold on end grid。 x(n+1)=10*normrnd(0,1,1)。 y(n+1)=10*normrnd(0,1,1)。 scatter(x(n+1),y(n+1),39。r39。)。 P=[x(n+1),y(n+1)] for i=1:n。 d(i)=sqrt((x(n+1)x(i))^2+(y(n+1)y(i))^2)。 end d(1,i)=d(i) for i=1:n。 內(nèi)蒙古工業(yè)大學本科畢業(yè)設計說明書 28 D(i)=d(i)+2*normrnd(0,1,1)。 end D(1,i)=D(i) for i=1:n1 M(i,1)=[2*(x(i)x(n))]。 M(i,2)=[2*(y(i)y(n))]。 end for i=1:n1。 N(i,1)=[x(i)^2x(n)^2+y(i)^2y(n)^2+D(n)^2D(i)^2]。 end T=M39。 X=inv(T*M)*T*N。 xp=X([1]) yp=X([2]) scatter(xp,yp,39。m39。)。 error=sqrt((xpx(n+1))^2+(ypy(n+1))^2) for i=1:3 x(i)=10*normrnd(0,1,1)。 % y(i)=10*normrnd(0,1,1)。 hold on scatter(x(i),y(i),39。b39。)。 end grid。 A=[x(1),y(1)] B=[x(2),y(2)] C=[x(3),y(3)] k = 1。 while (k==1) x(4)=10*normrnd(0,1,1)。 y(4)=10*normrnd(0,1,1)。 內(nèi)蒙古工業(yè)大學本科畢業(yè)設計說明書 29 a4=(x(4)x(3))/(x(2)x(3))(y(4)y(3))/(y(2)y(3))。 a1=(x(1)x(3))/(x(2)x(3))(y(1)y(3))/(y(2)y(3))。 b4=(x(4)x(1))/(x(3)x(1))(y(4)y(1))/(y(3)y(1))。 b2=(x(2)x(1))/(x(3)x(1))(y(2)y(1))/(y(3)y(1))。 c4=(x(4)x(1))/(x(2)x(1))(y(4)y(1))/(y(2)y(1))。 c3=(x(3)x(1))/(x(2)x(1))(y(3)y(1))/(y(2)y(1))。 if (a4*a10)amp。amp。(b4*b20)amp。amp。(c4*c30)。 k=0。 hold on scatter(x(4),y(4),39。b39。)。 P=[x(4),y(4)] end end for i=1:3。 d(i)=sqrt((x(4)x(i))^2+(y(4)y(i))^2)。 end d=[d(1),d(2),d(3)] for i=1:3。 D(i)=d(i)+2*normrnd(0,1,1)。 end D=[D(1),D(2),D(3)] M=[2*(x(1)x(3)),2*(y(1)y(3))。2*(x(2)x(3)),2*(y(2)y(3))]。 N=[x(1)^2x(3)^2+y(1)^2y(3)^2+D(3)^2D(1)^2。x(2)^2x(3)^2+y(2)^2y(3)^2+D(3)^2D(2)^2]。 T=M39。 X=inv(T*M)*T*N。 x4=X([1]) y4=X([2]) scatter(x4,y4,39。g39。) error=sqrt((x4x(4))^2+(y4y(4))^2) 算法誤差比較程序： 內(nèi)蒙古工業(yè)大學本科畢業(yè)設計說明書 30 close all clc clear disp(39。199。235。202。228。200。235。210。209。214。170。189。218。181。227。202。253。n(n3)。39。)。 n=input(39。input n=39。)。 if n3。 disp(39。202。228。200。235。189。218。181。227。202。253。180。237。39。)。 return。 end for m=1:100。 for i=1:n x(i)=10*normrnd(0,1,1,1)。 % y(i)=10*normrnd(0,1,1,1)。 end grid。 x(n+1)=10*normrnd(0,1,1,1)。 y(n+1)=10*normrnd(0,1,1,1)。 P=[x(n+1),y(n+1)]。 for i=1:n。 d(i)=sqrt((x(n+1)x(i))^2+(y(n+1)y(i))^2)。 end d(1,i)=d(i)。 for i=1:n。 D(i)=d(i)+2*normrnd(0,1,1)。 end D(1,i)=D(i)。 for i=1:n1 M(i,1)=[2*(x(i)x(n))]。 M(i,2)=[2*(y(i)y(n))]。 end 內(nèi)蒙古工業(yè)大學本科畢業(yè)設計說明書 31 for i=1:n1。 N(i,1)=[x(i)^2x(n)^2+y(i)^2y(n)^2+D(n)^2D(i)^2]。 end T=M39。 X=inv(T*M)*T*N。 xp=X([1])。 yp=X([2])。 error=sqrt((xpx(n+1))^2+(ypy(n+1))^2)。 errorList(m)=error。 end graNum=1:100。 plot(graNum,errorList)。 hold on for m = 1:100 for i=1:3 x(i)=10*normrnd(0,1,1,1)。 % y(i)=10*normrnd(0,1,1,1)。 end grid。 A=[x(1),y(1)]。 B=[x(2),y(2)]。 C=[x(3),y(3)]。 k = 1。 while (k==1) x(4)=10*normrnd(0,1,1,1)。 y(4)=10*normrnd(0,1,1,1)。 a4=(x(4)x(3))/(x(2)x(3))(y(4)y(3))/(y(2)y(3))。 a1=(x(1)x(3))/(x(2)x(3))(y(1)y(3))/(y(2)y(3))。 b4=(x(4)x(1))/(x(3)x(1))(y(4)y(1))/(y(3)y(1))。 b2=(x(2)x(1))/(x(3)x(1))(y(2)y(1))/(y(3)y(1))。 內(nèi)蒙古工業(yè)大學本科畢業(yè)設計說明書 32 c4=(x(4)x(1))/(x(2)x(1))(y(4)y(1))/(y(2)y(1))。 c3=(x(3)x(1))/(x(2)x(1))(y(3)y(1))/(y(2)y(1))。 if (a4*a10)amp。amp。(b4*b20)amp。amp。(c4*c30)。 k=0。 P=[x(4),y(4)]。 end end for i=1:3。 d(i)=sqrt((x(4)x(i))^2+(y(4)y(i))^2)。 end d=[d(1),d(2),d(3)]。 for i=1:3。 D(i)=d(i)+2*normrnd(0,1,1)。 end D=[D(1),D(2),D(3)]。 M=[2*(x(1)x(3)),2*(y(1)y(3))。2*(x(2)x(3)),2*(y(2)y(3))]。 N=[x(1)^2x(3)^2+y(1)^2y(3)^2+D(3)^2D(1)^2。x(2)^2x(3)^2+y(2)^2y(3)^2+D(3)^2D(2)^2]。 T=M39。 X=inv(T*M)*T*N。 x4=X([1])。 y4=X([2])。 error=sqrt((x4x(4))^2+(y4y(4))^2)。 errorList(m) = error。 end graNum = 1:100。 plot(graNum, errorList,39。r39。)