% filt.m
% 實現一個LCC濾波器

function[sys,xo]=filt(t,x,u,flag,Lr,Cr,Co,R)
% 由于不是離散系統，輸出參數沒有包括ts
% 此外需要用戶輸入4個新的參數
% 電感Lr，電感電阻R，與電感串聯的電容Cr，輸出并聯電容Co

% 初始化
if flag==0
   sys=[3 0 1 2 1 1];
   % 3個狀態，x1為電感電流；x2為串聯電容電壓；x3為輸出電壓
   % 1個輸出，sys(1)為輸出電壓
   % 2個輸入，u(1)為輸入電壓，u(2)為負載電流
   % 初始狀態為零
   xo=[0 0 0];
elseif flag==1
   % 微分方程描述
      sys(1)=[u(1)-x(2)-x(3)-R*x(1)]/Lr;
      sys(2)=x(1)/Cr;
      sys(3)=[x(1)-u(2)]/Co;
   elseif flag==3
         sys=x(3);
      else
         sys=[];
         end