design broadband matching networks for antennas -凯发k8网页登录

specify frequency and impedance

specify the center frequency, 350 mhz, and bandwidth, 110 mhz, of match to build a matching network with a bandpass response.

fc = 350e6;
bw = 110e6;

specify the source impedance, the reference impedance and the load resistance. in this example the load zl is modeled as a series r-l circuit. instead of calculating the load impedance, you could measure the impedance of the load.

zs = 50;                              % source impedance (ohm)              
z0 = 50;                              % reference impedance (ohm)
rl = 40;                              % load resistance (ohm)
l = 12e-8;                            % load inductance (henry)

define the number of frequency points to use for analysis and set up the frequency vector.

nfreq = 256;                          % number of frequency points
flower = fc - (bw/2);                 % lower band edge                 
fupper = fc   (bw/2);                 % upper band edge                
freq = linspace(flower,fupper,nfreq); % frequency array for analysis
w = 2*pi*freq;                        % frequency (radians/sec)

understand load behavior using reflection coefficient and power gain

use two simple expressions for calculating the load reflection coefficient and the power gain. this corresponds to directly connecting the source to the input terminals of an antenna i.e. in figure 1 there is no matching network.

xl = w*l;                                   % reactance (ohm)
zl = rl   1i*xl;                            % load impedance (ohm)
gammal = (zl - z0)./(zl   z0);              % load reflection coefficient
unmatchedgt = 10*log10(1 - abs(gammal).^2); % power delivered to load

use the function to plot the variation in the load reflection coefficient with frequency. an input reflection coefficient closer to center of the smith chart denotes a better matching performance.

figure
smithplot(freq,gammal,'legendlabels','#gamma load','linewidth',2,...
    'view','top-right');

this plot shows that the load reflection coefficient is far away from this point. therefore, there is an impedance mismatch. you can confirm this mismatch by plotting the transducer gain as a function of frequency.

figure
plot(freq.*1e-6,unmatchedgt,'r')
grid on;
title('power delivered to load - no matching network');
xlabel('frequency (mhz)');
ylabel('magnitude (decibels)');
legend('g_t','location','best');

as the plot shows, there is approximately 10 db power loss around the desired region of operation (295 - 405 mhz). as a result, the antenna needs a matching network that operates over a 110 mhz bandwidth centered at 350 mhz.

design matching network

the matching network must operate between 295 mhz and 405 mhz, therefore you choose a bandpass topology for the matching network shown below.

type - i: series lc first element followed by shunt lc

figure 2: matching network topology

the approach is to design an odd order 0.5 db chebyshev bandpass to obtain the initial design for the matching network shown in figure 2. this is a single match problem , i.e. the source is purely resistive while load is a combination of r and l, solution is you can begin by choosing a five-element prototype network.

n = 5; % order of matching network
filter = rffilter('filtertype',"chebyshev","filterorder",n, ...
    "implementation","lc tee","responsetype","bandpass",...
    "passbandfrequency",[flower fupper],"passbandattenuation",0.5);
lvals = filter.designdata.inductors;

use the object to build the bandpass tee matching network. note that the topology demands a bandpass tee prototype that begins with a series inductor. if the topology chosen is an lc bandpass pi then you would begin with shunt c for the lowpass prototype.

% create the matching network
matchingnw = lcladder(filter);
% copy initial values for comparison
l_initial = lvals;

optimize designed matching network

there are several factors to consider prior to the optimization.

the objective function used during the optimization process by fminsearch is shown here.

type('antennamatchobjectivefun.m')

function output = antennamatchobjectivefun(matchingnw,lvalues,freq,zl,z0)
%antennamatchobjectivefun is the objective function used by the example
% designing broadband matching networks for antennas.
%
% output = antennamatchobjectivefun(matchingnw,lvalues,freq,z0)
% returns the current value of the objective function stored in output
% evaluated after updating the inductor values in the object, matchingnw.
% the inductor values are stored in the variable lvalues.
%
% antennamatchobjectivefun is an objective function of rf toolbox demo:
% designing broadband matching networks (part i: antenna)
%   凯发官网入口首页 copyright 2008-2020 the mathworks, inc.
% ensure positive element values
if any(lvalues <= 0)                                                       
    output = inf;
    return
