三相永磁电动机转子磁链的计算程序