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

Aesthetic changes.

Browse Source
This commit is contained in:
Karthik 2025-02-20 22:43:15 +01:00
parent b2a5b172d8
commit 88c0d18066
3 changed files with 10 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
Data-Analyzer/plotImages.m
View File

@ -16,7 +16,7 @@ span = [50, 50];
fraction = [0.1, 0.1]; fraction = [0.1, 0.1];
pixel_size = 5.86e-6; pixel_size = 5.86e-6;
remove_fringes = false; removeFringes = true;
% Compute OD image, rotate and extract ROI for analysis % Compute OD image, rotate and extract ROI for analysis
% Get a list of all files in the folder with the desired file name pattern. % Get a list of all files in the folder with the desired file name pattern.
@ -42,7 +42,7 @@ for k = 1 : length(files)
end end
% Fringe removal % Fringe removal
if remove_fringes if removeFringes
optrefimages = removefringesInImage(absimages, refimages); optrefimages = removefringesInImage(absimages, refimages);
absimages_fringe_removed = absimages(:, :, :) - optrefimages(:, :, :); absimages_fringe_removed = absimages(:, :, :) - optrefimages(:, :, :);
@ -80,15 +80,15 @@ for k = 1 : length(od_imgs)
imagesc(x_range, y_range, od_imgs{k}) imagesc(x_range, y_range, od_imgs{k})
axis equal tight; axis equal tight;
hcb = colorbar; hcb = colorbar;
hL = ylabel(hcb, 'Optical Density', 'FontSize', 16); hL = ylabel(hcb, 'Normalised Optical Density', 'FontSize', 16);
set(hL,'Rotation',-90); set(hL,'Rotation',-90);
colormap jet; colormap jet;
set(gca,'CLim',[0 1.0]); set(gca,'CLim',[0 1.0]);
set(gca,'YDir','normal') set(gca,'YDir','normal')
set(gca, 'YTick', linspace(y_min, y_max, 5)); % Define y ticks set(gca, 'YTick', linspace(y_min, y_max, 5)); % Define y ticks
set(gca, 'YTickLabel', flip(linspace(y_min, y_max, 5))); % Flip only the labels set(gca, 'YTickLabel', flip(linspace(y_min, y_max, 5))); % Flip only the labels
xlabel('X', 'Interpreter', 'tex','FontSize',16); xlabel('Horizontal', 'Interpreter', 'tex','FontSize',16);
ylabel('Y', 'Interpreter', 'tex','FontSize',16); ylabel('Vertical', 'Interpreter', 'tex','FontSize',16);
drawnow drawnow
pause(0.5) pause(0.5)

6
Estimations/RotonInstability/ExtractingParameters.m
View File

@ -161,8 +161,8 @@ data_struct = struct;
wz_values = [500, 750, 1000, 2000]; wz_values = [500, 750, 1000, 2000];
kvec = linspace(0, 15e6, 1000); % Vector of magnitudes of k vector kvec = linspace(0, 15e6, 1000); % Vector of magnitudes of k vector
theta_values = 0:5:90; % Range of theta values theta_values = 0:1:7; % Range of theta values
phi = 90; % Azimuthal angle of momentum vector phi = 0; % Azimuthal angle of momentum vector
for mainloop_idx = 1:length(wz_values) for mainloop_idx = 1:length(wz_values)
format long format long
@ -216,7 +216,7 @@ for mainloop_idx = 1:length(wz_values)
end end
save('.\Results\ExtractingParameters_Result_FixedDensity_phi90.mat', 'data_struct'); save('.\Results\ExtractingParameters_Result_FixedDensity_phi0.mat', 'data_struct');
%% %%
function [Go,gamma4,Fka,Ukk] = computePotentialInMomentumSpace(k, gs, gdd, MeanWidth, theta, phi) function [Go,gamma4,Fka,Ukk] = computePotentialInMomentumSpace(k, gs, gdd, MeanWidth, theta, phi)
Go = sqrt(pi) * (k * MeanWidth/sqrt(2)) .* exp((k * MeanWidth/sqrt(2)).^2) .* erfc((k * MeanWidth/sqrt(2))); Go = sqrt(pi) * (k * MeanWidth/sqrt(2)) .* exp((k * MeanWidth/sqrt(2)).^2) .* erfc((k * MeanWidth/sqrt(2)));

4
Imaging-Response-Function-Extractor/extractIRF.m
View File

@ -20,7 +20,7 @@ lambda = 421e-9;
d = 0.61*lambda/NA; d = 0.61*lambda/NA;
k_cutoff = 2*NA/lambda/1e6; % in units of 1/µm) k_cutoff = 2*NA/lambda/1e6; % in units of 1/µm)
removeFringe = true; removeFringes = true;
%% Compute OD image, rotate and extract ROI for analysis %% Compute OD image, rotate and extract ROI for analysis
% Get a list of all files in the folder with the desired file name pattern. % Get a list of all files in the folder with the desired file name pattern.
@ -47,7 +47,7 @@ end
%% Fringe removal %% Fringe removal
if removeFringe if removeFringes
optrefimages = removefringesInImage(absimages, refimages); optrefimages = removefringesInImage(absimages, refimages);
absimages_fringe_removed = absimages(:, :, :) - optrefimages(:, :, :); absimages_fringe_removed = absimages(:, :, :) - optrefimages(:, :, :);