~ubuntu-branches/ubuntu/precise/yagiuda/precise

« back to all changes in this revision

Viewing changes to src/max_gain.c

Committer: Bazaar Package Importer
Author(s): Joop Stakenborg
Date: 2005-08-22 20:20:13 UTC
Revision ID: james.westby@ubuntu.com-20050822202013-mhhxp4xirdxrdfx1

Tags: upstream-1.19

Import upstream version 1.19

files added:

AUTHORS

BUGS

COPYING

ChangeLog

INSTALL

Makefile.am

Makefile.in

NEWS

README

TODO

aclocal.m4

config.guess

config.sub

configure

configure.ac

depcomp

doc/Makefile.am

doc/Makefile.in

doc/dipole.1.html

doc/first.1.html

doc/index.html

doc/input.1.html

doc/optimise.1.html

doc/output.1.html

doc/pattern.jpg

doc/yagi.1.html

doc/yagi.jpg

install-sh

man/Makefile.am

man/Makefile.in

man/man1

man/man1/Makefile.am

man/man1/Makefile.in

man/man1/dipole.1

man/man1/first.1

man/man1/input.1

man/man1/optimise.1

man/man1/output.1

man/man1/yagi.1

man/man5

man/man5/Makefile.am

man/man5/Makefile.in

man/man5/first.5

man/man5/input.5

man/man5/optimise.5

man/man5/output.5

man/man5/yagi.5

missing

mkinstalldirs

src/3-element-yagi

src/Makefile.am

src/Makefile.in

src/auto.c

src/better.c

src/change_max_percentage_changes.c

src/check_flags.c

src/ci.c

src/cin.c

src/cis_hack.c

src/com_hack.c

src/com_hack.h

src/copy_matrix.c

src/dipole.c

src/display_antenna_currents.c

src/dobetter.c

src/dynamic.c

src/end_stop.c

src/error_message.c

src/fill_v_vector.c

src/fill_z_matrix.c

src/first.c

src/fitness.c

src/gain.c

src/gaussian.c

src/genetic.c

src/genetic_algorithm_lib.c

src/get_command_line_options.c

src/get_number_of_elements.c

src/getfiles.c

src/getopt.c

src/globals.h

src/input.c

src/install-sh

src/linear.c

src/lub_hack.c

src/lud_hack.c

src/max_gain.c

src/max_side.c

src/missing

src/mkinstalldirs

src/mpr_hack.c

src/mutual.c

src/new_length.c

src/nr_hack.h

src/nrutil.c

src/nrutil.h

src/optimise.c

src/optimising_for.c

src/output.c

src/perform.c

src/power_input.c

src/print_relavent_performance_data.c

src/print_z_matrix.c

src/random.c

src/randomise.c

src/rantest.c

src/read_header.c

src/read_yagi_data.c

src/reflection_coefficient.c

src/self.c

src/self2.c

src/selftest.c

src/sens.c

src/set_mean_structure.c

src/set_performance_structures.c

src/show_all_optimise_parameters.c

src/solve.c

src/string.c

src/subtract_structures.c

src/test2.c

src/test_stop.c

src/usage_first.c

src/usage_input.c

src/usage_optimise.c

src/usage_output.c

src/usage_yagi.c

src/version.c

src/vswr.c

src/write_coordinates.c

src/write_gain.c

src/write_gnuplot_header.c

src/write_header_to_disk.c

src/write_input_data.c

src/yagi.c

src/yagi.h

src/z_input.c

src/zbr_hack.c

Show diffs side-by-side

added added

removed removed

src/max_gain.c

#ifdef HAVE_STDLIB_H

#include <stdlib.h>

#endif

#include <stdio.h>

#include <math.h>

#include <errno.h>

#include "yagi.h"

extern int errno;

/* The maximum gain of a Yagi is given in the ARRL handbook, 1991, pp33-17

as 12dB at one wavelength boom length, rising at 2.6 dB

each time boom lengtrh is doubled */

double determine_maximum_gain(double f, double boom_length_in_m)

{

double lambda, wavelengths, max_gain;

lambda=3e8/f;

wavelengths=boom_length_in_m/lambda; /* boom length in wavelengths */

if(wavelengths < 1)

max_gain=8.0;

else

max_gain = 12.0 + 2.6*log2(wavelengths);

#ifdef DEBUG

if(errno)

{

fprintf(stderr,"Errno =%d in determine_maximum_gain() of max_gain.c\n", errno);

exit(1);

}

#endif

return(max_gain);

}

/* This is another algorithm. Its not as accurate as above, but it gives the

same result form given number of elements. This allows more accurate comparison

of fitnesss'es for the GA and -W algorithm, which would otherwise not be possible. */

double determine_maximum_gain2(int elements)

{

double wavelengths, max_gain;

wavelengths=(elements-1)*0.25; /* boom length in wavelengths */

max_gain = 12.0 + 2.6*log2(wavelengths);

#ifdef DEBUG

if(errno)

{

fprintf(stderr,"Errno =%d in determine_maximum_gain2() of max_gain.c\n", errno);

exit(1);

}

#endif

return(max_gain);

}

double log2(double x)

{

return(log(x)/0.6931471);

}

Older »