数学建模社区-数学中国

标题: lingo嵌套for的分段函数求解答！！！！ [打印本页]

作者: 3325760473 时间: 2022-12-23 21:57
标题: lingo嵌套for的分段函数求解答！！！！
sets:
index/1..15/:PE,alpha,beta,rho;
index2/1..16/:P,phi;
endsets
data:
pi=3.141592653589;
l1=1;
l4=4;
enddata
min=@sum(index(i):(PE(i)-P(i))^2);
l2=x2;
l3=x3;
@for(index(i):@if((phi(i))#le#pi,
P(i)=pi-alpha(i)-beta(i),
P(i)=pi-alpha(i)+beta(i)));!问题出在这一块
@for(index(i):
      @cos(alpha(i))=(rho(i)^2+l3^2-l2^2)/(2*rho(i)*l3));
@for(index(i):
      @cos(beta(i))=(rho(i)^2+l4^2-l1^2)/(2*rho(i)*l4));
@for(index(i):
      rho(i)=(l1^2+l4^2-2*l1*l4*@cos(phi(i)))^(0.5));
@cos(phi(16))=((l1+l2)^2+l4^2-l3^2)/(2*l4*(l1+l2));
@cos(P(16))=((l1+l2)^2-l4^2-l3^2)/(2*l4*l3);
@for(index(i):
      PE(i)=P(16)+(phi(i)-phi(16))^2*2/(3*pi));
@for(index(i):
      phi(i)=phi(16)+pi*i/(2*15));
x2^2+x3^2+x2*x3*(2^0.5)-25<=0;
-x2^2-x3^2+x2*x3*(2^0.5)+9<=0;
end

想写一个if中嵌套for的分段函数，分段函数如下

作者: ilikenba 时间: 2022-12-24 07:26
函数的图片上传一下

作者: 3325760473 时间: 2022-12-24 14:08

ilikenba 发表于 2022-12-24 07:26
函数的图片上传一下

pu(i)=pai-alpha(i)-beta(i) (phi(i)<=pai)
pu(i)=pai-alpha(i)+beta(i) (phi(i)>pai)

作者: 3325760473 时间: 2022-12-24 14:16

ilikenba 发表于 2022-12-24 07:26
函数的图片上传一下

函数长这样，求助大佬怎么编写

EKS@JYA)W1@5SDK}CF]HY3T.png (7.24 KB, 下载次数: 51)

EKS@JYA)W1@5SDK}CF]HY3T.png

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