generate_plots.m

来自「implement basic routines in Matlab and O」· M 代码 · 共 87 行

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