程序写得没有层次,查错难。按这个程序y只计算出了一个点,原因是以分子除分母时,除号没有用点除(/),改为以下即可:
t=1:01:20;
y=(10(t^2 +100)(23^(1/2)(t^2 + 25)^(1/2) +
(3^(1/2)t^2(t^2 +25)^(1/2))/50 +
(35192^(1/2)(t^2 + 25)^(1/2)(2t^2 +25)^(1/2))/200 +
32^(1/2)(t^2 + 100)^(1/2)(2t^2 + 25)^(1/2) +
(6^(1/2)(t^2 + 100)^(1/2)(2t^2 + 25)^(1/2))/2))
/(9(100t(t^2 +25)^(1/2) + t^3(t^2 + 25)^(1/2) + 252^(1/2)t(t^2 +100)^(1/2)(2t^2 + 25)^(1/2)));
plot(t,y)
t=1900:10:1990;
p=[75995 91972 105711 123203 131669
150697 179323 203212 226505 249633];
a=interp1(t,p,1975);
plot(t,p,'-o');hold on
plot(1975,a,'+');hold off
text(1977,a,num2str(a));
用线性插值就差不多了
得到的a即使
1975年的人口估计
图中带有圈的线是原来数据画出的人口曲线
而+号表示估计的1975年人口,约214859万
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)