-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathdigitalNull.m
More file actions
13 lines (12 loc) · 845 Bytes
/
digitalNull.m
File metadata and controls
13 lines (12 loc) · 845 Bytes
1
2
3
4
5
6
7
8
9
10
11
12
13
function NulledZ = digitalNull(Needle, NoNeedle, AmpRef, stimFrequency, sampleRate)
%moves probe to zero DC and runs a long no signal
disp('nulling...')
AmpCpx = fftpoint_cpx3(AmpRef, sampleRate, stimFrequency);
NCpx = fftpoint_cpx3(Needle, sampleRate, stimFrequency);
NNCpx = fftpoint_cpx3(NoNeedle, sampleRate, stimFrequency);
outCpx = fftpoint_cpx3(data(:,1), BM.sampleRate,BM.stimFrequency);
inCpx = fftpoint_cpx3(data(:,2), BM.sampleRate,BM.stimFrequency);
digitalNulledSignal = NCpx + NNCpx; %they are out of phase from the probe thus adition negates
stackVelocity = -2i*pi* stimFrequency * AmpCpx(end); %a discrete frequency differentiation
NulledZ = digitalNulledSignal ./ stackVelocity;
disp('nulled!')