clear all
close all
clc
% coeficientes da curva de potência
A0=-1.278493e+05;
A1=2.195991e+02;
A2=-1.248258e-01;
A3=6.622842e-05;
A4=-2.172773e-08;
A5=3.559221e-12;
A6=-2.289778e-16;

NP=15;

nrpm=linspace(1000,4500,NP).';
PeW=A0+A1*nrpm+A2*nrpm.^2+A3*nrpm.^3+A4*nrpm.^4+A5*nrpm.^5+A6*nrpm.^6;
%% cálculo do torque (N.m)
T=PeW./(2*pi*nrpm/60);
%-------------------------------------------------------------------------
%% gráfico de torque e potência do motor
%-------------------------------------------------------------------------
figure('name','Torque e potência do motor')
[AX,H1,H2] = plotyy(nrpm,T,nrpm,PeW);
set(get(AX(1),'Ylabel'),'String','T (Nm)') 
set(get(AX(2),'Ylabel'),'String','Pe (W)') 
set(H1,'LineStyle',':')
set(H2,'LineStyle','--')
set(H1,'LineWidth',2)
set(H2,'LineWidth',2)
grid on
xlabel('n (rpm)')
