EMD（经验模态分解）算法 一

参考R语言提取时间序列的周期性成分应用EMD，小波滤波器，BAXTER过滤器等的方法进行应用。

第一次尝试

首先安装依赖包

# install.packages("reshape2")
# install.packages("plm")
# install.packages("ggplot2")

加载包

library(reshape2)
library(plm)
library(ggplot2)

读取数据

load<-data.frame(read.csv("安徽最高.csv"))

观察数据

load

重命名

names(load) <- c("Time","Load")

取对数

load[,"Load"] <- log(load[,"Load"]) 

观察对数负载随时间的变化

ggplot(load,aes(x=Time,y=Load)) + geom_line(size=.5) + theme_classic() + labs(title="对数负载")

形状如上，波动还是挺大，冬天夏天用电处于高峰，没有特别明显的线性或非线性趋势。
尝试去除线性趋势看看。

dat <- data.frame("Time"=load[,