?? rwg3.m
字號:
%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
%Load the data
load('mesh2');
%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');
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -