dipole.zip_HFSS API_hfss matlab_hfss matlab api_hfss matlab dip

% 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; false = 0; true = 1; % add paths to the required m-files. addpath('...\hfssapi/3dmodeler/'); addpath('...\hfssapi/general/'); addpath('...\hfssapi/analysis/'); addpath('...\hfssapi/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 = 'D:\yuzhu\tmpDipole.hfss'; tmpDataFile = 'D:\yuzhu\tmpData.m'; tmpScriptFile = 'D:\yuzhu\dipole_example.vbs'; % HFSS Executable Path. hfssExePath = 'D:\"Program Files"\Ansoft\HFSS10\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/');