U2

时序图

> yield=c(15.2,16.9,15.3,14.9,15.7,15.1,16.7)
> yield=ts(yield,start=1884)
> plot(yield)

点线结构参数

散点图

plot(yield,type="p")

点线图

plot(yield,type="o")

符号参数

> plot(yield,type="o",pch=17)

符号参数

连线类型参数

> plot(yield,lty=2)

线的宽度参数

> plot(yield,lwd=2)

颜色参数

> plot(yield,col=3)

添加文本

> plot(yield,main="1884-1890年英格兰和威尔士地区小麦平均亩产量",xlab="年份",ylab="亩产量")

指定坐标轴范围

#指定输出横轴范围
> plot(yield,xlim=c(1886,1890))

#指定输出纵轴范围
> plot(yield,ylim=c(15,16))

添加参考线

#添加一条垂线
> plot(yield)
> abline(v=1887,lty=2)

#添加多条垂直参考线
> plot(yield)
> abline(v=c(1885,1889),lty=2)

#添加水平线
> plot(yield)
> abline(h=c(15.5,16.5),lty=2)

绘制自相关图

> acf(yield)
    #acf(x,lag=)    x：变量名；lag：延迟阶数，若用户不特殊指定的话，系统会根据序列长度自动指定延迟阶数

虚线为自相关系数2倍标准差位置

例2.1 时序图

> sha=read.table("E:/例题data/file4.csv",sep=",",header=T)
> output=ts(sha$output,start=1964)
> plot(output)

自相关图

> sha=read.table("E:/例题data/file4.csv",sep=",",header=T)
> output=ts(sha$output,start=1964)
> plot(output)
> acf(output,lag=25)

习题2.5

https://wenku.baidu.com/view/984cefabfab069dc50220176.html

第一题第二问matlab程序
>> clear all;close all;
>> X_t=[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20];%时间序列值
>> k=6; %拟计算的自相关延迟期数值
>> rou_hat=zeros(k,1); %为拟计算的中阶自相关系数预留数组空间
>> Time_mean=mean(X_t); %计算时间序列值的样本均值
>> 
>> %计算自相关系数计算公式的分母部分
>> SST=0;
>> for t=1:20
SST=SST+(X_t(t)-Time_mean)^2;
end
>> 
>> %计算对应k阶自相关系数
>> for i=1:k
S=0;
for t=1:20-i
S=S+(X_t(t)-Time_mean)*(X_t(t+i)-Time_mean);%计算公式分子部分
end
rou_hat(i,1)=S/SST;
end
>> rou_hat=rou_hat %输出相关系数值

rou_hat =

    0.8500
    0.7015
    0.5560
    0.4150
    0.2801
    0.1526