www.pudn.com > loopantenna.rar > rwg3.m, change:2002-04-17,size:1737b

```%RWG3 Calculates the impedance matrix using function IMPMET
%   Uses the mesh file from RWG2, mesh2.mat, as an input.
%
%   The following parameters need to be specified prior to
%   calculations:
%
%   Frequency (Hz)                  f
%   Dielectric constant (SI)        epsilon_
%   Magnetic permeability (SI)      mu_
%
%   Copyright 2002 AEMM. Revision 2002/03/11
%   Chapter 2

clear all

%EM parameters (f=3e8 means that f=300 MHz)
%f=4.8e6 means that C=lambda/10 (loop with a =1 m)
f           =4.8e6;
epsilon_    =8.854e-012;
mu_         =1.257e-006;
%Speed of light
c_=1/sqrt(epsilon_*mu_);
%Free-space impedance
eta_=sqrt(mu_/epsilon_);

%Contemporary variables - introduced to speed up
%the impedance matrix calculation
omega       =2*pi*f;
k           =omega/c_;
K           =j*k;
Constant1   =mu_/(4*pi);
Constant2   =1/(j*4*pi*omega*epsilon_);
Factor      =1/9;

FactorA     =Factor*(j*omega*EdgeLength/4)*Constant1;
FactorFi    =Factor*EdgeLength*Constant2;

for m=1:EdgesTotal
RHO_P(:,:,m)=repmat(RHO_Plus(:,m),[1 9]);   %[3 9 EdgesTotal]
RHO_M(:,:,m)=repmat(RHO_Minus(:,m),[1 9]);  %[3 9 EdgesTotal]
end
FactorA=FactorA.';
FactorFi=FactorFi.';

%Impedance matrix Z
tic; %start timer

Z=  impmet( EdgesTotal,TrianglesTotal,...
EdgeLength,K,...
Center,Center_,...
TrianglePlus,TriangleMinus,...
RHO_P,RHO_M,...
RHO__Plus,RHO__Minus,...
FactorA,FactorFi);

toc %elapsed time

%Save result
FileName='impedance.mat';
save(FileName, 'f','omega','mu_','epsilon_','c_', 'eta_','Z');            ```