Andrew Wade awade

## randomUnitVectors.m
%% First a method using quaternions
qtn = randrot(1000,1);

pts = rotatepoint(qtn, [1 0 0]);

figure()
scatter3(pts(:,1), pts(:,2), pts(:,3))
axis equal


## RG58_loss.m
function dBLossPerM = RG58_loss(f)
% Function for estimate RG58 coaxial cable loss as a function of frequency
% per meter.  Data is drawn from https://www.w4rp.com/ref/coax.html and is
% not varified with real lab measurments. Take it with a
% grain of salt.
%
% Output interpolates from the spaced data points.  This is for estimating
% only.
%
% Author: Andrew Wade

## rms2.m
function [x_clipped, rms] = rms2(x, y, varargin)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%	find the rms value of y(x)
%	y(f) could be a spectral density
%	or y(t) could be a time series
%
%   optional var sets the cutoff of x to compute
%   from.  This will shorten the output vector

## cask-upgrade.sh
#!/bin/bash

# USE AT YOUR OWN RISK.

# This is a bash script for checking brew cask versions and reinstalling if
# there is a version mismatch. Obviously this may revert casks that have self
# updating to previous version, but then it is as easy as simply updating the
# repo with the cask-repair command. This way the community gets the updates
# too.

## PDH_shotNoise.py
def S_f_SN(Pcar, Finesse, Lcav, lambd=1064.0e-9):
    '''
    Compute PDH shot noise
    '''

    _h = 6.62607004e-34  # [J.s] Plank's constant
    _c = 2.99792458e8  # [m/s] Speed of light

    return np.sqrt(_h * _c**3) / 8 / (Finesse * Lcav * np.sqrt(lambd * Pcar))

## ZurichPythonAccess.ipynb

      
              1 file
            
          
              0 forks
            
          
              0 comments
            
          
              0 stars
            
          
                awade
                / ZurichPythonAccess.ipynb
            
            
              Created
              May 17, 2019 01:28
            
              
                Minimum working example of how to access Zurich HF2LI instrument via python API
              
          
        Loading

      Sorry, something went wrong. Reload?
      Sorry, we cannot display this file.
      Sorry, this file is invalid so it cannot be displayed.
      
          Viewer requires iframe.
      
    
## numpyfloatRandRound.py
def truncFloat(x,dp):
    '''Truncates a numpy float to given decimal places (dp)'''
    return np.trunc(x * 10 ** dp) / 10 ** dp

def probRound(x, dp):
    '''Rounds float to given decimal places (dp), with
       last digit round done on a probablistic basis.'''

    xTrunc = truncFloat(x, dp)  # Truncate to dp figures
    xrem = x - xTrunc  # find remainder of truncation

## rms.py
import numpy as np

'''
Find the rms value of y(x). Here y(f) could be a spectral density or y(t) could be a time series

usage:
    rms_of_y = rms(X,Y);
'''


## HostUp.py
def HostUp(hostname, waittime=1000):
    '''Function returns True if host IP returns a ping, else False'''
    assert isinstance(hostname, str), \
        "IP/hostname must be provided as a string."
    if os.system("ping -c 1 -W " + str(waittime) + " " +
                 hostname + " > /dev/null 2>&1") is 0:
        HOST_UP = True
    else:
        HOST_UP = False
    return HOST_UP

## roombaroo.sh
#!/bin/bash
#
# Cleans a mac computer directory when called from crontab
#
#Author: Andrew Wade
#Date: 20170125
#Last modified: yyyymmdd

TIMESTAMP=$(date +"%Y%m%d-%H%M%S")
	%% First a method using quaternions
	qtn = randrot(1000,1);

	pts = rotatepoint(qtn, [1 0 0]);

	figure()
	scatter3(pts(:,1), pts(:,2), pts(:,3))
	axis equal
	function dBLossPerM = RG58_loss(f)
	% Function for estimate RG58 coaxial cable loss as a function of frequency
	% per meter. Data is drawn from https://www.w4rp.com/ref/coax.html and is
	% not varified with real lab measurments. Take it with a
	% grain of salt.
	%
	% Output interpolates from the spaced data points. This is for estimating
	% only.
	%
	% Author: Andrew Wade
	function [x_clipped, rms] = rms2(x, y, varargin)

	%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	%
	% find the rms value of y(x)
	% y(f) could be a spectral density
	% or y(t) could be a time series
	%
	% optional var sets the cutoff of x to compute
	% from. This will shorten the output vector
	#!/bin/bash

	# USE AT YOUR OWN RISK.

	# This is a bash script for checking brew cask versions and reinstalling if
	# there is a version mismatch. Obviously this may revert casks that have self
	# updating to previous version, but then it is as easy as simply updating the
	# repo with the cask-repair command. This way the community gets the updates
	# too.
	def S_f_SN(Pcar, Finesse, Lcav, lambd=1064.0e-9):
	'''
	Compute PDH shot noise
	'''

	_h = 6.62607004e-34 # [J.s] Plank's constant
	_c = 2.99792458e8 # [m/s] Speed of light

	return np.sqrt(_h * _c*3) / 8 / (Finesse Lcav * np.sqrt(lambd * Pcar))
	def truncFloat(x,dp):
	'''Truncates a numpy float to given decimal places (dp)'''
	return np.trunc(x * 10 dp) / 10 dp

	def probRound(x, dp):
	'''Rounds float to given decimal places (dp), with
	last digit round done on a probablistic basis.'''

	xTrunc = truncFloat(x, dp) # Truncate to dp figures
	xrem = x - xTrunc # find remainder of truncation
	import numpy as np

	'''
	Find the rms value of y(x). Here y(f) could be a spectral density or y(t) could be a time series

	usage:
	rms_of_y = rms(X,Y);
	'''
	def HostUp(hostname, waittime=1000):
	'''Function returns True if host IP returns a ping, else False'''
	assert isinstance(hostname, str), \
	"IP/hostname must be provided as a string."
	if os.system("ping -c 1 -W " + str(waittime) + " " +
	hostname + " > /dev/null 2>&1") is 0:
	HOST_UP = True
	else:
	HOST_UP = False
	return HOST_UP
	#!/bin/bash
	#
	# Cleans a mac computer directory when called from crontab
	#
	#Author: Andrew Wade
	#Date: 20170125
	#Last modified: yyyymmdd

	TIMESTAMP=$(date +"%Y%m%d-%H%M%S")