William Matthews WillMatthews

## party-cam.sh
#!/bin/bash

# Creates a fake webcam pointing at a party parrot gif.
# Requires you to have a a party parrot somewhere - just make sure you set the path correctly

PARROT=~/party-parrot.gif

# Reload v4l2loopback
sudo modprobe -r v4l2loopback

## ulid.m
% MATLAB port of https://github.com/alizain/ulid
% inspiration from https://github.com/mdipierro/ulid
% License MIT

function data = ulid(varargin)

if nargin == 0
    data =  generate();
elseif nargin == 1
    %data = string(1,varargin(1));

## sendemail.m
function [] = sendemail(emailMessage)
setpref('Internet','SMTP_Server','smtp.gmail.com');
setpref('Internet','E_mail','your-email@gmail.com');
setpref('Internet','SMTP_Username','smtp-username-here');
setpref('Internet','SMTP_Password','p@$$w0rd-here');
props = java.lang.System.getProperties;
props.setProperty('mail.smtp.auth','true');
props.setProperty('mail.smtp.socketFactory.class', 'javax.net.ssl.SSLSocketFactory');
props.setProperty('mail.smtp.socketFactory.port','465');

## fourier-space-webcam.py
import cv2
import numpy as np

def fft_frame(frame):
    # b, g, r
    im = [np.array([[col[i] for col in row] for row in frame]) for i in range(3)]
    # shift to u,v space
    F = [np.fft.fft2(img) for img in im]

    # amplitude

## vis_psf.m
clear all; close all; clc;

xp = linspace(-2,2,1000);
yp = xp';

d = (xp.^2 + yp.^2).^0.5;
d2 = (abs(xp) + abs(yp));

% defines pupil
p = (d<0.1); % circle

## sine_fit.m
clear all; close all; clc;

% the below code fits a sine wave which has a zero DC component.
% the output is 'bestx' which is a 1x3 matrix with components
% [amplitude, omega, phase 0-2*pi]

% we also have a fitting function @fitfun which can be used y =
% fitfun(bestx,time_arry) to generate our fitted sine wave for some generic
% time array.

## serial-test.py
import serial # ONLY USE PYSERIAL!! pip install pyserial
import time
import sys
import glob

def serial_ports():
    """ Lists serial port names

        :raises EnvironmentError:
            On unsupported or unknown platforms

## flips.py
import matplotlib.pyplot as plt

flips = [0 for _ in range(100)]
flips[0] = 1
for i in range(1, len(flips)):
    flips[i] = flips[i-1]/2

total_mass = sum(flips)

prob = [f/total_mass for f in flips]

## playground.m
clear all; close all; clc;

% 3B13c Q1 Tutorial Exploration Sandbox
% Written by William Matthews 25/02/2021
% Q2, Q3, Q4, Q6 can be applied here,
% but is left as an exercise for the interested reader.

% For the interested reader, Q5 asks for the instantaneous frequency.
% This can be found with a signal processing technique,
% research 'Analytic signal' which is covered in 4th year's

## multi-detector-dot-prod.py
import matplotlib.pyplot as plt
import numpy as np
import math as m

def Rx(theta):
  return np.matrix([[ 1, 0           , 0           ],
                   [ 0, m.cos(theta),-m.sin(theta)],
                   [ 0, m.sin(theta), m.cos(theta)]])

def Ry(theta):
	#!/bin/bash

	# Creates a fake webcam pointing at a party parrot gif.
	# Requires you to have a a party parrot somewhere - just make sure you set the path correctly

	PARROT=~/party-parrot.gif

	# Reload v4l2loopback
	sudo modprobe -r v4l2loopback
	% MATLAB port of https://github.com/alizain/ulid
	% inspiration from https://github.com/mdipierro/ulid
	% License MIT

	function data = ulid(varargin)

	if nargin == 0
	data = generate();
	elseif nargin == 1
	%data = string(1,varargin(1));
	function [] = sendemail(emailMessage)
	setpref('Internet','SMTP_Server','smtp.gmail.com');
	setpref('Internet','E_mail','your-email@gmail.com');
	setpref('Internet','SMTP_Username','smtp-username-here');
	setpref('Internet','SMTP_Password','p@$$w0rd-here');
	props = java.lang.System.getProperties;
	props.setProperty('mail.smtp.auth','true');
	props.setProperty('mail.smtp.socketFactory.class', 'javax.net.ssl.SSLSocketFactory');
	props.setProperty('mail.smtp.socketFactory.port','465');
	import cv2
	import numpy as np

	def fft_frame(frame):
	# b, g, r
	im = [np.array([[col[i] for col in row] for row in frame]) for i in range(3)]
	# shift to u,v space
	F = [np.fft.fft2(img) for img in im]

	# amplitude
	clear all; close all; clc;

	xp = linspace(-2,2,1000);
	yp = xp';

	d = (xp.^2 + yp.^2).^0.5;
	d2 = (abs(xp) + abs(yp));

	% defines pupil
	p = (d<0.1); % circle
	clear all; close all; clc;

	% the below code fits a sine wave which has a zero DC component.
	% the output is 'bestx' which is a 1x3 matrix with components
	% [amplitude, omega, phase 0-2*pi]

	% we also have a fitting function @fitfun which can be used y =
	% fitfun(bestx,time_arry) to generate our fitted sine wave for some generic
	% time array.
	import serial # ONLY USE PYSERIAL!! pip install pyserial
	import time
	import sys
	import glob

	def serial_ports():
	""" Lists serial port names

	:raises EnvironmentError:
	On unsupported or unknown platforms
	import matplotlib.pyplot as plt

	flips = [0 for _ in range(100)]
	flips[0] = 1
	for i in range(1, len(flips)):
	flips[i] = flips[i-1]/2

	total_mass = sum(flips)

	prob = [f/total_mass for f in flips]
	clear all; close all; clc;

	% 3B13c Q1 Tutorial Exploration Sandbox
	% Written by William Matthews 25/02/2021
	% Q2, Q3, Q4, Q6 can be applied here,
	% but is left as an exercise for the interested reader.

	% For the interested reader, Q5 asks for the instantaneous frequency.
	% This can be found with a signal processing technique,
	% research 'Analytic signal' which is covered in 4th year's
	import matplotlib.pyplot as plt
	import numpy as np
	import math as m

	def Rx(theta):
	return np.matrix([[ 1, 0 , 0 ],
	[ 0, m.cos(theta),-m.sin(theta)],
	[ 0, m.sin(theta), m.cos(theta)]])

	def Ry(theta):