建模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

 

转载于:https://www.cnblogs.com/francischeng/p/9462379.html

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值