Infeasible appears after one iteration

Hi everyone,I really need your help. I have an iterative program that consists of a objective function and several constraints. In the iterative program, there are some initial points that have to be updated in each iteration. In each of the iterations,they are derived and used in next iteration as given points. my code is solved as follows. After the first iteration, lines 11 to 16 will be commented out.

N = 3;
w = 0.5;
No = 10^(-3);
p_max = 3;
D = [210^(2),310^(2),410^(2)];
c = [8,10,10];
apha = 10^(-19);
f_max = 3
10^(2);
B = 10^(2);
h = [0.6,1.8,1];
% beita_i = [0.4,0.4,0.1];
% p_i = [1,1.5,2.5];
% phai_i = No./h./p_i + 0.1;
% R_i = log2(1+1./phai_i) - 0.1;
% miut_i = [B./D.*beita_i.*R_i];
% Et_i = p_i./miut_i;

[optval,R_i,beita_i,phai_i,p_i,miut_i,Et_i] = interation(N,w,D,B,No,h,c,apha,f_max,p_max,R_i,beita_i,phai_i,p_i,miut_i,Et_i);

function [cvx_optval,R,beita,phai,p,miut,Et] = interation(N,w,D,B,No,h,c,apha,f_max,p_max,R_i,beita_i,phai_i,p_i,miut_i,Et_i)
cvx_begin
variable lamda(N)
variable p(N)
variable beita(N)
variable f(N)
variable tao
variable A(N)
variable E(N)
variable miut(N)
variable miuc(N)
variable Et(N)
variable Ec(N)
variable R(N)
variable phai(N)

minimize(tao)
subject to
for i = 1:N
    1*inv_pos(lamda(i)) + 1*inv_pos(miut(i)-lamda(i)) + 1*inv_pos(miuc(i)-lamda(i)) <= A(i);
    2/B*D(i)*miut(i) + R(i)^2 + beita(i)^2 - (R_i(i)+beita_i(i))^2 - 2*(R_i(i)+beita_i(i))*(R(i)+beita(i)-R_i(i)-beita_i(i)) <= 0;
    R(i) - log2(1+1/phai_i(i)) - (phai(i)-phai_i(i))/(log(2)*(phai_i(i)+phai_i(i)^2)) <= 0;
    2*No/h(i) + p(i)^2 + phai(i)^2 - (p_i(i)+phai_i(i))^2 - 2*(phai_i(i)+p_i(i))*(p(i)-p_i(i)+phai(i)-phai_i(i)) <= 0;
    miuc(i) - f(i)/(c(i)*D(i))  <= 0;
    2*p(i) + miut(i)^2 + Et(i)^2 - (miut_i(i)+Et_i(i))^2 - 2*(miut_i(i)+Et_i(i))*(miut(i)-miut_i(i)+Et(i)-Et_i(i)) <= 0;
    apha*c(i)*D(i)*f(i)*f(i) <= Ec(i);
end
w.*A + (1-w).*E <= tao;
Et + Ec <= E
lamda - miut <= 0;
lamda - miuc <= 0;
sum(beita) <= 1;
sum(f) <= f_max;
0 <= p < p_max;
0 <= lamda;
0 <= beita; 
0 <= f;    
cvx_end 

end

The first iteration is successful with the following information:
Calling SDPT3 4.0: 140 variables, 67 equality constraints
For improved efficiency, SDPT3 is solving the dual problem.

num. of constraints = 67
dim. of sdp var = 60, num. of sdp blk = 30
dim. of linear var = 50


SDPT3: Infeasible path-following algorithms


version predcorr gam expon scale_data
HKM 1 0.000 1 0
it pstep dstep pinfeas dinfeas gap prim-obj dual-obj cputime

