应用Elman神经网络的混沌时间序列预测-论文

应用Elman神经网络的混沌时间序列预测-论文
应用Elman神经网络的混沌时间序列预测-论文

华 东 理 工 大 学 学 报 Journal o f Eas t C hi na Un iv ersi ty of S cience and Tech no l ogy V o l.28增刊

2002-09

收稿日期35

作者简介张兴会(63),男,河北孟村人,副教授,博士生,研究方

向为智能控制和预测控制。

文章编号:1006-3080(2002)NN -S 0-0030-04

应用

Elman 神

经网络的混沌时间序列预测

张兴会1*, 刘 玲1, 陈增强2, 袁著祉2

(1.天津技术师范学院计算机系,天津300222; 2.南开大学自动化系,天津300071)

摘要:利用改进的Elma n 神经网络对3个典型的混沌时间序列在不同的噪声水平下进行预测,探讨了神经网络学习与泛化之间的关系,通过试凑法给出了Elman 最优的隐节点个数。并利用3种指标对预测结果进行了评估,结果显示Elma n 网络对混沌时间序列预测的良好特性。

关键词:混沌时间序列;预测;改进型Elman 神经网络;径向基函数神经网络

中图分类号:TP 18

文献标识码:A

Prediction of Chaos Time Series Using Elman Neural Networks

Z H AN G X ing -hui 1*

, LIU Ling 1

, C HEN Zeng-qia ng 2

, Y U AN Zhu-zhi

2

(1.Department of Computer ,Tia njin University of Technology and Educa tion ,Tianjin 300222,China ;

2.Depar tment of A uto mation,N ankai University,Tia njin 300071,China )Abstr act :This paper uses the improved Elman neural networks to predict three typical chaos time se-ries under different noise conditions .It also discusses the relationship between lear n ing and generalization of the neura l networks and gives the optimal number of Elma n 's hidden layer units.In addition,the pre-diction results a re evalua ted by three targets which show the perfect performance of Elman networks in the prediction of the chaos time series .

Key words :chaos time series;prediction ;improved Elman neural network;radial basis function neu-ral network

混沌现象是非线性系统研究的一个热点。混沌时间序列预测在工业、经济等复杂系统都有着重要的意义,其应用领域可见于经济预测和调整、金融分析、天气预报、地震预测等方面[1~2]

。神经网络技术已用于对混沌时间序列的预测研究,文献[3]指出该技术可以从被噪声污染的某些混沌系统中提取混沌吸引子。文献[4]探讨了BP 网络和回归网络对混沌系统的预测,仿真结果表明回归网络优于BP 网络,但是对其网络结构和隐含节点的选择未进行讨论,使得预测结果不如文献[5]的理想。文献[6~7]分别探讨了模糊神经网络、小波网络在混沌时间序列预测中的应用问题,但是它们算法较复杂,预测精度较低,且大多研究结果并未考虑观测噪声。文献[5]应

用径向基函数神经网络(RBF 网络)得到了较快的收敛速度、较高的预测精度,并分析了有噪声情况下的拟合误差和预测误差情况。

本文采用改进型回归记忆Elma n 网络,对Lo-gistic map 、Henon map 、Mackey-Glass 3种混沌时间序列进行了预测仿真研究。结果表明,改进型Elman 网络优于文献[4]提出的方法。在有噪声时,Elman 网络模型的拟合误差和预测误差均较小,明显优于文献[5]提出的方法,取得较好的效果。

1 改进型Elman 神经网络

网络是在B 网络基本结构的基础上,通过存储内部状态使其具备映射动态特性的功能,从而使系统具有适应时变特性的能力,在复杂过程的建模方面,有着广泛的应用前景。神经网络

30

D OI :10.14135/https://www.360docs.net/doc/8d4836417.html, k i .1006-3080.2002.s1.008

:2002-0-2:19-Elman P Elman

除输入层、隐层以及输出层单元外,还有一个独特的结构单元。与通常的多层前馈网络相同,输入层单元仅起信号传输作用,输出层单元起线性加权和作用,隐层单元可有线性或非线性激发函数。而结构单元则用来记忆隐层单元前一时刻的输出值,可以认为是一个时延算子。改进型Elman 神经网络[8]的结构单元具有一个固定增益T 的自反馈连接,改进型El-ma n 神经网络利用标准BP 学习算法辨识高阶动态系统,是解决高阶系统辨识问题的更好方案。改进型的Elman 神经网络结构如图1

所示。

图1 改进型Elma n 神经网络结构Fig .1

 Improved Elma n neura l network 设网络的外部输入为u (k -1)∈R r ,预测输出为y (k)∈R m

,若记隐层的输出为x (k)∈R n

,则有如下非线性状态空间表达式成立:

x (k )=f (W 1

x c (k )+W 2

u (k -1))

x c (k )=x (k -1)+Tx c (k -1)

y (k)=W 3x (k )

(1)

其中W 1

,W 2

,W 3分别是结构单元到隐层、输入层到隐层,以及隐层到输出层的连接权矩阵,T 为自连接反馈增益,f (W 1

x c (k)+W 2

u (k-1))为隐层单元激发函数所组成的非线性向量函数。

改进型Elman 网络可用标准BP 学习算法辨识

高阶系统,具体算法公式为:

Δw 3ij

=Z (y di (k )-y i (k))x j (k) i =1,2,…,m;

j =1,2,…,n

Δw 2

jq

=Z ∑m

i=1(y di (k)-y i (k ))w 3

i j f ′j (W 1

x c (k)+

W 2

u (k -1))u q (k -1)

j =1,2,…,n;q =1,2,…,r

(2)

Δw 1

jl =Z ∑m

=(y di (k )-y i (k))w 3

ij

x j (k )w 1j j =,,…,;=,,…, 其中Z 为网络学习率,y ()为y ()的预测值,而

x j (k)w 1

j l

=f ′j (W 1x c (k )+W 2

u (k -1))x l (k -1)+T

x j (k -1)

w 1jl

(3)

对于混沌系统建模,神经网络的输入变量个数的选择与估计的混沌吸引子的维数有关。文献[9]提出,对于含有噪声的混沌信号,利用有限的隐节点个数,可以得到最小预测误差,提高其泛化能力。就本文所选择的3个典型混沌系统,输入层节点个数为1,采用试凑法[10]确定隐层节点个数分别为:3,3,4。

2 仿真实验设计及结果

2.1 混沌模型的选择

我们研究以下3个混沌系统:

(1)Logistic map

