Your problem has nine constraints. Maybe your problem can be strongly simplified to five constraints if you make theta_i =exp(x_i) and replace theta_i by this new variable in all places. Therefore, you have eliminated the positivity constraints and the new problem depends on x_i (x_i is your new variable). Ok.... You find the optimal value of x_i and calculate theta_i=exp(x_i). This is a very common substitution in Econometrics when you are dealing for instance with variances or volatilities.
You can also try another substitution (I have not seem before, but it seems to work) to eliminate all lb or ub... Make y=exp(x)/(1+exp(x)) [logistic function]. Now your problem is much easier, since it has only one constraint (given by A and b) and follow the same procedure above.