快捷通道

QQ登录

只需要一步，快速开始

  注册地址  找回密码

数学建模社区-数学中国

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

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

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

5

主题

5

听众

58

积分

升级  55.79%

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

    • 签到天数: 8 天

    [LV.3]偶尔看看II

    自我介绍
    喜爱建模
    跳转到指定楼层
    1#
    发表于 2012-7-13 11:35 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta
    2体力
    下面这个程序为什么运行出错呢
    ' h$ G, M% I& U! xfunction r=cor_self(x,p,q)4 q" ~5 U. y% K% Q9 k
    %连续计算一个序列(列向量）的自相关系数，其中滞后时间k由p变化到q& A$ {/ E, l) n# U& J' T
    clc;close all;$ `3 v% [* H  I/ r: O$ p0 p2 l
    n=size(x,1);
    ' s8 w% r, X7 W% S5 b2 b4 n5 @; C5 `for k=p:q
    / N: q7 ~' `% v: T; j5 P    sum=0;' l! `: G% c/ Z( ~1 Z( }
        for i=1:n-k
    + x5 [# a4 z" s4 _, ?, ]- N) |        sum=sum+(x(i)-mean(x(1:n-k)))*(x(i+k)-mean(x(1+k:n)));
    ! ~) N2 ?3 V/ x# N2 h& ]6 _    end9 `2 W5 ~6 a4 g8 n0 p! U
        r(k+1)=sum/std(x(1:n-k))/std(x(1+k:n))/(n-k);
    - B$ t' [! _! Send: `9 P0 r7 I5 i/ n% B9 P: ^5 k
    k=[p:q];7 Z. y" x7 ^' D2 @: `
    plot(k,r,'-*');) Y* F# g0 v( j5 h) i7 P% P
    refline(0,0);
    9 Y1 `* t  B8 D3 @: c& D# Wend5 i, {. \3 t; H
    运行显示) m+ W8 B/ I' n0 _0 [5 q6 k) |
    Error using ==> plot
    4 X) a! d: \  N( p4 YVectors must be the same lengths.0 |2 M. a( \6 a4 f

    % O* q( q5 A! K+ y% c% e/ B+ P" wError in ==> cor_self at 13! p. b# L( A" k7 H4 [" ~
    plot(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等
    2#
    发表于 2012-7-13 16:31 |只看该作者 |招呼Ta 关注Ta
    k r 长度不一样   应该是r多了一位   可以进workspace里面看看
    回复

    使用道具 举报

    梓爱        

    5

    主题

    5

    听众

    58

    积分

    升级  55.79%

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

    • 签到天数: 8 天

    [LV.3]偶尔看看II

    自我介绍
    喜爱建模
    3#
    发表于 2012-7-15 10:13 |只看该作者 |招呼Ta 关注Ta
    zjqylcy 发表于 2012-7-13 16:31
    4 x& [. p% h" B, ?2 r, m- Mk r 长度不一样   应该是r多了一位   可以进workspace里面看看
    8 _& k& u8 v* _0 }3 v2 B( z- P* P9 U
    这个后来我又试了下，是矩阵的问题，x的输入应该是列向量的形式，不过k的范围是怎么确定的呀，我输入的是从0到12，可是出来的结果第十一和第十二的数显示为naf，帮忙解释下吧，谢谢了~~
    回复

    使用道具 举报

    zjqylcy 实名认证       

    7

    主题

    5

    听众

    1205

    积分

    宇宙无敌霹雳大头衔

    升级  20.5%

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

    • 签到天数: 101 天

    [LV.6]常住居民II

    自我介绍
    男～～

    新人进步奖

    群组: 杭州联盟

    群组: 小草的客厅

    群组: C 语言讨论组

    群组: 数学建模

    群组: 数学趣味、游戏、IQ等
    4#
    发表于 2012-7-16 10:49 |只看该作者 |招呼Ta 关注Ta
    我看不太明白的是   你前面K是一个变 后面K又成了一个矩阵   : y) V5 A! o% ]% E
      naf一般就是分母为0了
    回复

    使用道具 举报

    返回列表 发新帖
    高级模式
    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-6-14 21:00 , Processed in 0.470109 second(s), 68 queries .

    回顶部