数学建模社区-数学中国
标题:
灰色预测程序哪里出错
[打印本页]
作者:
一路追逐
时间:
2015-7-23 17:14
标题:
灰色预测程序哪里出错
a=[390.6,412,320,559.2, 380.8,542.4,553,310,561,...
, k( V B& \. a: a
300,632,540,406.2,313.8,576,587.6,318.5]';
' m9 k. X& U8 _8 q8 e7 O1 n
t0=find(a<=320);n=length(t0);
S% P$ ]% ?: E3 C6 ~3 |- R# u
t1=cumsum(t0); %累加运算
2 u/ M+ z; C- I+ H1 C
B=[-0.5*(t1(1:end-1)+t1(2:end)),ones(n-1,1)];Y=t0(2:end);
3 B. @( A: n) R$ i
r=B\Y
6 Z W) H) o% ~1 u) ]
y=dsolve('Dy+a*y=b','y(0)=y0');
6 a) J. T7 B5 B+ i
y=subs(y,{'a','b','y0'},{r(1),r(2),t1(1)});
) Q, T1 F, X4 h" b
yuce1=subs(y,'t',[0:n+1])
9 H6 D7 R2 ^( b' l; i
%为提高预测精度，先计算预测值，再显示微分方程的解
" o" e& s8 |3 U; G
y=vpa(y,6) %其中的6 表示显示6 位数字
5 i3 D6 I/ e2 D. L& P+ F$ i* K
yuce=diff(yuce1); %作差分运算，进行数据还原
, s% R7 r% h' w6 } [3 j! @) I% A* F/ o
yuce=[t0(1),yuce]
! {2 G; I3 M" I
yuce_new=yuce(n+1:end) %求得的两个预测值
4 P& e; Z- I% C, r$ @- [
7 }2 m* L. z. V0 J6 n. f) J
为什么在做差分运算时yuce=diff(yuce1）得到的结果是yuce =[ 0, 0, 0, 0, 0, 0, 0]
! d8 h8 o7 o' t% `: r2 o1 V
求赐教
Q# X: y! p! L: [/ U
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5