One approach —
K = 1; % Insert Correct Value
p1 = -3.02E4 + 2.52E4j;
p2 = conj(p1);
p3 = -1.16E5 + 1.13E6j;
p4 = conj(p3);
z1 = 2.35E5;
z2 = -3.33E5;
s = tf('s');
Gp = K*(s-z1)*(s-z2) / ((s+p1)*(s+p2)*(s+p3)*(s+p4))
figure
bodeplot(Gp)
figure
pzmap(Gp)
figure
stepplot(Gp)
The right-half-plane poles concern me. See the Step Response plot to understand the reason.
.
