已知电机参数为：
Rs=0.008Ω
Ld=200uH
Lq=400uH
Psif=0,04wb
p=4(电机极对数)
IsMAX=500A
UsMAX=400V
利用matlab求解或画出：
1.电流幅值为IsMAX不同的矢量角beta对应：磁阻转矩、励磁转矩、整体转矩；
2.最大转矩电流比（MTPA）曲线；
3.电流极限圆曲线；
4.不同扭矩对应的恒扭矩曲线簇；
5.电压极限椭圆曲线在不同速度下。
%% 电机参数 Rs=0.008; Ld=200e-6; Lq=400e-6; psif=0.04; p=4; ismax=500; usmax=400; %% 不同矢量角β对应:磁阻转矩、励磁转矩、整体转矩 beta=0:0.01*pi:pi; te=3/2*p*(psif*ismax*sin(beta) 1/2*(Ld-Lq)*ismax^2*sin(2*beta)); te_cizu=3/2*p*1/2*(Ld-Lq)*ismax^2*sin(2*beta); te_lici=3/2*p*psif*ismax*sin(beta); figure(1); plot(beta,te); hold on; plot(beta,te_lici); hold on; plot(beta,te_cizu); legend(合成转矩，励磁转矩，磁阻转矩)； title('te-β特征曲线； xlabel('β/rad'); ylabel('te/Nm');