% filt.m
% 實(shí)現(xiàn)一個(gè)LCC濾波器

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

% 初始化
if flag==0
   sys=[3 0 1 2 1 1];
   % 3個(gè)狀態(tài)，x1為電感電流；x2為串聯(lián)電容電壓；x3為輸出電壓
   % 1個(gè)輸出，sys(1)為輸出電壓
   % 2個(gè)輸入，u(1)為輸入電壓，u(2)為負(fù)載電流
   % 初始狀態(tài)為零
   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