سلام بیا این یک نمونه که قبلا خودم نوشتم مشابه سوال شماست
clear;
close all;
clc;
A=[0 , 1 , 0 , 0
0 , 0 ,-1 , 0
0 , 0 , 1 , 0
0 , 0 , 5 , 0]
B=[0
1
0
2];
C=[1 0 0 0];
D=0;
sys1=ss(A,B,C,D);
G=tf(sys1)
%step(G)
P=pole(sys1);
q=sum(real(P)>0);
if q>0
disp([' '...