Eatin520 发表于 2011-5-11 17:58

求各位帮帮忙。。急急急。。。

现在用matlab编一段代码。。可是老是得不到想要的答案。。不知哪里出错了。。可否请大家帮我看看。。指点一下。。谢谢。。信息有：
表
1
飞行记录数据
  飞机编号  横座标
x  纵座标
y  方向角（度）
  1  150  140  243
  2  85  85  236
  3  150  155  220.5
  4  145  50  159
  5  130  150  230
  新进入  0  0  52
http://www.madio.cn/forum.php?mod=image&aid=94872&size=300x300&key=503fxze95O1z32FgSweh%2B7QjbiSPhjdZJf%2FSmambDRK1%2F0vTLjHy8yYo&nocache=yes&type=fixnone。。


现在根据这个模型编程：


以下是我的代码：
clc;clear;
x0=;
y0=;
q=;

xy1=x0+i*y0;
xy2=exp(i*q*pi/180);
for m=1:6
  for n=1:6
      if n~=m
          b0(m,n)=angle((xy2(n)-xy2(m))/(xy1(m)-xy1(n)));

          d1=(x0(m)-x0(n))^2;
          d2=(y0(m)-y0(n))^2;
          d=d1+d2;
          r(m,n)=sqrt(d);
          a(m,n)=asin(8/r(m,n)) ;
      end
  end
end

C=';
t(m,n)=-(a(m,n)-b0(m,n));
A=[-1 -1 0 0 0 0
   -1 0 -1 0 0 0
   -1 0 0 -1 0 0
   -1 0 0 0 -1 0
   -1 0 0 0 0 -1
   0 -1 -1 0 0 0
   0 -1 0 -1 0 0
   0 -1 0 0 -1 0
   0 -1 0 0 0 -1
   0 0 -1 -1 0 0
   0 0 -1 0 -1 0
   0 0 -1 0 0 -1
   0 0 0 -1 -1 0
   0 0 0 -1 0 -1
   0 0 0 0 -1 -1];
B=;
VLB=-(pi/6)*eye(6,6);
VUB=(pi/6)*eye(6,6);

=linprog(C,A,B,[],[],VLB,VUB)


scarofsky 发表于 2011-5-11 18:59

可能不止一出错。楼主先把/ 和./    ^ 和.^弄清楚再说哦

alair004 发表于 2012-2-6 16:30

页: [1]
查看完整版本: 求各位帮帮忙。。急急急。。。