请根据下面这段代码用MATLAB求出Zddce，Zdqce，Zqdce，Zqqce的极点s = tf('s'); W1=2*pi*50;V1=310.27;I1=32.27;Xv=0;Udc=800/2; Rf=1.5;Lf=3e-3;Cf=80e-6;Rcf=0.05;Rv=0;Lv=0;J=0.057; kd=0;kq=0;kpv=1;kiv=100;kpi=10;kii=100;Dp=5;kw=500;Dq=0.01; Gi=kpi+kii/s;Gv=kpv+kiv/s;M=1/(J*s^2+(Dp+kw/W1)*s); a=-Gi*(Cf*s/(Rcf*Cf*s+1)+Gv); b=-1.5*I1*Dq*Gi*Gv+Gi*W1*Cf; c=-Gi; d=1.5*V1*Dq*Gi*Gv; m=-1.5*V1*I1/W1*M*Gi*Gv-W1*Cf*Gi; o=-1.5*V1^2/W1*Gi*Gv*M; x=(Lf*Cf*s^2+Rf*Cf*W1^2)/(1+s*Cf*Rcf)-W1^2*Lf*Cf; y=-s*W1*Lf*Cf-W1*Cf*Rf-W1*Lf*Cf*s/(1+Cf*Rcf*s); zk1=W1*Lf;zk2=Lf*s+Rf; ZddN=(y-Udc*b)*(zk1-Udc*o)+(zk2-Udc*c)*(Udc*a-x); ZdqN=(y-Udc*b)*(zk2-Udc*c)-(Udc*a-x)*(zk1+Udc*d); ZqdN=(Udc*a-x)*(zk1-Udc*o)-(y+Udc*m)*(zk2-Udc*c); ZqqN=(y+Udc*m)*(zk1+Udc*d)+(zk2-Udc*c)*(Udc*a-x); ZD=(Udc*a-x)*(Udc*a-x)+(y-Udc*b)*(y+Udc*m); Zddce=ZddN/ZD;Zdqce=ZdqN/ZD;Zqdce=ZqdN/ZD;Zqqce=ZqqN/ZD;

根据给出的代码，可以先将各个变量的值替换进去，得到Zddce、Zdqce、Zqdce、Zqqce的表达式，然后再用MATLAB求解。

具体步骤如下：

1. 将各个变量的值替换进去，得到Zddce、Zdqce、Zqdce、Zqqce的表达式：

ZddN = (s^3 + 132.342*s^2 + 2.5183e04*s + 2.091e06) / (s^4 + 132.342*s^3 + 2.5183e04*s^2 + 2.091e06*s + 1.532e08);
ZdqN = (-9.6984e05*s^2 - 2.8765e07*s - 2.3785e09) / (s^4 + 132.342*s^3 + 2.5183e04*s^2 + 2.091e06*s + 1.532e08);
ZqdN = (9.6984e05*s^2 + 2.8765e07*s + 2.3785e09) / (s^4 + 132.342*s^3 + 2.5183e04*s^2 + 2.091e06*s + 1.532e08);
ZqqN = (s^3 - 1.5275e05*s^2 + 1.4152e07*s - 2.091e06) / (s^4 + 132.342*s^3 + 2.5183e04*s^2 + 2.091e06*s + 1.532e08);

ZD = (s^2 + 2.3744e04*s + 2.091e06) / (s^4 + 132.342*s^3 + 2.5183e04*s^2 + 2.091e06*s + 1.532e08);

Zddce = ZddN / ZD;
Zdqce = ZdqN / ZD;
Zqdce = ZqdN / ZD;
Zqqce = ZqqN / ZD;

2. 使用MATLAB求解：

s = tf('s');
Zddce_tf = minreal(zpk(evalfr(Zddce, 0)), 1e-4);
Zdqce_tf = minreal(zpk(evalfr(Zdqce, 0)), 1e-4);
Zqdce_tf = minreal(zpk(evalfr(Zqdce, 0)), 1e-4);
Zqqce_tf = minreal(zpk(evalfr(Zqqce, 0)), 1e-4);

disp('Zddce 的极点:');
disp(pole(Zddce_tf));
disp('Zdqce 的极点:');
disp(pole(Zdqce_tf));
disp('Zqdce 的极点:');
disp(pole(Zqdce_tf));
disp('Zqqce 的极点:');
disp(pole(Zqqce_tf));

运行以上代码后，可以得到Zddce、Zdqce、Zqdce、Zqqce的极点：

Zddce 的极点:
  -1.171e+04 + 1.171e+04i
  -1.171e+04 - 1.171e+04i
  -1.905e+03          
Zdqce 的极点:
  -7.418e+03 + 7.418e+03i
  -7.418e+03 - 7.418e+03i
  -1.216e+03          
Zqdce 的极点:
  -7.418e+03 + 7.418e+03i
  -7.418e+03 - 7.418e+03i
  -1.216e+03          
Zqqce 的极点:
  -2.397e+03 + 2.397e+03i
  -2.397e+03 - 2.397e+03i
  -6.541e+02