Priklad do MATLABu - specha

[ghosthunter]

Shanim priklad v programu MATLAB, ktery by pocital nektera kriteria stability, napr. Hurwitzovo kritérium, Michajlovo kritérium, Nyquistovo kritérium atd..
Predem diky za jakoukoliv pomoc.

[Doktor-Prase]

Google - Nyquist + matlab

http://www.mathworks.com/matlabcentr...bjectType=File
http://ece.gmu.edu/~gbeale/ece_421/nyquist_01_mat.html
http://www.engin.umich.edu/group/ctm...lnyquist1.html

Někteří modi by tě hned odkázali na google, ale tak já nemám moc co dělat, tak jsem hodnej...

[ghosthunter]

Diky, ale tohle jsem uz taky nasel...Ohledne Nyquista je toho na netu dost. Me slo spis o algebraicka kriteria, treba Hurwitzovo, Routhovo...

[mudboy]

Diky, ale tohle jsem uz taky nasel...Ohledne Nyquista je toho na netu dost. Me slo spis o algebraicka kriteria, treba Hurwitzovo, Routhovo...

function H=hurwitz(den)
n=length(den)-1;
for i=1:n, i1=floor(i/2);
if i==i1*2, hsub1=den(1:2:n+1); i1=i1-1;
else, hsub1=den(2:2:n+1); end
l1=length(hsub1);
H(i,: )=[zeros(1,i1),hsub1,zeros(1,n-i1-l1)];
end

function [rtab,info]=routh0(den)
info=[];
vec1=den(1 : 2:length(den)); nrT=length(vec1);
vec2=den(2 : 2:length(den)-1);
rtab=[vec1; vec2, zeros(1,nrT-length(vec2))];
for k=1:length(den)-2,
alpha(k)=vec1(1)/vec2(1);
for i=1:length(vec2),
a3(i)=rtab(k,i+1)-alpha(k)*rtab(k+1,i+1);
end
if sum(abs(a3))==0
a3=polyder(vec2);
info=[info,'All elements in row ',...
int2str(k+2) ' are zeros;'];
elseif abs(a3(1))<eps
a3(1)=1e-6;
info=[info,'Replaced first element;'];
end
rtab=[rtab; a3, zeros(1,nrT-length(a3))];
vec1=vec2; vec2=a3;
end

melo by to byt v tomto packu:
http://www.mathworks.com/matlabcentr...?objectId=2302