LADRC_C代码

#ifndef _ADRC_H
#define _ADRC_H
#include “stm32f10x.h”
typedef struct ADRC
{
float v1,v2; //¸最速输出值
float r; //速度因子
float h; //步长
float z1,z2,z3; //观测器输出
float w0,wc,b0,u; //
}ADRC_NUM;
extern ADRC_NUM ADRC_IN1,ADRC_IN2,ADRC_IN3;
void ADRC_Loop(ADRC_NUM *ADRC_TYPE1);
#endif

#include “ADRC.h”
//
//函数名:void ADRC_Init(ADRC_NUM *ADRC_TYPE1)
//函数说明:LADRC参数初始化
//
void ADRC_Init(ADRC_NUM *ADRC_TYPE1)
{
ADRC_TYPE1->h=0.005; //定时时间及时间步长
ADRC_TYPE1->r=100; //跟踪速度参数
ADRC_TYPE1->w0=50; //观测器带宽
ADRC_TYPE1->wc=10; //状态误差反馈率带宽
}

//
//函数名:void ADRC_TD(ADRC_NUM *ADRC_TYPE1,float Expect)
//函数说明:LADRC跟踪微分部分
//入口参数,期望值Expect(v0)输出值v1,v2
//
void ADRC_TD(ADRC_NUM ADRC_TYPE1,float Expect)
{
float fh= -ADRC_TYPE1->rADRC_TYPE1->r*(ADRC_TYPE1->v1-Expect)-2ADRC_TYPE1->rADRC_TYPE1->v2;
ADRC_TYPE1->v1+=ADRC_TYPE1->v2ADRC_TYPE1->h;
ADRC_TYPE1->v2+=fhADRC_TYPE1->h;
}

void ADRC_ESO(ADRC_NUM ADRC_TYPE1,float FeedBack)
{
float Beita_01=3ADRC_TYPE1->w0;
float Beita_02=3ADRC_TYPE1->w0ADRC_TYPE1->w0;
float Beita_03=ADRC_TYPE1->w0ADRC_TYPE1->w0ADRC_TYPE1->w0;

float e= ADRC_TYPE1->z1-FeedBack;
ADRC_TYPE1->z1+= (ADRC_TYPE1->z2 - Beita_01e)ADRC_TYPE1->h;
ADRC_TYPE1->z2+= (ADRC_TYPE1->z3 - Beita_02e + ADRC_TYPE1->b0ADRC_TYPE1->u)ADRC_TYPE1->h;
ADRC_TYPE1->z3+=-Beita_03e*ADRC_TYPE1->h;
}

void ADRC_LF(ADRC_NUM ADRC_TYPE1)
{
float Kp=ADRC_TYPE1->wcADRC_TYPE1->wc;
float Kd=3ADRC_TYPE1->wc;
float e1=ADRC_TYPE1->v1-ADRC_TYPE1->z1;
float e2=ADRC_TYPE1->v2-ADRC_TYPE1->z2;
float u0=Kpe1+Kd*e2;
ADRC_TYPE1->u=(u0-ADRC_TYPE1->z3)/ADRC_TYPE1->b0;
}

void ADRC_Loop(ADRC_NUM *ADRC_TYPE1)
{
ADRC_TD(ADRC_TYPE1,Expect_Value);
ADRC_ESO(ADRC_TYPE1,Measure);
ADRC_LF(ADRC_TYPE1);
}

 

  • 4
    点赞
  • 54
    收藏
    觉得还不错? 一键收藏
  • 3
    评论
LADRC是一种控制算法,其C语言实现可以如下所示: ```c #define N 3 // 系统状态变量的数量 #define M 1 // 控制信号的数量 float x[N]; // 系统状态变量 float u[M]; // 控制信号 float r[N]; // 参考信号 float q[N][N]; // 状态误差加权系数矩阵 float p[N][N]; // 状态协方差矩阵 float k[N][M]; // 控制增益矩阵 float a[N][N]; // 系统矩阵 float b[N][M]; // 输入矩阵 float c[N][N]; // 观测矩阵 float d[M][M]; // 输入协方差矩阵 float e[N][N]; // 状态协方差矩阵 void ladrc_control(float *y, float *x, float *u) { float z[N]; // 观测向量 float v[M]; // 测量噪声 float w[N]; // 状态噪声 float kk[N][N]; // Kalman增益矩阵 // 计算观测向量 for (int i = 0; i < N; i++) { z[i] = c[i][0] * y[0]; for (int j = 1; j < N; j++) { z[i] += c[i][j] * x[j]; } } // 计算控制增益矩阵 for (int i = 0; i < N; i++) { for (int j = 0; j < M; j++) { k[i][j] = 0; for (int l = 0; l < N; l++) { k[i][j] += q[i][l] * c[l][j]; } k[i][j] /= d[j][j]; } } // 计算状态协方差矩阵 for (int i = 0; i < N; i++) { for (int j = 0; j < N; j++) { e[i][j] = 0; for (int l = 0; l < N; l++) { e[i][j] += a[i][l] * p[l][j]; } for (int l = 0; l < M; l++) { e[i][j] += b[i][l] * d[l][l] * b[j][l]; } } } // 计算Kalman增益矩阵 for (int i = 0; i < N; i++) { for (int j = 0; j < N; j++) { kk[i][j] = 0; for (int l = 0; l < N; l++) { kk[i][j] += e[i][l] * c[l][j]; } kk[i][j] /= (d[0][0] + kk[0][0]); } } // 计算状态估计值 for (int i = 0; i < N; i++) { x[i] = z[i]; for (int j = 0; j < N; j++) { x[i] -= kk[i][j] * z[j]; } } // 计算控制信号 for (int i = 0; i < M; i++) { u[i] = 0; for (int j = 0; j < N; j++) { u[i] += k[j][i] * (r[j] - x[j]); } } // 更新状态协方差矩阵 for (int i = 0; i < N; i++) { for (int j = 0; j < N; j++) { p[i][j] = e[i][j]; for (int l = 0; l < N; l++) { p[i][j] -= kk[i][l] * c[l][j] * e[0][0]; } } } } ``` 其中,`x`、`u`、`r`分别表示系统状态变量、控制信号和参考信号,`q`表示状态误差加权矩阵,`p`表示状态协方差矩阵,`k`表示控制增益矩阵,`a`、`b`、`c`、`d`分别表示系统矩阵、输入矩阵、观测矩阵和输入协方差矩阵,`e`表示状态协方差矩阵。函数`ladrc_control`表示LADRC控制算法的实现。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 3
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值