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

Eatin520

现在用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

。。


现在根据这个模型编程：


以下是我的代码：
clc;clear;
x0=[150 85 150 145 130 0];
y0=[140 85 155 50 150 0];
q=[243 236 220.5 159 230 52];

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=[1 1 1 1 1 1]';
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=[t(1,2);t(1,3);t(1,4);t(1,5);t(1,6);t(2,3);t(2,4);t(2,5);t(2,6);t(3,4);t(3,5);t(3,6);t(4,5);t(4,6);t(5,6)];
VLB=-(pi/6)*eye(6,6);
VUB=(pi/6)*eye(6,6);

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

scarofsky

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