建模2017A题 角度lingo代码
Posted francischeng
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了建模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
以上是关于建模2017A题 角度lingo代码的主要内容,如果未能解决你的问题,请参考以下文章