6504阅读
14回复

[资料共享]矩量法计算带状线特征阻抗Matlab代码

离线febi

UID ：2537

注册：2007-05-11
登录：2016-05-05
发帖：123
等级：退休版主

0楼发表于: 2008-01-15 01:39:48

% RJA - 2/23/07
%
% Translation of FORTRAN Code given in Fig. 5.9 to MATLAB. No attempt
% whatsoever was made to make the code efficient. It is a pure
% translation.

clear all; format compact;

%Compare to Table 5.2 page 313. These values were obtained with this
%program. There is obviously a factor of 2 error between this and the
%table but the source of the error is unknown. Several modification to the
%translation were made and noted in the code. Note that table 5.2 trends
%down while this program trends up as n increases.

% n Zo Zo/2
%-------------------------------
% 3 94.4473 47.2236
% 7 96.8223 48.4112
% 11 97.4035 48.7017
% 18 97.823 48.9115
% 39 98.23 49.115
% 59 98.3686 49.1843
% 100 98.49305 49.24652
% 1000 98.696565 49.348283

% *****************************************************************
% USING MOMENTS METHOD,
% THIS PROGRAM DETERMINES THE CHARACTERISTIC IMPEDANCE
% OF A STRIP TRANSMISSION LINE
% WITH CROSS-SECTION W X H
% THE STRIPS MAITAINED AT 1 VOLT AND -1 VOLT.
%
% ONE STRIP IS LOCATED ON THE Y = H/2 PLANE WHILE THE OTHER
% IS LOCATED ON THE Y = -H/2 PLANE.
%
% ALL DIMENSIONS ARE IS S.I. UNITS
%
% N IS THE NUMBER OF SUBSECTIONS INTO WHICH EACH STRIP IS DIVIDED
% *****************************************************************

% FIRST, SPECIFY THE PARAMETERS

NN = [3 7 11 18 39 59];

disp(' ')
disp('n Zo')

for N = NN
CL=3.0e8; % SPEED OF LIGHT IN FREE SPACE
ER=1.0;
EO=8.8541878176E-12;
H=2.0;
W=5.0;

NT=2*N;
DELTA = W/(N);

% SECOND, CALCULATE THE MATCH POINTS
% AND THE COEFFICIENT MATRIX [A]

for K=1:N
X(K) = DELTA*(K - .5);
Y(K) = -H/2.0;
X(K+N) = X(K);
Y(K+N) = H/2.0;
end

FACTOR = DELTA/(2.0*pi*EO);

for I=1:NT
for J=1:NT
if(I==J) %eqn (5.98)
A(I,J) = -(log(DELTA) - 1.5)*FACTOR;
else
R = sqrt( (X(I) - X(J))^2 + (Y(I) - Y(J))^2 );
A(I,J) = -log(R)*FACTOR;
end
end
end

% NOW DETERMINE THE MATRIX OF CONSTANT VECTOR [B]

for K=1:N
B(K)=1.0;
%The line below was commented!
B(K+N)= -1.0; %rja -this line needed to be un-commented
end

% INVERT MATRIX A(I,J) AND CALCULATE MATRIX RO(N)
% CONSISTING OF THE UNKNOWN ELEMENTS
% ALSO CALCULATE THE TOTAL CHARGE Q,
% THE CAPACITANCE C, AND THE CHARACTERISTIC IMPEDANCE Zo.

NIV=NT;
NMAX=100;
A = inv(A);
for I=1:NT
RO(I)=0.0;
%The line below had an error N instead of NT
for M=1:NT
RO(I)=RO(I) + A(I,M)*B(M);
end
end
SUM=0.0;
for I=1:N
SUM = SUM + RO(I);
end
Q=SUM*DELTA;
VD=2.0; & ..

未注册仅能浏览部分内容，查看全部内容及附件请先登录或注册

共1条评分

casey

rf币 +2

您热心助人+2

2008-01-15

磁砖

离线eloquent

UID ：1383