% Parameters
d1 = 258e-9; % meters
n1 = 1.45;
d2 = 163e-9; % meters
n2 = 2.3;
totalLength = 200e-6; % meters
layerCount = 500; % Approximate layers for 200 micrometers
frequencyRange = 100e12:1e12:200e12; % 100 to 200 THz, step by 1 THz for smooth plot
% Preallocate transmission array
transmission = zeros(length(frequencyRange), 1);
% Loop over each frequency
for i = 1:length(frequencyRange)
frequency = frequencyRange(i);
% Initialize total matrix as identity
M_total = eye(2);
% Generate random sequence of M1 and M2 multiplications
for j = 1:layerCount
if rand > 0.5
phi_1 = d1 * n1 * frequency;
M = [cos(phi_1), sin(phi_1); sin(phi_1), -cos(phi_1)];
else
phi_2 = d2 * n2 * frequency;
M = [cos(phi_2), sin(phi_2); sin(phi_2), -cos(phi_2)];
end
M_total = M_total * M;
end
% Calculate transmission
T = 1 / M_total(1,1);
transmission(i) = T;
end
% Plot Transmission vs Frequency
figure;
plot(frequencyRange * 1e-12, abs(transmission)); % Convert frequency to THz for plotting
xlabel('Frequency (THz)');
ylabel('Transmission');
title('Transmission vs Frequency');
% Plot Average Transmission vs Total Length
% Since total length is constant in this scenario, we show how to calculate the average
mean_T = mean(abs(transmission));
figure;
plot(totalLength * 1e6, mean_T, 'ro'); % Convert meters to micrometers for plotting
xlabel('Total Length (micrometers)');
ylabel('Average Transmission');
title('Average Transmission vs Total Length');
ylim([0, max(abs(transmission)) + 0.1]); % Adjust ylim to make the single point more visible
% Note: The calculation of average transmission vs total length doesn't vary with length in this script,
% as we've kept the length constant at 200 micrometers. This plot will just illustrate the average transmission for the given length.
-----------------------------------------------------------------------------------------------------------------------------------------------------
If you find the solution helpful and it resolves your issue, it would be greatly appreciated if you could accept the answer. Also, leaving an upvote and a comment are also wonderful ways to provide feedback.
It's important to note that the advice and code are based on limited information and meant for educational purposes. Users should verify and adapt the code to their specific needs, ensuring compatibility and adherence to ethical standards.
Professional Interests
- Technical Services and Consulting
- Embedded Systems | Firmware Developement | Simulations
- Electrical and Electronics Engineering
Feel free to contact me.
