NV_Matlab_Code/Wave_generator_Chase.m at main · cliulinnaeus/NV_Matlab_Code · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
function A = Wave_generator_Chase (Waveform, Sample_rate, Length)

% Sample_rate in GHz
% Length in ns

% Converts from human-readable waveform, to AWG-readable waveform

%Inputs to Wave_generator must be written in the form
% For Analog
%[t1,t2,freq1,phase1,amp1;...
% t3,t4,freq2,phase2,amp2;...
%...]
%
% For Digital
%[t1,t2,value1;...
%t3,t4,value2;...
%...]

% units in form
% t in ns
% Sample_rate in GHz
% Freq in GHz

% In Chase AWG, the vertical resolution is 2^12 = 4096
% 2048 = 0V
% The sum of Amplitude at each point must be < 1
% This program will automatically convert Amplitude to 2^12=4096 vertical resolution


i = 0;
Total = Length*Sample_rate; % Convert to program each points
% Totoal must be an integer number of 16
Total = ceil(Total/16)*16;

% result = 0;
A = zeros(1,Total,'single');
% A = [];
[dimx,dimy] = size(Waveform);  % dimy == 3 Digital waveform, dimy == 5 Analog wave
if (dimy == 5)
    while (i < dimx)
        i = i+1;
        % How many points do we have
        if Waveform(i,5) > 1
            error("Amplitude is out of bound.")
        end
        B = round(Waveform(i,1)*Sample_rate):(round(Waveform(i,2)*Sample_rate)-1);
        C = ceil(2047+...
            Waveform(i,5)*2047*sin(2 * pi * Waveform(i,3) / Sample_rate * B + Waveform(i,4)));
        if (i < dimx)
            Delay = ceil(2047*ones(1,round((Waveform(i+1,1)-Waveform(i,2))*Sample_rate)));
        elseif (i == dimx)
            Delay = ceil(2047*ones(1,Total-round(Waveform(i,2)*Sample_rate)));
        end
        StartPoint = 1+round(Waveform(i,1)*Sample_rate);
        EndPoint = StartPoint+length(C)+length(Delay)-1;

        A(StartPoint:EndPoint) = [C Delay];
    end

elseif (dimy == 6) % Gaussian pulse
    while (i < dimx)
        i = i+1;
        % How many points do we have
        B = round(Waveform(i,1)*Sample_rate):(round(Waveform(i,2)*Sample_rate)-1);
        Mid = (Waveform(i,1)+Waveform(i,2))*Sample_rate/2;
        Sigma = (Waveform(i,2)-Waveform(i,1))*Sample_rate/(2*Waveform(i,6));
        C = ceil(2047+2047*Waveform(i,5)*exp(-(B-Mid).^2/(2*Sigma^2)).*...
            sin(2 * pi * Waveform(i,3) / Sample_rate * B + Waveform(i,4)));
        if (i < dimx)
            Delay = ceil(2047*ones(1,round((Waveform(i+1,1)-Waveform(i,2))*Sample_rate)));
        elseif (i == dimx)
            Delay = ceil(2047*ones(1,Total-round(Waveform(i,2)*Sample_rate)));
        end
        StartPoint = 1+round(Waveform(i,1)*Sample_rate);
        EndPoint = StartPoint+length(C)+length(Delay)-1;

        A(StartPoint:EndPoint) = [C Delay];
    end


elseif (dimy == 8) % Adding of 2 frequencies
    while (i < dimx)
        i = i+1;
        if Waveform(i,5) + Waveform(i,8) > 1
            error("Amplitude is out of bound.")
        end
        B = round(Waveform(i,1)*Sample_rate):(round(Waveform(i,2)*Sample_rate)-1);
        C = ceil(2047+...
            2047*Waveform(i,5)*sin(2 * pi * Waveform(i,3) / Sample_rate * B + Waveform(i,4)) +  2047*Waveform(i,8)*sin(2 * pi * Waveform(i,6) / Sample_rate * B + Waveform(i,7)));
        if (i < dimx)
            Delay = ceil(2047*ones(1,round((Waveform(i+1,1)-Waveform(i,2))*Sample_rate)));
        elseif (i == dimx)
            Delay = ceil(2047*ones(1,Total-round(Waveform(i,2)*Sample_rate)));
        end
        StartPoint = 1+round(Waveform(i,1)*Sample_rate);
        EndPoint = StartPoint+length(C)+length(Delay)-1;

        A(StartPoint:EndPoint) = [C Delay];
    end


elseif (dimy == 3) % dimy == 3 Digital waveform, dimy == 5,6,8,9 Analog wave

    while (i <dimx)
        i = i+1;
        C = ceil(2047 ...
            + 2047 * Waveform(i,3) * ones(1,round(Waveform(i,2)*Sample_rate) - round(Waveform(i,1)*Sample_rate)));
        if (i < dimx)
            Delay = ceil(2047*ones(1,round((Waveform(i+1,1)-Waveform(i,2))*Sample_rate)));
            % Delay =
            % zeros(1,round((Waveform(i+1,1)-Waveform(i,2))*Sample_rate));
            % Seem to be Wrong
        elseif (i == dimx)
            Delay = ceil(2047*ones(1,Total-round(Waveform(i,2)*Sample_rate)));
        end

        A = [A, C, Delay];
    end
end

A = single(A);


end