分布式驱动电动汽车LQR控制与扭矩分配技术

1. 分布式驱动电动汽车控制架构解析

在分布式驱动电动汽车的稳定性控制系统中，我们采用分层控制架构来实现精确的横摆力矩控制。这种架构将复杂的控制问题分解为可管理的层级，每层专注于特定的控制目标。

1.1 上层控制：LQR直接横摆力矩控制

上层控制器采用线性二次型调节器（LQR）来计算所需的直接横摆力矩。其核心是基于车辆的二自由度动力学模型：

code复制m(v̇y + vxγ) = Fyf + Fyr
Izγ̇ = Fyf lf - Fyr lr + Mz

其中m为整车质量，vx和vy分别为纵向和侧向速度，γ为横摆角速度，Fyf和Fyr为前后轴侧向力，lf和lr为质心到前后轴距离，Iz为绕z轴的转动惯量，Mz为直接横摆力矩。

LQR控制器的设计关键在于状态方程构建和权重矩阵选择。我们选取的状态变量为：

code复制x = [β γ]T

其中β为质心侧偏角。控制输入为：

code复制u = [δ Mz]T

δ为前轮转角。状态方程可表示为：

code复制ẋ = Ax + Bu

1.2 下层控制：扭矩分配策略

下层控制器负责将上层计算的总横摆力矩和纵向力需求分配到四个驱动电机。我们采用基于二次规划的最优分配算法，目标是最小化轮胎附着利用率：

code复制min Σ(Fxi^2/(μiFzi)^2)

约束条件包括：

• 总纵向力等于需求：ΣFxi = Fx_des
• 总横摆力矩等于需求：Σ(-Fyi·xi + Fxi·yi) = Mz_des
• 各轮扭矩限制：Fxi_min ≤ Fxi ≤ Fxi_max

其中μi为各轮附着系数，Fzi为垂向载荷，xi和yi为轮胎位置坐标。

2. 系统实现与仿真验证

2.1 MATLAB/Simulink实现要点

在Simulink中实现该控制系统时，有几个关键模块需要特别注意：

1. 车辆模型更新模块：

matlab复制function [A, B] = update_model(vx, delta)
    % 参数定义
    m = 1500; lf = 1.2; lr = 1.5; 
    Iz = 2500; Cf = 80000; Cr = 100000;
    
    % 状态矩阵计算
    a11 = -(Cf + Cr)/(m*vx);
    a12 = -1 - (Cf*lf - Cr*lr)/(m*vx^2);
    a21 = -(Cf*lf - Cr*lr)/Iz;
    a22 = -(Cf*lf^2 + Cr*lr^2)/(Iz*vx);
    
    A = [a11 a12; a21 a22];
    B = [Cf/(m*vx) 0; Cf*lf/Iz 1/Iz];
end

2. LQR控制器模块：

matlab复制function [Mz_des] = lqr_controller(x_