In general, to create these types of curves, you need to find some variable that you can sweep, to generate the curve. In your case, that variable should be tied to either power or torque. For instance, you could sweep the calculations over the full range of torques, for a given alpha, and then also calculate the equivalent power since that is known if you know speed and torque. You can then iterate on this for each value of alpha. Unfortunately, your code is missing a lot of parameters and doesn't work, so you need to first get a set of equations that work for a single operating condition before worrying about trying to generate a curve.
Here is an example for a very simple system with basic drag and gravity effects on vehicle dynamics. These equations are too simplistic to be realistic, particularly at the higher speeds:
grade = [0 .04 .07 .14]; % percent
slope = atan(grade);
Mass = 1814; % 2 ton car in kg;
g = 9.8;
fx = -g*Mass*sin(slope); % reverse force from gravity
cd = 0.3; % average drag coefficient for a car
area = 7.33*0.0929; % area for an Acura CL
density = 1.21; % density of air
r = .48; % radius of tire to convert torque to force
final_drive_ratio = 4.2;
gear_ratio = 1.2;
TorqueSweep = 0:1:200; % sweep betweein 0 and 200 NM on the engine
TorqueTireSweep = TorqueSweep*final_drive_ratio*gear_ratio;
ForceSweep = TorqueTireSweep/r;
n = length(TorqueSweep);
for i = 1:length(fx)
dragForce = ForceSweep + fx(i); % need to solve for steady state operating condition
dragForce = max(dragForce',zeros(1,n)');
speed(i,:) = sqrt(dragForce/(0.5*cd*area*density))*2.236; % converstion to mph
end
figure(1)
plot(TorqueSweep,speed)
legend('0','4%','7%','14%')
xlabel('Torque (Nm)')
ylabel('Speed (mph)')
