matlab 预测方法的题目,54504《经济预测与决策技术及MATLAB实现（第2版）》杨德平、刘喜华（主要例题程序）教学资料-.doc...

经济预测与决策技术及MATLAB实现

(2016年第2版)

主要例题程序

【例1-11】

(1)在Matlab编辑器中建立一个任意命名的M文件(譬如a1.m)，其包含的程序为：

%输入矩阵A

A=[4.55 5.68 4.02

4.72 5.96 4.18

4.86 6.16 4.34

4.69 5.91 4.17

4.66 5.85 4.12

4.62 5.75 4.06

4.66 5.81 4.11

4.61 5.78 4.07

4.59 5.71 4.04

4.65 5.82 4.1

4.6 5.74 4.04

4.59 5.74 4.04

4.62 5.78 4.07

4.62 5.82 4.08

4.69 5.92 4.12

4.67 5.88 4.09

4.67 5.9 4.08

4.65 5.9 4.06

4.65 5.9 4.06

4.59 5.78 4.02

4.61 5.8 4.04

4.61 5.79 4.03

4.58 5.78 4.01];

save yinhang A %保存数据的命令

(2)运行a1.m文件，则将数据A保存在一个名为yinhang.dat的数据文件中。

(3)进行统计分析时，先用命令：load yinhang 调用数据文件yinhang中的数据。

例如，在命令行窗口输入命令可查看数据的属性：

>>load yinhang

>> [m,n]=size(A) %显示矩阵A的行和列

再如，创建一个坐标轴后，可以比较三支银行股票的收盘价曲线：

>>t=1:m;

>>x1=A(:,1); %取A的第一列数据，选取工商银行

>>x2=A(:,2); %取A的第二列数据，选取建设银行

>>x3=A(:,3); %取A的第三列数据，选取中国银行

>>plot(t,x1,'-ob',t,x2,'-+k',t,x3,'-^r') %绘制三条带有标识符的曲线

>>legend('工商银行','建设银行','中国银行')

>>xlabel('时间/天')

>>ylabel('收盘价/元')

【例1-12】 说明：本题演示数据与表1-3数据有几个数据做了改动！后面章节用到本题数据都以下面数据为准！

X=[16.84 16.77 16.16 16.65 16.75 17.29 17.12 18.35 18.01 17.88 ...

17.56 17.71 17.57 17.96 17.73 17.52 17.72 17.51 17.44 17.77 ...

17.73 17.96 18.19 18.06 17.85 18.34 18.58 18.25 18.22 18.69 ...

18.29 18.19 18.07 18.04 17.86 17.81 17.93 17.83 17.9 18.36 ...

18.57 18.68 19.08 18.64 18.44 18.27 18.14 18.1 18.14 18.16 ...

18.2 18.15 18.19 17.88 17.73 17.53 17.04 17.2 17.11 17.91 ...

17.78 17.25 17.39 17.42 17.39 17.26 16.92 17.06 17.03 16.96 ...

17.1 17.52 18.23 18.27 18.19 18.14 18.33 17.84 17.94 18.26];

%先画出直方图观察其服从的分布

hist(X,7)

%用正态分布拟合

[muhat,sigmahat,muci,sigmaci] = normfit(X,0.05)

%用最大似然估计法估计参数

[phat, pci]=mle('norm',X)

【例2-1】

x=[57346.5 57997.61 58844.99 57691.6458438.53 60259.53 ...

63040.51 72896.19 61949.81 60772.46 59075.9758604.26];

plot(x,'o-')

set(gca,'XTick',[1 2 3 4 5 6 7 8 9 10 11 1