منحنيات التحكم في سرعة الموتور الحثي عبر octave

قمت بإستخدام octave لعمل تقرير عن التحكم في سرعة الموتور الحثي عن طريق تغيير الفولت ، تغيير التردد ، تغيير المقاومة الخارجية ، تغيير الفولت والتردد معا بنسبة ثابتة

 
 
function t = Torque(V1, R1, R2, S, X1, X2, Ws)
t=(3*V1*V1*R2)/(S*Ws*( (R1+ (R2/S))^2 +(X1+X2)^2))
endfunction
 
 
global R1 = 0.641;               
global X1 = 0.750;                
global R2 = 0.300;                
global X2 = 0.500; 
global Ns = 1800;             
global Ws = 188.5;  
 
function voltagecontrol
 
 
        global R1
        global R2
        global X1
        global X2 
        global Ns
        global Ws
        
        S=1:-0.01:0.01;
 
        Nm=(1-S)*Ns;
        V1=(0:50:400);
        figure(1)
        for vi=1:8;
                for i=1:100
                        t(i)=Torque(V1(vi), R1, R2, S(i), X1, X2, Ws);
                end
                plot(Nm, t, 'k-', 'LineWidth', 2.0);
                hold on;
 
        end
        xlabel('NM');
        ylabel('T');
        title('Variable voltage control IM T-S c/cs');
 
endfunction 
 
 
function raddcontrol
         
 
 
        global R1
        global R2
        global X1
        global X2 
        global Ns
        global Ws
        V1 = 400 ;
 
        S=1:-0.01:0.01;
 
        Nm=(1-S)*Ns;
        radd=(0:.1:0.5);
        figure(2)
        for i=1:100
                t(i)=Torque(V1, R1, R2 ,S(i), X1, X2, Ws);
        end
        plot(Nm, t, 'k-', 'LineWidth', 1.0);
        hold on;
 
        for ri=1:5
                for i=1:100
                        tnew(i)=Torque(V1, R1, R2+radd(ri) ,S(i), X1, X2, Ws);
                end
                plot(Nm, tnew, 'k-', 'LineWidth', 1.0);
                hold on;
        end
        xlabel('NM');
        ylabel('T');
        title('Radd control IM T-S c/cs');
endfunction
 
 
function frequencycontrol
 
 
        global R1
        global R2
        global X1
        global X2 
        global Ns
        global Ws
        V1 = 380;
        F1=50;
 
        S=1:-0.01:0.01;
 
 
        s1=0.05
        s2=0.03
        Nm=(1-S)*Ns;
        
        for i=1:100
                t(i)=Torque(V1, R1, R2, S(i), X1, X2, Ws);
        end
        
        F=[50;60;70];
        figure(3)
 
                
        plot(Nm, t, 'k-', 'LineWidth', 2.0);
        hold on;
        for fi=1:3
                tnew=(t*F(fi)*s2)/(F1*s1);
                %new Nm
                Nm2=Nm*(sqrt(t(50)/tnew(50)));
                plot(Nm2, tnew,'k-', 'LineWidth', 2.0);
                hold on;
        end
        xlabel('NM');
        ylabel('T');
        title('Frequency control IM T-S c/cs');
 
endfunction 
 
 
 
function vfcontrol
 
 
 
        global R1
        global R2
        global X1
        global X2 
        global Ns
        global Ws
        Vph=200
        F1=50
 
        S=1:-0.01:0.01;
 
        Nm=(1-S)*Ns;
        %const ratio
        V=[200; 204; 208; 212;216;220];
        F=[50;51;52;53;54;55];
        
 
        s1=0.05
        s2=0.03
        figure(4)
        for i=1:100
                t(i)=Torque(Vph, R1, R2, S(i), X1, X2, Ws);
        end
        %plot(Nm, t, 'k-', 'LineWidth', 2.0);
        %hold on;
 
        
        for c=1:5
                vnew=V(c);
                fnew=F(c);
                Ns2=Ns*fnew/F1;
                Nm2=(1-S)*Ns2;
                tnew=t*vnew^2*F1*s1/(Vph^2*fnew*s2);
                %new NM
                plot(Nm2, tnew, 'k-', 'LineWidth', 2.0);
                hold on;    
        end
        xlabel('NM');
        ylabel('T');
        title('V/F control IM T-S c/cs');
endfunction