符号函数操作

Halcom

solve进行等式求解，然后采用subs进行赋值处理，最后进行double()数据类型转换

1. syms n1 q1 p1 A1 q2 p2 A2 B1 B2 n2
   
2. [p1,p2] = solve((1-n1)*q1+(p1-0.003*q1-11)*A1 + n1*((63-2*(q1+q2)-p1)*A1+(63-2*(q1+q2)-p2)*B1)==0,...
   
3. (1-n2)*q2+(p2-0.004*q2-15)*A2+n2*( (63-2*(q1+q2)-p2)*A2 + (63-2*(q1+q2)-p1)*B2 )==0 , p1, p2);
   
4. % p1 =
   
5. % -(11000*A1*A2 - 1000*A2*q1 - 63000*A1*A2*n1 - 11000*A1*A2*n2 - 48000*A2*B1*n1 + 3*A1*A2*q1 + 1000*A2*n1*q1 + 1000*A2*n2*q1 - 1000*B1*n1*q2 + 63000*A1*A2*n1*n2 - 63000*B1*B2*n1*n2 + 2000*A1*A2*n1*q1 + 2000*A1*A2*n1*q2 - 3*A1*A2*n2*q1 + 2000*A2*B1*n1*q1 + 2004*A2*B1*n1*q2 - 1000*A2*n1*n2*q1 + 1000*B1*n1*n2*q2 - 2000*A1*A2*n1*n2*q1 - 2000*A1*A2*n1*n2*q2 + 2000*B1*B2*n1*n2*q1 + 2000*B1*B2*n1*n2*q2)/(1000*(A1*A2*n1 - A1*A2 + A1*A2*n2 - A1*A2*n1*n2 + B1*B2*n1*n2))
   
6. % p2 =
   
7. % -(15000*A1*A2 - 1000*A1*q2 - 15000*A1*A2*n1 - 63000*A1*A2*n2 - 52000*A1*B2*n2 + 4*A1*A2*q2 + 1000*A1*n1*q2 + 1000*A1*n2*q2 - 1000*B2*n2*q1 + 63000*A1*A2*n1*n2 - 63000*B1*B2*n1*n2 - 4*A1*A2*n1*q2 + 2000*A1*A2*n2*q1 + 2000*A1*A2*n2*q2 + 2003*A1*B2*n2*q1 + 2000*A1*B2*n2*q2 - 1000*A1*n1*n2*q2 + 1000*B2*n1*n2*q1 - 2000*A1*A2*n1*n2*q1 - 2000*A1*A2*n1*n2*q2 + 2000*B1*B2*n1*n2*q1 + 2000*B1*B2*n1*n2*q2)/(1000*(A1*A2*n1 - A1*A2 + A1*A2*n2 - A1*A2*n1*n2 + B1*B2*n1*n2))
   
8. y = subs(p1, [A1,A2,B1,B2,n1,n2,q1,q2], [1, 1, 1, 1, 1, 1, 1, 1])
   
9. y1= double(y);

复制代码

Halcom

适应度函数写为：

1. function y = fun(x, p1, p2 )
   
2. %        q1,q2, p1,p2,e1,e2,n1,n2
   
3. syms n1 n2 q1 q2 A1 A2 B1 B2
   
4. q1_ = x(1)/500;
   
5. q2_ = x(2)/300;
   
6. e1_ = x(3);
   
7. e2_ = x(4);
   
8. n1_ = x(5);
   
9. n2_ = x(6);
   
10. if(e1_==0 && e2_==0)
    
11.     A1_ = -1/2;
    
12. elseif(e1_==0 && e2_~=0)
    
13.     A1_ = -1/2*(1+1/(0.002*e2_+eps));
    
14. elseif(e1_~=0 && e2_==0)
    
15.     A1_ = (e1_-1)/(0.003*e1_ + eps);
    
16. elseif(e1_~=0 && e2_~=0)
    
17.     A1_ = (1+1/(0.002*e2_+eps))*(e1_-1)/( 0.003*e1_*(1+1/(0.0015*e1_+eps) + 1/(0.002*e2_+eps)) +eps );
    
18. end
    
19. 
    
20. if(e1_==0 && e2_==0)
    
21.     B1_ = 0;
    
22. elseif(e1_==0 && e2_~=0)
    
23.     B1_ = 1/(0.004*e2+eps);
    
24. elseif(e1_~=0 && e2_==0)
    
25.     B1_ = (1-e1_)/(0.003*e1_+eps);
    
26. elseif(e1_~=0 && e2_~=0)
    
27.     B1_ = (1-e1_)/(0.002*e2_+eps)/( 0.003*e1_*(1+1/(0.0015*e1_+eps)+1/(0.002*e2_+eps)) + eps );
    
28. end
    
29. 
    
30. if(e1_==0 && e2_==0)
    
31.     A2_ = -1/2;
    
32. elseif(e1_==0 && e2_~=0)
    
33.     A2_ = (e2_-1)/(0.004*e2_+eps);
    
34. elseif(e1_~=0 && e2_==0)
    
35.     A2_ = -1/2*( 1+1/(0.0015*e1_+eps) );
    
36. elseif(e1_~=0 && e2_~=0)
    
37.     A2_ = (1+1/(0.0015*e1_+eps))*(e2_-1)/( 0.004*e2_*(1+1/(0.0015*e1_+eps) + 1/(0.002*e2_+eps)) +eps );
    
38. end
    
39. 
    
40. if(e1_==0 && e2_==0)
    
41.     B2_ = 0;
    
42. elseif(e1_==0 && e2_~=0)
    
43.     B2_ = (1-e2_)/(0.004*e2_+eps);
    
44. elseif(e1_~=0 && e2_==0)
    
45.     B2_ = (1)/(0.003*e1_+eps);
    
46. elseif(e1_~=0 && e2_~=0)
    
47.     B2_ = (1-e2_)/(0.0015*e1_+eps)/( 0.004*e2_*(1+1/(0.0015*e1_+eps)+1/(0.002*e2_+eps)) + eps );
    
48. end
    
49. p1_ = subs(p1, [A1,A2,B1,B2,n1,n2,q1,q2], [A1_,A2_,B1_,B2_,n1_,n2_,q1_,q2_]);
    
50. p1_ = double( p1_ );
    
51. p2_ = subs(p2, [A1,A2,B1,B2,n1,n2,q1,q2], [A1_,A2_,B1_,B2_,n1_,n2_,q1_,q2_]);
    
52. p2_ = double( p2_ );
    
53. % λ计算
    
54. if(n1~=0)
    
55.     lamda1 = 0;
    
56. end
    
57. if(n1~=1)
    
58.     lamda2 = 0;
    
59. end
    
60. if(n2~=0)
    
61.     lamda3 = 0;
    
62. end
    
63. if(n2~=1)
    
64.     lamda4 = 0;
    
65. end
    
66. lamda12 = (63-(q1_+q2_)-p1_)*q1_ + (63-(q1_+q2_)-p2_)*q2_;
    
67. lamda34 = (63-(q1_+q2_)-p1_)*q1_ + (63-(q1_+q2_)-p2_)*q2_;
    
68. %% 目标函数计算
    
69. if(e1_+e2_<=1)
    
70.     y = (63-(q1_+q2_)-p1_)*q1_ + (63-(q1_+q2_)-p2_)*q2_ + e1_*((p1_-11)*q1_-0.0015*q1_*q1_)+e2_*((p2_-15)*q2_-0.002*q2_*q2_);
    
71. else
    
72.     y=-1;
    
73. end
    

复制代码