ch2example11statefun.m

来自「清华大学出版社 邵玉斌编写的《通信系统建模与仿真实例分析》一书的所有MATLAB」· M 代码 · 共 15 行

function xdot=ch2example11statefun(t, x, flag, R,L,C,input,inputtimespan)
% 考虑RLC串联环路的状态方程函数
% 输入: t 当前计算时刻,  flag 此处不用
% x为2#1矩阵, x(1)为电感上的电流; x(2)为电容电压
% R: 电阻值
% L: 电感值
% C: 电容值
% input 为在时间点序列inputtimespan上给定输入信号
xdot=zeros(2,1);        % 状态变量矩阵初始化
A=[-R/L, -1/L; 1/C, 0]; % 状态方程系数矩阵
B=[1/L; 0]; 
f_t=interp1(inputtimespan, input, t); 
% 利用插值来计算任意时刻的信号值
% 因为求解器调用本函数时时间变量t是任意连续值
xdot=A*x+B*f_t;         % 状态方程