% rect.m
% 實現一個單相不控的全橋整流器

function[sys,xo]=rect(t,x,u,flag,R1,C,R)
% 由于不是離散系統，輸出參數沒有包括ts
% 此外需要用戶輸入3個新的參數
% 交流母線電阻R1，直流濾波電容C，負載電阻R
% 根據不控全橋整流電路的原理，列寫微分方程

a1=-(1/R1+1/R)/C;
a2=-1/R/C;
b1=1/R1/C;
b2=0;
if flag==0
   sys=[1 0 2 1 0 0];
   % 1個狀態，x為輸出直流電壓（直流濾波電容電壓）
   % 2個輸出，sys(1)為輸入交流母線電流；sys(2)為輸出直流電壓
   % 1個輸入，u為交流輸入電壓
   % 初始狀態為零
   xo=0;
elseif flag==1
      if abs(u)>=x
         sys=a1*x+b1*abs(u);
      else
         sys=a2*x;
         end
      elseif flag==3
            if abs(u)>=x
               sys(1)=sign(u)*(abs(u)-x)/R1;
               sys(2)=x;
            else
               sys(1)=0;
               sys(2)=x;
            end
         else
            sys=[];
            end