سلام. برنامه ی زیر صورت و مخرج حلقه بازو میگیره و با یافتن تابع تبدیل حلقه بسته پایداریو بررسی میکنه
n=input('enter num of G(s)=');
d=input('enter den of G(s)=');
na=length
![Thumbs down (n) (n)](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
;
nb=length(d);
D=d+[zeros(1,nb-na),n]; % yaftane moadele moshakhase
r=roots(D);
k=0;
for i=1:length(r)
if real(r(i))>0
k=k+1;
end
end
if k>0
disp('the system is unstable' );
else
disp('the system is stable' );
end