A+20004005+ 周旭、赵基宇、毛建锋(6)

Qp(:,6+i)=b;

Xp(6+i,:)=f3(Xp(6+i-1,:),s0,b); end

n=length(Xp(:,1)); bmnr1; for i=1:n-6 if i==1

Kp(i)=(X(i+5,:)*s0*K(6)/1000+a*X(i+6,1))/sum(X(i+6,:)); else

Kp(i)=(Xp(i+6,:)*s0*Kp(i-1)/1000+bmnr1(i)*Xp(i+6,1))/sum(Xp(i+6,:)) ; end end

function f12(p,k,str)

%画图，不同力度下总人口的预测； load matlab.mat str='r.b0yxg+k*'; for i=1:length(p)

[qp,xp,kp]=f11(X,Q,p(i),s0,1.185,k,K); sw=sum(xp');

kp=[(K(1:6))',kp]; figure(1);

plot(1:length(kp),kp,'*'); hold on;

sall=sw.*(kp+1); figure(2);

plot(1:k+6,sall*1000,'.'); hold on; end

24