openofdm/verilog/complex_to_mag.v
2017-04-07 11:36:21 -04:00

55 lines
1.0 KiB
Verilog

module complex_to_mag
#(
parameter DATA_WIDTH = 16
)
(
input clock,
input enable,
input reset,
input signed [DATA_WIDTH-1:0] i,
input signed [DATA_WIDTH-1:0] q,
input input_strobe,
output reg [DATA_WIDTH-1:0] mag,
output mag_stb
);
reg [DATA_WIDTH-1:0] abs_i;
reg [DATA_WIDTH-1:0] abs_q;
reg [DATA_WIDTH-1:0] max;
reg[ DATA_WIDTH-1:0] min;
delayT #(.DATA_WIDTH(1), .DELAY(3)) stb_delay_inst (
.clock(clock),
.reset(reset),
.data_in(input_strobe),
.data_out(mag_stb)
);
// http://dspguru.com/dsp/tricks/magnitude-estimator
// alpha = 1, beta = 1/4
// avg err 0.006
always @(posedge clock) begin
if (reset) begin
mag <= 0;
abs_i <= 0;
abs_q <= 0;
max <= 0;
min <= 0;
end else if (enable) begin
abs_i <= i[DATA_WIDTH-1]? (~i+1): i;
abs_q <= q[DATA_WIDTH-1]? (~q+1): q;
max <= abs_i > abs_q? abs_i: abs_q;
min <= abs_i > abs_q? abs_q: abs_i;
mag <= max + (min>>2);
end
end
endmodule