MATLAB Fundamentals>>>Smoothing Electricity Data

本文介绍了如何在MATLAB中使用`smoothdata`函数对电力数据进行1年和2年的移动平均平滑,以及处理不均匀采样点的方法。作者展示了如何计算住宅、商业和工业部门的用电量移动平均,并探讨了`SamplePoints`参数在处理非等间距时间序列数据时的应用。
摘要由CSDN通过智能技术生成

MATLAB Fundamentals>Common Data Analysis Techniques>Smoothing Data> (3/5) Smoothing Electricity Data

数据准备:This code sets up the activity.

load electricityData
whos
total = usage(:,end);
sectorUsage = usage(:,1:3);
sectors = sectors(1:3);
plot(dates,total,".-")

说明1:

The table edata contains monthly electricity usage, in MWh, for the U.S., separated by sector. The matrix sectorUsage contains the consumption for three sectors (residential, commercial, and industrial). The vector total contains the total consumption. The months are stored in the datetime vector dates.

You can smooth the data using the smooth data function.

dataSm = smoothdata(data,"movmean",n)


When the data are evenly sampled, you do not need to specify the sample points.

任务1:

Calculate the 1-year moving average of total. Add this to the existing plot using point markers and a solid line. Repeat for the 2-year moving average.

解答1:

total12 = smoothdata(total,"movmean",12)
hold on 
plot(dates,total12,".-")
total24 = smoothdata(total,"movmean",24)
plot(dates,total24,".-")
hold off

结果1:


说明2:

The smoothdata function works column-wise on matrices. The sectorUsage variable contains usage data from three sectors in its three columns.

任务2:

Calculate the 2-year moving average of all the sectors in sectorUsage. Create a new plot of the result against dates. Use solid lines and no markers.

解答2:

sectorUsage24 = smoothdata(sectorUsage,"movmean",24)
plot(dates,sectorUsage24,"-")

结果2:


附加练习:

The dates are spaced one calendar month apart. However, the months have different lengths, so the usage data is spaced slightly unevenly. To use the dates as sample points, you need to specify the window length as a duration.
 

nDur = years(2)
useSp = smoothdata(sectorUsage,"movmean",...
    nDur,"SamplePoints",dates)
plot(dates,useSp)


The data look almost identical, but a closer look at the difference reveals that they are not.

stem(dates,sectorUsage24-useSp)
legend(sectors)


Here, sectorUsage24 represents the 2-year moving average of all the sectors in sectorUsage, but you may have given this variable from the last task a different name.

附加代码:

nDur = years(2)
useSp = smoothdata(sectorUsage,"movmean",nDur,"SamplePoints",dates)
plot(dates,useSp)

stem(dates,sectorUsage24-useSp)
legend(sectors)

附加结果:


笔记:附加部分的练习到底有何不同,”SamplePoints“到底该怎么用?

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值