FrFT DPW processing validade in terms of dimensions and code generation

codegen_fracFdpw/TBc_fracFdpw.m

@@ -0,0 +1,20 @@
+a = single(1);
+
+[Achirp,AchirpOut,H,scale] = fracF_init(a);
+
+X = complex(randn(1024,1024,'single'), ...
+            randn(1024,1024,'single'));
+
+Fdpw = fracF_dpw(X,...
+                 Achirp,...
+                 AchirpOut,...
+                 H,...
+                 scale);
+
+Fref = complex(zeros(512,64,'single'));
+
+for k = 1:1024
+    Fref(:,k) = fracF_cg(X(:,k),a);
+end
+
+maxErr = max(abs(Fdpw(:)-Fref(:)))

codegen_fracFdpw/fracF_dpw.m

@@ -0,0 +1,51 @@
+function F = fracF_dpw(f,...
+                       Achirp,...
+                       AchirpOut,...
+                       H,...
+                       scale)
+%#codegen
+%
+% fracF_dpw
+%
+% Matrix FrFT processing for an entire DPW.
+%
+% INPUT:
+%   f          [1024 x Nframes] complex(single)
+%   Achirp     [1024 x 1]
+%   AchirpOut  [512  x 1]
+%   H          [2048 x 1]
+%   scale      scalar
+%
+% OUTPUT:
+%   F          [512 x Nframes]
+
+N    = 1024;
+Nfft = 2048;
+
+Nframes = size(f,2);
+
+%% First chirp multiplication
+
+g = f .* Achirp;
+
+%% Zero-padding
+
+g_pad = complex(zeros(Nfft,Nframes,'single'));
+
+g_pad(1:N,:) = g;
+
+%% FFT convolution
+
+Gfft = fft(g_pad);
+
+G = ifft(Gfft .* H);
+
+%% Extract valid region and decimate
+
+G_valid = G(N+1:2:end,:);
+
+%% Final chirp multiplication
+
+F = scale .* G_valid .* AchirpOut;
+
+end

codegen_fracFdpw/fracF_init.m

