MATLAB: I have a non-closed figure but i need to link the limits

pi controller

this is the curve
and this is the code
w=0:0.001:1.57;
kp=zeros(1,length(w));
ki=zeros(1,length(w));
mu=1;
teta=0.5;
phi=0
A=1;
lamda=1;
kd=0;
for i=1:length(w)  
  kp(i)=-((1/A)*((w(i).^(lamda+1.5)*sin(w(i)*teta+phi+((lamda+1.5)*pi/2))))+(kd*w(i).^(lamda+mu)*sin((lamda+mu)*pi/2)))/(w(i).^(lamda)*sin((lamda*pi/2)));
  ki(i)=((1/A)*(w(i).^(2*lamda+1.5)*sin(teta*w(i)+phi+2.36))+(kd*w(i).^(2*lamda+mu)*sin((mu)*pi/2)))/(w(i).^(lamda)*sin((lamda*pi/2)));
end
plot(kp,ki)

Best Answer

Maybe try area() instead of plot(). If not, then describe what closed means to you.
w=0:0.001:1.57;
kp=zeros(1,length(w));
ki=zeros(1,length(w));
mu=1;
teta=0.5;
phi=0
A=1;
lamda=1;
kd=0;
for i=1:length(w)  
  kp(i)=-((1/A)*((w(i).^(lamda+1.5)*sin(w(i)*teta+phi+((lamda+1.5)*pi/2))))+(kd*w(i).^(lamda+mu)*sin((lamda+mu)*pi/2)))/(w(i).^(lamda)*sin((lamda*pi/2)));
  ki(i)=((1/A)*(w(i).^(2*lamda+1.5)*sin(teta*w(i)+phi+2.36))+(kd*w(i).^(2*lamda+mu)*sin((mu)*pi/2)))/(w(i).^(lamda)*sin((lamda*pi/2)));
end
area(kp,ki)
xlabel('kp', 'FontSize', 20);
ylabel('ki', 'FontSize', 20);
grid on;
Related Question