x t =λx t-1(1-x t-1)

(4)

当λ∈(3.5669,4]时,Logistic ma p 呈现混沌

状态,本文取λ=4。

(2)Henon map

x t =1-a x 2

t-1+z t-1

z t =bx t-1

(5)

当a = 1.4,b =0.3时,Henon ma p 系统将处于混段状态。

(3)Mackey-Glass 方程的离散形式

x t =x t-1+a x t-s

1+(x t -s )

c -bx t -1

(6)

其中a =0.2,b=0.1,c=10.0,s=17。

仿真研究将在无任何噪声和有附加噪声两种情况下进行。对于无噪声时,观测数据为:

y t =x t

(7)把y t 看作目标变量。对于有附加噪声的情形,

y t =x t +v t (8)

其中v t 为噪声部分。

v t =U e u t

(9)

式中U 分别取0.01,0.1,0.25,0.5,U 是x t 的

样本标准差,u t 是标准正态分布的随机变量。对于每种模型,取m 个样本,其中前n 个样本用于建模,后m-n 个样本用于预测。对于Logistic

map 模型和Henon map 模型,取m=200,n=180。对于Mackey -Glass 方程,取m =1000,n =900。当存在附加噪声时,对于U 的每一个取值,都取组数据集分别进行建模和预测。

 评价指标

对于每种模型的每次仿真实验,都采用3个统

31

增刊张兴会等:应用El man 神经网络的混沌时间序列预测 i 1

l

12n l 12n

i k di k 202.2

计量来评价模型的拟合程度和预测效果,这3个统计量是RMSE 、RMSP E 和RMSP E /SD 。

为比较训练后网络的拟合能力,计算出如下的误差指标。

RMSE =

1

n ∑n

t=1

(y t -y t )2

(10)

其中y t 是y t 的估计值。

为评价预测的效果,计算一步预测的均方误差的平方根,即

RMSPE =

1

m -n ∑m

t =n +1

(y t -y t )2

(11)

其中y t 是y t 的一步预测值。

还可以通过计算RMSP E 与SD 的比值(RM-SP E /SD )来做进一步的比较。对于无噪声情况,SD

由下式计算:

SD =

1

m -n ∑m

t =n +1

(y t -y -)2

(12)

这里,y -是y t 的样本均值,也就是说,SD 是预

测阶段的样本标准差。RMSP E /SD 的值越接近于零,则估计模型的预测效果越好。

对于存在噪声的情况,SD 由下式计算

SD =

1m -n ∑m

t =n+1

(v t -v -)2(13)

这里,v -是在预测阶段v t 的样本均值。因v 是由式(9)产生的,SD 的值接近于v t 的标准差,即为

U e 。当估计模型的预测效果较好时,估计模型的RMSP E 值应该接近于SD 值,即RMSP E 与SD 的比值应该接近1。2.3 仿真实验结果

采用改进型Elma n 神经网络分别对Logistic map 、Henon map 和Mackey-Glass 方程进行学习,计算出不同噪声水平下(U=0.01,0.1,0.25,0.5)的3个统计量,即RMSE 、RMSPE 和RMSP E /SD ,并与相同噪声条件下对3种混沌时间序列利用

RBF 网络的仿真结果[5]

进行了比较,如表1、表2、表3所示。

表1 Logistic ma p 的仿真结果Table 1 Simula tion results of Logistic map

U

R MSE

Elman

RBF RM SP E

Elman RB F RM SP E /SD

El man RBF 0.01

1.2415×10-40.0033 6.4514×10-40.00361.06071.06250.1

2.0842×10-3

0.0335 2.8044×10-3

0.03601.18121.06210.250.08320.08550.08270.09121.10391.02070.5

0.1637

0.1673

0.1694

0.1802

1.2140

1.0626

表2 Henon map 的仿真结果Table 2 Simula tion r esults of Henon map

U

R MSE

Elman

RBF RM SP E

Elman RB F RM SP E /SD

El man RBF 0.01

1.3611×10-40.0057 1.7856×10-40.00831.13091.18130.1

2.7598×10-3

0.0560 4.2000×10-3

0.08481.12961.19290.250.12140.14110.12660.21051.12871.19780.5

0.2243

0.2805

0.2287

0.4232

1.0220

1.2067

表3 Mackey-Glass 方程的仿真结果Table 3 Sim ulation results of Mackey-Glass ma p

U

R MSE

Elman

RBF RM SP E

Elman

RB F

RM SP E /SD

El man RBF 0.01

6.6726×10-40.0022 2.6223×10-4

0.00231.02361.03250.11.1191×10-3

0.02200.01010.02251.00071.030955555633535

 

32 华 东 理 工 大 学 学 报第28卷

0.20.04010.000.02820.01.081.0120.0.08880.1108

0.0749

0.1118

1.0001

1.0210

从表中可以看出,对于3个模型在不同噪声水平下,采用改进型Elman 网络时的3个评价指标(RMSE 、RMSPE 和RMSP E /SD )小于采用RBF 网络时的相同值,表明改进型Elman 网络在混沌系统建模与预测方面效果比RBF 网络要好,说明该网络有很好的抗噪能力。

3 讨 论

混沌是确定性系统内在的随机性,其内部有着确定的规律性。这种规律性用Elman 神经网络得到

了较好的表达,取得了比其他神经网络方法更好的效果。如何设计网络结构使其性能(学习速度、精度、泛化能力、结构复杂性、计算复杂性)最优,挖掘回归神经网络的运行机理,给出严格的数学证明,证明该类型网络具有最好的预测精度与泛化能力等问题有待进一步工作。参考文献:

[1] 郭 刚,王丽敏,付大鹏.基于混沌理论进行股票市场的多步

预测[J ].信息与控制,2000,29(2):177-180.

[2] 梁志珊,史忠科,戴冠中.应用混沌理论的电力系统短期负荷

预测[J ].控制与决策,1998,13(1):87-90.

[3] Albano A M,Pass aman te A,Hed i ger T,et https://www.360docs.net/doc/8d4836417.html,ing neural

net s t o l o ok for chaos [J].Ph ysica D,1992,58:1-9.

[4] Gencay R ,Li u T .Nonli near mo d eli ng and p rediction wi th

feedfoward and recu rren t networks [J ].Ph ysica D,1997,108:119-134.

[5] 李冬梅,王正欧.基于RBF 神经网络的混沌时间序列预测

[J ].模式识别与人工智能,2001,14(2):231-234.

[6] Magui re L P ,Roche B ,Mcginni ty T M ,et al .Pred i cting a

