필터 지우기
필터 지우기

How to plot multiple curves in the same graph window?

조회 수: 9 (최근 30일)
Aishwarya Govekar
Aishwarya Govekar 2020년 3월 15일
댓글: Star Strider 2020년 3월 16일
I want to plot graphs for multiple Kp values for a proportional controller. I tried using for loop to perform this. But I didn't get the output.Below, I attach the code which I have written for your reference. Please guide me on this. If you can tell me if there's any other way to do it, then that would be great.
Kindly reply asap.
Thank you.
  댓글 수: 5
dpb
dpb 2020년 3월 15일
Just paste short demo code as text, use the code button to format ... much easier than having to download a file and open it externally.
For figures, SaveAs a .jpg and then attach the immage with the picture icon at the INSERT group.
I'll guess just from the words you didn't use "hold on" after the first plot to add subsequent to the same axes...
Aishwarya Govekar
Aishwarya Govekar 2020년 3월 15일
clc;
clear all;
close all;
J=20.48*10e-6;
b=0.702;
kt=14.70;
kb=0.48;
R=1.6;
L=0.02*0.001;
num=kt;
den=[(J*L) ((J*R)+(L*b)) ((b*R)+(kb*kt))];
for kp=0:50:300 %plotting graph for P controller with Kp values in steps of 50
ki=0;
kd=0;
numc=[kd,kp,ki];
denc=[0 1];
numa=conv(num,numc);
dena=conv(den,denc);
[numac,denac]=cloop(numa,dena);
step(numac,denac)
hold on;
xlabel('time(sec)'), ylabel('velocity(rad/sec)')
title ('pid control');
grid;
end

댓글을 달려면 로그인하십시오.

채택된 답변

Star Strider
Star Strider 2020년 3월 15일
Try this:
J=20.48*10e-6;
b=0.702;
kt=14.70;
kb=0.48;
R=1.6;
L=0.02*0.001;
num=kt;
den=[(J*L) ((J*R)+(L*b)) ((b*R)+(kb*kt))];
kp=0:50:300; %plotting graph for P controller with Kp values in steps of 50
for k = 1:numel(kp)
ki=0;
kd=0;
numc=[kd,kp(k),ki];
denc=[0 1];
numa=conv(num,numc);
dena=conv(den,denc);
[numac,denac]=cloop(numa,dena);
sys = tf(numac,denac)
[y{k},t{k}] = step(sys);
% QL(k,:) = [min(y{k}), max(y{k})] % Information
end
figure
hold on
for k = 2:numel(y)
plot(t{k}, y{k})
end
hold off
xlabel('time(sec)'), ylabel('velocity(rad/sec)')
title ('pid control');
grid
It is necessary to create a system object to use the step funciton. The first step output is uniformly zero, and that causes problems with the plot. Start with the second one instead.
Also, use feedback rather than cloop (that has been deprecated).
  댓글 수: 8
Aishwarya Govekar
Aishwarya Govekar 2020년 3월 16일
편집: Aishwarya Govekar 2020년 3월 16일
How to change the range of x and y axes of the graph which is already plotted?
I want to obtain the graph similar to the graph I have attached below there. Scale of the graph seems to be different on execution of this code. How do we change it? I tried to use linspace() function for this before giving range for kp values. But it didn't work.
Please help.
Thanks.
J=0.01;
b=0.1;
p=0.01;
R=1;
L=0.5;
num=p;
den=[ (J*L) ((J*R)+(L*b)) ((b*R) + p^2)] ;
%%%%% Different values of kp , ki and kd
%%%%%%%%%
kp=0:30:300;
for k = 1:numel(kp)
ki=0;
kd=0;
numc=[kd,kp(k),ki];
denc=[0 1];
numa=conv(num,numc);
dena=conv(den,denc);
[numac,denac]=cloop(numa,dena);
sys = tf(numac,denac)
[y{k},t{k}] = step(sys);
% QL(k,:) = [min(y{k}), max(y{k})] % Information
end
figure
hold on
for k = 2:numel(y)
plot(t{k}, y{k})
end
hold off
xlabel('time(sec)'), ylabel('velocity(rad/sec)')
title ('pid control');
Star Strider
Star Strider 2020년 3월 16일
To change the axis limits use the xlim, ylim or axis functions. To set limits on axes that are already plotted, see the documentation on Axes Properties.

댓글을 달려면 로그인하십시오.

추가 답변 (0개)

카테고리

Help CenterFile Exchange에서 Graphics Performance에 대해 자세히 알아보기

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by