dipole_example_no.m

对电磁波入射的源为偶极子进行仿真

% Description:
% ------------
% A simple example to demonstrate the HFSS-MATLAB-API.  This script optimizes
% a dipole antenna design to resonate at a specified frequency. The initial
% length of the dipole is taken to be half the wavelength and is optimized 
% so that the simulated resonance frequency and the desired resonance 
% frequency are close.

clear all;

% add paths to the required m-files.
addpath('../3dmodeler/');
addpath('../general/');
addpath('../analysis/');
addpath('../boundary/');

% Antenna Parameters.
fC = 150e6;		% Frequency of Interest.
Wv = 3e8/fC;	% Wavelength.
L = Wv/2; 		% Antenna Length.
gapL = 5e-2;	% Antenna Gap.
aRad = 2e-2;    % Antenna Radius.

% Simulation Parameters.
fLow = 100e6;
fHigh = 200e6;
nPoints = 201;

% AirBox Parameters.
AirX = Wv/2 + L;	% Include the antenna length.
AirY = Wv/2; 
AirZ = Wv/2;

% Temporary Files. These files can be deleted after the optimization
% is complete. We have to specify the complete path for all of them.
tmpPrjFile = 'C:\temp\tmpDipole.hfss';
tmpDataFile = 'C:\temp\tmpData.m';
tmpScriptFile = 'C:\temp\dipole_example.vbs';

% HFSS Executable Path.
hfssExePath = 'C:\"Program Files"\Ansoft\HFSS9\hfss.exe';

% Plot Colors.
pltCols = ['b', 'r', 'k', 'g', 'm', 'c', 'y'];
nCols = length(pltCols);

% Optimization stop conditions.
maxIters = 15;		% max # of iterations.
Accuracy = 0.01;	% accuracy required (1%).
hasConverged = false;

disp(sprintf('The Initial Dipole Length is %.2f meter ...', L));
for iIters = 1:maxIters,
	disp(sprintf('Running iteration #%d ...', iIters));
	disp('Creating the Script File ...');
	
	% Create a new temporary HFSS script file.
	fid = fopen(tmpScriptFile, 'wt');

	% Create a new HFSS Project and insert a new design.
	hfssNewProject(fid);
	hfssInsertDesign(fid, 'without_balun');

	% Create the Dipole.
	hfssDipole(fid, 'Dipole', 'X', [0, 0, 0], L, 2*aRad, gapL, 'meter');

	% Assign PE boundary to the antenna elements.
	hfssAssignPE(fid, 'Antennas',  {'Dipole1', 'Dipole2'});

	% Create a Lumped Gap Source (a rectangle normal to the Y-axis)
	hfssRectangle(fid, 'GapSource', 'Y', [-gapL/2, 0, -aRad], 2*aRad, gapL, ...
    	          'meter');
	hfssAssignLumpedPort(fid, 'LumpedPort', 'GapSource', [-gapL/2, 0, 0], ...
    	                 [gapL/2, 0, 0], 'meter');

	% Add an AirBox.
	hfssBox(fid, 'AirBox', [-AirX, -AirY, -AirZ]/2, [AirX, AirY, AirZ], ...
	        'meter');
	hfssAssignRadiation(fid, 'ABC', 'AirBox');

	% Add a Solution Setup.
	hfssInsertSolution(fid, 'Setup150MHz', fC/1e9);
	hfssInterpolatingSweep(fid, 'Sweep100to200MHz', 'Setup150MHz', ...
	                       fLow/1e9, fHigh/1e9, nPoints);

	% Save the project to a temporary file and solve it.
	hfssSaveProject(fid, tmpPrjFile, true);
	hfssSolveSetup(fid, 'Setup150MHz');

	% Export the Network data as an m-file.
	hfssExportNetworkData(fid, tmpDataFile, 'Setup150MHz', 'Sweep100to200MHz');

	% Close the HFSS Script File.
	fclose(fid);

	% Execute the Script by starting HFSS.
	disp('Solving using HFSS ..');
	hfssExecuteScript(hfssExePath, tmpScriptFile);

	% Load the data by running the exported matlab file.
	run(tmpDataFile);

	% The data items are in the f, S, Z variables now. 
	% Plot the data.
	disp('Solution Completed. Plotting Results for this Iteration ...');
	figure(1);
	hold on;
	plot(f/1e6, 20*log10(abs(S)), pltCols(mod(iIters, nCols) + 1)); 
    hold on;
	xlabel('Frequency (MHz) ->');
	ylabel('S_{11} (dB) ->');
	axis([fLow/1e6, fHigh/1e6, -20, 0]);
	
	% Find the Resonance Frequency.
	[Smin, iMin] = min(S);
	fActual = f(iMin);
	disp(sprintf('Simulated Resonance Frequency: %.2f MHz', fActual/1e6));

	% Check if the required accuracy is met.
	if (abs((fC - fActual)/fC) < Accuracy)
		disp('Required Accuracy is met !!');
		disp(sprintf('Optimized Antenna Length is %.2f meter.', L));
		hasConverged = true;
		break;
	end;
	
	% Adjust the antenna length in accordance with the discrepancy between
	% the estimated and desired resonance frequency.
	L = L*fActual/fC;

	% Loop all over again ...
	disp('Required accuracy not yet met ...');
	disp(sprintf('The new estimate for the dipole length is %.2f meter', L));
end;

if (~hasConverged)
	disp('Max Iterations exceeded. Optimization did NOT converge ...');
end;

% remove all the added paths.
rmpath('../3dmodeler/');
rmpath('../general/');
rmpath('../analysis/');
rmpath('../boundary/');