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.
Printable View
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.
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...
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)Citace:
Původně odeslal ghosthunter
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