„Szabályozástechnika - Soros kompenzátorok tervezése” változatai közötti eltérés

David14 (vitalap | szerkesztései)
David14 (vitalap | szerkesztései)
281. sor: 281. sor:
% Ez NEM egy általános célú függvény lesz, mivel a T1, T2, Phit és umax  
% Ez NEM egy általános célú függvény lesz, mivel a T1, T2, Phit és umax  
% változókat nem paraméterként kapja meg, hanem az egyszerűség kedvéért
% változókat nem paraméterként kapja meg, hanem az egyszerűség kedvéért
% fixen beleírtam a kódba!!!  
% fixen beleírtam a kódba!!!
 
% I. egyenlet kifejtése (pólus zérus kiejtést kapásból elvégezve):
 
%                                          Ap                  Aplant
% W0(j*wc)= Wc(j*wc)*Wp(j*wc)= --------------------------- * -------------
%                              Ti * (j*wc) * (1 + j*wc*Tc)  (1 + j*wc*T3)
 
%                              Ap * Aplant
% |W0(j*wc)|= ------------------------------------------------
%              Ti * wc * sqrt{ (1+wc^2*Tc^2) * (1+wc^2*T3^2) }
 
% Továbbá tudjuk az 1. egyenletből, hogy Ti = T1+T2-Tc
 


function f=myPID(x)
function f=myPID(x)
296. sor: 309. sor:
umax=10;
umax=10;


f1=Ap*Aplant/(T1+T2-Tc)/(wc*sqrt((1+wc^2*T3^2)*(1+wc^2*Tc^2)))-1;
f1=Ap*Aplant/(T1+T2-Tc)/(wc*sqrt((1+wc^2*Tc^2)*(1+wc^2*T3^2)))-1;
f2=pi-pi/2-atan(wc*T3)-atan(wc*Tc)-Phit;
f2=pi-pi/2-atan(wc*T3)-atan(wc*Tc)-Phit;
f3=Ap*T1*T2/(T1+T2-Tc)/Tc-umax;
f3=Ap*T1*T2/(T1+T2-Tc)/Tc-umax;