how can get rid of "ValueError: all the input array dimensions except for the concatenation axis must match exactly" error when use cvxopt function?

Question

first i define some matrix and vector in proper shape .

initialization

I=np.eye(24)
Z=np.zeros((24,24))
a=0.012
b=1.1
gamma1=0.9/80
gamma2=1.1/80
MM=np.eye(24)
for i in range (22):
    MM[i+1,i]=-1
    MM[0,23]=-1
M=random.randint(200,300, size=(24,1))
max_pch=5
max_pdch=5
ppp=random.randint(150,200, size=(24,))

define matrix of objective function

Q and C is matrix and vector of objective function 1/2 x^T Q x +C^T x , respectively.

Q=np.asarray(np.bmat([[a*I,Z,Z,Z],[Z,a*I,Z,Z],[Z,Z,Z,Z],[Z,Z,Z,Z]  ]))
C=np.asarray(np.bmat([[b*np.ones(24),b*np.ones(24),0*np.ones(24),ppp]]))

##create equal subject

In my problem, I have just equal constraint and upper bound and lower bound that define blew.

Aeq=np.asarray(np.bmat([[-I,I,Z,I], [-gamma1*I, gamma2*I,MM,Z],[np.zeros((48,96))]]))
beq=np.asarray(np.bmat([[M],[np.zeros((72,1))]]))

##create upper and lower bound in shape (1,96)

lb=np.asarray(np.bmat([[0*np.ones(24),0*np.ones(24),[0.1],0.1*np.ones(22), 
    [0.1],100*np.ones(24)]]))

ub=np.asarray(np.bmat([[max_pch*np.ones(24),max_pdch*np.ones(24),[0.1],0.9*np.ones(22), 
    [0.9],500*np.ones(24)]]))

x = solve_qp(P=matrix(Q), q=C.T,
         G=None,h=None, A=matrix(Aeq), b=beq, lb=lb.T, ub=ub.T,solver='quadprog')

##error

---------------------------------------------------------------------------
ValueError                                Traceback (most recent call last)
<ipython-input-48-111d9695d5a8> in <module>
     25 
     26 x = solve_qp(P=matrix(Q), q=C.T,
---> 27              G=None,h=None, A=matrix(Aeq), b=beq, lb=lb.T, ub=ub.T,solver='quadprog')
     28 
     29 

~\Anaconda3\lib\site-packages\qpsolvers\__init__.py in solve_qp(P, q, G, h, A, b, lb, ub, 
solver, initvals, sym_proj, verbose, **kwargs)
    271     kwargs["verbose"] = verbose
    272     try:
--> 273         return __solve_function__[solver](*args, **kwargs)
    274     except KeyError:
    275         raise SolverNotFound(f"solver '{solver}' is not available")

~\Anaconda3\lib\site-packages\qpsolvers\quadprog_.py in quadprog_solve_qp(P, q, G, h, A, b, 
initvals, verbose)
     85         else:
     86             qp_C = -vstack([A, G]).T
---> 87             qp_b = -hstack([b, h])
     88         meq = A.shape[0]
     89     else:  # no equality constraint

~\Anaconda3\lib\site-packages\numpy\core\shape_base.py in hstack(tup)
    338         return _nx.concatenate(arrs, 0)
    339     else:
--> 340         return _nx.concatenate(arrs, 1)
    341 
    342 

ValueError: all the input array dimensions except for the concatenation axis must match 
exactly

If anybody can help me, I am glad.