该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

N=length(data);

max_d=200;%the maximum value of the time delay

sigma=std(data);%calcute standard deviation s_d

for t=1:max_d

s_t=0;

delt_s_s=0;

for m=2:5

s_t1=0;

for j=1:4

r=sigma*j/2;

data_d=disjoint(data,N,t);%将时间序列分解成t个不相交的时间序列

N_d=size(data_d);

s_t3=0;

for i=1:t

i;

Y=data_d(i,:);

C_1(i)=correlation_integral1(Y,N_d,r);%计算C(1,N_d,r,t)

X=reconstitution(Y,N_d,m,t);%相空间重构

N_r=N_d-(m-1)*t;

C_I(i)=correlation_integral(X,N_r,r);%计算C(m,N_r,r,t)

s_t3=s_t3+(C_I(i)-C_1(i)^m);%对t个不相关的时间序列求和

end

s_t2(j)=s_t3/t;

s_t1=s_t1+s_t2(j);%对rj求和

end

delt_s_m(m)=max(s_t2)-min(s_t2);%求delt S(m,t)

delt_s_s=delt_s_s+delt_s_m(m);%delt S(m,t)对m求和

s_t0(m)=s_t1;

s_t=s_t+s_t0(m);%S对m求和

end

s(t)=s_t/16;

delt_s(t)=delt_s_s/4;

s_cor(t)=delt_s(t)+abs(s(t));

end

fid=fopen('result.txt','w');

fprintf(fid,'%f %f %f %f/n',t,s(t),delt_s(t),s_cor(t));

fclose(fid);

t=1:max_d;

plot(t,s,t,delt_s,'.',t,s_cor,'*')

function C_I=correlation_integral(Y,M,r)

% 此函数用于计算关联积分,取无穷范数

% Y为重构的相空间

% M为相空间中点的个数

% r为搜索半径

% Y为输出,是M*m维矩阵

C=0;

for i=1:M-1

for j=i+1:M

d=norm((Y(:,i)-Y(:,j)),inf);      % 计算状态空间中每两点之间的距离,取无穷范数

if r-d>=0

C=C+1;

end

end

end

C_I=2*C/(M*(M-1));

%%%%%%%%%%%%%%%%%%%%%%%%%

??? Error using ==> mtimes

Inner matrix dimensions must agree.

Error in ==> correlation_integral1 at 17

C_I=2*C/(M*(M-1));

Error in ==> c_c at 18

C_1(i)=correlation_integral1(Y,N_d,r);%计算C(1,N_d,r,t)

这是怎么回事啊,望各位指点

matlab时间延迟怎么写,大家帮我看一下这个c_c方法时间延迟的程序运行结果??...相关推荐

  1. 一拳超人—写给码农看的数据库优化方法:everything is a file

    Everything is a file是UNIX世界中的名言,指的是将系统中各种资源都认为是文件,通过字符串读写方式实现统一接口调用. 起这个标题仅为蹭一下名言热度,我们真正要聊的是数据库世界中的& ...

  2. 程序员:为什么几个月后我自己写的代码也看不懂了?

    写在开始的"注释很重要","注释很重要","注释很重要",重要的事情说三遍. 程序员们大多都会有过这样的经历,就是要看别人写的代码.比如说公 ...

  3. 程序员必看,学习的方法

    程序员必看,学习的方法 0.引言 在学习编程语言的路上,发现初学者(包括俺)问 得最多的问题莫过于"如何学习编程?Java该如何学习?".类似的问题回答多了,难免会感觉厌烦,就萌生 ...

  4. 写给后端开发看的安卓入门

    写给后端开发看的安卓入门 写给后端开发看的安卓入门 可以看也可以不看的序 安卓的一些入门知识 ADT SDK API Intent Android的开发IDE Eclipse IDEA Android ...

  5. 小娴的男友小旭不幸患了一种怪病,这种怪病吞噬了他的大部分记忆,同时让他突然间不会书写符合正确语序的英文。神奇的是,虽然他写出的句子看上去杂乱无章,不过经过仔细分析可以发现,如果把单词的顺序倒过来,语法

    题目描述 小娴的男友小旭不幸患了一种怪病,这种怪病吞噬了他的大部分记忆,同时让他突然间不会书写符合正确语序的英文.神奇的是,虽然他写出的句子看上去杂乱无章,不过经过仔细分析可以发现,如果把单词的顺序倒 ...

  6. 写给理工科人看的乐理(五)调性

    原文地址为: 写给理工科人看的乐理(五)调性 这一讲是乐理中有关音调知识的终曲,不过我们还是从历史讲起. 一.乐理发展史 · 之四 前面的历史课我们了解到"确定相对音高"的律制系统 ...

  7. 报错信息找不到原因,麻烦帮我看一下

    今天辛辛苦苦写了一上午的代码,突然报错了,找不到原因,麻烦帮我看一下,以下是报错信息:java.lang.NoMoneyException: KFC CrazyThursday need $50; 挺 ...

  8. 三本书帮你看淡元宇宙

    前言 有的人说了,这是在凑热点,其实并不是.不是我去找的元宇宙,而是元宇宙来找上的我们.那元宇宙到底靠不靠谱呢? 元宇宙(Metaverse)一词起源于29年前科幻小说<雪崩>,小说中将人 ...

  9. php输出后以密文显示,加密解密,_求大神帮看一下,如何根据下面的程序,把明文生成密文谢谢,加密解密 - phpStudy...

    求大神帮看一下,如何根据下面的程序,把明文生成密文谢谢 求大神帮看一下,如何根据下面的程序,把明文生成密文谢谢 def parserLicense(self,encCode): data = base ...

最新文章

  1. Linux控制删除深度!除find ./ -maxdepth 1 -name *.c | xargs rm -rf
  2. UUID介绍与生成的方法
  3. 用D触发器构造边沿触发器
  4. VB 字节数组和字符串的转换问题 (StringByte)
  5. mysql重启电脑提示145_MYSQL错误:Can’t open file: ‘#215;#215;#215;.MYI’ (errno: 145)修复方法...
  6. python有序队列_【python】collections模块(有序字典,计数器,双向队列)
  7. 贪心算法-Best cow line-字典序问题
  8. 腾讯视频如何退出登陆
  9. C/C++混淆点-转义字符
  10. allegro设置快捷键
  11. 【教学类-10-01】20221025《空心图案4*2-有重复》( 随机图案拼贴)(大班主题《动物花花衣》)
  12. 机械制图--零件图--过渡线
  13. UCHome全面大解析【二】----基本体系结构
  14. AppWeb认证绕过漏洞(CVE-2018-8715)
  15. js 主动触发 a 标签 href 链接写法,
  16. 厉害了,Python也能使用动态链接库
  17. Kolmogorov–Smirnov test
  18. JVM —— Java 对象占用空间大小计算
  19. 软件设计师--系统开发与运行
  20. 约瑟夫环问题的几种解法

热门文章

  1. 知道两点经纬度,计算双方距离与方向的原理
  2. python 英文关键词提取_python 利用jieba.analyse进行 关键词提取
  3. 不评价别人的生活,是一个人最基本的修养
  4. 助力低碳出行 | 基于ACM32 MCU的电动滑板车方案
  5. python中x%y是什么意思_python是什么
  6. 10部关于数学的顶级纪录片
  7. 那些年我吃过的豆子【记录ing】
  8. 『教师节』程序猿用文心大模型带你一键加速祝福,祝老师们节日快乐
  9. 关于c++中引用的基本用法
  10. java校验中文姓名合法性(需要支持生僻字与分隔符·)