数学建模社区-数学中国

标题: 求各位帮帮忙。。急急急。。。 [打印本页]

作者: 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

。。

现在根据这个模型编程：

以下是我的代码：
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)

2.jpg (3.76 KB, 下载次数: 147)

作者: scarofsky 时间: 2011-5-11 18:59
可能不止一出错。楼主先把/ 和./ ^ 和.^弄清楚再说哦

作者: alair004 时间: 2012-2-6 16:30
楼主分享的很好。。。5060873639047146

欢迎光临数学建模社区-数学中国 (http://www.madio.net/)