Homework 1: Word transmiter

Hi! fot this first homework I had to make a python program that generates words, the main idea is that the program determinates the percentage of success that a word is transmitted with correctly. The words are binary.

The parameters of the program are:

- Word length

- Frequemcy of zeros

- Percentage of correct zeros

- Percentage of correct ones 

- Iterations

The scrip of bash just take the values that are in the code and send a the script of python then take the results of python and save in the resultados.txt ready to plot with script gnuplot

The script of gnuplot just plot the values that we send

Here the scripts, there are in this order:

• awk
• gnuplot
• bash
• python  

	set term png
	set output 'algo.png'
	set xrange [0:9] #rangos
	set yrange [0.1:0.9]
	set zrange [0.1:0.9]
	set xlabel "Potencias"
	set ylabel "Frecuencia"
	set zlabel "Exitos"
	set pm3d implicit at s
	splot 'resultados.txt'

view raw algo.plot hosted with ❤ by GitHub

	import sys, math, random
	def Words(chain, freczero): #function to create the words
	word = []
	contzero = 0
	contone = 0
	limite0 = int (chain * freczero)
	limite1 = chain - limite0
	for i in range(chain):
	nums = random.randint (0,1)
	if (nums == 0):
	contzero = contzero + 1
	if (contzero > limite0):
	word.append = (1)
	else:
	word.append = (0)
	else:
	contone = contone + 1
	if (contone >= limite1):
	word.append(0)
	else:
	word.append(1)
	return word
	def Transmision(message, probzero, probone): # function to send the word and create a list with the message
	msg = []
	for i in message:
	x = random()
	if i == 0:
	if x < probzero:
	msg.append = (0)
	else:
	msg.append = (1)
	else:
	if x < probone:
	msg.append = (1)
	else:
	msg.append = (0)
	if message == msg:
	return 1
	else:
	return 0
	def main():
	chain = int(argv[1]) # lenght of the word
	probzero = float(argv[2]) # frequency to transmit a zero
	probone = float(argv[3]) # frequency to transmit a one
	freczero = float(argv[4])
	iteration = int(argv[5])
	wins = 0
	for i in range(iteration):
	message = Words(chain, freczer)
	restrans = Transmision(message, probzero, probone)
	if (restrans == 1):
	wins = wins + 1
	goods = float(wins)/float(iteration)
	print goods

view raw canal.py hosted with ❤ by GitHub

	BEGIN {
	count = 0;
	sumall = 0;
	}
	{
	sumall = sumall + $1;
	count = count + 1;
	}
	END {
	Resfinal = sumall/count;
	print Resfinall;
	}

view raw Obed.awk hosted with ❤ by GitHub

	#!/bin/bash
	# Las siguientes funciones son para checar
	# si existe el documento, para que no haya
	# problemas.
	dat=resultados.txt
	if [ -e $dat ]; then
	rm resultados.txt
	fi
	dat2=awk.txt
	if [ -e $dat2 ]; then
	rm awk.txt
	fi
	touch resultados.txt
	for chain in 20 15 10 5; do
	for probzero in 0.5 0.3 0.7 0.9; do
	for probone in 0.6 0.8 0.3 0.4; do
	for freczero in 0.8 0.9 0.6 0.3; do
	touch awk.txt
	for iteration in 1 2 3 6; do
	goods=`python canal.py $chain $probzero $probone $freczero $iteration`
	echo $goods >> awk.txt
	done
	resgoods=`awk -f Obed.awk awk.txt` #instruccion para tomar todas las impresion del archivo awk
	echo $chain $probzero $probone $freczero $resgoods >> resultados.txt #pasamos los datos
	rm awk.txt
	done
	done
	done
	done
	gnuplot algo.plot #creamos el plot
	eog algo.png

view raw script.sh hosted with ❤ by GitHub

Here a part of resultados.txt:

the order is:

- Word length

- Percentage of correct zeros

- Percentage of correct ones

- success

Here the result: this result is not that i want