动态电路的状态变量分析

动态电路的状态变量分析StateVariableAnalysisofDynamicCircuitsContents课程目录状态变量与状态方程分析，从理论基础到工程实践的完整知识脉络。01状态变量基本概念02状态方程的建立03状态方程的求解方法04工程应用与案例分析CHAPTER01状态变量基本概念从系统状态到状态变量的理论演进CIRCUITDYNAMICS动态系统的定义与特征动态系统是指系统行为不仅取决于当前输入，还依赖于系统内部储能状态的历史积累。在电路中，电容和电感作为储能元件赋予系统"记忆"能力，使电路响应呈现时间延续性和状态依赖性，这是动态电路区别于纯电阻电路的本质特征。RLC电路实验板·电阻、电感、电容储能元件记忆性·系统当前输出由当前输入和内部状态共同决定，与纯电阻电路的即时响应有本质区别储能元件·电容C储存电场能、电感L储存磁场能，是动态行为的物理来源，能量不能突变数学描述·动态系统需要微分方程或差分方程，方程阶数等于独立储能元件的个数典型系统·RLC电路、机械振动、热传导等，都遵循能量守恒和状态连续演化规律StateDefinition状态的定义与物理意义系统的"状态"是能够完全确定系统未来演化行为的最小信息集合——用最少的变量捕捉系统的全部"记忆"，实现从历史信息到未来预测的完整映射。01严格定义初始时刻t₀的状态x(t₀)与输入u(t)已知时，系统在[t₀,t)上任何瞬时的行为都被唯一确定。02物理直觉质点：位置与速度；理想气体：温度、压力、体积；电路：电容电压与电感电流。03充分必要性状态变量既不冗余——个数最少，也不遗漏——能完全描述系统行为。04内部特性突破状态是系统内部特性的集中体现，突破了传统输入输出法只关注外部行为的局限。相空间中的系统状态轨迹·状态变量的几何表达DynamicCircuitAnalysis电路中的状态变量选择在动态电路分析中，状态变量通常选择独立的电容电压u_C和电感电流i_L。选择依据是它们直接对应储能元件的能量状态，且变量之间线性无关。但需注意特殊情况：纯电容回路、纯电感割集会导致独立状态变量数目少于储能元件总数。标准选择电容电压u_C和电感电流i_L，直接表征电场能量(½Cu²)和磁场能量(½Li²)的存储状态。这两个变量能够完整描述动态电路的瞬态行为。½Cu²·½Li²独立性要求状态变量必须线性无关，任意一个不能由其他变量的线性组合表示，确保变量组的最小性。这是状态空间方程可解性的基本前提。线性无关纯电容回路由KVL约束，回路中各电容电压不独立，独立状态变量数减少一个。此时需选择回路中部分电容电压作为独立变量。KVL约束纯电感割集由KCL约束，割集中各电感电流不独立，独立状态变量数同样减少一个。需选择割集中部分电感电流作为独立变量。KCL约束STATEVARIABLES·CASESTUDY状态变量确定实例：RLC串联电路以RLC串联电路为例，系统包含一个电感和一个电容两个独立储能元件。选择电感电流iL和电容电压uC作为状态变量，可以完全描述电路的动态行为。RLC串联电路实验测量场景01电路结构：电压源us(t)与R、L、C串联，包含两个独立储能元件，需要两个状态变量描述02状态变量选择：x₁=iL（电感电流），x₂=uC（电容电压），二者线性无关且数目最少03完整性验证：uR=R·iL，i=iL，iC=C·duC/dt=iL，所有响应均可由状态变量表示04状态数确定：独立状态变量数=储能元件数−纯电容回路数−纯电感割集数=2−0−0=2方法论对比状态变量分析法vs传统输入输出法状态变量分析法通过描述系统内部状态演化，突破了传统输入输出法只能处理单输入单输出线性时不变系统的局限。该方法能够处理多输入多输出、非线性、时变等复杂系统，且天然适合计算机数值求解，是现代控制理论和复杂电路分析的基础工具。传统输入输出法Input–Output01仅关注输入与输出的外部关系，将系统视为"黑箱"，不揭示内部结构和状态演化02适用于单输入单输出(SISO)线性时不变系统，对MIMO、非线性、时变系统几乎无能为力03数学工具为传递函数或频率响应，物理直觉清晰但计算复杂度随系统阶数急剧增加状态变量分析法State-Space01通过状态变量描述系统内部特性，全面反映系统结构和能量流动，实现"白箱"分析02天然适用于多输入多输出(MIMO)、非线性和时变系统，应用范围显著扩展03数学形式为一阶微分方程组，易于转化为矩阵形式并借助计算机进行数值求解STATESPACEANALYSIS状态空间与状态轨迹状态空间以n个状态变量为坐标轴，系统状态随时间描绘出轨迹，将微分方程转化为几何图像。01状态空间定义—由n个状态变量x₁,x₂,…,xₙ构成的n维向量空间，记为状态向量x=[x₁,x₂,…,xₙ]ᵀ02状态轨迹—系统状态随时间演化在状态空间中描绘的曲线，反映从初始状态到终态的完整动态过程03几何直觉—二维状态空间可直接可视化（如RLC电路i_L–u_C平面），高维空间需借助投影或数值方法04平衡点分析—状态导数为零的点对应系统稳态工作点，是稳定性分析的关键状态空间相平面图：系统状态轨迹的几何可视化Summary核心概念小结状态变量分析法的核心概念体系包括动态系统、状态、状态变量、状态空间和状态轨迹等相互关联的概念。理解这些概念的物理意义和数学定义，是掌握状态方程建立和求解方法的前提。状态变量的完全性、最小性和独立性是贯穿整个分析过程的三条基本原则。状态变量分析核心概念一览概念定义电路实例动态系统输出取决于当前输入和内部历史状态的系统含电容或电感的电路状态完全确定系统未来行为的最小信息集合某时刻所有u_C和i_L的值状态变量描述状态的线性无关最小变量组独立的u_C和i_L状态空间以状态变量为坐标的n维向量空间i_L-u_C平面（二阶电路）状态轨迹状态点随时间在状态空间中的运动路径阻尼振荡的螺旋线五个核心概念构成状态变量分析的完整理论框架Chapter02状态方程的建立从物理电路到标准矩阵形式的系统化方法STATESPACE·EQUATIONS状态方程的标准形式状态空间表达式由状态方程和输出方程组成，标准矩阵形式为ẋ=Ax+Bu和y=Cx+Du，四组矩阵分别刻画系统动态、输入耦合与输出映射。状态方程ẋ(t)=Ax(t)+Bu(t)，描述状态变量一阶导数与当前状态和输入的关系，是n个一阶微分方程的矩阵形式。ẋ=Ax+Bu输出方程y(t)=Cx(t)+Du(t)，描述系统输出与状态和输入的代数关系，通常不含微分运算，是系统可观测部分的数学表达。y=Cx+Du矩阵维度A:n×n，B:n×m，C:p×n，D:p×m；n为状态数，m为输入数，p为输出数。矩阵维度匹配保证运算合法性。n×n·p×m物理意义A矩阵对角元素反映各状态自衰减特性，非对角元素反映状态间的耦合关系，B矩阵表征输入对状态的影响程度。自衰减·耦合METHODOLOGY状态方程建立的步骤状态方程的建立遵循系统化的四步流程：选择状态变量、列写基本方程、代入元件关系、整理标准形式。关键在于正确处理储能元件的伏安关系，并通过代数消元将所有非状态变量用状态变量和输入表示，最终得到仅含状态变量及其导数的标准一阶微分方程组。01选择状态变量确定电路中独立的电容电压uC和电感电流iL，排除纯电容回路和纯电感割集带来的冗余变量。uC·iL02列写基本方程对含电容节点列KCL，对含电感回路列KVL，同时列写必要的辅助方程（如其他节点或回路方程）。KCL·KVL03代入元件关系利用iC=C(duC/dt)和uL=L(diL/dt)将电容电流和电感电压转化为状态变量的导数。iC·uL04整理标准形式通过代数消元去除非状态变量，将方程整理为ẋ=Ax+Bu的矩阵形式，提取A、B矩阵。ẋ=Ax+BuSTATEEQUATION·DERIVATION实例推导：RC串联电路的状态方程以RC串联电路为例，通过完整的四步推导建立状态方程。该一阶系统只有一个状态变量uC，最终得到的状态方程清晰展示了系统矩阵与输入矩阵的物理意义。01电路描述电压源us(t)与电阻R、电容C串联，一个储能元件，阶数n=102选择状态变量x=uC（电容电压），输入u=us（电压源），输出y=uC03列方程KVL给出us=Ri+uC，元件关系i=C(duC/dt)，联立得RC(duC/dt)+uC=us04标准形式ẋ=(−1/RC)x+(1/RC)u，即A=−1/RC，B=1/RC，C=1，D=0，τ=RCRC电路在面包板上的实际搭建State-SpaceModeling实例推导：RLC串联电路的状态方程RLC串联电路作为典型的二阶系统，需要两个状态变量iL和uC来描述，状态方程呈现2×2系统矩阵，R/L项反映阻尼作用，1/L和1/C项反映能量交换耦合。01状态变量选择x₁=iL（电感电流），x₂=uC（电容电压），输入u=us，阶数n=2n=202KVL方程us=RiL+L(diL/dt)+uC，整理得diL/dt=−(R/L)iL−(1/L)uC+(1/L)usKVL03电容伏安关系iC=C(duC/dt)=iL，整理得duC/dt=(1/C)iLiC=iL04矩阵形式A=[−R/L,−1/L;1/C,0]，B=[1/L;0]，特征方程s²+(R/L)s+1/(LC)=02×2STATE-SPACEMODELING复杂电路的状态方程建立受控源、多储能元件及互感耦合是复杂电路的三大特殊因素，系统化处理是保证状态方程正确性的关键。含受控源电路将受控源的控制量（电压或电流）用状态变量和输入表示后代入方程受控源可能在A矩阵中引入额外的非对角元素，反映状态间的受控耦合受控耦合多储能元件电路独立状态变量数=储能元件总数−纯电容回路数−纯电感割集数需逐一检查拓扑约束，确定哪些电容电压或电感电流可被其他变量表示拓扑约束含互感电路互感电压uM=M(di/dt)在状态方程中引入耦合导数项，需额外代数处理可将互感等效为受控源，或直接列写含互感的KVL方程后消元交叉耦合State-SpaceModel输出方程的建立输出方程y=Cx+Du描述系统输出与状态和输入的代数关系，是状态空间表达式的必要组成部分。01📐输出方程形式y(t)=Cx(t)+Du(t)是纯代数方程，不含状态变量的导数项，直接反映输出与状态、输入的瞬时关系。02🔧C矩阵确定将输出y用状态变量线性表示，各项系数构成C矩阵行向量。例如y=u_R=R·i_L，则C=[R,0]。03⚡D矩阵确定若输出直接包含输入分量则D≠0；大多数物理系统D=0，表示输出不依赖输入的瞬时直通路径。04📊多输出情况p个输出对应p行C矩阵和p行D矩阵，每个输出独立确定其对应行向量，各行之间互不影响。STATESPACE状态变量的非唯一性与相似变换状态变量不唯一，经非奇异变换z=Tx关联；矩阵元素改变而特征值与传递函数不变。非唯一性原因状态空间中任意非奇异线性变换z=Tx都可产生新的有效状态变量组非奇异变换相似变换公式A'=TAT⁻¹，B'=TB，C'=CT⁻¹，D'=D，其中T为任意n×n非奇异矩阵TAT⁻¹不变量特征多项式、特征值、传递函数矩阵、可控性与可观性在相似变换下保持不变特征值不变实用意义选择合适的T将A化为对角型或约当标准型，简化后续分析和计算约当标准型CHAPTER03状态方程的求解方法时域解析与复频域变换的双重求解路径LinearTime-InvariantSystems状态方程解的结构ẋ=Ax+Bu的通解由零输入响应与零状态响应叠加：x(t)=e^(At)x(0)+∫₀ᵗe^(A(t-τ))Bu(τ)dτ，体现线性系统的叠加原理。01零输入响应x_zi(t)=Φ(t)·x(0)，初始状态在无外部激励下的自由演化过程Φ(t)·x(0)02零状态响应x_zs(t)=∫₀ᵗΦ(t−τ)Bu(τ)dτ，外部输入通过B矩阵对状态的累积驱动∫ΦBudτ03完整解x(t)=x_zi(t)+x_zs(t)，符合叠加原理，两部分可独立计算后相加叠加原理04输出响应将x(t)代入输出方程y=Cx+Du即可得到系统的完整输出响应y=Cx+DuLinearSystemTheory状态转移矩阵的性质状态转移矩阵Φ(t)=eAt是描述系统自由演化规律的核心数学对象，完全由系统矩阵A的特征值决定，是连接系统结构与动态响应的桥梁。PROPERTY01初值性质Φ(0)=eA·0=I（单位矩阵），零时刻状态保持不变，是解的初始条件验证依据Φ(0)=IPROPERTY02半群性质Φ(t₁+t₂)=Φ(t₁)·Φ(t₂)，从0到t₁+t₂的转移等于先到t₁再到t₂的复合Φ(t₁+t₂)PROPERTY03可逆性质Φ⁻¹(t)=Φ(-t)=e-At，状态转移过程可逆，可从当前状态回溯历史Φ⁻¹(t)PROPERTY04导数性质dΦ(t)/dt=AΦ(t)=Φ(t)A，转移矩阵导数等于系统矩阵左乘或右乘转移矩阵dΦ/dt=AΦSTATESPACEANALYSIS矩阵指数e(At)的计算方法矩阵指数e(At)的计算是状态方程时域求解的核心步骤。主要方法包括泰勒级数展开法、拉普拉斯变换法和对角化方法。其中拉氏变换法e(At)=L⁻¹{(sI-A)⁻¹}最为通用，适用于任意阶系统；对角化方法计算最简洁但要求A可对角化；级数法主要用于理论推导和数值近似。01拉普拉斯变换法计算(sI-A)⁻¹，即预解矩阵，然后对每个元素做拉氏反变换得到e(At)适用性最广，无需判断A是否可对角化，是解析求解的首选方法(sI-A)⁻¹02对角化方法若A有n个不同特征值，A=PΛP⁻¹，则e(At)=P·diag(eλ₁t,...,eλₙt)·P⁻¹计算简洁直观，但需要A可对角化；有重特征值时需改用约当标准型PΛP⁻¹03凯莱-哈密顿定理法利用特征方程将e(At)表示为I,A,A²,...,Aⁿ⁻¹的有限多项式，系数由特征值确定特别适合手工计算低阶系统（n≤3），避免了矩阵求逆的复杂运算n≤3LaplaceTransformMethod复频域求解方法对状态方程取拉普拉斯变换，将微分方程转化为代数方程：sX(s)−x(0)=AX(s)+BU(s)。整理得X(s)=(sI−A)⁻¹[x(0)+BU(s)]，其中(sI−A)⁻¹为预解矩阵。变换过程对ẋ=Ax+Bu两边取拉氏变换，利用L{ẋ}=sX(s)−x(0)，将微分方程转化为代数方程sX(s)−x(0)求解公式X(s)=(sI−A)⁻¹x(0)+(sI−A)⁻¹BU(s)，分解为零输入响应和零状态响应两部分零输入+零状态预解矩阵Φ(s)=(sI−A)⁻¹，其极点det(sI−A)=0即为系统特征方程，决定模态与稳定性(sI−A)⁻¹与时域法的联系L⁻¹{(sI−A)⁻¹}=eAt=Φ(t)，复频域与时域两种方法在数学上完全等价eAt=Φ(t)CASESTUDY求解实例：RC电路阶跃响应以RC串联电路的阶跃响应为例，完整演示拉普拉斯变换法求解状态方程的过程。该结果与经典三要素法一致，验证了状态变量法的正确性。01问题设定RC串联电路，状态方程ẋ=-(1/RC)x+(1/RC)u，输入u=U₀（阶跃），初始条件x(0)=002拉氏变换sX(s)-0=(-1/RC)X(s)+(1/RC)(U₀/s)，整理得X(s)=U₀/[RCs(s+1/RC)]03部分分式展开X(s)=U₀[1/s-1/(s+1/RC)]，分解为稳态项和暂态项04反变换结果x(t)=uC(t)=U₀(1-e-t/RC)，t→∞时趋于稳态值U₀，时间常数τ=RCStabilityAnalysis基于状态方程的稳定性分析线性时不变系统的稳定性完全由系统矩阵A的特征值分布决定。特征值实部全为负则渐近稳定，存在正实部则不稳定，纯虚根对应临界稳定的等幅振荡。渐近稳定A的所有特征值λᵢ满足Re(λᵢ)<0，状态转移矩阵e^(At)→0（t→∞），系统自由响应衰减至零。Re(λ)<0不稳定存在至少一个特征值满足Re(λᵢ)>0，对应e^(λᵢt)模态发散，系统响应无限增长。Re(λ)>0临界稳定存在Re(λᵢ)=0的纯虚根且无重根，对应等幅振荡模态，如LC无阻尼振荡电路。Re(λ)=0物理对应RLC电路中R>0保证特征值负实部（阻尼耗散），R=0时特征值在虚轴上（无阻尼振荡）。RLC电路ControlTheory传递函数矩阵与状态空间的转换H(s)=C(sI-A)⁻¹B+D联通状态空间与传递函数，极点即A的特征值，零点由系统矩阵共同决定。转换公式H(s)=C(sI-A)⁻¹B+D，将状态空间四元组(A,B,C,D)转化为频域传递函数矩阵(A,B,C,D)→H(s)极点确定传递函数极点由det(sI-A)=0给出，即矩阵A的特征值，决定系统自然频率det(sI-A)=0SISO特殊情况单输入单输出时H(s)为标量，分子阶次≤分母阶次，相对阶次≥0ScalarH(s)零极点对消状态空间实现可能含不可控或不可观模态，使传递函数阶次低于状态变量数OrderReductionChapter04工程应用与案例分析从理论方法到工程实践的应用拓展CaseStudy应用案例：电力系统暂态分析电力系统暂态分析是状态变量法的典型工程应用。通过将发电机转子运动方程、励磁系统、调速器等建模为状态方程组，可以分析系统在扰动后的动态响应和稳定性。电力系统变电站与输电线路实景01状态变量选择：发电机转子角度δ、转速ω、暂态电动势E'q等，每台发电机贡献3–7个状态变量02非线性特性：转子运动方程d²δ/dt²包含sin(δ)非线性项，需在运行点线性化后建立增量状态方程03暂态稳定性：求解故障后系统状态轨迹，判断发电机能否保持同步运行，即δ是否趋于新稳态值04数值求解：大规模系统状态方程阶数可达数百，采用Runge-Kutta等数值积分方法在计算机上求解APPLICATION·ACTIVEFILTERS应用案例：状态变量有源滤波器状态变量滤波器是一种经典的有源滤波器拓扑，使用三个运放同时提供低通、带通和高通三种输出。其设计过程直接基于状态空间模型：通过选择电路参数配置系统矩阵A的特征值，实现期望的中心频率和品质因数。该滤波器的优势在于三个输出独立可调，参数灵敏度高，广泛应用于音频处理和通信系统。01电路结构：三个运放构成积分-求和环路，两个积分器对应二阶系统的两个状态变量02三种输出：低通输出LP、带通输出BP、高通输出HP同时可用，无需额外电路03参数设计：中心频率ω₀=1/(RC)，品质因数Q由反馈电阻比决定，三者可独立调节04状态空间实现：直接将传递函数H(s)转化为状态方程ẋ=Ax+Bu的电路实现，理论与电路一一对应含有运算放大器的有源滤波器电路板State-SpaceModeling应用案例：DC-DC变换器建模DC-DC变换器通过开关管周期性切换实现电压变换，经状态空间平均法得到统一模型，是设计反馈控制器、确保输出稳定的基础。Buck变换器状态方程ONL(diL/dt)=Vin−uC，C(duC/dt)=iL−uC/ROFFL(diL/dt)=−uC，C(duC/dt)=iL−uC/R2拓扑状态空间平均法WEIGHTED两种拓扑的A、B矩阵按占空比d加权平均：Aavg=dA₁+(1−d)A₂LINEAR平均模型为连续线性系统，可用经典控制理论设计补偿器Aavg小信号建模LINEARIZE在稳态工作点附近对平均模型做小信号扰动线性化，得到小信号状态方程TRANSFER推导控制到输出的传递函数，为PID或状态反馈控制器设计提供依据PID数值方法·COMPUTATIONALMETHODS状态方程的数值求解方法对于无法获得解析解的复杂状态方程，数值求解是工程主要手段。欧拉法简单但精度低，龙格-库塔法精度高但计算量大，MATLAB的ode45是标准工具。欧拉法（一阶）递推公式x(t+h)=x(t)+h·f(x,t)，计算简单但截断误差为O(h²)，需要很小步长适用于快速原型验证、实时仿真等对精度要求不高但对速度敏感的场景O(h²)四阶龙格-库塔法每步计算4次函数值并加权平均，截断误差O(h⁵)，精度远高于欧拉法工程标准方法，MATLAB的ode45采用其变步长版本，自动调节步长O(h⁵)刚性系统处理A矩阵特征值实部差异悬殊时系统呈刚性，显式方法需极小步长才能稳定采用隐式方法（如ode15s）或专用刚性求解器，牺牲单步效率换取稳定性隐式求解ComparativeAnalysis状态变量分析法的优势与局限状态变量分析法在描述系统内部特性、处理多变量复杂系统方面具有不可替代的优势，同时也存在建模复杂度高、物理直觉不如频域法等