你是一名专业的算法大师，请帮我使用matlab编写代码求出电机1和电机2的最大转矩电流比(MTPA)曲线和最大转矩电压比(MTPV)曲线，计算不考虑饱和特性 已知：电机1参数分别如下:44kW PMSM:定子相电阻R=4m2定子链相绕组永磁磁链峰值;=0055Wb电机极对数p=4，转动惯量J=0.048kg·m，定子绕组相漏电感L=2pHd轴励磁电感L;=88pH，q轴励磁电感L=300pH，额定转速为 4000r/min，峰值转矩temax=210N·m。不考虑饱和特性。电机2参数分别如下：2100kW PMSM:定子相电阻R=83m定子链相绕组永磁磁链峰值=0.071Wb，电机极对数p=4，转动惯量J=0.1kgm，定子绕组相漏电感L=30uHd轴励磁电感L;=174uHq轴励磁电感 L=293pH，额定转速为4700r/min峰值转矩 temax256N·m。

好的，我会尽力回答您的问题。根据您提供的电机参数，可以使用以下Matlab代码计算电机1和电机2的MTPA和MTPV曲线：

% 电机1参数
P_1 = 44e3; % 额定功率（W）
R_1 = 4e-2; % 定子相电阻（Ω）
Ld_1 = 88e-6; % d轴励磁电感（H）
Lq_1 = 300e-6; % q轴励磁电感（H）
lambda_m_1 = 0.055; % 定子链相绕组永磁磁链峰值（Wb）
J_1 = 0.048; % 转动惯量（kg·m²）
p_1 = 4; % 极对数
temax_1 = 210; % 峰值转矩（N·m）
n_rated_1 = 4000; % 额定转速（r/min）

% 电机2参数
P_2 = 2100e3; % 额定功率（W）
R_2 = 83e-3; % 定子相电阻（Ω）
Ld_2 = 174e-6; % d轴励磁电感（H）
Lq_2 = 293e-6; % q轴励磁电感（H）
lambda_m_2 = 0.071; % 定子链相绕组永磁磁链峰值（Wb）
J_2 = 0.1; % 转动惯量（kg·m²）
p_2 = 4; % 极对数
temax_2 = 256; % 峰值转矩（N·m）
n_rated_2 = 4700; % 额定转速（r/min）

% 计算MTPA曲线
n1 = 0:100:n_rated_1; % 电机1转速范围
i1_max = sqrt((P_1/R_1)./((n1*2*pi/60/p_1).^2.*(J_1*(Ld_1+Lq_1)-Lq_1^2)));
n2 = 0:100:n_rated_2; % 电机2转速范围
i2_max = sqrt((P_2/R_2)./((n2*2*pi/60/p_2).^2.*(J_2*(Ld_2+Lq_2)-Lq_2^2)));

% 计算MTPV曲线
v1_max = sqrt(R_1*temax_1/lambda_m_1); % 电机1最大转矩电压比
v2_max = sqrt(R_2*temax_2/lambda_m_2); % 电机2最大转矩电压比
v1 = 0:0.01:v1_max; % 电机1电压范围
i1_mtpv = sqrt((P_1/R_1)./((v1/sqrt(2)/lambda_m_1).^2.*(J_1*(Ld_1+Lq_1)-Lq_1^2)));
v2 = 0:0.01:v2_max; % 电机2电压范围
i2_mtpv = sqrt((P_2/R_2)./((v2/sqrt(2)/lambda_m_2).^2.*(J_2*(Ld_2+Lq_2)-Lq_2^2)));

% 绘制MTPA曲线
figure;
plot(n1, i1_max, 'b');
hold on;
plot(n2, i2_max, 'r');
xlabel('转速（r/min）');
ylabel('最大转矩电流（A）');
legend('电机1', '电机2');
title('MTPA曲线');

% 绘制MTPV曲线
figure;
plot(v1, i1_mtpv, 'b');
hold on;
plot(v2, i2_mtpv, 'r');
xlabel('电压（V）');
ylabel('最大转矩电流（A）');
legend('电机1', '电机2');
title('MTPV曲线');

运行以上代码可以得到MTPA和MTPV曲线的图像。需要注意的是，上述代码不考虑饱和特性，因此只适用于小转速范围。

希望能够帮到您！