数学建模社区-数学中国

标题: 请高手帮我看看这个GM（1：N）程序错在哪啊？谢谢了 [打印本页]

作者: 武器159951 时间: 2010-9-12 14:26
标题: 请高手帮我看看这个GM（1：N）程序错在哪啊？谢谢了

请高手帮我看看这个GM（1：N）程序错在哪啊？我看了好久都不知道哪错了，改了一些还是错的，请大家帮帮忙，谢谢。
[img][/img]
??? Error using ==> horzcat
All matrices on a row in the bracketed expression must have the
same number of rows.
这些是他给出的错误提示。软件用的是Matlab的。
程序的目的是要求解出时间相应函数和对以后的一个预测

clc,clear
x10=[272.81 286.64 311.89 324.76 337.07 351.81 390.85 466.75 490.83 545.46 648.3 696.54 781.66 893.77 1114.32 1519.23 1990.86 2499.43 2957.55 3438.79 3801.09 4188.73 4771.17 5210.12 5741.03
];
x20=[211.05 221.21 236.1 244.34 249.32 255.32 275.37 325.63 336.02 364.38 433.05 466.18 505.6 550.64 677.39 902.38 1148.45 1419.41 1596.72 1774.02 1871.89 1984.64 2207.63 2403.18 2622.45
];
x30=[50.76 53.83 65.69 69.84 74.44 82.97 98.22 121.59 135.12 159.48 187.89 200.73 241.82 309.07 402.77 579.03 794.8 1020.2 1292.11 1592.74 1855.36 2129.6 2486.86 2728.94 3038.9
];
n=length(x10);
x11=cumsum(x10)
x21=cumsum(x20)
x31=cumsum(x30)
for i=2:n
z11(i)=0.5*(x11(i)+x11(i-1));
end
B=[-z11(2:n)',x21(2:n)',x31(2:n)];
Y=x10(2:n)';
u=B\Y
x=dsolve('Dx+a*x=b2*x2+b3*x3','x(0)=x0');%求解符号常微分方程
x=subs(x,{'a','b2',''b3','x0','x2','x3'},{u(1),u(2),u(3),x10(1),'x21','x31'});
digits(6),x=vpa(x);%确定符号数值的有效数字位数和给出数值型符号结果
x=**(x)%运用各种指令化简符号表达式
x=subs(x,{'t','x21','x31'},{[0:n-1],x21(1:n),x31(1:n)})
xhat=[x(1),diff(x)]%求导数、差分和近似微分
epsilon=x10-xhat
delta=abs(epsilon./x10)%模

作者: zhou625315237 时间: 2010-9-12 14:32
好复杂呀，不懂

作者: 雨文diudiu 时间: 2010-9-12 14:47
这是什么？

作者: tom1990 时间: 2010-9-12 14:53
qinagdaldkj;ld

作者: lizejun77 时间: 2010-9-12 15:00
我基本上是采用看英语文章的办法，先泛读，再精读，再一句一句看，最后再提纲挈领，总算是明白一点了，当然，也可能还是领悟错了。最后要说的一句话是：楼主，你很牛叉，希望你不是真的有病。

作者: jimmy900210 时间: 2010-9-12 15:03
万恶过热鬼地方快乐到死见风使舵看风景；韩国附近

作者: 转身孤单 时间: 2010-9-12 16:11
灰色预测，恩……

作者: jimmy900210 时间: 2010-9-12 18:40

作者: jcyou 时间: 2010-9-12 20:00
我基本上是采用看英语文章的办法，先泛读，再精读，再一句一句看，最后再提纲挈领，总算是明白一点了，当然，也可能还是领悟错了。最后要说的一句话是：楼主，你很牛叉，希望你不是真的有病。

作者: hbchina99 时间: 2010-9-13 08:00
鉴定完毕！

作者: ask4kevin 时间: 2010-9-13 12:00
看起来好~~像啊~~~~~

作者: kkly 时间: 2010-9-13 15:00
声明一下：本人看贴和回贴的规则，好贴必看，精华贴必回。

作者: amocilin 时间: 2010-9-13 20:00
强烈支持。楼主万岁

作者: pankg 时间: 2010-9-14 08:00
我基本上是采用看英语文章的办法，先泛读，再精读，再一句一句看，最后再提纲挈领，总算是明白一点了，当然，也可能还是领悟错了。最后要说的一句话是：楼主，你很牛叉，希望你不是真的有病。

作者: franklb888 时间: 2010-9-14 12:00
声明一下：本人看贴和回贴的规则，好贴必看，精华贴必回。

作者: millerwyj 时间: 2010-9-14 15:00
我来了~~~~~~~~~ 闪人~~~~~~~~~~~~~~~~

作者: zr1052310 时间: 2010-9-14 20:00
我来了~~~~~~~~~ 闪人~~~~~~~~~~~~~~~~

作者: chunxiaozhim 时间: 2010-9-15 08:00
留个脚印```````

作者: hehuan_hh 时间: 2010-9-15 12:00
来报道!!!!!!!!!!!

作者: wzzh 时间: 2010-9-15 15:00
鉴定完毕！

作者: vicky9269 时间: 2010-9-15 20:00
哦～～

作者: mjw3918 时间: 2010-9-15 20:00
我基本上是采用看英语文章的办法，先泛读，再精读，再一句一句看，最后再提纲挈领，总算是明白一点了，当然，也可能还是领悟错了。最后要说的一句话是：楼主，你很牛叉，希望你不是真的有病。

作者: lzz50648483 时间: 2010-9-16 08:00
顶顶更健康，越顶吃的越香。

作者: qian7981 时间: 2010-9-16 12:00
鉴定完毕！

作者: leuse 时间: 2010-9-16 15:00
楼主，你写得实在是太好了。我惟一能做的，就只有把这个帖子顶上去这件事了

作者: tengfei110 时间: 2010-9-16 20:00
留个脚印```````

作者: alair001 时间: 2012-1-13 17:42
没有体力啦，资料能发给我一份吗？我的邮箱是18633525948圈163邮箱，谢啦

欢迎光临数学建模社区-数学中国 (http://www.madio.net/)