您现在的位置:首页 >> 技术文章 >> MATLAB技术 >> 内容

MATLAB代做|FPGA代做--直达波

时间:2018-10-6 15:53:43 点击:

  核心提示:MATLAB代做|FPGA代做--直达波...

%模型直达波距离R=800,ros模型半径r=100,基站坐标为原点(0,0),定位目标真实坐标为(800,0)

for j=1:4   %4条多径
    xs=700+200*rand(1);
    l(j)=sqrt(xs^2+100^2-(800-xs)^2)+100;%nlos的toa
   
end
lmax=l(1);
lmin=l(1);
for j=2:4
  if  l(j)>lmax
      lmax=l(j);%最大toa
  elseif l(j)<lmin
      lmin=l(j);%最小toa
  end
end
R=700;  %初始化R
r=200;
% for i=1:100
%  pr=@(r)(-(2/(3*R-lmax))*(l(1)-r)/(R*r*sqrt(1-((R-l(1))*(R+l(1)-2*r)/(2*R*r)+1)^2))*(l(2)-r)/(R*r*sqrt(1-((R-l(2))*(R+l(2)-2*r)/(2*R*r)+1)^2))*(l(3)-r)/(R*r*sqrt(1-((R-l(3))*(R+l(3)-2*r)/(2*R*r)+1)^2))*(l(4)-r)/(R*r*sqrt(1-((R-l(4))*(R+l(4)-2*r)/(2*R*r)+1)^2)));
%     [r,m]=fminsearch(pr,200);
% pR=@(r)(-(1/(lmin-lmax+2*r))*(l(1)-r)/(R*r*sqrt(1-((R-l(1))*(R+l(1)-2*r)/(2*R*r)+1)^2))*(l(2)-r)/(R*r*sqrt(1-((R-l(2))*(R+l(2)-2*r)/(2*R*r)+1)^2))*(l(3)-r)/(R*r*sqrt(1-((R-l(3))*(R+l(3)-2*r)/(2*R*r)+1)^2))*(l(4)-r)/(R*r*sqrt(1-((R-l(4))*(R+l(4)-2*r)/(2*R*r)+1)^2)));
%      [R,n]=fminsearch(pR,700);
% end
for i=1:100     %em循环100次
%     k1=int(r*(l(1)-r)/(R*r*sqrt(1-((R-l(1))*(R+l(1)-2*r)/(2*R*r)+1)^2)),r,(l(1)-R)/2,R);
%     pr=@(r)(l(1)-r)/(R*r*sqrt(1-((R-l(1))*(R+l(1)-2*r)/(2*R*r)+1)^2))*k1*(2/(3*R-lmax));
    for j=1:4
        k1(j)=int('r*(l(j)-r)/(R*r*sqrt(1-((R-l(j))*(R+l(j)-2*r)/(2*R*r)+1)^2))','r',(l(j)-R)/2,R);
    end
    %pr为r的条件概率密度函数,用以求r
    pr=@(r)(-(2/(3*R-lmax))*(l(1)-r)/(R*r*sqrt(1-((R-l(1))*(R+l(1)-2*r)/(2*R*r)+1)^2))*k1(1)*(l(2)-r)/(R*r*sqrt(1-((R-l(2))*(R+l(2)-2*r)/(2*R*r)+1)^2))*k1(2)*(l(3)-r)/(R*r*sqrt(1-((R-l(3))*(R+l(3)-2*r)/(2*R*r)+1)^2))*k1(3)*(l(4)-r)/(R*r*sqrt(1-((R-l(4))*(R+l(4)-2*r)/(2*R*r)+1)^2))*k1(4));
    [r,m]=fminsearch(pr,200);%求r
    
     for j=1:4
        k2(j)=int('r*(l(j)-r)/(R*r*sqrt(1-((R-l(j))*(R+l(j)-2*r)/(2*R*r)+1)^2))','R',l(j)-2*r,l(j));
     end
     %pR为R的条件概率密度函数,用以求R
     pR=@(r)(-(1/(lmin-lmax+2*r))*(l(1)-r)/(R*r*sqrt(1-((R-l(1))*(R+l(1)-2*r)/(2*R*r)+1)^2))*k2(1)*(l(2)-r)/(R*r*sqrt(1-((R-l(2))*(R+l(2)-2*r)/(2*R*r)+1)^2))*k2(2)*(l(3)-r)/(R*r*sqrt(1-((R-l(3))*(R+l(3)-2*r)/(2*R*r)+1)^2))*k2(3)*(l(4)-r)/(R*r*sqrt(1-((R-l(4))*(R+l(4)-2*r)/(2*R*r)+1)^2))*k2(4));
     [R,n]=fminsearch(pR,700);%求R
end


Tags:直达波 
作者:直达波 来源:直达波
  • 您是如何找到本站的?
  • 百度搜索
  • Google搜索
  • 查阅资料过程中
  • 论坛发现
  • 百度贴吧发现
  • 朋友介绍
本站最新成功开发工程项目案例
相关文章
  • 没有相关文章
相关评论
发表我的评论
  • 大名:
  • 内容:
  • matlab代做|matlab专业代做|matlab淘宝代做(www.hslogic.com) © 2018 版权所有 All Rights Reserved.
  • Email:highspeed_logic@163.com 站长QQ: 1224848052

    专业代做/代写/承接、MATLAB、SIMULINK、FPGA项目、博士/硕士/本科毕业设计、课题设计、论文,毕业论文,Coursework、Eassy、Assignment