MATLAB解宋健人口模型

1.%求出p(t,r)函数symsstr;u=0.00001*s^4-0.001*s^3+0.032*s^2-34*s+2.2;%U(r)为起始年份即t=0时的各个年龄段的人口死亡率I=simple(int(u/1000,s,r-t,r));p=(123.06*(r-t)^3-23601*(r-t)^2+1022300*(r-t)+9794300)*exp(-I)%p(r,0)=123.06*(r-t)^3-23601*(r-t)^2+1022300*(r-t)+9794300为t=0时的人口密度2.%代入p(t,r)函数，定义M文件functionyuce=ppp(t,r)t=t;p=((44665600*r)/53-(44665600*t)/53-(1042640*(r-t)^2)/53+(132252*(r-t)^3)/1325+11232000)/exp(((((r-t)*(t/500000000-r/500000000+1/4000000)-1/93750)*(r-t)+17/100000)*(r-t)-11/5000)*(r-t)+r*(r*(r*(r*(r/500000000-1/4000000)+1/93750)-17/100000)+11/5000));3.运行程序，求解不同年份a时的人口数symsria;a=;%令t=0时为2005年,代入不同的年份a,从而求其总人口t=a-2005;b=1.5;s1=0;k=4.3969*10^(-8)*r^3-8.3066*10^(-6)*r^2+0.00034379*r+0.0043419;h=(r-15)^4.5*exp(-(r-15)/2)/1595;fori=1:ts1(i+1)=s1(i)+(double(symsum(b*k*prpr(i,r)*h,15,49)));ends2=double(symsum(prpr(t,r),t,100));s=s1+s2;