这个代码:
alpha=0.3;
a=0.15;
b=0.07;
c=0.05;
m=1.12;
u1=26.5;
u2=31.6;
g=0.15;
n=0.07;
h=0.05;
K0=1000
H0=500
L0=750
K1(1:2,1)=K0;
L1(1:2,1)=L0;
H1(1:2,1)=H0;
K2(1:2,1)=K0;
L2(1:2,1)=L0;
H2(1:2,1)=H0;
MPK1(:,1)=alpha.*u1.*K1(:,1).^(alpha-1).*L1(:,1).^(1-alpha).*H1(:,1).^m;
MPK2(:,1)=alpha.*u2.*K2(:,1).^(alpha-1).*L2(:,1).^(1-alpha).*H2(:,1).^m;
for t=1:100
K1(:,t)=((1+g).*(1-a)).^(t-1).*K0;
L1(:,t)=(1+n).^(t-1).*L0.*(t>=0&t<50)+(1+n).^(t-1).*(1-c).^(t-50).*L0.*(t>=50);
H1(:,t)=((1+h).*(1-b)).^(t-1).*H0;
MPK1(:,t)=alpha.*u1.*K1(:,t).^(alpha-1).*L1(:,t).^(1-alpha).*H1(:,t).^m;
K2(:,t)=(1+g).^(t-1).*(1+(1-a).^(t-1)).*K0;
L2(:,t)=(1+n).^(t-1).*(1+(1-c).^(t-1)).*L0;
H2(:,t)=(1+h).^(t-1).*(1+(1-b).^(t-1)).*H0;
MPK2(:,t)=alpha.*u2.*K2(:,t).^(alpha-1).*L2(:,t).^(1-alpha).*H2(:,t).^m;
end
plot(t,MPK1,t,MPK2);
想吧MPK1和MPK2画在一张图上(横坐标都是t),怎么画捏。。。
求大神告知Orz
alpha=0.3;
a=0.15;
b=0.07;
c=0.05;
m=1.12;
u1=26.5;
u2=31.6;
g=0.15;
n=0.07;
h=0.05;
K0=1000
H0=500
L0=750
K1(1:2,1)=K0;
L1(1:2,1)=L0;
H1(1:2,1)=H0;
K2(1:2,1)=K0;
L2(1:2,1)=L0;
H2(1:2,1)=H0;
MPK1(:,1)=alpha.*u1.*K1(:,1).^(alpha-1).*L1(:,1).^(1-alpha).*H1(:,1).^m;
MPK2(:,1)=alpha.*u2.*K2(:,1).^(alpha-1).*L2(:,1).^(1-alpha).*H2(:,1).^m;
for t=1:100
K1(:,t)=((1+g).*(1-a)).^(t-1).*K0;
L1(:,t)=(1+n).^(t-1).*L0.*(t>=0&t<50)+(1+n).^(t-1).*(1-c).^(t-50).*L0.*(t>=50);
H1(:,t)=((1+h).*(1-b)).^(t-1).*H0;
MPK1(:,t)=alpha.*u1.*K1(:,t).^(alpha-1).*L1(:,t).^(1-alpha).*H1(:,t).^m;
K2(:,t)=(1+g).^(t-1).*(1+(1-a).^(t-1)).*K0;
L2(:,t)=(1+n).^(t-1).*(1+(1-c).^(t-1)).*L0;
H2(:,t)=(1+h).^(t-1).*(1+(1-b).^(t-1)).*H0;
MPK2(:,t)=alpha.*u2.*K2(:,t).^(alpha-1).*L2(:,t).^(1-alpha).*H2(:,t).^m;
end
plot(t,MPK1,t,MPK2);
想吧MPK1和MPK2画在一张图上(横坐标都是t),怎么画捏。。。
求大神告知Orz