24小时热门版块排行榜

返回列表

二哥哥33

木虫 (职业作家)

应助: 1 (幼儿园)
金币: 6834.5
散金: 4990
红花: 17
帖子: 3312
在线: 911.2小时
虫号: 1074962
注册: 2010-08-15
专业: 泛函分析

[交流] 帮我看一段Matlab程序哪里出错了已有1人参与

代入2001年～2005年的数据，预测一下我国未来50年的人口变化情况，2001年-2005年的数据为12627 128453 129227 130000 130756.
程序如下：
clc
i=input('你想预测今后(从2005年开始)多少年的中国人口数据？请输入参数：');
while i<0
disp('历史还需要预测吗？请您重新输入大于0的数字：')
i=input('');
end
if i==0
plot(2005,d1(16,2),'*')
else
y(2)=d1(12,2);
y(3)=d1(13,2);
y(4)=d1(14,2);
y(5)=d1(15,2);
n=d1(12:end,2);
for k=5:i+5
      x=k-4:k;
      y=[y(2),y(3),y(4),y(5),n(end)];
      p=polyfit(x,y,2);
      n(k+1)=polyval(p,k+1);
end
t=2005:2005+i;
plot(t,n(6:end),'*-')
str1='2005年到';
str2=num2str(i+2005);
str1=strcat(str1,str2);
str2='年(总共';
str1=strcat(str1,str2);
str2=num2str(i);
str1=strcat(str1,str2);
str2='年)中国人口总数预测曲线';
tstring=strcat(str1,str2);
title(tstring)
xlabel('(时间单位：年)');
ylabel('(人数单位：万)');
end

[ 来自科研家族皇家数理科学协会 ]

回复此楼

» 猜你喜欢

之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？已经有5人回复
博士读完未来一定会好吗已经有29人回复
博士申请都是内定的吗？已经有5人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有12人回复
投稿精细化工已经有4人回复
高职单位投计算机相关的北核或SCI四区期刊推荐，求支招！已经有4人回复
导师想让我从独立一作变成了共一第一已经有9人回复
读博已经有4人回复
JMPT 期刊投稿流程已经有4人回复
心脉受损已经有5人回复

» 本主题相关价值贴推荐，对您同样有帮助:

1楼2013-05-12 14:34:21

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

hubery.zhu

金虫 (正式写手)

应助: 8 (幼儿园)
金币: 1269.9
散金: 232
红花: 22
帖子: 671
在线: 796小时
虫号: 2421584
注册: 2013-04-17
性别: GG
专业: 人工智能与知识工程

★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
xzhdty: 金币+1, 谢谢参与 2013-05-12 19:04:04

d1()是一个函数吗，如果是的话，你得把d1()这个函数的代码也复制上来。下面是我根据你给的程序改出来的程序：
clear
clc

% d1=zeros(55,1);
d1(1)=12627;
d1(2)=128453;
d1(3)=129227;
d1(4)=130000;
d1(5)=130756;
hold on
plot(2001:2005,d1(1:5));
plot(2001:2005,d1(1:5),'ro')
i=input('你想预测今后(从2005年开始)多少年的中国人口数据？请输入参数：');
while i<0
disp('历史还需要预测吗？请您重新输入大于0的数字：')
i=input('');
end
if i==0
plot(2005,d1(5),'*')
else
y(2)=d1(2);
y(3)=d1(3);
y(4)=d1(4);
y(5)=d1(5);
n=d1(1:end);
for k=5:i+5
      x=k-4:k;
      y=[y(2),y(3),y(4),y(5),n(end)];
      p=polyfit(x,y,2);
      n(k+1)=polyval(p,k+1);
end
t=2005:2005+i;
plot(t,n(6:end),'*-')
str1='2005年到';
str2=num2str(i+2005);
str1=strcat(str1,str2);
str2='年(总共';
str1=strcat(str1,str2);
str2=num2str(i);
str1=strcat(str1,str2);
str2='年)中国人口总数预测曲线';
tstring=strcat(str1,str2);
title(tstring)
xlabel('(时间单位：年)');
ylabel('(人数单位：万)');
end

赞一下

回复此楼

耐得住寂寞，抵的住诱惑，拥得了繁华!

2楼2013-05-12 15:59:31

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主二哥哥33 的主题更新

返回列表