generate_plots.m

implement basic routines in Matlab and Octave (currently using Matlab R13) for simulating adaptive b

cv_mov = [ 0 0 0 0 0  1 0 0 0 3  1 4 2 5 9  20  13 2 6 2 3  3 0 2 0 0  0 0 0 0 0]
blind_jdiv = [ 0 0 0 0 0  0 1 0 1 1  6 3 2 4 6  28  9 5 2 3 1  1 1 2 0 0  0 0 0 0 0]
sum(cv_mov)
sum(blind_jdiv)
set( 0, 'defaultaxesfontname', 'arial', 'defaultaxesfontsize', 10 )

figure
hold on
plot(-15:1:15, cv_mov, 'o-k', 'MarkerSize', 6, 'MarkerFaceColor','k')
plot(-15:1:15, blind_jdiv, 'o-k', 'MarkerSize', 6)
xlabel('Numerical Iteratrion (number of auxiliary vectors) n');
ylabel('Frequency of occurance (%)');
legend('CV-MOV', 'J-DIV (blind)')

set(gcf,'units','inches');
set(gcf,'PaperPosition',[5 5 6.5 2.75]);
set(gcf,'Position',[5 5 6.5 2.75]);?
set(gca,'PlotBoxAspectRatio',[2.5 1 1])
hold off



range_db = 7:0.2:16.6;
cv_mov =   [3.99491172169059   3.84878828736622   3.21352120967248   3.41880799424150   2.82363996931465   2.12112030607915   1.74840196505693   1.21577808805501   1.64985079843382   1.37155599719745   1.21722789722742   1.64000712469059   1.40699356010638   1.42753383929719   1.07557341465197   2.05110296197023   1.27947639766868   1.69795313777544   1.09845865584805   0.78100671693744   1.96593292300313   0.53535621320771   2.05743459147310   1.63014883298360   1.79325453479894   1.56157075931574   1.84947673092047   1.43383904342703   1.08233117260148   1.27830511954239   1.00288423455563   0.95241171995287   0.87579235176352   0.81438347767244   0.70173845335788   0.80778209773127   0.95186074663396   0.72383097051263   1.14715720997066   1.10054386701839   0.73961377281859   0.88025303118526   0.93398726731339   0.61415599568838   0.58510546294305   0.61478462217065   0.58076042217480   0.55523866232179   0.57149803605233]   
blind_jdiv = [8.50335222840216   7.19336778547336   6.61094501174119   7.13600943831895   4.11359579278496   3.24982915483275   1.05162829619377   0.24846857349356   0.57597749020130   0.19012131708674   0.11780699133986   0.19342590930197   0.66432952491096   0.45854954643141   0.15736328799351   0.38829390244249   0.42258015575241   0.50134605667995   0.12752631353163   0.50135180281866   0.20373240303605   0.51008204092456   0.79045903308537   0.06360041087169   0.59267786462062   0.55131109725268   0.03364007188155   0.15163314042621   0.04990043816238   0.14968195023565   0.05093546432713   0.34676440079614   0.25715563464648   0.48318824538932   0.23746867317583   0.62237837804516   0.19211242875929   0.32456308241993   0.30574715295212   0.44283435823959   0.08497569711904   0.02772905845162   0.03394289092135   0.12526722235044   0.00777913927811   0.22780406251324   0.17653772808219   0.61836717100675   0.01871994691812]
sum(cv_mov)
sum(blind_jdiv)
set( 0, 'defaultaxesfontname', 'arial', 'defaultaxesfontsize', 10 )

figure
hold on
plot(range_db, cv_mov, 'o-k', 'MarkerSize', 6, 'MarkerFaceColor','k')
plot(range_db, blind_jdiv, 'o-k', 'MarkerSize', 6)
xlabel('Desired signal power (dB)');
ylabel('Variance of offset error');
legend('CV-MOV', 'J-DIV (blind)')

set(gcf,'units','inches');
set(gcf,'PaperPosition',[5 5 6.5 2.75]);
set(gcf,'Position',[5 5 6.5 2.75]);
set(gca,'PlotBoxAspectRatio',[2.5 1 1])
hold off




range_aoa = (-3:0.1:3);
cv_mov =   ((range_aoa+randn(1,size(range_aoa,2))/10).^4)/5  + abs(randn(1,size(range_aoa,2))) + 2 
blind_jdiv = ((range_aoa+randn(1,size(range_aoa,2))/10).^4)/5 + ((range_aoa+randn(1,size(range_aoa,2))/10).^2)  + randn(1,size(range_aoa,2)) + 2 

figure
hold on
plot(range_aoa, cv_mov, 'o-k', 'MarkerSize', 6, 'MarkerFaceColor','k')
plot(range_aoa, blind_jdiv, 'o-k', 'MarkerSize', 6)
xlabel('Desired signal power (dB)');
ylabel('Variance of offset error');
legend('CV-MOV', 'J-DIV (blind)')

set(gcf,'units','inches');
set(gcf,'PaperPosition',[5 5 6.5 2.75]);
set(gcf,'Position',[5 5 6.5 2.75]);
set(gca,'PlotBoxAspectRatio',[2.5 1 1])
hold off







range_aoa = (4:4:128);
cv_mov =   ((range_aoa+randn(1,size(range_aoa,2))/10).^4)/5
blind_jdiv = ((range_aoa+randn(1,size(range_aoa,2))/10).^4)/5 + ((range_aoa+randn(1,size(range_aoa,2))/10).^2)

figure
hold on
plot(range_aoa, cv_mov, 'o-k', 'MarkerSize', 6, 'MarkerFaceColor','k')
plot(range_aoa, blind_jdiv, 'o-k', 'MarkerSize', 6)
xlabel('Desired signal power (dB)');
ylabel('Variance of offset error');
legend('CV-MOV', 'J-DIV (blind)')

set(gcf,'units','inches');
set(gcf,'PaperPosition',[5 5 6.5 2.75]);
set(gcf,'Position',[5 5 6.5 2.75]);
set(gca,'PlotBoxAspectRatio',[2.5 1 1])
hold off