建模2017A题 角度lingo代码

 1 model:
 2 
 3 sets:
 4 set/1..37/:mm,t,l,x,y,delta,m;
 5 endsets
 6 
 7 data:
 8 mm = 
 9 6.5414     
10 15.1468     
11 20.3436     
12 24.4151     
13 27.8592     
14 30.8865     
15 33.6098     
16 36.0980     
17 38.3967     
18 40.5383     
19 42.5464     
20 44.4393     
21 46.2309     
22 47.9328     
23 49.5540     
24 51.1024     
25 52.5843     
26 54.0053     
27 55.3699     
28 56.6824     
29 57.9461     
30 59.1643     
31 60.3397     
32 61.4747     
33 62.5716     
34 63.6324     
35 64.6587     
36 65.6522     
37 66.6145     
38 67.5468     
39 68.4503     
40 69.3262     
41 70.1754     
42 70.9991     
43 71.7980     
44 72.5730     
45 73.3248     
46 ;
47 
48 enddata
49 theta>29*3.1415/180;
50 theta<30*3.1415/180;
51 d=0.2768;
52 x1<-60;
53 x1>-70;
54 a = 15;
55 b = 40;
56 y1 = @sin(theta)/@cos(theta)*x1;
57 k = -1*@cos(theta)/@sin(theta);
58 u = 1.7724;
59 r = 189;
60 @for(set(i):x(i)=x1+(i+r-1)*@cos(theta)*d);
61 @for(set(i):y(i)=y1+(i+r-1)*@sin(theta)*d);
62 @for(set(i):m(i) = y(i)-k*x(i));
63 @for(set(i):delta(i)=4*k^2*m(i)^2*a^4-4*(a^2*k^2+b^2)*(a^2*m(i)^2-a^2*b^2));
64 @for(set(i):delta(i)>=0);
65 
66 min = @sum(set(i):(mm(i)-(@sqrt((1+k^2)*(delta(i)))/(a^2*k^2+b^2))*u)^2);
67 
68 
69 end