chaotic t ime seri es using a fuzzy neu ral net work [J ].Infor-mati on sciences,1998,112:125-136.

[7] Cao L Y,Hong Y G,Fang H P,et al.

Predi cti ng chaotic

tim e s eri es wi th wavel et net wo rk s [J ].Ph ysica D,1995,85:225-238.

[8] 孙增圻,张再兴,邓志东.智能控制理论与技术[M ].北京:清

华大学出版社,2000.

[9] Henry L,Ti t us L O,Wang S.Predicti on of noising chaotic

tim e seri es usi ng an opti mal radial basis function neural Net -work [J ].IEEE Trans on Neural Networks ,2001,12(5):1163-1171.

[10] 何述东,瞿 坦,黄献青,等.多层前项神经网络结构的研究进

展[J ].控制理论与应用,1998,15(3):313-319.

(上接第29页

)

图4 模型的泛化能力

Fig .4

 Generalization a bility of the model ○—Desi red ou tpu t;+—Act ual out put

参考文献:

[1] 张 彤,王宏伟,王子才.变尺度混沌优化方法及其应用[J ].

控制与决策,1999,14(3):285-288.

[2] 李 兵,蒋慰孙.混沌优化方法及其应用[J ].控制理论与应

用,1997,14(4):613-615.

[3] 童树鸿,沈 毅,刘志言.遗传算法在模糊系统优化设计中的

应用研究[J ].系统工程与电子技术,2001,23(1):73-76.

33

增刊张兴会等:应用El man 神经网络的混沌时间序列预测

小波神经网络的时间序列预测短时交通流量预测.doc

%% 清空环境变量 clc clear %% 网络参数配置 load traffic_flux input output input_test output_test M=size(input,2); %输入节点个数 N=size(output,2); %输出节点个数 n=6; %隐形节点个数 lr1=0.01; %学习概率 lr2=0.001; %学习概率 maxgen=100; %迭代次数 %权值初始化 Wjk=randn(n,M);Wjk_1=Wjk;Wjk_2=Wjk_1; Wij=randn(N,n);Wij_1=Wij;Wij_2=Wij_1; a=randn(1,n);a_1=a;a_2=a_1; b=randn(1,n);b_1=b;b_2=b_1; %节点初始化 y=zeros(1,N); net=zeros(1,n); net_ab=zeros(1,n); %权值学习增量初始化 d_Wjk=zeros(n,M); d_Wij=zeros(N,n); d_a=zeros(1,n);

d_b=zeros(1,n); %% 输入输出数据归一化 [inputn,inputps]=mapminmax(input'); [outputn,outputps]=mapminmax(output'); inputn=inputn'; outputn=outputn'; %% 网络训练 for i=1:maxgen %误差累计 error(i)=0; % 循环训练 for kk=1:size(input,1) x=inputn(kk,:); yqw=outputn(kk,:); for j=1:n for k=1:M net(j)=net(j)+Wjk(j,k)*x(k); net_ab(j)=(net(j)-b(j))/a(j); end temp=mymorlet(net_ab(j)); for k=1:N y=y+Wij(k,j)*temp; %小波函数 end end

《时间序列分析》第二章 时间序列预处理习题解答

《时间序列分析》习题解答?0?2习题2.3?0?21考虑时间序列10判断该时间序列是否 平稳计算该序列的样本自相关系数 kρ∧绘制该样本自相关图并解释该图形. ?0?2解根据时序图可以看出该时间序列有明显的递增趋势所以它一定不是平稳序列?0?2即可判断该时间序是非平稳序列其时序图程序见后。?0?2 时间序描述程序data example1 input number timeintnxyear01jan1980d _n_-1 format time date. cards 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 proc gplot dataexample1 plot numbertime1 symbol1 cblack vstar ijoin run?0?2?0?2?0?2当延迟期数即k本题取值1 2 3 4 5 6远小于样本容量n本题为20时自相关系数kρ∧计算公式为 number1234567891011121314151617181920time01JAN8001J AN8101JAN8201JAN8301JAN8401JAN8501JAN8601JAN870 1JAN8801JAN8901JAN9001JAN9101JAN9201JAN9301JAN9 401JAN9501JAN9601JAN9701JAN9801JAN99121nkttktknttX XXXXXρ?6?1∧?6?1?6?1≈?6?1∑∑ 0kn4.9895?0?2 注20.05125.226χ接受原假设认为该序列为纯随机序列。?0?2解法三、Q统计量法计算Q统计量即12214.57kkQnρ∑?0?2?0?2?0?2?0?2?0?2?0?2?0?2?0?2?0?2?0?2查表得210.051221.0261χ?6?1由于Q统

最新时间序列分析期末考试B

精品文档 浙江农林大学 2009 - 2010 学年第 二 学期考试卷(A 卷) 课程名称: 应用时间序列分析 课程类别: 必修 考试方式: 闭卷 注意事项:1、本试卷满分100分。 2、考试时间 120分钟。 一、单项选择题(在每小题的四个备选答案中,选出一个正确 答案,并将正确答案的选项填在题后的括号内。每小题2分,共12分) 1. 关于严平稳与(宽)平稳的关系,不正确的为 。 ( ) A. 严平稳序列一定是宽平稳序列 B. 当序列服从正态分布时,两种平稳性等价 C. 二阶矩存在的严平稳序列一定为宽平稳的 D. MA(p)模型一定是宽平稳的 2. 下图为某时间序列的相关检验图,图1为自相关函数图,图2为偏自相关函数图,请选择模型 。 ( ) 图1 图2 学院: 专业班级: 姓名: 学号: 装 订 线 内 不 要 答 题

A. AR(1) B. AR(2) C. MA(1) D. MA(2) 3. 下图中,图3为某序列一阶差分后的自相关函数图,图4为某序列一阶差分后的 偏自相关函数图,请对原序列选择模型。( ) 图3 图4

