ExperimentA77.m

load('BaseParams.mat')

expnum = 'A77';
time = 71222:7:73413; % days after January 1, 1800

max_iter = 12;
numModes = 5000;

dictfn = @(numModes) NOAAComboDictEven(numModes,time);

dataname = 'NOAAdata';

numTrials = 1;
errors = zeros(numTrials, 1);
corrs = zeros(numTrials, 1);
sparsities = zeros(numTrials, 1);

for t = 1:numTrials
    [errors(t), corrs(t), sparsities(t)] = BaseExperiment(decompfn, dictfn, dataname, initfnA, ...
    initfnB, initfnZ, noisefn, numModes, Zthresholdfn, max_iter, lowtau, ...
    hightau, perturbAfn, perturbBfn, perturbZfn, expnum, t, ...
    lenSqn, slopeThreshold, diffThreshold);
end

csvwrite(sprintf('./exp%s/summary.csv',expnum),[errors, corrs, sparsities])