LLT1 QP EMI Simulation Study
LLT1 QP EMI Simulation Study
(a) Consider a transmission line with characteristic impedance Z0 = 50 Ω and 7.5 4 CO4 3
various load impedances ZL= 25 Ω, 75 Ω and 1000 Ω. Write a MATLAB code and
estimate the reflection coefficient and VSWR for each load. Assuming the
operating frequency to be 1 GHz and Transmission line length to be 50 cm, plot the
2. Voltage standing wave pattern for each case and estimate VSWR from the plots.
Compare the VSWR obtained analytically and from the graphs.
(b) Estimate the power transfer efficiency for each case and show that the power 7.5 3 CO4 2
transfer is maximum when the characteristic impedance matches with the load
impedance.
Course Outcome (CO) and Bloom’s level (BL) Coverage in Questions
BL Coverage (%)
CO Coverage (%)
100 100
80 80
60 60
40 40
20 20
0 0
CO1 CO2 CO3 CO4 CO5 BL1 BL2 BL3 BL4
Evaluation Sheet
Consolidated Marks:
CO Maximum Marks PO Maximum Marks
Marks Obtained Marks Obtained
3 15 2 15
4 15 3 15
Total 30 Total 30
Students can refer to the following MATLAB codes for answering the Simulation study
questions. They also need to print and attach their MATLAB codes along with the
submission.
1. MATLAB Code for Question 1(a) and 1(b) (Students need to modify the code as per
their requirement)
% Parameters
a = ??; % Waveguide width (in meters)
b = ??; % Waveguide height (in meters)
lambda = ??; % Operating wavelength (in meters)
c = ??; % Speed of light in vacuum
% Magnetic field (H_z) and electric field (E_x, E_y) for TE10 mode
Hz_TE10 = sin(pi * X / a); % H_z for TE10
Ex_TE10 = -(pi / b) * sin(pi * X / a); % E_x for TE10
Ey_TE10 = zeros(size(X)); % E_y for TE10 is zero
% Magnetic field (H_z) and electric field (E_x, E_y) for TE01 mode
Hz_TE01 = sin(pi * Y / b); % H_z for TE01
Ex_TE01 = zeros(size(X)); % E_x for TE01 is zero
Ey_TE01 = -(pi / a) * sin(pi * Y / b); % E_y for TE01
2. MATLAB Code for Question 2(a) and 2(b) (Students need to modify the code as per
their requirement)
% Given values
Z0 = ??; % Characteristic impedance in Ohms
ZL_values = [??]; % Load impedances
frequency = ??; % Operating frequency (1 GHz)
c = ??; % Speed of light (m/s)
lambda = c / frequency; % Wavelength
line_length = ??; % Transmission line length (in m)
z = linspace(0, line_length, 500); % Points along the line
% Reflection Coefficient
Gamma = (ZL - Z0) / (ZL + Z0);
reflection_coefficients(i) = Gamma;
% VSWR
VSWR = (1 + abs(Gamma)) / (1 - abs(Gamma));
VSWR_values(i) = VSWR;
% Display results
fprintf('For Z_L = %d Ohms:\n', ZL);
fprintf('Reflection Coefficient (Gamma) = %.3f\n', Gamma);
fprintf('VSWR = %.3f\n', VSWR);
fprintf('Power Transfer Efficiency = %.3f%%\n\n', efficiency * 100);
end
% Plot
figure;
plot(z, V_magnitude, 'LineWidth', 2);
xlabel('Position along the line (m)');
ylabel('Voltage Magnitude |V(z)|');
title('Standing Wave Pattern for Z_L = ?? Ohms');
grid on;