QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 2512|回复: 0
打印 上一主题 下一主题

[代码资源] 灰色预测程序哪里出错

[复制链接]
字体大小: 正常 放大

7

主题

11

听众

48

积分

升级  45.26%

  • TA的每日心情
    慵懒
    2016-4-15 21:13
  • 签到天数: 23 天

    [LV.4]偶尔看看III

    社区QQ达人

    跳转到指定楼层
    1#
    发表于 2015-7-23 17:14 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta
    a=[390.6,412,320,559.2, 380.8,542.4,553,310,561,...; l: k& D* q, T/ a( d5 B
    300,632,540,406.2,313.8,576,587.6,318.5]';$ p% h8 n6 A7 a  A* I& f( J
    t0=find(a<=320);n=length(t0);
    " h& |. O% ]1 q1 \& Z, W( D4 | t1=cumsum(t0); %累加运算
    9 n2 }% G6 W- G5 k, y: U: a9 jB=[-0.5*(t1(1:end-1)+t1(2:end)),ones(n-1,1)];Y=t0(2:end);4 Z+ j* M) Z, d. z4 _
    r=B\Y
    3 I3 y1 n! I# W. }$ a3 o' E9 R' v y=dsolve('Dy+a*y=b','y(0)=y0');$ [9 R! \9 \8 H3 |2 |
    y=subs(y,{'a','b','y0'},{r(1),r(2),t1(1)});
    4 e8 U# ^" W  }/ ^" E yuce1=subs(y,'t',[0:n+1])
    ! i3 P4 W. t, u0 o% J %为提高预测精度,先计算预测值,再显示微分方程的解) P: w, {( L3 Z
    y=vpa(y,6) %其中的6 表示显示6 位数字
    ' r; @4 S# T( j5 K0 c4 a5 o0 s* [yuce=diff(yuce1); %作差分运算,进行数据还原; I8 T( \0 W7 C
    yuce=[t0(1),yuce]2 Q2 b& c; X1 S- r) }& N! I
    yuce_new=yuce(n+1:end) %求得的两个预测值: H2 u* C/ `2 @' a: L! V
    4 ]! f/ [; _+ C2 `& F
    为什么在做差分运算时yuce=diff(yuce1)得到的结果是yuce =[ 0, 0, 0, 0, 0, 0, 0]1 O! W  R# i# `. {4 ~, e
    求赐教
    8 I6 y$ R0 I6 o
    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-4-11 13:15 , Processed in 0.445497 second(s), 54 queries .

    回顶部