0|0.000|0.000|9.5e+02|3.1e+00|1.8e+05|-5.857139e+03 0.000000e+00| 0:0:00| chol 1 1
1|0.354|0.313|6.1e+02|2.1e+00|1.5e+05|-2.387643e+03 -3.641937e+02| 0:0:00| chol 1 1
2|0.800|0.492|1.2e+02|1.1e+00|8.7e+04| 2.655878e+03 -8.345807e+02| 0:0:00| chol 1 1
3|0.651|0.715|4.3e+01|3.1e-01|3.6e+04| 3.807633e+03 -1.532060e+03| 0:0:00| chol 1 1
4|0.860|0.565|6.0e+00|1.3e-01|2.2e+04| 4.289365e+03 -1.812670e+03| 0:0:00| chol 1 1
5|0.887|0.834|6.8e-01|2.2e-02|6.7e+03| 2.504203e+03 -1.309745e+03| 0:0:00| chol 1 1
6|0.775|0.750|1.5e-01|5.6e-03|2.5e+03| 1.110197e+03 -5.958138e+02| 0:0:00| chol 1 1
7|0.902|0.862|1.5e-02|7.6e-04|8.8e+02| 5.081015e+02 -2.498136e+02| 0:0:00| chol 1 1
8|0.674|0.511|4.8e-03|4.4e-04|6.0e+02| 2.844856e+02 -1.869016e+02| 0:0:00| chol 1 1
9|0.808|0.762|9.3e-04|1.5e-04|2.5e+02| 8.490869e+01 -1.094776e+02| 0:0:00| chol 1 1
10|0.859|0.882|1.3e-04|3.0e-05|6.4e+01|-2.374245e+01 -7.573400e+01| 0:0:00|# chol 1 1
11|0.851|0.964|2.0e-05|4.1e-06|1.5e+01|-5.599873e+01 -6.859277e+01| 0:0:00|# chol 1 1
12|1.000|1.000|6.8e-07|1.4e-06|7.2e+00|-6.148946e+01 -6.782088e+01| 0:0:00|# chol 1 1
13|0.940|1.000|2.2e-07|1.4e-07|2.1e+00|-6.439544e+01 -6.634956e+01| 0:0:00|# chol 1 1
14|1.000|1.000|7.8e-08|4.3e-08|7.8e-01|-6.513145e+01 -6.586961e+01| 0:0:00|# chol 1 1
15|0.957|1.000|2.2e-08|1.6e-08|2.1e-01|-6.542955e+01 -6.562190e+01| 0:0:00|# chol 1 1
16|1.000|1.000|3.4e-09|4.3e-09|3.5e-02|-6.551603e+01 -6.554817e+01| 0:0:00|# chol 1 1
17|0.986|0.960|5.1e-11|3.9e-10|1.1e-03|-6.553188e+01 -6.553275e+01| 0:0:00|# chol 1 1
18|0.966|0.931|7.7e-12|3.7e-11|9.0e-05|-6.553212e+01 -6.553219e+01| 0:0:00|# chol 1 1
19|0.993|1.000|1.9e-11|1.5e-12|9.6e-06|-6.553214e+01 -6.553215e+01| 0:0:00|# chol 2 2
20|1.000|0.995|4.4e-13|4.0e-14|1.6e-07|-6.553214e+01 -6.553214e+01| 0:0:00|
stop: max(relative gap, infeasibilities) < 1.49e-08

number of iterations = 20
primal objective value = -6.55321431e+01
dual objective value = -6.55321432e+01
gap := trace(XZ) = 1.60e-07
relative gap = 1.21e-09
actual relative gap = 1.04e-09
rel. primal infeas (scaled problem) = 4.39e-13
rel. dual " " " = 4.04e-14
rel. primal infeas (unscaled problem) = 0.00e+00
rel. dual " " " = 0.00e+00
norm(X), norm(y), norm(Z) = 1.7e+03, 1.6e+05, 1.6e+05
norm(A), norm(b), norm© = 1.9e+02, 2.0e+00, 9.4e+02
Total CPU time (secs) = 0.14
CPU time per iteration = 0.01
termination code = 0
DIMACS: 4.4e-13 0.0e+00 4.3e-14 0.0e+00 1.0e-09 1.2e-09


Status: Solved
Optimal value (cvx_optval): +65.5321

And the second iteration shows infeasible:
Calling SDPT3 4.0: 140 variables, 67 equality constraints
For improved efficiency, SDPT3 is solving the dual problem.

num. of constraints = 67
dim. of sdp var = 60, num. of sdp blk = 30
dim. of linear var = 50


SDPT3: Infeasible path-following algorithms


version predcorr gam expon scale_data
HKM 1 0.000 1 0
it pstep dstep pinfeas dinfeas gap prim-obj dual-obj cputime

