Iron_man. 发表于 2014-4-2 11:46

数学建模lingo 程序报错 求大神指导啊

model:
sets:
row/1..6/;
col/1..6/:da,x,y,a;
matrix(row,col):b,c,r;
endsets
min=@sum(col:(da)^2);
@for(col:@bnd(-0.0833333,da,0.0833333));
@for(row(i):@for(col(j)|i#gt#j:
@if(b(i,j)#gt#0,2*b(i,j)+da(i)+da(j),-(2*b(i,j)+da(i)+da(j)))>2*c(i,j)));
@for(row(i):@for(col(j)|i#GT#j:r(i,j)=@sqrt((y(j)-y(i))^2+(x(j)-x(i))^2)));
@for(row(i):@for(col(j)|i#GT#j:c(i,j)=@asin(8/r(i,j))));
@for(row(i):@for(col(j)|i#GT#j:b(i,j)=@atan((@sin(a(j))-@sin(a(i)))/(@cos(a(j))-@cos(a(i))))-@atan((y(i)-y(j))/(x(i)-x(j)))));
@for(row(i):@for(col(j):@free(c(i,j))));
@for(row(i):@for(col(j):@free(b(i,j))));
@for(col:@free(da));
init:
da=0 0 0 0 0 0;
endinit
data:
x=150, 85, 150, 145, 130, 0;
y=140, 85, 155, 50, 150, 0;
a=0.675, 0.6555, 0.6125, 0.4417, 0.6389, 0.1444;
enddata
end






数学建模飞行管理问题的一个老题,报错内容为 Arithmetic error in constraint: 48
lingo初学者,感激不尽。

Iron_man. 发表于 2014-4-2 11:47

- -怎么冒号都变成愤怒的表情了

madio 发表于 2014-4-2 12:19

需要禁用表情,表情符号和代码一样的

Iron_man. 发表于 2014-4-2 13:13

madio 发表于 2014-4-2 12:19 static/image/common/back.gif
需要禁用表情,表情符号和代码一样的

嗯嗯谢谢!

Iron_man. 发表于 2014-4-2 13:19

madio 发表于 2014-4-2 12:19 static/image/common/back.gif
需要禁用表情,表情符号和代码一样的

Arithmetic error in constraint: 48   
算术错误意味着什么呢?
我写了三条双层的for循环赋值语句,可我已经考虑到了不能出现r(i,i)的情况(分母不能为零)并且限制了j>i
为什么还报错

Iron_man. 发表于 2014-4-2 13:20

madio 发表于 2014-4-2 12:19 static/image/common/back.gif
需要禁用表情,表情符号和代码一样的

Arithmetic error in constraint: 48   
算术错误意味着什么呢?
我写了三条双层的for循环赋值语句,可我已经考虑到了不能出现r(i,i)的情况(分母不能为零)并且限制了j>i
为什么还报错

madio 发表于 2014-4-2 14:43

@for(row(i):@for(col(j)|i#GT#j:b(i,j)=@atan((@sin(a(j))-@sin(a(i)))/(@cos(a(j))-@cos(a(i))))-@atan((y(i)-y(j))/(x(i)-x(j)))));就是这句代码的问题,你具体说说要写一个什么样子的约束?

Iron_man. 发表于 2014-4-2 18:20

  就是这样一个式子。

Iron_man. 发表于 2014-4-2 18:21

madio 发表于 2014-4-2 14:43 static/image/common/back.gif
就是这句代码的问题,你具体说说要写一个什么样子的约束?

就是这样一个式子。

madio 发表于 2014-4-3 20:14

Iron_man. 发表于 2014-4-2 18:21 static/image/common/back.gif
就是这样一个式子。

你这个错误就在分母上,因为x的输入值有相同的值,就是两个150,这样一做减法就可能出现分母为零的情况,估计是你的输入数据有误.
页: [1] 2
查看完整版本: 数学建模lingo 程序报错 求大神指导啊