算法：LPI-TL

输入： $\left(B_I^s,P_I^s\right),I=1,2,\cdots ,t_1$ ， $\left(B_J^t,P_J^t\right),J=1,2,\cdots ,t_2$ ，

${\lambda }_1,{\lambda }_2,C_1,C_2,\epsilon$

输出： $w_1,w_2,b_1,b_2$

方法：

a)：For i=1:(t₁+t₂)

b): $y_i=-\log \left(\frac{1}{p_i}-1\right)$

c): End

d): 调用Matlab的CVX包求解公式6，解得拉格朗日乘子a

e): 将a代入公式(7)-(9)求得w₁，w₂

f): 将w₁，w₂，y_i代入公式3和公式4可解得b₁，b₂