19-11-16 04:36 PM D:\Mis documentos\MATERI...\RiveroMTLB.
m 1 of 1
%PROGRAMACION GRAFICADOR
%NOMBRE: Rivero Albornoz Carlos Erasmo
%CARRERA: Ing. Electromecanica
close all;
clc;
clear;
E_o=input('Ingrese el valor de Eo: ');
Q_u=input('Ingrese el valor de Qu: ');
R=-Q_u;
P_v=input('Ingrese el valor de Pv: ');
Aux=E_o*R;
p1=[Q_u P_v];
p2=[1 P_v];
plot([p1(1) p2(1)],[p1(2) p2(2)],'g','linewidth',3);
hold on;
p3=[Q_u 0];
p4=[Q_u E_o];
plot([p3(1) p4(1)],[p3(2) p4(2)],'c','linewidth',3);
tita=pi/2:0.01:pi;
x_1=Aux*cos(tita);
y_1=Aux*sin(tita);
plot(Q_u-x_1,y_1,'r','linewidth',3);
hold on;
alfa=0:0.01:pi;
x_2=cos(alfa);
y_2=sin(alfa);
plot(x_2,y_2,'b','linewidth',3);
hold on;
x1=Q_u:0.01:-0.15;
y1=(tan(75*pi/180))*(x1-Q_u);
plot(x1,y1);
text(0.5,0.5,'Lim Pot. Turbina');
text(0.1,1.3,'zona sobreexitado');
text(-0.5,1.3,'zona subexitado');
text(-0.3,0.1,'75°');
text(0.5,0.97,'Fdp=0,8');
plot([0 0],[0 1.5],'black','linewidth',2);
grid on;
axis square;
axis equal;
axis([-1 1 0 1.5]);
caxis('auto')
axis on;