This plot shows two indices of solar activity:
The data is taken from the SpaceWeather.com website once per day.
It also shows the correlation that has been observed between radio flux and sunspot number using the equations from the Bureau of Meteorology in Australia. These equations are only valid on a statistical (ie, average) basis, so the correlation is plotted as a smoothed 30 day moving average. The smoothed Sunspot Number (also a 30 day moving average) is plotted for comparison.
% Sunspot number % time vector is: time2088 % data vector is: data2088 plot(time2088,data2088,'.','MarkerSize',7,'Color','magenta'); datetick hold all % Solar radio flux % time vector is: time2089 % data vector is: data2089 plot(time2089,data2089,'.','MarkerSize',7,'Color','cyan'); datetick hold all % Calculate correlation between flux and sunspot number F = (data2089 - 67.0); R = (1.61*F) - (0.0733*F).^2 + (0.0240*F).^3; span=30; correlation = smooth(R,span); smoothsunspot = smooth(data2088,span); plot(time2089,correlation,'--','Color','blue'); hold all plot(time2088,smoothsunspot,'--','Color','red'); datetick hleg1 = legend('Sunspot Number','Solar Flux at 10.7 cm','Smoothed Correlation','Smoothed Sunspot Number'); set(hleg1,'Location','Northwest')