快捷通道

QQ登录

只需要一步，快速开始

  注册地址  找回密码

数学建模社区-数学中国

留学>>
返回列表 发新帖
查看: 8640|回复: 3
打印 上一主题 下一主题

[代码资源] 求：用MATLAB求自相关系数

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

5

主题

5

听众

58

积分

升级  55.79%

  • TA的每日心情
    无聊
    2012-8-28 02:10

    • 签到天数: 8 天

    [LV.3]偶尔看看II

    自我介绍
    喜爱建模
    跳转到指定楼层
    #
    发表于 2012-7-13 11:35 |只看该作者 |正序浏览
    |招呼Ta 关注Ta
    2体力
    下面这个程序为什么运行出错呢7 s, T) z! w( U4 f) l) H. e' f
    function r=cor_self(x,p,q)
    / b( y" T& l* L5 I& y%连续计算一个序列(列向量）的自相关系数，其中滞后时间k由p变化到q
    , u+ d! |! m; N! ~  G: D/ tclc;close all;
    ; G- V3 `4 i7 y3 N4 dn=size(x,1);
    0 i! A2 Q% A% P8 i; W2 w. S. A6 pfor k=p:q  |3 N# r, @% o2 y  i
        sum=0;
    : u% V4 m; ^4 \4 K; ?! V    for i=1:n-k* g( m& V% A& `& l- r& ]# Y+ \
            sum=sum+(x(i)-mean(x(1:n-k)))*(x(i+k)-mean(x(1+k:n)));
    $ G1 L, n$ H$ q3 H/ a    end% a! e0 Y' v! A* j6 }
        r(k+1)=sum/std(x(1:n-k))/std(x(1+k:n))/(n-k);( B- b4 j) T+ m. x3 A* c. t' \
    end
    % M3 U: `8 n2 X1 v. W" nk=[p:q];
    ! D( [+ m  ~3 l: _( ~plot(k,r,'-*');) |- s6 [* y; Z9 ^) [& ~
    refline(0,0);
    2 |% y1 r/ O0 E" I& O8 U$ Jend7 v/ S* f; R& I0 D
    运行显示8 x" U: |# O1 b0 i" N4 N  `
    Error using ==> plot
    9 y2 T8 S* C' L2 r+ u7 UVectors must be the same lengths.( }/ O- Q6 X& z
    + A8 `8 L5 L2 Y0 a/ T
    Error in ==> cor_self at 13
    7 M0 M9 w# g3 |+ f% E3 eplot(k,r,'-*');

    zan
    function, 系数, close, Error, 程序
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信

    相关帖子
    回复

    使用道具 举报

    zjqylcy 实名认证       

    7

    主题

    5

    听众

    1205

    积分

    宇宙无敌霹雳大头衔

    升级  20.5%

  • TA的每日心情
    奋斗
    2013-5-18 15:14

    • 签到天数: 101 天

    [LV.6]常住居民II

    自我介绍
    男～～

    新人进步奖

    群组: 杭州联盟

    群组: 小草的客厅

    群组: C 语言讨论组

    群组: 数学建模

    群组: 数学趣味、游戏、IQ等
    3#
    发表于 2012-7-16 10:49 |只看该作者 |招呼Ta 关注Ta
    我看不太明白的是   你前面K是一个变 后面K又成了一个矩阵   1 z: {' ?# E6 I4 ?# S- T3 m$ Y
      naf一般就是分母为0了
    回复

    使用道具 举报

    梓爱        

    5

    主题

    5

    听众

    58

    积分

    升级  55.79%

  • TA的每日心情
    无聊
    2012-8-28 02:10

    • 签到天数: 8 天

    [LV.3]偶尔看看II

    自我介绍
    喜爱建模
    2#
    发表于 2012-7-15 10:13 |只看该作者 |招呼Ta 关注Ta
    zjqylcy 发表于 2012-7-13 16:31
    ! g& S- ?1 V& ?k r 长度不一样   应该是r多了一位   可以进workspace里面看看
    ) v: l2 _$ J. ^4 P/ W
    这个后来我又试了下，是矩阵的问题，x的输入应该是列向量的形式，不过k的范围是怎么确定的呀，我输入的是从0到12，可是出来的结果第十一和第十二的数显示为naf，帮忙解释下吧，谢谢了~~
    回复

    使用道具 举报

    zjqylcy 实名认证       

    7

    主题

    5

    听众

    1205

    积分

    宇宙无敌霹雳大头衔

    升级  20.5%

  • TA的每日心情
    奋斗
    2013-5-18 15:14

    • 签到天数: 101 天

    [LV.6]常住居民II

    自我介绍
    男～～

    新人进步奖

    群组: 杭州联盟

    群组: 小草的客厅

    群组: C 语言讨论组

    群组: 数学建模

    群组: 数学趣味、游戏、IQ等
    1#
    发表于 2012-7-13 16:31 |只看该作者 |招呼Ta 关注Ta
    k r 长度不一样   应该是r多了一位   可以进workspace里面看看
    回复

    使用道具 举报

    返回列表 发新帖
    高级模式
    B Color Image Link Quote Code Smilies
    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩

    • 电话咨询

    • 04714969085

    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2026-10-9 03:58 , Processed in 0.425509 second(s), 71 queries .

    回顶部