BLDC速度控制仿真请教

doubushe

在主程序中采用龙格-库塔法计算相电流，但是计算相电流切线斜率的函数如何定义，seq_ias()如何定义，请指教。

%4阶龙格-库塔法计算ias
if(ias_det==1)
phase=1;
grad_1=seq_ias(tt,ias);
grad_2=seq_ias(tt+dt/2,ias+dt*grad_1/2);
grad_3=seq_ias(tt+dt/2,ias+dt*grad_2/2);
grad_4=seq_ias(tt+dt,ias+dt*grad_3);
ias=ias+(grad_1+2*grad_2+2*grad_3+grad_4)*dt/6;
end