karthik/Calculations
1
0
Fork 0
Code Issues Pull Requests repo.project_board Releases Wiki Activity

Minor modifications

Browse Source
This commit is contained in:
Karthik 2025-06-10 18:42:37 +02:00
parent 8ff8fbae58
commit b5169b7e97
3 changed files with 13 additions and 13 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
Data-Analyzer/analyzePhaseTransitionSimulation.m
View File

@ -1,6 +1,6 @@
%% Extract Images %% Extract Images
baseDir = 'D:/Results - Numerics/Data_Full3D/PhaseTransition/DTS/'; baseDir = 'D:/Results - Numerics/Data_Full3D/PhaseTransition/DTS/Point2';
JobNumber = 0; JobNumber = 0;
runFolder = sprintf('Run_%03d', JobNumber); runFolder = sprintf('Run_%03d', JobNumber);
movieFileName = 'DropletsToStripes.mp4'; % Output file name movieFileName = 'DropletsToStripes.mp4'; % Output file name
@ -10,7 +10,7 @@ reverseOrder = false; % Set this to true to reverse the theta ordering
TitleString = 'Change across transition: Droplets To Stripes'; TitleString = 'Change across transition: Droplets To Stripes';
%% %%
baseDir = 'D:/Results - Numerics/Data_Full3D/PhaseTransition/STD/'; baseDir = 'D:/Results - Numerics/Data_Full3D/PhaseTransition/STD/Point2';
JobNumber = 0; JobNumber = 0;
runFolder = sprintf('Run_%03d', JobNumber); runFolder = sprintf('Run_%03d', JobNumber);
movieFileName = 'StripesToDroplets.mp4'; % Output file name movieFileName = 'StripesToDroplets.mp4'; % Output file name
@ -184,6 +184,7 @@ for k = 1:nimgs
x_min = min(x); x_min = min(x);
x_max = max(x); x_max = max(x);
% Display the cropped OD image % Display the cropped OD image
ax1 = nexttile; ax1 = nexttile;
imagesc(x, y, IMG') imagesc(x, y, IMG')
@ -234,8 +235,7 @@ for k = 1:nimgs
% Plot the angular distribution % Plot the angular distribution
nexttile nexttile
spectral_contrast(k) = computeSpectralContrast(fft_imgs{k}, 10, 25, Threshold); spectral_contrast(k) = computeSpectralContrast(fft_imgs{k}, 10, 35, Threshold);
[theta_vals, S_theta] = computeNormalizedAngularSpectralDistribution(fft_imgs{k}, 10, 25, N_bins, Threshold, Sigma);
spectral_weight(k) = trapz(theta_vals, S_theta); spectral_weight(k) = trapz(theta_vals, S_theta);
plot(theta_vals/pi, S_theta,'Linewidth',2); plot(theta_vals/pi, S_theta,'Linewidth',2);
axis square; axis square;

2
Data-Analyzer/matchTFRadii.m
View File

@ -1,7 +1,7 @@
% ---------------------------- % ----------------------------
% Experimental data % Experimental data
% ---------------------------- % ----------------------------
PixelSize = 5.86; % microns PixelSize = 4.6; % microns
AtomNumbers = [9.14950197, 8.6845907 , 8.82521245, 8.5899089 , 8.21675841, ... AtomNumbers = [9.14950197, 8.6845907 , 8.82521245, 8.5899089 , 8.21675841, ...
8.96234044, 8.8636914 , 8.70332154, 8.82930706, 8.91869919, ... 8.96234044, 8.8636914 , 8.70332154, 8.82930706, 8.91869919, ...

16
Estimations/LinearizeScatteringLengthScan.m
View File

@ -1,5 +1,5 @@
% --- User setup --- % --- User setup ---
resIndex = 3; % resonance index (1-based) resIndex = 1; % resonance index (1-based)
duration = 0.5; % total ramp time [s] duration = 0.5; % total ramp time [s]
N_points = 300; % number of time points N_points = 300; % number of time points
FR_choice = 1; % resonance data choice (1=new, 0=old) FR_choice = 1; % resonance data choice (1=new, 0=old)
@ -16,14 +16,14 @@ doPlot = true; % show resonance and ramp region
figure; figure;
subplot(2,1,1); subplot(2,1,1);
plot(t, B_t, 'b-', 'LineWidth', 1.5); plot(t, B_t, 'b-', 'LineWidth', 1.5);
ylabel('B (G)'); ylabel('B (G)', 'Interpreter', 'tex');
title('Magnetic field B(t)'); title('Magnetic field B(t)', 'Interpreter', 'tex');
subplot(2,1,2); subplot(2,1,2);
plot(t, a_t, 'r-', 'LineWidth', 1.5); plot(t, a_t, 'r-', 'LineWidth', 1.5);
ylabel('Scattering length a (a_0)'); ylabel('Scattering length a ($a_0$)', 'Interpreter', 'latex');
xlabel('Time (s)'); xlabel('Time (s)', 'Interpreter', 'tex');
title('Linear ramp of a(t)'); title('Linear ramp of a(t)', 'Interpreter', 'tex');
%% Function: makeLinearScatteringRamp %% Function: makeLinearScatteringRamp
@ -134,8 +134,8 @@ if nargin >= 8 && doPlot
hold on; hold on;
plot(B_ramp, a_ramp, 'r-', 'LineWidth', 2, 'DisplayName', 'Ramp a(t)'); plot(B_ramp, a_ramp, 'r-', 'LineWidth', 2, 'DisplayName', 'Ramp a(t)');
xlabel('B (G)'); xlabel('B (G)');
ylabel('a (a_0)'); ylabel('a ($a_0$)', 'Interpreter', 'latex');
title(sprintf('Feshbach Resonance #%d at B0=%.3f G (%s side)', resIndex, B0, side)); title(sprintf('Feshbach Resonance #%d at B0=%.3f G (%s side)', resIndex, B0, side), 'Interpreter', 'tex');
legend('Location','best'); legend('Location','best');
grid on; grid on;
end end