Prechádzať zdrojové kódy

add polar plots, only print high frequency of cardinal plots, add options (not automatic) for color and monochromatic printing

Fernando Pablo Lopez-Lezcano 10 rokov pred
rodič
commit
c394a961ee
1 zmenil súbory, kde vykonal 18 pridanie a 18 odobranie
  1. 18 18
      make_figures_aes.m

+ 18 - 18
make_figures_aes.m

@@ -4,7 +4,7 @@
 %% Copyright Fernando Lopez-Lezcano 2016
 %% nando@ccrma.stanford.edu
 
-function make_figures_aes(col)
+function [FIR] = make_figures_aes(col)
 
   %% start from scratch and load common functions
   close all;
@@ -13,28 +13,23 @@ function make_figures_aes(col)
 
   global AESDIR = "../papers/AES_SFCconf2016/AES_SFCconf2016_template_abstractOnly/";
 
-  %% select color or mono
-%  if (! exist(col))
-%    col = 0;
-%  end
-%  if (col == 1)
-    %% color figures
-%    cext = "_color";
-%    copt = "epsc";
-%  else
-    %% mono figures
-    cext = "_mono";
-    copt = "eps";
-%  end
-  
+  %% color figures
+  %% cext = "_color";
+  %% copt = "epsc";
+  %% mono figures
+  cext = "_mono";
+  copt = "eps";
+    
   %% frequency response with just A2B matrices
   [F FE Fs A2BF] = sphear_a2b_all((1:16), 300, 20000, 3, 256, 32);
+  %% deliberately use nyquist to extend this to high frequencies
   FIR = sphear_fir_design(F, A2BF, 256, 20, Fs, 24000, copt);
   sphear_fir_render(FIR, 300, 20000, 256, 32, 8192, copt);
   close all;
 
   system(["cp -p figures/sphear_bf_azimuth_000" cext ".eps " AESDIR "sphear_bf_azimuth_000_a2b_all.eps"]);
   system(["cp -p figures/sphear_bf_azimuth_045" cext ".eps " AESDIR "sphear_bf_azimuth_045_a2b_all.eps"]);
+  system(["cp -p figures/wxyz_filter_lf" cext ".eps " AESDIR "wxyz_filter_a2b_only.eps"]);
 
   %% frequency response with just the A2B matrices and limit frequency
   FIR = sphear_fir_design(F, A2BF, 256, 20, Fs, 4000, copt);
@@ -52,7 +47,7 @@ function make_figures_aes(col)
   [W X Y Z] = sphear_fir_extract(B, 0, 1);
   FIR = sphear_fir_design(F, A2BF, 256, 20, Fs, 4000, copt, W, X, Y);
   sphear_fir_render(FIR, 300, 20000, 256, 32, 8192, copt);
-  sphear_fir_render_directions(FIR, 300, 20000, 256, 32, 4096, copt);
+  sphear_fir_render_directions(FIR, 3000, 20000, 256, 32, 4096, copt);
   close all;
 
   system(["cp -p figures/sphear_bf_azimuth_000" cext ".eps " AESDIR "sphear_bf_azimuth_000.eps"]);
@@ -67,12 +62,17 @@ function make_figures_aes(col)
   [W X Y Z] = sphear_fir_extract(B, 1, 1);
   FIR = sphear_fir_design(F, A2BF, 256, 20, Fs, 4000, copt, W, X, Y);
   sphear_fir_render(FIR, 300, 20000, 256, 32, 8192, copt);
-  sphear_fir_render_directions(FIR, 300, 20000, 256, 32, 4096, copt);
+  sphear_fir_render_directions(FIR, 3000, 20000, 256, 32, 4096, copt);
+  [B] = sphear_a2b_fir_all(300, 20000, 3, 8192, 48000, 1);
   close all;
 
   system(["cp -p figures/sphear_bf_azimuth_000" cext ".eps " AESDIR "sphear_bf_azimuth_000_all.eps"]);
   system(["cp -p figures/sphear_bf_azimuth_045" cext ".eps " AESDIR "sphear_bf_azimuth_045_all.eps"]);
   system(["cp -p figures/wxyz_filter_ff" cext ".eps " AESDIR "wxyz_filter_ff_all.eps"]);
   system(["cp -p figures/sphear_directions.eps " AESDIR "sphear_directions_all.eps"]);
-
+  system(["cp -p figures/sphear_polar_multi_00604_00762_mono.eps " AESDIR "sphear_polar_multi_600.eps"]);
+  system(["cp -p figures/sphear_polar_multi_01216_01536_mono.eps " AESDIR "sphear_polar_multi_1200.eps"]);
+  system(["cp -p figures/sphear_polar_multi_04932_06228_mono.eps " AESDIR "sphear_polar_multi_5000.eps"]);
+  system(["cp -p figures/sphear_polar_multi_09932_12542_mono.eps " AESDIR "sphear_polar_multi_10000.eps"]);
+  
 end