{"id":7551,"date":"2021-01-21T13:37:08","date_gmt":"2021-01-21T18:37:08","guid":{"rendered":"https:\/\/blogs.mathworks.com\/community\/?p=7551"},"modified":"2021-02-08T16:07:39","modified_gmt":"2021-02-08T21:07:39","slug":"matlab-drive-and-british-coal","status":"publish","type":"post","link":"https:\/\/blogs.mathworks.com\/community\/2021\/01\/21\/matlab-drive-and-british-coal\/","title":{"rendered":"MATLAB Drive and British Coal"},"content":{"rendered":"
\r\n
The other day I received this email from my friend Jos.<\/div>\r\n
\"\"<\/div>\r\n
Very thoughtful of him! After a recent meeting, we were discussing the UK's notable decline in coal consumption. Jos had found some data from a government website<\/a> and shared it with me via MATLAB Drive. Please note that all the numbers shown here have to do with the energy used to generate electricity (as opposed to, say, heat or transportation).<\/div>\r\n
Because we both use MATLAB Drive<\/a>, Jos can leave his code and data in place and just issue the invitation. Once I accept it, I can work in the same exact folder, thereby avoiding a lot of back and forth messaging with attachments.<\/div>\r\n
I accepted the invitation and saw a folder with a spreadsheet, two import routines, and a file called runme.m.<\/div>\r\n
\"\"<\/div>\r\n
Here's that runme.m file from Jos.<\/div>\r\n
\"\"<\/div>\r\n
And when I run it, this is the result.<\/div>\r\n
\"\"<\/div>\r\n
Sure enough, the usage of coal is now close to zero! But since Jos invited me to collaborate with this shared data, I thought I would fill out the analysis a bit more. Let's dig a little deeper.<\/div>\r\n
Step one was to move his code into a Live Script (much later, I would take this Live Script and publish it directly to WordPress, so you're effectively looking at the Live Script I created).<\/div>\r\n
Next, I loaded the data and cleaned up the table a bit more.<\/div>\r\n
\r\n
\r\n
filename = \"UK-energy-fuels.xls\"<\/span>;<\/span><\/div>\r\n<\/div>\r\n
\r\n
data = importdata(filename);<\/span><\/div>\r\n<\/div>\r\n
\r\n
headers = importheaders(filename);<\/span><\/div>\r\n<\/div>\r\n
\r\n
<\/div>\r\n<\/div>\r\n
\r\n
% Clean up some headers<\/span><\/span><\/div>\r\n<\/div>\r\n
\r\n
headers(5) = \"Hydro\"<\/span>;<\/span><\/div>\r\n<\/div>\r\n
\r\n
headers(7) = \"Wind_Onshore\"<\/span>;<\/span><\/div>\r\n<\/div>\r\n
\r\n
headers(8) = \"Wind_Offshore\"<\/span>;<\/span><\/div>\r\n<\/div>\r\n
\r\n
headers(end) = \"Total\"<\/span>;<\/span><\/div>\r\n<\/div>\r\n
\r\n
<\/div>\r\n<\/div>\r\n
\r\n
% Missing values in data (NaN should be zero)<\/span><\/span><\/div>\r\n<\/div>\r\n
\r\n
data(isnan(data)) = 0;<\/span><\/div>\r\n<\/div>\r\n
\r\n
<\/div>\r\n<\/div>\r\n
\r\n
t = array2table(data','VariableNames'<\/span>,headers);<\/span><\/div>\r\n<\/div>\r\n
\r\n
% Throw out some unused columns<\/span><\/span><\/div>\r\n<\/div>\r\n
\r\n
t(:,[9 12 13 14]) = [];<\/span><\/div>\r\n<\/div>\r\n<\/div>\r\n
Let's switch to a timetable. Timetables are awesome.<\/div>\r\n
\r\n
\r\n
t.Time = datetime(1998:2019,1,1)';<\/span><\/div>\r\n<\/div>\r\n
\r\n
t = t(:,[end 1:(end-1)]);<\/span><\/div>\r\n<\/div>\r\n
\r\n
t = table2timetable(t)<\/span><\/div>\r\n
\r\n
\r\n
\r\n
t = 22\u00d711 timetable <\/span><\/div>\r\n
\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n
<\/th>\r\nTime<\/th>\r\nCoal<\/th>\r\nOil<\/th>\r\nGas<\/th>\r\nNuclear<\/th>\r\nHydro<\/th>\r\nWind<\/th>\r\nWind_Onshore<\/th>\r\nWind_Offshore<\/th>\r\nSolar<\/th>\r\nBioenergy<\/th>\r\nTotal<\/th>\r\n<\/tr>\r\n<\/thead>\r\n
1<\/th>\r\n01-Jan-1998<\/td>\r\n117.0348<\/td>\r\n6.8340<\/td>\r\n116.2895<\/td>\r\n90.5900<\/td>\r\n5.0939<\/td>\r\n0.8769<\/td>\r\n0<\/td>\r\n0<\/td>\r\n0<\/td>\r\n2.9637<\/td>\r\n355.1683<\/td>\r\n<\/tr>\r\n
2<\/th>\r\n01-Jan-1999<\/td>\r\n101.2574<\/td>\r\n6.0717<\/td>\r\n139.6720<\/td>\r\n87.6721<\/td>\r\n5.3034<\/td>\r\n0.8510<\/td>\r\n0<\/td>\r\n0<\/td>\r\n0<\/td>\r\n3.7010<\/td>\r\n361.9159<\/td>\r\n<\/tr>\r\n
3<\/th>\r\n01-Jan-2000<\/td>\r\n114.7365<\/td>\r\n5.9284<\/td>\r\n144.8912<\/td>\r\n78.3339<\/td>\r\n5.0586<\/td>\r\n0.9470<\/td>\r\n0<\/td>\r\n0<\/td>\r\n0<\/td>\r\n4.0673<\/td>\r\n371.4396<\/td>\r\n<\/tr>\r\n
4<\/th>\r\n01-Jan-2001<\/td>\r\n125.4019<\/td>\r\n4.7766<\/td>\r\n138.7174<\/td>\r\n82.9850<\/td>\r\n4.0324<\/td>\r\n0.9650<\/td>\r\n0<\/td>\r\n0<\/td>\r\n0<\/td>\r\n4.7554<\/td>\r\n374.5751<\/td>\r\n<\/tr>\r\n
5<\/th>\r\n01-Jan-2002<\/td>\r\n118.4754<\/td>\r\n4.2168<\/td>\r\n148.8707<\/td>\r\n81.0903<\/td>\r\n4.7628<\/td>\r\n1.2590<\/td>\r\n0<\/td>\r\n0<\/td>\r\n0<\/td>\r\n5.3066<\/td>\r\n375.0735<\/td>\r\n<\/tr>\r\n
6<\/th>\r\n01-Jan-2003<\/td>\r\n131.7604<\/td>\r\n4.1707<\/td>\r\n145.1360<\/td>\r\n81.9114<\/td>\r\n3.2117<\/td>\r\n1.2880<\/td>\r\n0<\/td>\r\n0<\/td>\r\n0<\/td>\r\n6.2900<\/td>\r\n378.6893<\/td>\r\n<\/tr>\r\n
7<\/th>\r\n01-Jan-2004<\/td>\r\n125.6896<\/td>\r\n4.0933<\/td>\r\n153.7336<\/td>\r\n73.6816<\/td>\r\n4.8205<\/td>\r\n1.9390<\/td>\r\n0<\/td>\r\n0<\/td>\r\n0<\/td>\r\n7.4275<\/td>\r\n380.8878<\/td>\r\n<\/tr>\r\n
8<\/th>\r\n01-Jan-2005<\/td>\r\n128.5131<\/td>\r\n4.6509<\/td>\r\n149.2124<\/td>\r\n75.1728<\/td>\r\n4.7505<\/td>\r\n2.9120<\/td>\r\n0<\/td>\r\n0<\/td>\r\n0<\/td>\r\n8.9544<\/td>\r\n385.0995<\/td>\r\n<\/tr>\r\n
9<\/th>\r\n01-Jan-2006<\/td>\r\n141.4950<\/td>\r\n5.4071<\/td>\r\n137.7543<\/td>\r\n69.2372<\/td>\r\n4.5655<\/td>\r\n4.2360<\/td>\r\n0<\/td>\r\n0<\/td>\r\n0<\/td>\r\n9.1093<\/td>\r\n381.3778<\/td>\r\n<\/tr>\r\n
10<\/th>\r\n01-Jan-2007<\/td>\r\n129.0301<\/td>\r\n4.4649<\/td>\r\n162.3896<\/td>\r\n57.2489<\/td>\r\n5.0317<\/td>\r\n5.2879<\/td>\r\n0<\/td>\r\n0<\/td>\r\n0<\/td>\r\n8.5343<\/td>\r\n379.2790<\/td>\r\n<\/tr>\r\n
11<\/th>\r\n01-Jan-2008<\/td>\r\n118.0531<\/td>\r\n5.8852<\/td>\r\n172.9876<\/td>\r\n47.6731<\/td>\r\n5.1109<\/td>\r\n7.1385<\/td>\r\n0<\/td>\r\n0<\/td>\r\n0<\/td>\r\n8.6219<\/td>\r\n378.2265<\/td>\r\n<\/tr>\r\n
12<\/th>\r\n01-Jan-2009<\/td>\r\n97.7979<\/td>\r\n5.3651<\/td>\r\n163.4552<\/td>\r\n62.7617<\/td>\r\n5.1968<\/td>\r\n9.3016<\/td>\r\n0<\/td>\r\n0<\/td>\r\n0<\/td>\r\n9.6011<\/td>\r\n358.2008<\/td>\r\n<\/tr>\r\n
13<\/th>\r\n01-Jan-2010<\/td>\r\n102.1764<\/td>\r\n4.3082<\/td>\r\n172.4532<\/td>\r\n56.4417<\/td>\r\n3.5658<\/td>\r\n10.3278<\/td>\r\n0<\/td>\r\n0<\/td>\r\n0<\/td>\r\n11.1318<\/td>\r\n364.3963<\/td>\r\n<\/tr>\r\n
14<\/th>\r\n01-Jan-2011<\/td>\r\n103.0063<\/td>\r\n2.8074<\/td>\r\n143.8056<\/td>\r\n62.6552<\/td>\r\n5.6546<\/td>\r\n16.2075<\/td>\r\n0<\/td>\r\n0<\/td>\r\n0<\/td>\r\n11.8370<\/td>\r\n353.9106<\/td>\r\n<\/tr>\r\n
\u00ff<\/th>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n<\/div>\r\n<\/div>\r\n<\/div>\r\n<\/div>\r\n<\/div>\r\n<\/div>\r\n
A quick stacked plot is a fun way to get our bearings. The nice thing about stacked plots is that they know all about timetables, so there's not much you have to worry about. Just one line of code, and no messing around with labels and limits.<\/div>\r\n
\r\n
\r\n
stackedplot(t(:,[1 2 3 4 5 6 9 10]))<\/span><\/div>\r\n
\r\n
\r\n
<\/div>\r\n<\/div>\r\n<\/div>\r\n<\/div>\r\n<\/div>\r\n
Since the y-scales are all different, the stacked plot is mostly useful for trends. Notice that Wind and Solar are moving steadily upward, while Hydro hasn't varied much. And there's Coal doing its disappearing act.<\/div>\r\n
To understand the relative contributions of these energy sources, let's do an area plot.<\/div>\r\n
Tables make this so pleasant!<\/div>\r\n
\r\n
\r\n
h = area(t.Time,[t.Coal t.Oil t.Gas t.Nuclear t.Hydro t.Wind t.Solar t.Bioenergy]);<\/span><\/div>\r\n<\/div>\r\n<\/div>\r\n
Fun fact: you can now set colors in MATLAB plots with CSS-style hex color designations. Let's find some appropriate colors for the chart. I picked these colors with this tool: https:\/\/www.w3schools.com\/colors\/colors_picker.asp.<\/a><\/div>\r\n
\r\n
\r\n
sources = [ ...<\/span><\/span><\/div>\r\n<\/div>\r\n
\r\n
\"#4d4d4d\"<\/span>, \"Coal\"<\/span>;<\/span><\/div>\r\n<\/div>\r\n
\r\n
\"#663300\"<\/span>, \"Oil\"<\/span>;<\/span><\/div>\r\n<\/div>\r\n
\r\n
\"#999966\"<\/span>, \"Gas\"<\/span>;<\/span><\/div>\r\n<\/div>\r\n
\r\n
\"#aa80ff\"<\/span>, \"Nuclear\"<\/span>;<\/span><\/div>\r\n<\/div>\r\n
\r\n
\"#0066ff\"<\/span>, \"Hydro\"<\/span>;<\/span><\/div>\r\n<\/div>\r\n
\r\n
\"#66ccff\"<\/span>, \"Wind\"<\/span>;<\/span><\/div>\r\n<\/div>\r\n
\r\n
\"#ffff66\"<\/span>, \"Solar\"<\/span>;<\/span><\/div>\r\n<\/div>\r\n
\r\n
\"#33cc33\"<\/span>, \"Bioenergy\"<\/span>]; <\/span><\/div>\r\n<\/div>\r\n
\r\n
<\/div>\r\n<\/div>\r\n
\r\n
for <\/span>i = 1:length(h)<\/span><\/div>\r\n<\/div>\r\n
\r\n
h(i).FaceColor = sources(i,1);<\/span><\/div>\r\n<\/div>\r\n
\r\n
end<\/span><\/span><\/div>\r\n<\/div>\r\n
\r\n
<\/div>\r\n<\/div>\r\n
\r\n
legend(sources(:,2), ...<\/span><\/span><\/div>\r\n<\/div>\r\n
\r\n
\"Location\"<\/span>,\"northeastoutside\"<\/span>)<\/span><\/div>\r\n<\/div>\r\n
\r\n
xlabel(\"Date\"<\/span>)<\/span><\/div>\r\n<\/div>\r\n
\r\n
ylabel(\"Energy (TWh)\"<\/span>)<\/span><\/div>\r\n<\/div>\r\n
\r\n
grid on<\/span><\/span><\/div>\r\n
\r\n
\r\n
<\/div>\r\n<\/div>\r\n<\/div>\r\n<\/div>\r\n<\/div>\r\n
Here's Jos's Coal plot, but normalized against the total, so we're looking at energy share as a percentage.<\/div>\r\n
\r\n
\r\n
plot(t.Time, 100*t.Coal.\/t.Total, \"o-\"<\/span>,\"LineWidth\"<\/span>,3)<\/span><\/div>\r\n<\/div>\r\n
\r\n
grid on<\/span><\/span><\/div>\r\n<\/div>\r\n
\r\n
xlabel(\"Date\"<\/span>)<\/span><\/div>\r\n<\/div>\r\n
\r\n
ylabel(\"Percentage\"<\/span>)<\/span><\/div>\r\n<\/div>\r\n
\r\n
title(\"Percentage of energy derived from Coal in the UK\"<\/span>)<\/span><\/div>\r\n
\r\n
\r\n
<\/div>\r\n<\/div>\r\n<\/div>\r\n<\/div>\r\n<\/div>\r\n
Who, at the beginning of this century, would have predicted this outcome? There is a certain poetic symmetry to the fact that Britain was the first country to make massive use of coal for industrial power, and it may well be the first to let it go.<\/div>\r\n
There's an interesting side story here about the horse race between onshore and offshore wind. The UK has been a pioneer in offshore wind. You might say they're the Saudi Arabia of blustery sealanes.<\/div>\r\n
\r\n
\r\n
plot(t.Time,[t.Wind_Onshore t.Wind_Offshore],\"o-\"<\/span>,\"LineWidth\"<\/span>,3)<\/span><\/div>\r\n<\/div>\r\n
\r\n
xlim(datetime([2012 2020],1,1))<\/span><\/div>\r\n<\/div>\r\n
\r\n
legend({'Onshore'<\/span>,'Offshore'<\/span>},'Location'<\/span>,\"best\"<\/span>)<\/span><\/div>\r\n<\/div>\r\n
\r\n
xlabel(\"Date\"<\/span>)<\/span><\/div>\r\n<\/div>\r\n
\r\n
ylabel(\"Energy (TWh)\"<\/span>)<\/span><\/div>\r\n<\/div>\r\n
\r\n
grid on<\/span><\/span><\/div>\r\n
\r\n
\r\n
<\/div>\r\n<\/div>\r\n<\/div>\r\n<\/div>\r\n<\/div>\r\n
From a standing start at the beginning of the century, wind now generates 20% of the UK's power. And between the two flavors, onshore and offshore, as of 2019 it's a dead heat! Given the trends, you'd have to put your money on offshore wind. There's a lot more ocean than island. And fewer crinkly bits to get in the way.<\/div>\r\n
As I was doing these plots, I couldn't help but notice that the total amount of energy consumed in Britain is shrinking over time. Which is all to the good. After all, it's much easier to use less coal when you need less overall energy in the first place.<\/div>\r\n
\r\n
\r\n
plot(t.Time, t.Total,\"LineWidth\"<\/span>,3)<\/span><\/div>\r\n<\/div>\r\n
\r\n
grid on<\/span><\/span><\/div>\r\n<\/div>\r\n
\r\n
xlabel(\"Date\"<\/span>)<\/span><\/div>\r\n<\/div>\r\n
\r\n
ylabel(\"Energy (TWh)\"<\/span>)<\/span><\/div>\r\n<\/div>\r\n
\r\n
title(\"Overall power consumption in the UK\"<\/span>)<\/span><\/div>\r\n
\r\n
\r\n
<\/div>\r\n<\/div>\r\n<\/div>\r\n<\/div>\r\n<\/div>\r\n
Since I am a speculative person, and since MATLAB is the master of curve fitting, I naturally wondered: suppose this downward trend keeps up? What then?<\/div>\r\n
\r\n
\r\n
tm = year(t.Time(8:end));<\/span><\/div>\r\n<\/div>\r\n
\r\n
dt = t.Total(8:end);<\/span><\/div>\r\n<\/div>\r\n
\r\n
<\/div>\r\n<\/div>\r\n
\r\n
% p is the polynomial equation for the linear fit to the data<\/span><\/span><\/div>\r\n<\/div>\r\n
\r\n
p = polyfit(tm,dt,1);<\/span><\/div>\r\n<\/div>\r\n
\r\n
<\/div>\r\n<\/div>\r\n
\r\n
% ti is where we intersect the x axis<\/span><\/span><\/div>\r\n<\/div>\r\n
\r\n
ti = -p(2)\/p(1);<\/span><\/div>\r\n<\/div>\r\n
\r\n
<\/div>\r\n<\/div>\r\n
\r\n
t0 = 2005;<\/span><\/div>\r\n<\/div>\r\n
\r\n
tf = [t0 ti];<\/span><\/div>\r\n<\/div>\r\n
\r\n
df = polyval(p,tf);<\/span><\/div>\r\n<\/div>\r\n
\r\n
hold on<\/span><\/span><\/div>\r\n<\/div>\r\n
\r\n
plot(datetime(round(tf),1,1),df,\"o-\"<\/span>,\"Color\"<\/span>,\"red\"<\/span>)<\/span><\/div>\r\n<\/div>\r\n
\r\n
hold off<\/span><\/span><\/div>\r\n<\/div>\r\n
\r\n
xlim(datetime([2000 2100],1,1))<\/span><\/div>\r\n<\/div>\r\n
\r\n
xlabel(\"Date\"<\/span>)<\/span><\/div>\r\n<\/div>\r\n
\r\n
ylabel(\"Energy (TWh)\"<\/span>)<\/span><\/div>\r\n
\r\n
\r\n
<\/div>\r\n<\/div>\r\n<\/div>\r\n<\/div>\r\n<\/div>\r\n
Based on current trends, we can confidently assert that the UK will need no power of any kind by 2096.<\/div>\r\n
So thanks for this data, Jos! Sharing data and code on MATLAB Drive made this investigation easy. And now that I've done my bit, I'll see what Jos has to add. I'm curious to discover, in 80 years time, how he'll manage with no electricity.<\/div>\r\n<\/div>\r\n","protected":false},"excerpt":{"rendered":"
<\/div>

\r\nThe other day I received this email from my friend Jos.\r\n\r\nVery thoughtful of him! After a recent meeting, we were discussing the UK's notable decline in coal consumption. Jos had found some data... read more >><\/a><\/p>","protected":false},"author":69,"featured_media":0,"comment_status":"open","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":[],"categories":[1],"tags":[],"_links":{"self":[{"href":"https:\/\/blogs.mathworks.com\/community\/wp-json\/wp\/v2\/posts\/7551"}],"collection":[{"href":"https:\/\/blogs.mathworks.com\/community\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/blogs.mathworks.com\/community\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/blogs.mathworks.com\/community\/wp-json\/wp\/v2\/users\/69"}],"replies":[{"embeddable":true,"href":"https:\/\/blogs.mathworks.com\/community\/wp-json\/wp\/v2\/comments?post=7551"}],"version-history":[{"count":4,"href":"https:\/\/blogs.mathworks.com\/community\/wp-json\/wp\/v2\/posts\/7551\/revisions"}],"predecessor-version":[{"id":7561,"href":"https:\/\/blogs.mathworks.com\/community\/wp-json\/wp\/v2\/posts\/7551\/revisions\/7561"}],"wp:attachment":[{"href":"https:\/\/blogs.mathworks.com\/community\/wp-json\/wp\/v2\/media?parent=7551"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/blogs.mathworks.com\/community\/wp-json\/wp\/v2\/categories?post=7551"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/blogs.mathworks.com\/community\/wp-json\/wp\/v2\/tags?post=7551"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}