脉振方波高频注入代码+增强型滑膜esmo代码,永磁同步电机高频注入程序 资料为C代码一份,大厂代码,可运行,经典流传; 配套一篇代码对应的说明文档,详细算法说明; 脉振方波注入方法相对于脉振正弦信号注入的形式,信号处理的过程少了一些滤波器 ,计算更简单,并且由于信号频段较高,可以实现更高的动态响应能力。 增强滑膜控制;
//高频信号生成模块 void HF_SquareWave_Inject(float amplitude, uint32_t freq) { static uint8_t hf_cycle = 0; hf_signal = (hf_cycle < HF_CYCLE/2) ? amplitude : -amplitude; hf_cycle = (hf_cycle + 1) % HF_CYCLE; //叠加到基波电压指令 Vd_inj = Vd_ref + hf_signal; Vq_inj = Vq_ref; }这代码妙在哪儿?相比正弦注入方案,直接采用方波省掉了正弦查表和滤波器。注意看那个HF_CYCLE参数,它控制着方波的切换频率。实际调试时遇到过坑——当载波频率和注入频率不成整数倍关系时,会在电流频谱上产生混叠。解决办法是在PWM中断里做计数周期对齐。
接着看增强型滑模观测器的实现重点:
//ESMO核心计算函数 void SMO_Enhanced(float i_alpha, float i_beta) { //滑模面计算 s_alpha = i_alpha_hat - i_alpha; s_beta = i_beta_hat - i_beta; //改进型切换函数 sign_alpha = (s_alpha > 0) ? 1 : -1; sign_beta = (s_beta > 0) ? 1 : -1; //边界层处理 if(fabs(s_alpha) < epsilon) sign_alpha = s_alpha / epsilon; if(fabs(s_beta) < epsilon) sign_beta = s_beta / epsilon; //参数更新 i_alpha_hat += Ts*( -Rs/Ls*i_alpha + ... + K_slide*sign_alpha ); i_beta_hat += Ts*( -Rs/Ls*i_beta + ... + K_slide*sign_beta ); }传统滑模的抖振问题在这里被两个技巧化解:一是用饱和函数替代符号函数(边界层处理),二是动态调整滑模增益Kslide。实测发现当电机转速超过基速的30%时,将Kslide降低30%能显著抑制高频噪声。
信号解调环节有个骚操作:
//高频响应提取 void HF_Response_Extract(void) { static float hf_buffer[4] = {0}; //移动平均滤波 hf_buffer[3] = hf_buffer[2]; hf_buffer[2] = hf_buffer[1]; hf_buffer[1] = hf_buffer[0]; hf_buffer[0] = I_alpha_meas * hf_signal; //相关运算 //差分处理 delta_hf = (hf_buffer[0] - hf_buffer[2])/2; //位置误差计算 position_error = sign_hf * delta_hf * HF_GAIN; }这里用移动平均代替FIR滤波器,牺牲了一点滤波效果但节省了30%的运算时间。注意sign_hf变量需要与注入信号严格同步,之前调试时因为时序错位导致估算位置漂移的坑,后来在注入信号翻转前插入5us的延时才解决。
整套方案在零低速下的实测表现:30rpm时位置估算误差小于0.2弧度,载波频率18kHz下CPU占用率仅12%。但要注意电机参数敏感性——当电感值偏差超过20%时需要进行在线补偿,否则高频电流幅值会异常增大。
最后说个实战经验:高频注入方案对逆变器死区特别敏感。曾经在某个量产项目中,电机反转时出现周期性震荡,后来发现是死区补偿没做好。解决办法是在位置估算环节加入死区电压前馈,代码里加了两行:
Vd_comp += DeadTime_Compensation(qei.speed); //死区电压补偿 Vq_comp += Cross_Coupling_Comp(); //交叉耦合补偿这种细节处理才是大厂代码的精髓,算法文档里可不会写这些实战技巧。