@@ -0,0 +1,70 @@
+function [Achirp,AchirpOut,H,scale] = fracF_init(a)
+%#codegen
+%
+% fracF_init
+%
+% Precompute FrFT coefficients for DPW processing.
+%
+% INPUT:
+%   a  - FrFT order (single)
+%
+% OUTPUTS:
+%   Achirp     [1024 x 1]
+%   AchirpOut  [512  x 1]
+%   H          [2048 x 1]
+%   scale      scalar
+
+N    = 1024;
+
+pi_s = single(pi);
+
+phi = a * (pi_s/2);
+
+tan_half_phi = tan(phi/2);
+sin_phi      = sin(phi);
+cos_phi      = cos(phi);
+
+csc_phi = 1/sin_phi;
+cot_phi = cos_phi/sin_phi;
+
+twoDelta = 2*sqrt(single(N)/2);
+
+%% Chirp A
+
+n = single((-N/2:N/2-1).') / twoDelta;
+
+Achirp = exp(-1j*pi_s*(n.^2)*tan_half_phi);
+
+%% Chirp B
+
+m = single((-N:N-1).') / twoDelta;
+
+Bchirp = exp(1j*pi_s*csc_phi*(m.^2));
+
+%% FFT of Chirp B
+
+H = fft(Bchirp);
+
+%% Output chirp
+
+AchirpOut = Achirp(1:2:end);
+
+%% Scale factor
+
+scale = sqrt(1 - 1j*cot_phi) / twoDelta;
+
+%% Force single precision complex
+
+Achirp    = complex(single(real(Achirp)), ...
+                    single(imag(Achirp)));
+
+AchirpOut = complex(single(real(AchirpOut)), ...
+                    single(imag(AchirpOut)));
+
+H         = complex(single(real(H)), ...
+                    single(imag(H)));
+
+scale     = complex(single(real(scale)), ...
+                    single(imag(scale)));
+
+end

resources/project/EEtUlUb-dLAdf0KpMVivaUlztwA/0W-QdLjPl5JXWY5WJJ8XuZdPRxsd.xml

@@ -0,0 +1,2 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<Info Ref="codegen_fracFdpw" Type="Relative"/>

resources/project/EEtUlUb-dLAdf0KpMVivaUlztwA/0W-QdLjPl5JXWY5WJJ8XuZdPRxsp.xml

@@ -0,0 +1,2 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<Info location="1617913e-1a81-4c05-b67f-42a35e5cd27b" type="Reference"/>

resources/project/EEtUlUb-dLAdf0KpMVivaUlztwA/IYUg9GI8QpDmyAvmDmvwktdmhjEd.xml

@@ -1,2 +0,0 @@
-<?xml version="1.0" encoding="UTF-8"?>
-<Info Ref="codegen_frft/codegen/mex/fracF_cg/interface" Type="Relative"/>

resources/project/EEtUlUb-dLAdf0KpMVivaUlztwA/IYUg9GI8QpDmyAvmDmvwktdmhjEp.xml

@@ -1,2 +0,0 @@
-<?xml version="1.0" encoding="UTF-8"?>
-<Info location="4a7e46c8-b608-4097-9e1f-6b08f40b724c" type="Reference"/>

resources/project/EEtUlUb-dLAdf0KpMVivaUlztwA/N9n74IL2KrtysYy1MJ3GXg0B9cwd.xml

@@ -1,2 +0,0 @@
-<?xml version="1.0" encoding="UTF-8"?>
-<Info Ref="codegen_frft/codegen" Type="Relative"/>

resources/project/EEtUlUb-dLAdf0KpMVivaUlztwA/N9n74IL2KrtysYy1MJ3GXg0B9cwp.xml

@@ -1,2 +0,0 @@
-<?xml version="1.0" encoding="UTF-8"?>
-<Info location="7593a6a4-6687-4364-8723-1bc8794d0177" type="Reference"/>

resources/project/EEtUlUb-dLAdf0KpMVivaUlztwA/SmoDhotSC9DRW6RAWIZ3HhfJ1Gkd.xml

@@ -1,2 +0,0 @@
-<?xml version="1.0" encoding="UTF-8"?>
-<Info Ref="codegen_frft/codegen/mex" Type="Relative"/>

resources/project/EEtUlUb-dLAdf0KpMVivaUlztwA/SmoDhotSC9DRW6RAWIZ3HhfJ1Gkp.xml

@@ -1,2 +0,0 @@
-<?xml version="1.0" encoding="UTF-8"?>
-<Info location="732585ee-100e-439f-988e-c97de5008a4c" type="Reference"/>

resources/project/EEtUlUb-dLAdf0KpMVivaUlztwA/XNALo81LXKoPpl_sWox0XlMZwrgp.xml

@@ -1,2 +0,0 @@
-<?xml version="1.0" encoding="UTF-8"?>
-<Info location="a92ac691-f104-48a3-8517-c0b00eec410f" type="Reference"/>

resources/project/EEtUlUb-dLAdf0KpMVivaUlztwA/yHR9AGbyWFGhYHIDxoJYqebuMa8p.xml

@@ -0,0 +1,2 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<Info location="c2875307-deb1-4e40-b64d-d77c1eb908cb" type="Reference"/>

resources/project/EEtUlUb-dLAdf0KpMVivaUlztwA/yPxdZXCe9dRoqAjp-aSxSvfCy5Ed.xml

@@ -1,2 +0,0 @@
-<?xml version="1.0" encoding="UTF-8"?>
-<Info Ref="codegen_frft/codegen/mex/fracF_cg" Type="Relative"/>

resources/project/EEtUlUb-dLAdf0KpMVivaUlztwA/yPxdZXCe9dRoqAjp-aSxSvfCy5Ep.xml

@@ -1,2 +0,0 @@
-<?xml version="1.0" encoding="UTF-8"?>
-<Info location="81080468-b05c-494e-b13d-b9983b9ef350" type="Reference"/>

resources/project/SDMLmaWG7GHkKoV_kJ80kV-luys/1X5PSVWg8jp3h1SyNiJtJmD4jnQd.xml

@@ -0,0 +1,6 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<Info>
+    <Category UUID="FileClassCategory">
+        <Label UUID="design"/>
+    </Category>
+</Info>

resources/project/SDMLmaWG7GHkKoV_kJ80kV-luys/1X5PSVWg8jp3h1SyNiJtJmD4jnQp.xml

@@ -0,0 +1,2 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<Info location="fracF_dpw.m" type="File"/>

resources/project/SDMLmaWG7GHkKoV_kJ80kV-luys/GIFBHwpY7Mj_mw96S-Wz-4P1djYd.xml

@@ -0,0 +1,6 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<Info>
+    <Category UUID="FileClassCategory">
+        <Label UUID="design"/>
+    </Category>
+</Info>

resources/project/SDMLmaWG7GHkKoV_kJ80kV-luys/GIFBHwpY7Mj_mw96S-Wz-4P1djYp.xml

@@ -0,0 +1,2 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<Info location="TBm_fracFdpw.slx" type="File"/>

resources/project/SDMLmaWG7GHkKoV_kJ80kV-luys/VdNularF-Bizb3hmjaa7F00asGMd.xml

@@ -0,0 +1,2 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<Info/>

resources/project/SDMLmaWG7GHkKoV_kJ80kV-luys/VdNularF-Bizb3hmjaa7F00asGMp.xml

@@ -0,0 +1,2 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<Info location="1" type="DIR_SIGNIFIER"/>

resources/project/SDMLmaWG7GHkKoV_kJ80kV-luys/_GBXjYOPGNzpaGXZHNU9eBKs2y4d.xml

@@ -0,0 +1,6 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<Info>
+    <Category UUID="FileClassCategory">
+        <Label UUID="design"/>
+    </Category>
+</Info>

resources/project/SDMLmaWG7GHkKoV_kJ80kV-luys/_GBXjYOPGNzpaGXZHNU9eBKs2y4p.xml

@@ -0,0 +1,2 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<Info location="fracF_init.m" type="File"/>

resources/project/SDMLmaWG7GHkKoV_kJ80kV-luys/rFgbE7J0q1IfyjFFVSXP14Sbh30d.xml

@@ -0,0 +1,6 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<Info>
+    <Category UUID="FileClassCategory">
+        <Label UUID="design"/>
+    </Category>
+</Info>

resources/project/SDMLmaWG7GHkKoV_kJ80kV-luys/rFgbE7J0q1IfyjFFVSXP14Sbh30p.xml

@@ -0,0 +1,2 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<Info location="TBc_fracFdpw.m" type="File"/>

resources/project/qaw0eS1zuuY1ar9TdPn1GMfrjbQ/SDMLmaWG7GHkKoV_kJ80kV-luysd.xml

@@ -0,0 +1,2 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<Info/>

resources/project/qaw0eS1zuuY1ar9TdPn1GMfrjbQ/SDMLmaWG7GHkKoV_kJ80kV-luysp.xml

@@ -0,0 +1,2 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<Info location="codegen_fracFdpw" type="File"/>

utilities/soc_rfsoc_init.m

@@ -25,7 +25,7 @@ NCOCountIncDT = numerictype(1,NCOAccumWL*2,NCOAccumWL);
 
 % Pulse start/end frequencies
 pulseCentFreq = 0e6;
-pulseBw = 40e6; % Pulse bandwidth
+pulseBw = 32e6; % Pulse bandwidth
 
 % Number of pulses
 numPulses = 10;
@@ -37,7 +37,7 @@ PRI = 1/PRF;
 
 % Pulse time duration
 %pulseT = 10; % use very long pulse help emulate CW
-pulseT = 20e-6;
+pulseT = 32e-6;
 
 % CW mode (bypass pulse generation)
 %CwMode = false;
@@ -83,8 +83,17 @@ channelizerCoeffs = channelizer.coeffs.Numerator;
 %Starting frequency for each channel
 %chanFStart = chanBW/2:chanBW:(fs/2-chanBW/2);
 
-%Number of frames in the DPW
+%% Frame and DPW capture parameters 
-nFrames = 1024;%nChan/SamplesPerCycle;
+samplesFrame = 512; % number of samples in one frame
+
+TFrame = samplesFrame*Ts_eff; % time duration of a frame
+
+%Number of frames in the DPW
+nFrames = 1024;
+
+%% FrFT tests
+pulseBeta = pulseBw/pulseT; % chirp-rate in Hz/s
+
+aMatch = -(2/pi)*(atan(fs_eff/(pulseBeta*TFrame)));
+%aMatch = -(2/pi)*(atan(fs_eff/(pulseBw*pulseT)));
 
-% Frame size after serializing x2
-%frameSize = SamplesPerCycle/2;

utilities/soc_rfsoc_prj_startup.m

@@ -1,8 +1,24 @@
 %% Startup Tasks
 %
-% Configure HDL Coder to use Xilinx Vivado for HDL workflows.
+% Configure HDL Coder to use Xilinx Vivado for HDL workflows
+% only when Vivado is installed on this machine.
 %
-hdlsetuptoolpath('ToolName','Xilinx Vivado', ...
-                 'ToolPath','/tools/Xilinx/Vivado/2024.1/bin/vivado');
 
-%% 
+vivadoPath = '/tools/Xilinx/Vivado/2024.1/bin/vivado';
+
+if exist(vivadoPath, 'file') == 2
+
+    hdlsetuptoolpath( ...
+        'ToolName', 'Xilinx Vivado', ...
+        'ToolPath', vivadoPath);
+
+    fprintf('[FrFT Project] Vivado detected. HDL workflow enabled.\n');
+
+else
+
+    warning([ ...
+        '[FrFT Project] Vivado not found.\n' ...
+        'HDL synthesis / bitstream generation disabled.\n' ...
+        'Project running in MATLAB/Simulink modeling-only mode.']);
+
+end