求MATLAB的高手，怎么把两张仿真曲线合并显示在一张图上？

关注:194  答案:5  手机版

解决时间 2021-02-10 14:48

提问者终究是陌生了

2021-02-09 16:38

程序1%using  pole  placement-Ackermann%Modellinga=[0 1.0000 0 0;0 -0.1818 2.6727 0;0 0 0 1.0000;0 -0.4545 31.1818 0];b=[0;1.8182;0;4.5455];c=[1 0 0 0;0 0 1 0];d=[0;0];%Check the controllabilitydisp('The rank of controllability matrix')rc=rank(ctrb(a,b))%Designp=[-7.4527+9.666j,-7.4527-9.666j ,-3.1538+1.8334j, -3.1538-1.8334j]k=acker(a,b,p),x0=[0.2,0,0,0]%The close loop state system is denoted as (ac,b,c,d)ac=[a-b*k ];t=0:0.005:3;figureinitial(ac,b,c,d,x0,t)title(' initial respond of inverted pendulum system ')xlabel('Sec')ylabel('Output y=x3')figure[y,x,t]=initial(ac,b,c,d,x0,t);plot(t,x,'y')figure,x1=[1 0 0 0]*x',plot(t,x1);grid;title('x1的响应曲线') 程序2a=[0 1.0000 0 0;0 -0.1818 2.6727 0;0 0 0 1.0000;0 -0.4545 31.1818 0];b=[0;1.8182;0;4.5455];c=[1 0 0 0;0 0 1 0];d=[0;0];eig(a)Qc=ctrb(a,b);Rc=rank(Qc)if (Rc==4)    disp('状态完全能控')else    disp('状态不完全能控')endQo=obsv(a,c);Ro=rank(Qo)if( Ro==4)    disp('状态完全能观测')else    disp('状态不完全能观测')endQ=[25 0 0 0;0 0 0 0;0 0 81 0;0 0 0 0];R=1;                               K=lqr(a,b,Q,R);x0=[0.2,0,0,0];ac=[(a-b*K)];bc=[b];cc=[c];dc=[d];[K,P,e]=lqr(a,b,Q,R),t=0:0.005:7;figure      initial(ac,b,c,d,x0)title('step respond of inverted pendulum system ')xlabel('Sec')ylabel('Output y=x3')figure[y,x,t]=initial(ac,bc,cc,dc,x0,t);plot(t,x,'y')figure,x1=[1 0 0 0]*x',plot(t,x1);grid;title('x1的响应曲线') 怎么样才能把这两个程序的XI 显示在一张图上

最佳答案

二级知识专家暖心欧巴

2021-02-09 17:11

既然小编提到我了，我就好好看了一下程序，因为不确定画的其他图像是否有用，所以没删。最后一张图满足要求，修改如下：

a=[0 1.0000 0 0;0 -0.1818 2.6727 0;0 0 0 1.0000;0 -0.4545 31.1818 0];

b=[0;1.8182;0;4.5455];c=[1 0 0 0;0 0 1 0];d=[0;0];

%Check the controllability

disp('The rank of controllability matrix')

rc=rank(ctrb(a,b))

%Design

p=[-7.4527+9.666j,-7.4527-9.666j ,-3.1538+1.8334j, -3.1538-1.8334j]

k=acker(a,b,p),x0=[0.2,0,0,0]

%The close loop state system is denoted as (ac,b,c,d)

ac=[a-b*k ];t=0:0.005:3;

figure

initial(ac,b,c,d,x0,t)

title(' initial respond of inverted pendulum system ')

xlabel('Sec')

ylabel('Output y=x3')

figure

[y,x,t]=initial(ac,b,c,d,x0,t);

plot(t,x,'y')

figure,

x1=[1 0 0 0]*x',

plot(t,x1);

grid;

title('x1的响应曲线')

T=t;

X1=x1; % 后面用

a=[0 1.0000 0 0;0 -0.1818 2.6727 0;0 0 0 1.0000;0 -0.4545 31.1818 0];

b=[0;1.8182;0;4.5455];c=[1 0 0 0;0 0 1 0];d=[0;0];

eig(a)

Qc=ctrb(a,b);

Rc=rank(Qc)

if (Rc==4)    disp('状态完全能控')

else    disp('状态不完全能控')

end

Qo=obsv(a,c);Ro=rank(Qo)

if( Ro==4)    disp('状态完全能观测')

else    disp('状态不完全能观测')

end

Q=[25 0 0 0;0 0 0 0;0 0 81 0;0 0 0 0];

R=1;

K=lqr(a,b,Q,R);

x0=[0.2,0,0,0];

ac=[(a-b*K)];

bc=[b];cc=[c];dc=[d];

[K,P,e]=lqr(a,b,Q,R),t=0:0.005:7;

figure

initial(ac,b,c,d,x0)

title('step respond of inverted pendulum system ')

xlabel('Sec')

ylabel('Output y=x3')

figure

[y,x,t]=initial(ac,bc,cc,dc,x0,t);

plot(t,x,'y')

figure,

x1=[1 0 0 0]*x',

hold on

plot(T,X1,'r');

plot(t,x1);

grid;

title('x1的响应曲线')

legend   程序1的响应曲线  程序2的响应曲线两张图时间轴不一般长，但能画。需要的就是最后一张图，如下：

全部回答

1楼一个狠有粪量的人

2021-02-09 19:26

使用hold on 命令

2楼何必执着

2021-02-09 18:44

很简单啊，用figure

plot(a,b);

figure;

plot(x,y);

3楼雨落轻尘

2021-02-09 17:27

hold on，在两句画图之间加一句这个代码。

4楼清酒孤欢

2021-02-09 17:22

很简单，

plot(y1);

hold on;

plot(y2);

这样就能显示在一张图上了，赶快试试吧。

我要举报

如以上信息为低俗、色情、不良、暴力、侵权、涉及违法等信息，可以点下面链接进行举报！

推荐资讯

大家都在看