def main():
    loaddir = "../../../data/"
    probdir = "prob_trans.inputs.sigma2_u1.100.0.sigma2_u2.1.0.sigma2_v.1.0.m.20.nofu1s.5.N.25000" 
    utils.setsizes()
    
    params, data = dataload(dir=loaddir+"lms/"+probdir, paramfile="stepsizes.txt")
    ymax = data[params==0.001]
    plot(params, data/ymax, label="LMS")

    #params, data = dataload(dir=loaddir+"bmp/"+probdir, paramfile="metastepsizes.txt")
    #plot(params, data/ymax, label="K-Y")

    params, data = dataload(dir=loaddir+"idbd/"+probdir, paramfile="metastepsizes.txt")
    plot(params, data/ymax, label="IDBD")

    params, data = dataload(dir=loaddir+"autostep/"+probdir, paramfile="metastepsizes.txt")
    plot(params, data/ymax, label="Autostep")

    params, data = dataload(dir=loaddir+"rls/"+probdir, paramfile="one_minus_gammas.txt")
    plot(params, data/ymax, label="RLS")

    legend(loc='lower left')
    xlim([10**-4, 1])
    ylim([0, 1.2])
    xscale("log")
    xlabel("tunable parameter", fontsize=utils.getxlabelsize(), labelpad=utils.getxlabelpadsize())
    ylabel("$S$ (MSE relative to best LMS)", fontsize=utils.getxlabelsize(), labelpad=utils.getxlabelpadsize())
def main():
    loaddir = "../../../data/"
    probdir = "prob.s2u.1.0.a.0.9998.s2v.1.0.s2o1.100.0.s2o1.1.0.m.20.o1s.5.N.50000" 
    utils.setsizes()
    
    params, data = dataload(dir=loaddir+"lms/"+probdir, paramfile="stepsizes.txt")
    ymax = 1
    plot(params, data/ymax, label="LMS")

#    params, data = dataload(dir=loaddir+"bmp/"+probdir, paramfile="metastepsizes.txt")
#    plot(params, data/ymax, label="BMP")

    params, data = dataload(dir=loaddir+"idbd/"+probdir, paramfile="metastepsizes.txt")
    plot(params, data/ymax, label="IDBD")
    params, data = dataload(dir=loaddir+"autostep/"+probdir, paramfile="metastepsizes.txt")
    plot(params, data/ymax, label="Autostep")

    params, data = dataload(dir=loaddir+"rls/"+probdir, paramfile="one_minus_gammas.txt")
    plot(params, data/ymax, label="RLS")

    legend(loc='lower left')
    xlim([10**-11, 1])
    ylim([0, 140])
    xscale("log")
    xlabel("tunable parameter", fontsize=utils.getxlabelsize(), labelpad=utils.getxlabelpadsize())
    ylabel("RMSE", fontsize=utils.getxlabelsize(), labelpad=utils.getxlabelpadsize())
Exemplo n.º 3
0
def main():
    utils.setsizes()
    mpl.rcParams['lines.markeredgewidth'] = utils.geterrorbarlinewidth()
    loaddir = "../../../data/"
    probdir1 = "prob.s2u.1.0.a.0.9998.s2v.1.0.s2o1.0.1.s2o1.0.0.m.20.o1s.5.N.50000"
    probdir2 = "prob.s2u.1.0.a.0.9998.s2v.1.0.s2o1.100.0.s2o1.0.0.m.20.o1s.5.N.50000" 
    
    params, data1, data2 = jointdataload(dir1=loaddir+"lms/"+probdir1, dir2=loaddir+"lms/"+probdir2, paramfile="stepsizes.txt")
    ymax1 = data1[params==0.05]
    ymax2 = data2[params==0.05]
    #plot(params, (data1/ymax1 + data2/ymax2)/2.0)

    params, data1, data2 = jointdataload(dir1=loaddir+"lms/"+probdir1, dir2=loaddir+"lms/"+probdir2, paramfile="stepsizes.txt")
    plot(params, (data1/ymax1 + data2/ymax2)/2.0, label="LMS")

    params, data1, data2 = jointdataload(dir1=loaddir+"idbd/"+probdir1, dir2=loaddir+"idbd/"+probdir2, paramfile="metastepsizes.txt")
    plot(params, (data1/ymax1 + data2/ymax2)/2.0, label="IDBD")

    params, data1, data2 = jointdataload(dir1=loaddir+"autostep/"+probdir1, dir2=loaddir+"autostep/"+probdir2, paramfile="metastepsizes.txt")
    plot(params, (data1/ymax1 + data2/ymax2)/2.0, label="Autostep")

    #params, data1, data2 = jointdataload(dir1=loaddir+"bmp/"+probdir1, dir2=loaddir+"bmp/"+probdir2, paramfile="metastepsizes.txt")
    #plot(params, (data1/ymax1 + data2/ymax2)/2.0, label="K-Y")

    params, data1, data2 = jointdataload(dir1=loaddir+"rls/"+probdir1, dir2=loaddir+"rls/"+probdir2, paramfile="one_minus_gammas.txt")
    plot(params, (data1/ymax1 + data2/ymax2)/2.0, label="RLS")

    legend(loc='lower left')
    xlim([10**-9, 1])
    ylim([0.6, 1.2])
    xscale("log")
    xlabel("tunable parameter", fontsize=utils.getxlabelsize())
    ylabel("$S$ (RMSE relative to best LMS)", fontsize=utils.getxlabelsize(), labelpad=utils.getxlabelpadsize())