tphase.m

Matlab实现无刷直流电机控制系统建模仿真.zip

%This programme solves various formulaes for three phase half wave star rectifier
%This function runs independently as well as with "power_electronics"
%In order to execute this function type ''tphase'' on command prompt
function b=tphase(varargin)
fprintf('Welcome From Worldhitter')
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    fprintf('Provide The Following Data')
    Vm=input('\nEnter The Value of Peak Voltage(Volts) :')
    phi=input('Enter The Value Of Phase Angle(Degrees) :')
    Is1=input('Enter The Value Of Component Of Secondary Current(Amperes) :')
    Isp=input('Enter The Value Of Peak Current(Amperes) :')
    R=input('Enter The Value Of Resistence(Ohms): ')
     pause(1)
     %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    o=menu('Performance Parameters(Plz Select)' ,' Dc Voltage(Vdc)','Dc Current(Idc)','Root Mean Square Voltage(Vrms)','Root Mean Square Current(Irms)','Ac Voltage(Vac)','Dc Power(Pdc)','Ac power(Pac)','Secondary Voltage(Vs)','Peak Current(Im)','Seconadary current(Is)','Form Factor(FF)','Ripple Factor(RF)','Transformer Utilization Factor(TUF)','Displacement Factor(DF)','Harmonic Factor','Power Factor','Crest Factor(CF)');
    fprintf('Wait For 2 Seconds')
    pause(2)
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
    if(o == 1)
     fprintf('\nThe Magnitude Of Dc Voltage Is As Under')
      Vdc=3*sqrt(3)*Vm/6.28
      %^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
      Data=1:64;Data=(Data'*Data)/64;
      CreateStruct.WindowStyle='replace';
      CreateStruct.Interpreter='tex';
       h=msgbox('Vdc=3*sqrt(3)*Vm/6.28','Graph Equation','custom',Data,hsv(64),CreateStruct);
       pause(1)
       closereq
       msgbox('Thanx For Using My Programme ','','custom',Data,cool(64))
      pause(1)
      closereq   
      thanks(1,1,1,1,1,1) %Function call
    %$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
elseif(o == 2)
        fprintf('\nThe Magnitude Of Dc Current Is As Under')
        Idc=3*sqrt(3)*Vm/(6.28*R)
        fprintf('Proceeding Towards The Graph')
        %^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
        Data=1:64;Data=(Data'*Data)/64;
        CreateStruct.WindowStyle='replace';
        CreateStruct.Interpreter='tex';
        h=msgbox('Idc=3*sqrt(3)*Vm/(6.28*R)','Graph Equation','custom',Data,bone(64),CreateStruct);
        pause(1)
        closereq
        msgbox('Thanx For Using My Programme ','','custom',Data,cool(64))
        pause(1)
        closereq
        thanks(1,1,1,1,1,1) %Function call
       %$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
   elseif(o == 3)
        fprintf('\nThe Magnitude Of Rms Voltage Is As Under')
        Vrms=sqrt((3/pi)*Vm^2*(pi/3 +sqrt(3)/4 ))
        fprintf('Proceeding Towards The Graph')
        %^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
        Data=1:64;Data=(Data'*Data)/64;
        CreateStruct.WindowStyle='replace';
        CreateStruct.Interpreter='tex';
        h=msgbox('Vrms=sqrt((3/pi)*Vm^2*(pi/3 +sqrt(3)/4 ))','Graph Equation','custom',Data,winter(64),CreateStruct);
        pause(1)
        closereq
        msgbox('Thanx For Using My Programme ','','custom',Data,cool(64))
        pause(1)
        closereq   
       thanks(1,1,1,1,1,1) %Function call
       %$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
   elseif(o == 4)
    
        fprintf('\nThe Magnitude Of Rms Current Is As Under')
        Irms=sqrt((3/pi)*Vm^2*(pi/3 +sqrt(3)/4 ))/R;
        fprintf('Proceeding Towards The Graph')
        %^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
        Data=1:64;Data=(Data'*Data)/64;
        CreateStruct.WindowStyle='replace';
        CreateStruct.Interpreter='tex';
        h=msgbox('Irms=sqrt((3/pi)*Vm^2*(pi/3 +sqrt(3)/4 ))/R','Graph Equation','custom',Data,copper(64),CreateStruct);
       pause(1)
       closereq
       msgbox('Thanx For Using My Programme ','','custom',Data,cool(64))
       pause(1)
       closereq   
       thanks(1,1,1,1,1,1) %Function call
       %$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
   elseif(o == 5)
        fprintf('\nThe Magnitude Of Ac Voltage Is As Under')
        Vrms=sqrt((3/pi)*Vm^2*(pi/3 +sqrt(3)/4 ));
        Vdc=3*sqrt(3)*Vm/6.28;
        Vac=sqrt(Vrms^2 - Vdc^2)
         %^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
        Data=1:64;Data=(Data'*Data)/64;
        CreateStruct.WindowStyle='replace';
        CreateStruct.Interpreter='tex';
        h=msgbox('Vac=sqrt(Vrms^2 - Vdc^2)','Equation','custom',Data,hot(64),CreateStruct);
        pause(1)
        closereq
        msgbox('Thanx For Using My Programme ','','custom',Data,cool(64))
        pause(1)
        closereq   
        thanks(1,1,1,1,1,1) %Function call
        %$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
   elseif(o == 6)
       
       fprintf('\nThe Magnitude Of Dc Power Is As Under')
        Vdc=3*sqrt(3)*Vm/6.28;
        Idc=3*sqrt(3)*Vm/(6.28*R);
        Pdc=Vdc*Idc
        %^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
        Data=1:64;Data=(Data'*Data)/64;
        CreateStruct.WindowStyle='replace';
        CreateStruct.Interpreter='tex';
        h=msgbox('Pdc=Vdc*Idc','Equation','custom',Data,hot(64),CreateStruct);
        pause(1)
        closereq
        msgbox('Thanx For Using My Programme ','','custom',Data,cool(64))
        pause(1)
        closereq   
        thanks(1,1,1,1,1,1) %Function call
        %$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
     elseif(o == 7)
        fprintf('\nThe Magnitude Of Ac Power Pac As Under')
        Vrms=sqrt((3/pi)*Vm^2*(pi/3 +sqrt(3)/4 ));
        Irms=sqrt((3/pi)*Vm^2*(pi/3 +sqrt(3)/4 ))/R;
        Pac=Vrms*Irms
         %^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
        Data=1:64;Data=(Data'*Data)/64;
        CreateStruct.WindowStyle='replace';
        CreateStruct.Interpreter='tex';
        h=msgbox('Pac=Vrms*Irms','Equation','custom',Data,hot(64),CreateStruct);
        pause(1)
        closereq
        msgbox('Thanx For Using My Programme ','','custom',Data,cool(64))
        pause(1)
        closereq
        thanks(1,1,1,1,1,1) %Function call
         %$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
   elseif(o == 8)    
       fprintf('\nThe Magnitude Of Secondary Voltage Is As Under')
       Vs=Vm/sqrt(2)
         %^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
        Data=1:64;Data=(Data'*Data)/64;
        CreateStruct.WindowStyle='replace';
        CreateStruct.Interpreter='tex';
        h=msgbox('Vs=Vm/sqrt(2)','Equation','custom',Data,hot(64),CreateStruct);
        pause(1)
        closereq
        msgbox('Thanx For Using My Programme ','','custom',Data,cool(64))
        pause(1)
        closereq   
       thanks(1,1,1,1,1,1) %Function call       
       %$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
   elseif(o == 9)
       fprintf('\nThe Magnitude Of Peak(Maximum) Current Is As Under')
       Im=sqrt(3)*Vm/R
         %^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
        Data=1:64;Data=(Data'*Data)/64;
        CreateStruct.WindowStyle='replace';
        CreateStruct.Interpreter='tex';
        h=msgbox(' Im=sqrt(3)*Vm/R','Equation','custom',Data,hot(64),CreateStruct);
        pause(1)
        closereq
        msgbox('Thanx For Using My Programme ','','custom',Data,cool(64))
        pause(1)
        closereq   
       thanks(1,1,1,1,1,1) %Function call
       %$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
   elseif(o == 10)
       fprintf('\nThe Magnitude Of Seconadary Current Is As Under')
       Im=sqrt(3)*Vm/R;
       Is=.7804*Im
         %^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
        Data=1:64;Data=(Data'*Data)/64;
        CreateStruct.WindowStyle='replace';
        CreateStruct.Interpreter='tex';
        h=msgbox('  Is=.7804*Im)','Equation','custom',Data,hot(64),CreateStruct);
        pause(1)
        closereq
        msgbox('Thanx For Using My Programme ','','custom',Data,cool(64))
        pause(1)
        closereq   
       thanks(1,1,1,1,1,1) %Function call
      %$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
     elseif(o == 11)
      fprintf('\nThe Form Factor Is As Under')
      Vrms=sqrt((3/pi)*Vm^2*(pi/3 +sqrt(3)/4 ));
      Vdc=3*sqrt(3)*Vm/6.28;
      FF=Vrms/Vdc
     %^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
        Data=1:64;Data=(Data'*Data)/64;
        CreateStruct.WindowStyle='replace';
        CreateStruct.Interpreter='tex';
        h=msgbox('FF=Vrms/Vdc','Equation','custom',Data,summer(64),CreateStruct);
        pause(1)
        closereq
        msgbox('Thanx For Using My Programme ','','custom',Data,gray(64))
        pause(1)
        closereq
        thanks(1,1,1,1,1,1) %Function call
        %$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
    elseif(o == 12)
        fprintf('\nThe Ripple Factor Is As Under')
        Vrms=sqrt((3/pi)*Vm^2*(pi/3 +sqrt(3)/4 ));
        Vdc=3*sqrt(3)*Vm/6.28;
        RF=sqrt(Vrms^2 - Vdc^2)/Vdc
         %^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
        Data=1:64;Data=(Data'*Data)/64;
        CreateStruct.WindowStyle='replace';
        CreateStruct.Interpreter='tex';
        h=msgbox('RF=sqrt(Vrms^2 - Vdc^2)/Vdc','Equation','custom',Data,flag(64),CreateStruct);
        pause(1)
        closereq
        msgbox('Thanx For Using My Programme ','','custom',Data,spring(64))
        pause(1)
        closereq
        thanks(1,1,1,1,1,1) %Function call
        %$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
    elseif(o == 13)
        fprintf('\nThe Transformer Utilization  Factor Is As Under')
        Idc=3*sqrt(3)*Vm/(6.28*R);
        Vdc=3*sqrt(3)*Vm/6.28;
        Pdc=Vdc*Idc;
        Vs=Vm/sqrt(2);
        TUF=Pdc/(2*Vs*Is)
         %^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
        Data=1:64;Data=(Data'*Data)/64;
        CreateStruct.WindowStyle='replace';
        CreateStruct.Interpreter='tex';
        h=msgbox('TUF=Pdc/(2*Vs*Is)','Equation','custom',Data,autumn(64),CreateStruct);
        pause(1)
        closereq
        msgbox('Thanx For Using My Programme ','','custom',Data,cool(64))
        pause(1)
        closereq
        thanks(1,1,1,1,1,1) %Function call
        %$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
    elseif(o == 14)
        fprintf('\nThe Displacement  Factor Is As Under')
        DF=cos(phi)
         %^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
        Data=1:64;Data=(Data'*Data)/64;
        CreateStruct.WindowStyle='replace';
        CreateStruct.Interpreter='tex';
        h=msgbox('DF=cos(phi)','Equation','custom',Data,lines(64),CreateStruct);
        pause(1)
        closereq
        msgbox('Thanx For Using My Programme ','','custom',Data,cool(64))
        pause(1)
        closereq
        thanks(1,1,1,1,1,1) %Function call
     %$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
    elseif(o == 15)
      fprintf('\nThe Harmonic Factor Is As Under')
      Im=sqrt(3)*Vm/R;
      Is=.7804*Im
      HF=sqrt((Is/Is1)^2 -1)
        
     %^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
        Data=1:64;Data=(Data'*Data)/64;
        CreateStruct.WindowStyle='replace';
        CreateStruct.Interpreter='tex';
        h=msgbox('HF=sqrt((Is/Is1)^2 -1)','Equation','custom',Data,prism(64),CreateStruct);
        pause(1)
        closereq
        msgbox('Thanx For Using My Programme ','','custom',Data,cool(64))
        pause(1)
        closereq 
        thanks(1,1,1,1,1,1) %Function call
        %$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
    elseif(o == 16)
        fprintf('\nThe  Power Factor Is As Under')
        Im=sqrt(3)*Vm/R;
       Is=.7804*Im
        PF=Is1/Is*cos(phi)
        %^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
        Data=1:64;Data=(Data'*Data)/64;
        CreateStruct.WindowStyle='replace';
        CreateStruct.Interpreter='tex';
        h=msgbox('PF=Is1/Is*cos(phi))','Equation','custom',Data,white(128),CreateStruct);
        pause(1)
        closereq
        msgbox('Thanx For Using My Programme ','','custom',Data,cool(64))
        pause(1)
        closereq
        thanks(1,1,1,1,1,1) %Function call
        %$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
    elseif(o == 17)
        fprintf('\nThe  Crest Factor Is As Under')
        Im=sqrt(3)*Vm/R;
        Is=.7804*Im
        CF=Isp/Is
        %^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
        Data=1:64;Data=(Data'*Data)/64;
        CreateStruct.WindowStyle='replace';
        CreateStruct.Interpreter='tex';
        h=msgbox('CF=Isp/Is)','Equation','custom',Data,hsv(128),CreateStruct);
        pause(1)
        closereq
        msgbox('Thanx For Using My Programme ','','custom',Data,cool(64))
        pause(1)
        closereq 
        thanks(1,1,1,1,1,1) %Function call
    end