順着昨天的思路,明顯可以看到速率偏低,於是運行代碼,發現第四項的結果非常大,仔細一看發現距離產生的有問題,修改距離的產生
d_out_ture=sqrt(2^2+d_out.^2-2*2*d_out.*cos(theta_out));
forth_2=sum(pu*z*d_out_ture(1,:).^r)-pu*z*d_out_ture(1,j)^r;
forth_3=sum(pu*z*d_out_ture(2,:).^r)-pu*z*d_out_ture(2,j)^r;
forth=forth_1+forth_2+forth_3;
發現結果能稍微好一點。
但是還是無法確定速率的高低,所以無法確定寫法的對錯。基於這種情況,將這種算法應用於時間同步的情況中,和paper中的結果進行比較,進而確定算法的正確性,再回過頭來應用於空間劃分的情況。