end
% update the element values in the matching network
matchingnw.inductances(1) = lvalues(1);
matchingnw.inductances(end) = lvalues(end);
% perform analysis on tuned matching network
s = sparameters(matchingnw,freq,z0);
% calculate input reflection coefficient 'gammain'
gin = gammain(s,zl);
% cost function
output = mean(abs(gin));
% other possible choices for objective function could be : -
% output =  max(abs(gin));
% output = -1*mean(gt_pass);
% animate
smithplot(freq,gin);
drawnow

there are several ways to choose the cost function and some options are shown within the objective function above (in comments). the optimization variables are the first and last inductors, l1 and l5 respectively. the element values are stored in the variable l_optimized.

niter = 125;
options = optimset('display','iter','maxiter',niter);  % set options structure
l_optimized = [lvals(1) lvals(end)];
l_optimized = ...
    fminsearch(@(l_optimized)antennamatchobjectivefun(matchingnw, ...
    l_optimized,freq,zl,z0),l_optimized,options);

 
 iteration   func-count         f(x)         procedure
     0            1         0.933982         
     1            3         0.933982         initial simplex
     2            5         0.920323         expand
     3            7         0.911353         expand
     4            9         0.853255         expand
     5           11         0.730444         expand
     6           13         0.526448         reflect
     7           15         0.526448         contract inside
     8           17         0.421103         reflect
     9           19         0.421103         contract inside
    10           20         0.421103         reflect
    11           22         0.421103         contract inside
    12           24         0.421103         contract inside
    13           26         0.339935         expand
    14           27         0.339935         reflect
    15           29          0.28528         reflect
    16           31          0.28528         contract inside
    17           32          0.28528         reflect
    18           34         0.283527         reflect
    19           36         0.283527         contract inside
    20           38         0.278939         contract inside
    21           40         0.278123         reflect
    22           41         0.278123         reflect
    23           43          0.27636         contract inside
    24           45         0.275782         contract inside
    25           47         0.275637         contract inside
    26           49         0.275498         reflect
    27           51         0.275282         contract inside
    28           52         0.275282         reflect
    29           54         0.275282         contract inside

    30           56         0.275282         contract inside
 
optimization terminated:
 the current x satisfies the termination criteria using options.tolx of 1.000000e-04 
 and f(x) satisfies the convergence criteria using options.tolfun of 1.000000e-04 

update matching network elements with optimal values

when the optimization routine stops, the optimized element values are stored in l_optimized. the following code updates the input and output matching network with these values.

matchingnw.inductances(1) = l_optimized(1);     % update the matching network inductor l1
matchingnw.inductances(end) = l_optimized(end); % update the matching network inductor l5

analyze and display optimization results

compare and plot the input reflection coefficient of the matched and unmatched results.

s = sparameters(matchingnw,freq,z0);
gin = gammain(s,zl);
smithplot(freq,[gin transpose(gammal)],'legendlabels',...
    {'#gamma in (matched)','#gamma in (unmatched)'})

the optimized matching network improves the performance of the circuit. in the passband (295 - 405 mhz), the input reflection coefficient is closer to the center of the smith chart. plot the power delivered to load for both the matched and unmatched system.

matchedgt = powergain(s,zs,zl,'gt');
figure;
plot(freq*1e-6,matchedgt)
hold all;
plot(freq*1e-6,unmatchedgt,'r')
grid on;
hold off;
title('power delivered to load');
legend('optimized network','no matching network','location','best');

the power delivered to the load is approximately 1 db down for the optimized matching network.

display optimized element values

the following code shows the initial and optimized values for inductors l1 and l5.

l1_initial = l_initial(1)

l1_initial = 1.2340e-07

l1_optimized = l_optimized(1)

l1_optimized = 1.2111e-07

l5_initial = l_initial(end)

l5_initial = 1.2340e-07

l5_optimized = l_optimized(end)

l5_optimized = 1.7557e-09

there are few things to consider when setting up an optimization:

references

see also

design broadband matching networks for amplifier.