0|0.000|0.000|3.5e+02|3.6e+00|1.7e+05|-3.641329e+03 0.000000e+00| 0:0:00| chol 1 1
1|0.593|0.416|1.4e+02|2.1e+00|1.2e+05| 5.717214e+01 -4.524804e+02| 0:0:00| chol 1 1
2|0.610|0.694|5.6e+01|6.5e-01|5.1e+04| 1.681368e+03 -1.028308e+03| 0:0:00| chol 1 1
3|0.827|0.480|9.6e+00|3.4e-01|3.3e+04| 3.489181e+03 -1.263383e+03| 0:0:00| chol 1 1
4|0.942|0.767|5.6e-01|7.9e-02|1.2e+04| 2.823071e+03 -1.421320e+03| 0:0:00| chol 1 1
5|0.799|0.677|1.1e-01|2.5e-02|5.1e+03| 1.133171e+03 -1.002391e+03| 0:0:00| chol 1 1
6|0.625|0.473|4.2e-02|1.3e-02|3.5e+03|-1.172505e+02 -7.946819e+02| 0:0:00| chol 1 1
7|0.714|0.337|1.2e-02|9.0e-03|3.4e+03|-3.627430e+03 -6.797546e+02| 0:0:00| chol 1 1
8|0.231|0.192|9.3e-03|7.3e-03|3.7e+03|-1.902427e+04 -7.257486e+02| 0:0:00| chol 1 1
9|0.028|0.042|9.1e-03|7.0e-03|7.1e+03|-9.612565e+04 -1.161559e+03| 0:0:00| chol 1 1
10|0.012|0.022|9.0e-03|6.8e-03|2.0e+04|-5.392080e+05 -2.630191e+03| 0:0:00| chol 1 1
11|0.006|0.014|8.9e-03|6.7e-03|7.8e+04|-3.445642e+06 -7.929886e+03| 0:0:00| chol 1 1
12|0.003|0.008|8.9e-03|6.7e-03|3.6e+05|-2.475795e+07 -2.748553e+04| 0:0:00| chol 1 2
13|0.001|0.006|8.9e-03|6.7e-03|1.7e+06|-1.914023e+08 -1.131707e+05| 0:0:00| chol 1 2
14|0.001|0.003|8.9e-03|6.6e-03|9.5e+06|-1.562514e+09 -4.497484e+05| 0:0:00| chol 2 2
15|0.000|0.002|8.9e-03|6.6e-03|6.4e+07|-1.313695e+10 -1.415724e+06| 0:0:00| chol 2 2
16|0.000|0.001|8.9e-03|6.6e-03|4.1e+08|-1.115148e+11 -5.295921e+06| 0:0:00| chol 2 2
17|0.000|0.000|8.9e-03|6.6e-03|3.5e+09|-9.585540e+11 -9.997203e+06| 0:0:00| chol 2 2
18|0.000|0.000|8.9e-03|6.6e-03|3.2e+10|-8.137121e+12 -1.375726e+07| 0:0:00| chol 2 2
19|0.000|0.000|1.1e-02|6.6e-03|2.8e+11|-6.628284e+13 -1.562868e+07| 0:0:00| chol 2 2
20|0.000|0.000|3.4e-02|6.6e-03|2.2e+12|-4.835603e+14 -3.454159e+07| 0:0:00| chol 2 2
21|0.000|0.000|2.7e-01|6.6e-03|1.4e+13|-2.850740e+15 -4.661043e+07| 0:0:00| chol 2 2
22|0.000|0.000|5.3e-01|6.6e-03|3.5e+13|-7.106397e+15 -2.085241e+08| 0:0:00| chol 2 2
23|0.000|0.000|4.0e+00|6.6e-03|1.4e+14|-2.897465e+16 -2.824141e+08| 0:0:00| chol 2 2
24|0.000|0.000|7.9e+00|6.6e-03|3.9e+14|-8.026409e+16 -1.404260e+09| 0:0:00| chol 2 2
25|0.000|0.000|3.2e+01|6.6e-03|2.4e+15|-4.978376e+17 -1.919086e+09| 0:0:00|
sqlp stop: dual problem is suspected of being infeasible

number of iterations = 25
residual of dual infeasibility
certificate X = 1.28e-16
reldist to infeas. <= 1.32e-18
Total CPU time (secs) = 0.14
CPU time per iteration = 0.01
termination code = 2
DIMACS: 3.2e+01 0.0e+00 7.3e-03 0.0e+00 -1.0e+00 4.8e-03


Status: Infeasible
Optimal value (cvx_optval): +Inf

check http://ask.cvxr.com/t/how-to-write-log-1-ax-a-are-two-constant/9474/21 for similar advice.