声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 990|回复: 9

[编程技巧] 求助:一个关于自相关的问题

[复制链接]
发表于 2009-1-6 14:21 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
有一组二进制的随机数,以txt格式储存,大小有1.9M,怎样能求自相关并作图?
我用下面的方法试验过,但是运行3个小时左右才能出来(我的机器的cpu是core2duo P8400 2.26GHz  2.27GHz)

d=dlmread('ransuu2.txt');  %ransuu2.txt是储存随机数的文件
x=d(0:1:2000);                 %2000以上无法运行
[a,b]=xcorr(x,'unbiased');
plot(b,a)

望给与指教,谢谢.
回复
分享到:

使用道具 举报

发表于 2009-1-6 16:21 | 显示全部楼层

回复 楼主 takagigi 的帖子

x=d(0:1:2000);  ???
index可以0 ??
发表于 2009-1-6 16:36 | 显示全部楼层
没你的资料, 用rand试, 没问题!?
x=rand(1,100000); [a,b]=xcorr(x,'unbiased'); plot(b,a)

评分

1

查看全部评分

 楼主| 发表于 2009-1-7 12:29 | 显示全部楼层
非常感谢您的指教,是否应该这样做:
d=dlmread('ransuu2.txt');  
x=rand(1,100000);
[a,b]=xcorr(x,'unbiased');
plot(b,a)
目前正在运算中.
发表于 2009-1-7 14:25 | 显示全部楼层
d=dlmread('ransuu2.txt');  
x=d(1:10000);  %% depend on you
[a,b]=xcorr(x,'unbiased');
plot(b,a)

[ 本帖最后由 ChaChing 于 2009-1-7 14:34 编辑 ]
 楼主| 发表于 2009-1-8 12:59 | 显示全部楼层
反应慢的原因找到了,是我的系统(vista)的问题,用朋友的机器(xp)就很快出来了。谢谢您的指导。
发表于 2009-1-9 12:45 | 显示全部楼层

回复 板凳 ChaChing 的帖子

麻烦解释一下为啥噪声的自相关中间会有一个较强的分量
发表于 2009-1-9 14:22 | 显示全部楼层

回复 7楼 wwbeyondww 的帖子

千万勿指名回覆! 个人懂的很少啊!:@L

[ 本帖最后由 ChaChing 于 2009-1-9 14:40 编辑 ]
发表于 2009-1-9 14:24 | 显示全部楼层

回复 8楼 ChaChing 的帖子

晕,以为你知道哇。
哎,最近正被类似问题困扰着
发表于 2009-1-9 14:48 | 显示全部楼层

回复 9楼 wwbeyondww 的帖子

即然是自相关又是噪声资料, 应该就是lag t=0相关性最高吧!
我是从自相关的定义猜的! 不知是否对!?
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-9-22 21:28 , Processed in 0.053234 second(s), 19 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表