% seti contrast Fresnel Synth % is used in tests/fresnelSynth.m % -- setGrid seti.dim = 2; seti.rCD = 0.2; seti.nCD = 256; % original % -- fresnel data (loadData.m) seti.expData = 'fresnel'; seti.fresnelFreq = 5*1E9; % standard 3 GHz (4 GHz) seti.fresnelFile = 'inexpdata/fresnel_opus_1/dielTM_dec8f.exp'; % 1--8 GHz, test: single cylinder, 4 GHz %seti.fresnelFile = 'inexpdata/fresnel_opus_1/twodielTM_8f.exp'; % 1--8 GHz filename = seti.fresnelFile; seti.nuMax = 7; seti.ampCalc = 1; % default (A1) %seti.nuMax = 7; seti.ampCalc = 2; % test (A2) %seti.nuMax = 7; seti.ampCalc = 3; % test (A3) seti.nuMax = 10; seti.ampCalc = 1; % default (B1) %seti.nuMax = 10; seti.ampCalc = 2; % test (B2) %seti.nuMax = 10; seti.ampCalc = 3; % test (B3) % 3 GHz: single dielectricum % B1: 1.3 & 14.9 & 15.3 % B2: 1.8 & 15.0 & 15.7 % B3 (with 2E3): 1.7 & 15.0 & 15.7 % B3 (with 1E5): 1.6 & 14.9 & 15.7 % B3 (with 1E8): 1.5 & 14.9 & 15.5 % 5 GHz: single dielectricum % B1: 3.4 & 20.1 & 22.4 % B2: 21.2 & 20.5 & 26.9 % B3 (with 2E3): 18.9 & 19.6 & 24.7 % B3 (with 1E5): 15.2 & 19.1 & 23.1 % single dielectricum with method (B1) %1\,GHz & 0.9 & 29.8 & 39.7\\ %2\,GHz & 0.6 & 18.0 & 14.9\\ %3\,GHz & 1.3 & 14.9 & 15.3\\ %4\,GHz & 0.8 & 25.2 & 19.7\\ %5\,GHz & 3.4 & 20.1 & 22.4\\ %6\,GHz & 4.2 & 39.7 & 32.7\\ %7\,GHz & 16.4 & 22.2 & 28.1\\ %8\,GHz & 18.1 & 44.3 & 39.1 % two dielectrics with method (B1) %1\,GHz & 1.0 & 11.5 & 64.5\\ %2\,GHz & 0.6 & 46.7 & 8.1\\ %3\,GHz & 1.3 & 13.4 & 21.6\\ %4\,GHz & 0.8 & 47.9 & 15.6\\ %5\,GHz & 3.3 & 20.3 & 27.4\\ %6\,GHz & 4.0 & 59.5 & 36.8\\ %7\,GHz & 16.1 & 25.2 & 42.9\\ %8\,GHz & 17.8 & 72.2 & 41.1 %----- % -- setKernel cLight = 2.99792458E8; % light velocity seti.k = 2*pi*seti.fresnelFreq/cLight; seti.model = 'helmholtz'; % -- setMeas (measurements) % simulate fresnel_opus_1: seti.incType = 'pointSource'; % used in op. 1 and 2 seti.measType = 'nearField'; % used in op. 1 and 2 seti.incNb = 36; seti.measNb = 72; seti.radSrc = 0.72; seti.radMeas = 0.76; % -- inversion process (regularization) %seti.invNo = 3; seti.invNo = 6; seti.delta = 0; % to test %seti = setiAuxRec(seti,usefindpdatau,usevaralpha); seti.alpha = 1E8; % old: 2E4 corresponding 0.05 (in old code) seti.beta = 0; % because we expect no sharp edges... seti.physBounds = [-500,500,0,inf]; seti.nOut = 30; % -- minimization with pda seti.pdaN = 50;