function xdot=ch2example11statefun(t, x, flag, R,L,C,input,inputtimespan)
% 考慮RLC串聯環(huán)路的狀態(tài)方程函數
% 輸入: t 當前計算時刻,  flag 此處不用
% x為2#1矩陣, x(1)為電感上的電流; x(2)為電容電壓
% R: 電阻值
% L: 電感值
% C: 電容值
% input 為在時間點序列inputtimespan上給定輸入信號
xdot=zeros(2,1);        % 狀態(tài)變量矩陣初始化
A=[-R/L, -1/L; 1/C, 0]; % 狀態(tài)方程系數矩陣
B=[1/L; 0]; 
f_t=interp1(inputtimespan, input, t); 
% 利用插值來計算任意時刻的信號值
% 因為求解器調用本函數時時間變量t是任意連續(xù)值
xdot=A*x+B*f_t;         % 狀態(tài)方程