MATLAB之Impedance

%***********************************************************************
% IMPEDANCE.m
%**********************************************************************
%
% SELF AND MUTUAL IMPEDANCES
%******************************************************************
% THIS IS A MATLAB BASED PROGRAM THAT COMPUTES:
%
% I. SELF IMPEDANCE OF ANY LENGTH DIPOLE
% II. MUTUAL IMPEDANCE BETWEEN TWO IDENTICAL LINEAR DIPOLES
%
% BASED ON THE INDUCED EMF METHOD AND THE IDEAL CURRENT DISTRIBUTION
% OF EQUATION (4-56).
%
% I. SELF IMPEDANCE (INPUT IMPEDANCE)
%
% 1. BASED ON CURRENT AT THE INPUT
% Zin = Rin + jXin
% Rin = INPUT RESISTANCE [EQU. (8-60a) & (8-61a)]
% Xin = INPUT REACTANCE [EQU. (8-60b) & (8-61b)]
% 2. BASED ON CURRENT MAXIMUM
%
% Zinm = Rinm + jXinm
% Rinm = SELF RESISTANCE [EQU. (8-60a)]
% Xinm = SELF REACTANCE [EQU. (8-60b)]
%
% II. MUTUAL IMPEDANCE
%
% 1. BASED ON CURRENT AT THE INPUT
% Z21i = R21i + jX21i [EQU. (8-68)]
% R21i = INPUT MUTUAL RESISTANCE
% Xin = INPUT MUTUAL REACTANCE
% 2. BASED ON CURRENT MAXIMUM
% Z21m = R21m + jX21m [EQU. (8-70)]
% R21m = SELF MUTUAL RESISTANCE
% X21m = SELF MUTUAL REACTANCE
%
% THE DIPOLES FOR MUTUAL IMPEDANCE COMPUTATIONS MUST BE IDENTICAL
% WITH LENGTH OF ODD MULTIPLES OF HALF WAVELENGTH.
%
% OPTION I. SELF IMPEDANCE (INPUT IMPEDANCE)
%
% ** ARRAY INPUT PARAMETERS
%
% 1. LENGTH OF DIPOLE (IN WAVELENGTHS)
% 2. RADIUS OF DIPOLE (IN WAVELENGTHS)
%
% OPTION II. MUTUAL IMPEDANCE
%
% CHOICE A: SIDE-BY-SIDE [Fig. 8.19(a)]
%
% ** ARRAY INPUT PARAMETERS:
%
% 1. LENGTH OF THE DIPOLES (IN WAVELENGTHS)
% 2. HORIZONTAL DISPLACEMENT OF DIPOLES (IN WAVELENGTHS)
%
% CHOICE B: COLLINEAR [Fig. 8.19(b)]
%
% ** ARRAY INPUT PARAMETERS:
%
% 1. LENGTH OF THE DIPOLES (IN WAVELENGTHS)
% 2. VERTICAL DISPLACEMENT OF DIPOLES (IN WAVELENGTHS)
%
% CHOICE C. PARALLEL-IN-ECHELON [Fig. 8.19©]
%
% ** ARRAY INPUT PARAMETERS:
%
% 1. LENGTH OF THE DIPOLES (IN WAVELENGTHS)
% 2. VERTICAL DISPLACEMENT OF DIPOLES (IN WAVELENGTHS)
% 3. HORIZONTAL DISPLACEMENT OF DIPOLES (IN WAVELENGTHS)
%
%** NOTE: ALL THE INPUT PARAMETERS ARE IN WAVELENGTHS.
%******************************************************************
% Written by: Panayiotis Anastasiou, Arizona State University
%******************************************************************
clear all;
close all;
k = 2pi;
ETA = 120pi;
C=0.5772156649;

% Choose Output Device —> Screen or Output File

device_option = sprintf(‘OUTPUT DEVICE OPTION\n’);
device_option = sprintf(’%s OPTION (1): SCREEN\n’,device_option);
device_option = sprintf(’%s OPTION (2): OUTPUT FILE\n\n’,device_option);
device_option = sprintf(’%sOUTPUT DEVICE = ',device_option);

output_device = input(device_option);
if output_device == 2
filename = sprintf(’\nINPUT THE DESIRED OUTPUT FILENAME (in single quotes) = ');
filename = input(filename);
end

impedance_option = sprintf('\nSELF OR MUTUAL IMPEDANCE OPTION\n');
impedance_option = sprintf('%s   OPTION (1): SELF IMPEDANCE\n',impedance_option);
impedance_option = sprintf('%s   OPTION (2): MUTUAL IMPEDANCE\n\n',impedance_