Since you are new, let us show you how Trendy works.

Close-button

Solar Activity

This plot is broken.

Undefined function 'smooth' for input arguments of type 'double'.

Did you mean: correlation = smooth3(R,span); ?
  • Created by: Michael Umbricht
  • Latest result: Undefined function 'smooth' for input arguments of type 'double'. Did you mean: correlation = smooth3(R,span); ?
  • Created on: 08 Sep 2013

This plot shows two indices of solar activity:

  • the Boulder Index of the observed sunspot number
  • the solar radio flux measured at a wavelength of 10.7 cm

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.

Plot Image
% 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.0;
correlation = R;
% 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')
Tags:

Add Tags


0 comments