A.ARIMA(4,1,0) B. ARIMA(0,2,1) C. ARIMA(0,1,2) D.ARI MA(0,1,4) 4. 记B 为延迟算子,则下列不正确的是 。 ( ) A. 0 1B = B. (1)k t t k t X X B X --=- C. 12t t BX X --= D. 11()t t t t B X Y X Y --±=± 5.对于平稳时间序列,下列错误的是 ( ) A.)(212εσεE = B.),(),(k t t k t t y y Cov y y Cov -+= C.k k -=ρρ D.)(?)1(?1k y k y t t +=+ 6.下图为对某时间序列的拟合模型进行显著性水平0.05α=的显著性检验,请选择 该序列的拟合模型 。 ( )

混沌时间序列分析

第四章
混沌时间序列分析
一.相空间重建 二.相关维数 三.最优延迟时间 四. Lyapunov特征指数 五.应用举例
一、相空间重建
1

2

Embedding
Φ
A M System with dynamics f has an attractor A ? M
Z ?d
A is transformed into a set Z ? ?d such that the all the important geometric characteristics of A will be preserved. Lets also assume Φ is invertible.
Ruelle(1981),法国科学家
对m维动力系统:
? x1 = f1 ( x1 , x 2 , .... x m ) ? ? x 2 = f 2 ( x1 , x 2 , ... x m ) ? ? ................... ? ? x m = f m ( x1 , x 2 , ... x m )
( x1 , x2 ,.....xm )
是状态空间坐标
x(t ), x(t + τ), x(t + 2τ),......., x [t + (m ? 1)τ]
3

Phase Space Reconstruction
一个单变量时间序列:x0 , x1 , x2 ,...
?τ = 1.?t ? ?n = 3 ( x0 , x1 , x2 ) ( x1 , x2 , x3 ) ( x2 , x3 , x4 ) ............... ( xn ?1 , xn , xn +1 )
4

(完整版)基于MATLAB的混沌序列图像加密程序

设计题目:基于MATLAB的混沌序列图像加密程序 一.设计目的 图像信息生动形象,它已成为人类表达信息的重要手段之一,网络上的图像数据很多是要求发送方和接受都要进行加密通信,信息的安全与保密显得尤为重 要,因此我想运用异或运算将数据进行隐藏,连续使用同一数据对图像数据两次异或运算图像的数据不发生改变,利用这一特性对图像信息进行加密保护。 熟练使用matlab运用matlab进行编程,使用matlab语言进行数据的隐藏加密,确保数字图像信息的安全,混沌序列具有容易生成,对初始条件和混沌参数敏感等特点,近年来在图像加密领域得到了广泛的应用。使用必要的算法将信息进行加解密,实现信息的保护。 .设计内容和要求 使用混沌序列图像加密技术对图像进行处理使加密后的图像 使用matlab将图像信息隐藏,实现信息加密。 三.设计思路 1. 基于混沌的图像置乱加密算法 本文提出的基于混沌的图像置乱加密算法示意图如图1所示 加密算法如下:首先,数字图像B大小为MX N( M是图像B的行像素数,N是图像B的列像素数),将A的第j行连接到j-1行后面(j=2,3, A,M,形成长度为MX N的序列C。其次,用Logistic混沌映射产生一个长度为的混沌序列{k1,k2,A,kMX N},并构造等差序列D: {1,2,3, A,MX N-1,MX N}。再次,将所

产生的混沌序列{kl, k2. A, kMX N}的M N个值由小到大排序,形成有序序列{k1', k2'. A' kMX N' },确定序列{k1, k2, A, kMX N}中的每个ki在有序序列{k1', k2', A , kMX N' }中的编号,形成置换地址集合 {t1 , t2 , A, tM X N},其中ti为集合{1 , 2, A, MX N}中的一个;按置换地址集合{t1 , t2 , A, tM X N}对序列C进行置换,将其第i个像素置换至第ti列, i=1 , 2, A, MX N,得到C'。将等差序列D做相同置换,得到D'。 最后,B'是一个MX N 的矩阵,B' (i ,j)=C ' ((i-1) X M+j),其中i=1 , 2, A, M j=i=1 , 2, A, N,则B'就是加密后的图像文件。 解密算法与加密算法相似,不同之处在于第3步中,以序列C'代替随机序列{k1, k2, A, kMX N},即可实现图像的解密。 2. 用MATLAB勺实现基于混沌的图像置乱加密算法 本文借助MATLAB^件平台,使用MATLAB!供的文本编辑器进行编程实现加密功能。根据前面加密的思路,把加密算法的编程分为三个主要模块:首先,构造一个与原图a等高等宽的矩阵b加在图像矩阵a后面形成复合矩阵c: b=zeros(m1, n1); ifm1>=n1 ifm1> n1 fore=1: n1 b=(e,e); end else fore=1: n1 end fore=1:( n1-m1) b((m1+e-1),e)=m1+e-1 end end c=zeros(m1*2, n1); c=zeros(m1*2,1); c=[b,a]; 然后,用Logitic映射产生混沌序列:

时间序列分析期末考试

时间序列分析期末考试 TYYGROUP system office room 【TYYUA16H-TYY-

诚信应考,考试作弊将带来严重后果! 湖南大学课程考试试卷 课程名称: 时间序列分析 ;课程编码: 试卷编号: A ;考试时间: 一、 简答题(每小题5分,共计20分) 1、 说明平稳序列建模的主要步骤。 2、 ADF 检验与PP 检验的主要区别是什么? 3、 如何进行两变量的协整检验? 4、 简述指数平滑法的基本思想。 二、 填空题(每小题2分,共计20分) 1. 对平稳序列,在下列表中填上选择的的模型类别 2. 时间序列模型建立后,将要对模型进行显着性检验,那么检验的对象为___________,检验的原假设是___________。 3. 时间序列预处理常进行两种检验,即为_______检验和_______检验。 4. 根据下表,利用AIC 和BIC 准则评判两个模型的相对优劣,你认为______模型优 于______模型。 5. 设ARMA(2, 1):1210.50.1t t t t t X X aX εε---=++-,当a 满足_________时,模型平稳。 6. 设ARMA (2, 1):

则所对应的特征方程为_______________________。 7. 简单季节差分模型的模型结构为: ______________________。 8、对于时间序列{}t X ,如果___________________,则()~2t X I 。 9. 设时间序列{}t X 为来自GARCH(p, q)模型,则其模型结构可写为_____________。 10. k 步差分的定义为k t X ?=___________________________。 三、 (15分)设{}t ε为正态白噪声序列,()()2t t 0,E Var εεσ==,时间序列}{t X 来自 试检验模型的平稳性与可逆性。

基于BP神经网络的时序预测及其应用

目录 摘要 (1) 前言 (2) 第一章时间序列的预测函数及其评价指标 (4) 第一节预测函数 (5) 第二节评价预测的数量指标 (5) 第二章 BP神经网络 (6) 第一节 BP神经网络的结构 (6) 第二节 BP神经网络算法及公式推导 (7) 第三节 BP神经网络算法的步骤 (9) 第三章基于BP神经网络的时间序列预测及其应用 (11) 第四章结论 (14) 总结与体会 (15) 致谢词 (15) 参考文献 (15) 附录 (16)

摘要 首先,本文介绍了时间序列的含义和时间序列预测在国内外的研究情况,列举了两个时间序列预测的实际例子。文中阐述了时间序列预测及其评价指标,比较了各评价指标之间的长处和短处。其次, 本文阐述了BP神经网络算法及其公式推导。给出了BP神经网络算法的流程图。最后,本文从实用出发,列出了1993年至2006年我国GDP的数据,此组数据呈现出增长趋势,这种增长趋势反映了近十几年我国经济的快速增长。用BP神经网络预测出我国2007年的GDP是200790亿元, 这表明今后我国经济有减缓的迹象,这也说明我国近几年宏观经济调控获得了一定的成果。 【关键词】时间序列神经网络预测 GDP Abstract This grade paper, times series, and the development of times series forecast are introduced at first, and then the practical examples of times series forecast are enumerated. The function of times series forecast and its evaluative index are given. We compare the advantage and disadvantage of these evaluative indexes. Secondly, The principles of BP neural network and BP neural network’s algorithm are presented. Finally, we particularize our country GDP statistics, which it increases, which it indicates economy’s fast increasing, year by year, from 1993 to 2006. We also study BP neural network’s forecast algorithm. Our country GDP in 2007,wiche it is about 200790 hundred millions is forecasted by BP neural network, and it shows that the Chinese macro-economy policy in ten years are succeed. Keywords time series neural network prediction GDP

时间序列分析期末论文 (1)

课程论文时间序列分析 题目时间序列模型在人口增长中的应用学院数学与统计学院 专业统计学 班级统计(二)班 学生殷婷 2010101217 指导教师翠霞 职称 2012 年10 月29 日

引言 人口问题是一个世界各国普遍关注的问题。人作为一种资源,主要体现在人既是生产者,又是消费者。作为生产者,人能够发挥主观能动性,加速科技进步,促进社会经济的发展;作为消费者,面对有限的自然资源,人在发展的同时却又不得不考虑人口数量的问题。我国是一个人口大国,人口数量多,增长快,人口素质低;由于人口众多,不仅造成人均资源的数量很少,而且造成住房、教育、就业等方面的很大压力。所以人口数量是社会最为关注的问题,每年新增加的国民生产总值有相当一部分被新增加的人口所抵消,从而造成社会再生产投入不足,严重影响了国民经济的可持续发展。因此,认真分析研究我国目前的人口发展现状和特点,采取切实可行的措施控制人口的高速增长,已经成为我国目前经济发展中需要解决的首要问题。 本文通过时间序列模型对人口的增长进行预测,国家制定未来人口发展目标和生育政策等有关人口政策的基础,对于国民经济计划的制定和社会战略目标的决策具有重要参考价值。人口的预测,作为经济、社会研究的需要,应用越来越广泛,也越来越受到人们的重视。在描绘未来小康社会的蓝图时,首先应要考虑的是未来中国的人口数量、结构、分布、劳动力、负担系数等等,而这又必须通过人口的预测来一一显示。人口数量在时间上的变化,可以用时间序列模型来预测其继后期的数量。 本文通过时间序列分析的方法对人口增长建立模型,取得了较好

的预测结果。时间序列分析是研究动态数据的动态结构和发展变化规律的统计方法。以1990年至2008年中国人口总数为例,用时间序列分析Eviews软件建立模型,并对人口的增长进行预测,研究时间序列模型在人口增长中的应用。 基本假设 (1) 在预测中国人口的增长趋势时,假设全国人口数量的变化是封闭的即人口的出生率和死亡率是自然变化的,而不考虑与其他国家的迁移状况; (2)在预测的年限,不会出现意外事件使人口发生很大的波动,如战争,疾病; (3) 题目数据能够代表全国的整体人数。。 问题分析 根据抽样的基本原理,预测人口增长趋势最直接的方法就是预测出人口总数的增长量,因此我们运用中华人民国国家统计局得到的1990年到2008年度总人口数据。考虑到迁移率、死亡率、出生率、年龄结构等多个因素对人口数量的影响,求解人口增长趋势的关键是如何在我们的模型中充分的利用这些影响因素从而使我们的预测结果具有较高的精确性。 研究数据:

MATLAB动态神经网络在时间序列预测中的应用

MATLAB动态神经网络在时间序列预测中的应用 摘要:本文在介绍了Matlab神经网络工具箱的基础上,主要对时间序列预测工具箱的使用作了说明,并用实例仿真说明如何进行时间序列预测的调用实现,通过不断的调整参数,最后使训练的模型比较理想,满足实际的需求,表明了直接使用时间序列预测的有效性,并为Matlab神经网络工具箱的使用提供了新的方法。 关键词:Matlab;神经网络;时间序列;预测 引言 时间序列是根据时间顺序得到跟时间相关的变量或者参数的观测数据[1]。对时间序列的研究主要是挖掘其中有价值的信息,找到其中变化的内在规律[2]。时间序列预测是时间序列分析研究的主要内容,是指根据现有的和历史的时间序列的数据,建立能反映时间序列中所包含的动态依存关系的数学模型[3],从而能对序列未来的趋势做出合理的预测。简单的说,时间序列预测就是用已有的数据预测下一个时间段的值。目前,时间序列预测已经广泛应用在自然界、经济、化学、科学工程等各个领域。 随着Matlab版本的不断更新,神经网络工具箱不断的完善,使得仿真的实现日益简单,R2010b后的版本对时间序列预测的实现不需要手动写代码,网络训练完毕,从Simple Script可看到网络代码,并可对代码进行编辑、改编,因此,只要调用就可应用在各个领域。本文结合时间序列预测的特点,将Matlab神经网络工具箱中的时间序列预测应用到温度预测的实例中,通过快速的仿真及不断的调整参数,从而形成较理想的数学模型,为后期进行温度的预测奠定了基础。 1Matlab神经网络工具箱简介 神经网络分为静态和动态两类。静态神经网络是无反馈、无记忆的,输出仅依赖于当前的输入,例如BP神经网络和RBF神经网络。动态神经网络是有记忆的神经网络,其输出依赖于当前和以前的输入。动态神经网络又分为有反馈和无反馈,有反馈指输出依赖于当前输入和前一个输入输出,无反馈指输出依赖于当前和之前的输入。因此,动态神经网络比静态神经网络功能强,本文选择动态神经网络进行时间序列预测。 Matlab神经网络工具箱提供了一系列用于模型训练的工具,包括曲线拟合工具箱、模式识别工具箱、聚类工具箱和时间序列工具箱,利用这些工具箱可进行快速的调整参数,通过仿真得到直观的结果。另外,Matlab神经网络工具箱还提供人机交互界面,可根据提示一步一步的完成模型的训练,并对仿真的结果进行分析,直到满足要求为止。 选择时间序列工具箱或者直接在命令窗口中输入ntstool,可打开时间序列预测工具箱界面,根据数据选择符合哪种情况,根据人机交互界面的提示,将数据

混沌序列产生及其在图像

混沌序列产生及其在图像、视频加密中 的应用研究 左建政 【摘要】:伴随着科技的不断进步,信息技术已经渗透进人们生活的各个方面,信息安全问题已经引起越来越多的关注,因此如何加强信息的保密性就成为了一个急需解决的难题。混沌信号所固有的非周期, 宽频谱和对初始值非常敏感等突出特征,使得其在信息加密中有着良好的应用前景。而要想提高混沌在信息加密中的保密性以及实用性,需要做的工作主要是以下两个方面:一方面,提高混沌自身的性能;另一方面,提高加密系统的性能。本文以此为背景,分别从上述两方面着手进行研究,取得了一系列新的结果。本文的主要工作和创新体现在以下几个方面: (1)系统地研究了各种混沌序列的产生方式。混沌序列的产生是混沌应用到信息加密的前提也是一项关键技术。从最初的模拟电路到现在数字系统,在不断地进步,当然也会产生一系列新的问题。本文总结以前所做的研究,系统地介绍了各种混沌序列的产生方式,包括模拟电路、FPGA、LabVIEW、DSP等。通 过分析和比较,为以后的继续深入研究发挥重要的参考作用。 (2)设计了一个新的混沌系统。从混沌信息加密工程的观点考虑,构造一个庞大的混沌函数库是必要的。为了设计性能良好的混沌系统,在研究Sprott系统的基础上,构造了一个新的混沌系统。对构造的混沌系统进行了动力学分析,其中包括分岔特性以及Ly apunov指数等特性分析,同时设计了相应的模拟电路,通过电路实验获得了与 仿真相符的混沌吸引子,验证了混沌系统的特性。 (3)设计了一个新的分数阶混沌系统。并且介绍了两种分数阶微积分的分析方法,时域求解法对其进行数值仿真,时频域转换法对其进行电路仿真。数值仿真结果表明系统存在混沌的最低阶数是2.31。设计了该系统的分数阶混沌振荡电路,电路仿真与数值仿真结果相符,证实了该分数阶混沌振荡电路的可行性。分数阶混沌系统更能反映系统呈现的工程物理现象,一个确定的分数阶混沌系统随着其阶数即分数值的不同而呈现不同的状态,因而这种系统具有更大的密钥空间,不容易被复制。 (4)首次利用数字系统产生分数阶混沌序列。对分数阶混沌系统的广泛研究开始于最近十几年,目前的研究大多处于理论阶段。本文通过利用LabVIEW等数字系统,获得了真实的分数阶混沌序列。通过LabVIEW与MATLAB的接口,首先利用MATLAB编程计算混沌状态方程,然后再在LABVIEW平台设计前面板,调节参数,最后通过数据采集卡即可实现分数阶混沌序列输出。数字系统可以做到参数相同,并且精度可控,容易控制。 (5)借助DSP平台,利用分数阶混沌序列,成功实现了图像、视频的实时加密、解密。利用前面得到的分数阶混沌序列作为图像、视频加密的密钥流,对图像、视频中每一帧图像的像素点进行异或加密。分数阶混沌系统密钥空间大,因此安全性高,需要考虑的主要就是实时性的问题,而借助于运算速度非常快的DSP芯片,能很好地满足实时性的要求。这种加密方式突破了传统软件加密时,加密速度慢、容易被破译的缺点,具有广阔的应用前景

时间序列分析考试卷及答案

考核课程 时间序列分析(B 卷) 考核方式 闭卷 考核时间 120 分钟 注:B 为延迟算子,使得1-=t t Y BY ;?为差分算子,。 一、单项选择题(每小题3 分,共24 分。) 1. 若零均值平稳序列{}t X ,其样本ACF 和样本PACF 都呈现拖尾性,则对{}t X 可能建立( B )模型。 A. MA(2) B.ARMA(1,1) C.AR(2) D.MA(1) 2.下图是某时间序列的样本偏自相关函数图,则恰当的模型是( B )。 A. )1(MA B.)1(AR C.)1,1(ARMA D.)2(MA 3. 考虑MA(2)模型212.09.0--+-=t t t t e e e Y ,则其MA 特征方程的根是( C )。 (A )5.0,4.021==λλ (B )5.0,4.021-=-=λλ (C )5.2221==λλ, (D ) 5.2221=-=λλ, 4. 设有模型112111)1(----=++-t t t t t e e X X X θφφ,其中11<φ,则该模型属于( B )。 A.ARMA(2,1) B.ARIMA(1,1,1) C.ARIMA(0,1,1) D.ARIMA(1,2,1) 5. AR(2)模型t t t t e Y Y Y +-=--215.04.0,其中64.0)(=t e Var ,则=)(t t e Y E ( B )。 A.0 B.64.0 C. 1 6.0 D. 2.0 6.对于一阶滑动平均模型MA(1): 15.0--=t t t e e Y ,则其一阶自相关函数为( C )。 A.5.0- B. 25.0 C. 4.0- D. 8.0 7. 若零均值平稳序列{}t X ?,其样本ACF 呈现二阶截尾性,其样本PACF 呈现拖尾性,则可初步认为对{}t X 应该建立( B )模型。 A. MA(2) B.)2,1(IMA C.)1,2(ARI D.ARIMA(2,1,2) 8. 记?为差分算子,则下列不正确的是( C )。 A. 12-?-?=?t t t Y Y Y B. 212 2--+-=?t t t t Y Y Y Y C. k t t t k Y Y Y --=? D. t t t t Y X Y X ?+?=+?) ( 二、填空题(每题3分,共24分);

基于神经网络的Mackey-Glass时间序列预测

目录 1引言 (1) 2MG时间序列 (1) 2.1MG时间序列简介 (1) 2.2利用dde23函数求解MG时间序列 (1) 3BP神经网络 (3) 3.1神经网络总体思路 (3) 3.2MATLAB中的newff函数 (3) 3.3BP神经网络的训练 (4) 3.4构建输入输出矩阵 (6) 3.5对MG时间序列未来值预测 (6) 4参考文献 (7) 5附录 (8)

1 引言 本文选用的神经网络的是BP 神经网络,利用MATLAB 编程实现。首先通过求解Mackey-Glass 方程得到具有513个数据的Mackey-Glass 时间序列,其中一半用于训练神经网络,一半用于检测预测值。BP 神经网络输入层神经元个数为4,隐含层为8,输出层为1。利用BP 神经网络工具箱构建神经网络并对其进行训练,然后利用训练好的神经网络对未来值进行预测,画出比较图。 2 MG 时间序列 2.1 MG 时间序列简介 Mackey-Glass 混沌系统一类非常典型的混沌系统,混沌系统模型由以下的时滞微分方程来描述: )() (1) ()(t x t x t x dt t dx βτταγ--+-= 其中 α =0.2,β =0.1,γ =10,τ是可调参数,x(t)是在t 时刻的时间序列的值。MG 方程表现出了某种周期性与混沌特性,在τ<16.8时,表现出周期性,在 τ>16.8时,则表现出混沌特性。 2.2 利用dde23函数求解MG 时间序列 本课程设计中取τ=10,也就是说MG 时间序列会表现为周期性。可以利用MATLAB 求解MG 方程,MG 方程是一个时滞微分方程,其中一种求解方法是利用MATLAB 的dde23函数。具体求解方法是:首先建立MG .m 函数文件,代码如下 function y = MG(t,x,z) %UNTITLED Summary of this function goes here % Detailed explanation goes here

时间序列分析期末考试2010B

浙江农林大学 2009 - 2010 学年第 二 学期考试卷(A 卷) 课程名称: 应用时间序列分析 课程类别: 必修 考试方式: 闭卷 注意事项:1、本试卷满分100分。 2、考试时间 120分钟。 一、单项选择题(在每小题的四个备选答案中,选出一个正确 答案,并将正确答案的选项填在题后的括号内。每小题2分,共12分) 1. 关于严平稳与(宽)平稳的关系,不正确的为 。 ( ) A. 严平稳序列一定是宽平稳序列 B. 当序列服从正态分布时,两种平稳性等价 C. 二阶矩存在的严平稳序列一定为宽平稳的 D. MA(p)模型一定是宽平稳的 2. 下图为某时间序列的相关检验图,图1为自相关函数图,图2为偏自相关函数图,请选择模型 。 ( ) 图1 图2 学院: 专业班级: 姓名: 学号: 装 订 线 内 不 要 答 题

A. AR(1) B. AR(2) C. MA(1) D. MA(2) 3. 下图中,图3为某序列一阶差分后的自相关函数图,图4为某序列一阶差分后的 偏自相关函数图,请对原序列选择模型。( ) 图3 图4

A.ARIMA(4,1,0) B. ARIMA(0,2,1) C. ARIMA(0,1,2) D.ARI MA(0,1,4) 4. 记B 为延迟算子,则下列不正确的是 。 ( ) A. 0 1B = B. (1)k t t k t X X B X --=- C. 12t t BX X --= D. 11()t t t t B X Y X Y --±=± 5.对于平稳时间序列,下列错误的是 ( ) A.)(212εσεE = B.),(),(k t t k t t y y Cov y y Cov -+= C.k k -=ρρ D.)(?)1(?1k y k y t t +=+ 6.下图为对某时间序列的拟合模型进行显著性水平0.05α=的显著性检验,请选择 该序列的拟合模型 。 ( )

神经网络预测时间序列

神经网络预测时间序列 如何作预测?理想方法是利用已知数据建立一系列准则,用于一般条件下预测,实际上由于系统的复杂性而不太可能,如股票市场预测。另一种途径是假设一次观测中过去、未来值之间存在联系。其中一种选择是发现一个函数,当过去观测值作为输入时,给出未来值作为输出。这个模型是由神经网络来实现的。 1.2 神经网络预测时间序列 (1) 简单描述 在时间序列预测中,前馈网络是最常使用的网络。在这种情形下,从数学角度看,网络成为输入输出的非线性函数。记一个时间序列为}{n x ,进行其预测可用下式描述: ),,(1+-1-+=m n n n k n x x x f x (1) 时间序列预测方法即是用神经网络来拟合函数)(?f ,然后预测未来值。 (2) 网络参数和网络大小 用于预测的神经网络性质与网络参数和大小均有关。网络结构包括神经元数目、隐含层数目与连接方式等,对一个给定结构来说, 训练过程就是调整参数以获得近似基本联系,误差定义为均方根误差,训练过程可视为一个优化问题。 在大多数的神经网络研究中,决定多少输入与隐层单元数的定量规则问题目前尚未有好的进展,近有的是一些通用指导:首先, 为使网络成为一个完全通用的映射,必须至少有一个隐层。1989年证明一个隐层的网可逼近闭区间内任意一个连续函数。其次,网络结构要尽可能紧致,即满足要求的最小网络最好。实际上,通常从小网络开始。逐步增加隐层数目。同样输入元数目也是类似处理。 (3) 数据和预测精度 通常把可用的时间序列数据分为两部分:训练数据和检验数据。训练数据一般多于检验数据两倍。检验过程有三种方式: 短期预测精度的检验。用检验数据作为输入,输出与下一个时间序列点作比较,误差统计估计了其精度。 长期预测中迭代一步预测。以一个矢量作为输入,输出作为下一个输入矢量的一部分,递归向前传播。 直接多步预测。即用1+-1-m n n n x x x ,,直接进行预测,输出k n x +的预测值,其中 1>k 。

基于混沌序列的密钥生成新方法

第36卷第12期2006年12月数学的实践与认识M A TH EM A T I CS I N PRA CT I CE AND TH EO R Y V o l 136 N o 112  D ecem.,2006  基于混沌序列的密钥生成新方法杨文安, 袁德明 (徐州建筑职业技术学院计算机技术工程系,江苏徐州 221008) 摘要: 设计了一种从混沌序列生成密钥的新方法.其基本原理是从混沌序列依次取若干数据构成实值序列,将其按非线性规则映射成二值序列,再用实值序列和任意指定序列分别置乱这个二值序列,被置乱后的二值序列即为所生成密钥.实验表明,在混沌密码体制研究中,这种密钥较一般序列密钥更具有独立性、均匀性和不可预测性. 关键词: 混沌序列;密码体制;密钥 1 引 言 收稿日期:2006209216 基金项目:徐州建筑职业技术学院自然科学基金(JYA 30409) 随着网络技术发展,网络安全问题成为当今社会的焦点.人们也在不断寻求解决网络安全问题的有效方法[1].混沌密码体制[2,3]就是适应网络安全的需要应运而生的. 混沌是确定性系统中具有理论意义上的完全随机运动,而不是通常所用的伪随机[4]运动.混沌系统具有确定性、有界性、初值敏感性、拓朴传递性与混合性、宽带性、快速衰退的自相关性和长期不可预测性[5]等.它所具有的基本特性恰好满足了信息保密通信和密码学的基本要求(即Shannon 提出的密码系统设计的基本原则:扩散原则和混淆原则,以及加 解密过程中的可靠性).在著名的L ogistic 混沌映射的基础上,作者设计了一种新型密码体制.该密码体制利用非线性变换从双混沌序列中产生两个二值序列,并对其进行置乱得到两组密钥,由这两组密钥构成密钥矩阵作用于明 密文信息实现信息的加 解密操作,进而达到信息保密的目的.有关新型密码体制的详细内容将另文介绍,本文重点介绍利用非线性变换从混沌序列中生成一组密钥的方法.利用文中方法生成的密钥较一般序列密钥更具有独立性、均匀性和不可预测性. 2 L og istic 混沌映射 美国普林斯顿大学的生态学家R .M ay 在研究昆虫群体繁殖规律时提出的L ogistic 混沌模型的离散形式为: x i +1=Λx i (1-x i ), 1<Λ<4, 0

时间序列分析期末考试

时间序列分析期末考试 Prepared on 22 November 2020

诚信应考,考试作弊将带来严重后果! 湖南大学课程考试试卷 课程名称:时间序列分析;课程编码:试卷编号: A ;考试时间:120分 题号一二三四五六七八九十总分 应得分20 20 15 15 20 10 100 实得分 评卷人 一、简答题(每小题5分,共计20分) 1、说明平稳序列建模的主要步骤。 2、ADF检验与PP检验的主要区别是什么 3、如何进行两变量的协整检验 4、简述指数平滑法的基本思想。 二、填空题(每小题2分,共计20分) 1.对平稳序列,在下列表中填上选择的的模型类别 ____年___月___日 考试用

2. 时间序列模型建立后,将要对模型进行显着性检验,那么检验的对象为___________,检验的原假设是___________。 3. 时间序列预处理常进行两种检验,即为_______检验和_______检验。 4. 根据下表,利用AIC 和BIC 准则评判两个模型的相对优劣,你认为______模型优 于______模型。 5. 设ARMA(2, 1):1210.50.1t t t t t X X aX εε---=++-,当a 满足_________时,模型平稳。 6. 设ARMA (2, 1): 1210.50.40.3t t t t t X X X εε---=++- 则所对应的特征方程为_______________________。 7. 简单季节差分模型的模型结构为: ______________________。 8、对于时间序列{}t X ,如果___________________,则()~2t X I 。 9. 设时间序列{}t X 为来自GARCH(p, q)模型,则其模型结构可写为_____________。 10. k 步差分的定义为k t X ?=___________________________。

12-13时间序列分析期末试卷

诚信应考 考出水平 考出风格 浙江大学城市学院 2012— 2013学年第二学期期末考试试卷 《时间序列分析》 开课单位:计算学院 ;考试形式:闭卷;考试时间:2013年7月7日; 所需时间:120分钟 一.简答和计算题(本大题共9题,第1到5题每题5分,第6到9题每题7分,共53分。) 1. 写出(,,)ARIMA p d q 模型的结构。 2. 写出(,)ARMA p q 模型的传递形式和格林函数的递推式。 3. 写出(,)ARMA p q 模型的逆转形式和逆函数的递推式。 第1页共5页

4.计算模型120.5t t t t x x x ε--=--+的偏自相关系数。 5.判断模型121 0.80.5 1.1t t t t t x x x εε---=-++-的平稳性与可逆性。 6. 对于(1)AR 模型: 11()t t t x x μφμε--=-+,根据t 个历史观察值数据: ,10.1,9,6,已求 出?10μ=,1?0.3φ=,29εσ=,求: (1)之后3期的预测值及95%置信区间。 (2)假定获得新的观察值数据为110.5 t x +=,求之后2期的预测值及95%置信区间。 第2页共5页

7.已知某地区每年常住人口数量近似服从(3)MA 模型(单位:万人): 21231000.80.60.2,25t t t t t x εεεεεσ---=+-+-= 最近3年的常住人口数量及一步预测数量如下: 年份 统计人数 预测人数 2002 104 110 2003 108 100 2004 105 109 请预测未来5年该地区常住人口的95%置信区间。 8. 使用指数平滑法得到 ?5t x =, 2? 5.26t x +=,已知序列观察值 5.25 t x =, 1 5.5 t x +=,求指数 平滑系数α。 9. 某一10期观察值序列为5.43, 6.19, 6.63, 7.18, 8.95, 10.14, 11.74, 12.60, 17.26, 21.07 (1)使用6期移动平均法预测12?x 。 (2)使用指数平滑法确定12?x ,其中平滑系数为0.4α= 第3页共5页

用小波神经网络来对时间序列进行预测

/* Note:Your choice is C IDE */ #include"stdio.h" void main() { }/*用小波神经网络来对时间序列进行预测 */ /*%File name : nprogram.m %Description : This file reads the data from %its source into their respective matrices prior to % performing wavelet decomposition. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %Clear command screen and variables */ clc; clear; /*%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % user desired resolution level (Tested: resolution = 2 is best)*/ level = menu('Enter desired resolution level: ', '1',... '2 (Select this for testing)', '3', '4'); switch level case 1, resolution = 1; case 2, resolution = 2; case 3, resolution = 3; case 4, resolution = 4; end msg = ['Resolution level to be used is ', num2str(resolution)]; disp(msg); /*%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % user desired amount of data to use */ data = menu('Choose amount of data to use: ', '1 day', '2 days', '3 days', '4 days',... '5 days', '6 days', '1 week (Select this for testing)'); switch data case 1, dataPoints = 48; /*%1 day = 48 points */ case 2, dataPoints = 96; /* %2 days = 96 points */ case 3, dataPoints = 144; /*%3 days = 144 points */ case 4, dataPoints = 192; /*%4 days = 192 points */ case 5, dataPoints = 240; /* %5 days = 240 points */

相关文档
最新文档