The code:
clc;clear;
t=0:0.01:10;
wm=188.5.*(1-exp(-0.2.*t));
Tind=10.*exp(-0.2.*t);
p=Tind.*wm;
plot(t,wm)
title('Rotational speed Vs. Time');
xlabel('Time (s)');
ylabel('\it{\omega}_m (rad/s)');
figure, plot(t,Tind)
title('Induced torque Vs. Time');
xlabel('Time (s)');
ylabel('\it{\tau}_{ind} (N.m)');
figure, plot(t,p)
title('Power Vs. Time');
xlabel('Time (s)');
ylabel('P (Watt)');
[maxp i]=max(p);
fprintf('The maximum is %f and its at %f seconds \n',maxp,i);
the output in the Command Window:
The maximum is 471.249658 and its at 348.000000 seconds
No comments:
Post a Comment