Verilog 找出任意六个数中的最大、次最大和第三最大值以及三个数对应的序号
思路：1、把六个数看成是两组的3个数
2、对每组的三个数进行从大到小的排序
3、找出有序的两组数中最大、次最大和第三最大值
//模块1：对三个数按照从大到小排序
//例：若 i_t_1st_max=20 i_t_2nd_max=5 i_t_3th_max=30
i_Num_t_1st_max=1 i_Num_t_2nd_max=2 i_Num_t_3th_max=3
则 o_all_1st_max=30 o_all_2nd _max=20 o_all_3th_max=5
o_index_1stmax=3 o_index_2ndmax=1 o_index_3thmax=2
module rank3data(
input i_rst_n,
i_start
//输入的三个数
input signed[`EXTENDWORDSIZE:0] i_t_1st_max,
i_t_2nd_max,
i_t_3th_max,
//输入的三个数对应的序列号
input [`INDEX:0] i_Num_t_1st_max,
i_Num_t_2nd_max,
i_Num_t_3th_max,
//排序后输出的三个数对应的序号
output reg [`INDEX:0] o_index_1stmax,
o_index_2ndmax,
o_index_3thmax,
//排序后输出的三个数
output reg signed[`EXTENDWORDSIZE:0] o_all_1st_max,
o_all_2nd_max,
o_all_3th_max
);
always @(i_rst_n or i_start or i_t_1st_max or i_t_2nd_max or i_t_3th_max or i_Num_t_1st_max or i_Num_t_2nd_max or i_Num_t_3th_max) begin
if(!i_rst_n)begin
o_all_1st_max=0;
o_all_2nd_max=0;
o_all_3th_max=0;
o_index_1stmax=0;
o_index_2ndmax=0;
o_index_3thmax=0;
end
else if(i_start)begin
if(i_t_1st_max>=i_t_2nd_max)begin
if(i_t_1st_max>=i_t_3th_max) begin
if(i_t_2nd_max>=i_t_3th_max) begin
o_all_1st_max=i_t_1st_max;
o_all_2nd_max=i_t_2nd_max;
o_all_3th_max=i_t_3th_max;
o_index_1stmax=i_Num_t_1st_max;
o_index_2ndmax=i_Num_t_2nd_max;
o_index_3thmax=i_Num_t_3th_max;
end
else begin
o_all_1st_max=i_t_1st_max;
o_all_2nd_max=i_t_3th_max;
o_all_3th_max=i_t_2nd_max;
o_index_1stmax=i_Num_t_1st_max;
o_index_2ndmax=i_Num_t_3th_max;
o_index_3thmax=i_Num_t_2nd_max;
end
end
else begin
o_all_1st_max=i_t_3th_max;
o_all_2nd_max=i_t_1st_max;
o_all_3th_max=i_t_2nd_max;
o_index_1stmax=i_Num_t_3th_max;
o_index_2ndmax=i_Num_t_1st_max;
o_index_3thmax=i_Num_t_2nd_max;
end
end
else begin
if(i_t_1st_max
if(i_t_2nd_max>=i_t_3th_max) begin
o_all_1st_max=i_t_2nd_max;
o_all_2nd_max=i_t_3th_max;
o_all_3th_max=i_t_1st_max;
o_index_1stmax=i_Num_t_2nd_max;
o_index_2ndmax=i_Num_t_3th_max;
o_index_3thmax=i_Num_t_1st_max;
end
else begin
o_all_1st_max=i_t_3th_max;
o_all_2nd_max=i_t_2nd_max;
o_all_3th_max=i_t_1st_max;
o_index_1stmax=i_Num_t_3th_max;
o_index_2ndmax=i_Num_t_2nd_max;
o_index_3thmax=i_Num_t_1st_max;
end
end
else begin
o_all_1st_max=i_t_2nd_max;
o_all_2nd_max=i_t_1st_max;
o_all_3th_max=i_t_3th_max;
o_index_1stmax=i_Num_t_2nd_max;
o_index_2ndmax=i_Num_t_1st_max;
o_index_3thmax=i_Num_t_3th_max;
end
end
end
end
endmodule
//模块2：两组3个有序的数比较得到最大、次最大和第三最大值已经对应序号
//例：若i_t_1st_max=30 i_t_2nd_max=20 i_t_3th_max=5
i_b_1st_max=10 i_b_2nd_max=8 i_b_3th_max= -30
i_Num_t_1st_max=1 i_Num_t_2nd_max=2 i_Num_t_3th_max=3
i_Num_t_1st_max=4 i_Num_t_2nd_max=5 i_Num_t_3th_max=6
则o_all_1st_max=30 o_all_2nd _max=20 o_all_3th_max=10
o_index_1stmax=1 o_index_2ndmax=2 o_index_3thmax=4
module two3rank_compare(
input i_rst_n,
i_start,
input signed[`EXTENDWORDSIZE:0] i_t_1st_max,
i_t_2nd_max,
i_t_3th_max,
i_b_1st_max,
i_b_2nd_max,
i_b_3th_max,
input [`INDEX:0] i_Num_t_1st_max,
i_Num_t_2nd_max,
i_Num_t_3th_max,
i_Num_b_1st_max,
i_Num_b_2nd_max,
i_Num_b_3th_max,
output reg [`INDEX:0] o_index_1stmax,
o_index_2ndmax,
o_index_3thmax,
output reg signed[`EXTENDWORDSIZE:0] o_all_1st_max,
o_all_2nd_max,
o_all_3th_max
);
always @(i_rst_n or i_start or i_t_1st_max or i_t_2nd_max or i_t_3th_max or i_b_1st_max or i_b_2nd_max or i_b_3th_max or i_Num_t_1st_max or i_Num_t_2nd_max or i_Num_t_3th_max or i_Num_b_1st_max or i_Num_b_2nd_max or i_Num_b_3th_max)begin
if(!i_rst_n) begin
o_all_1st_max=0;
o_all_2nd_max=0;
o_all_3th_max=0;
o_index_1stmax=0;
o_index_2ndmax=0;
o_index_3thmax=0;
end
else if(i_start)begin
if(i_t_3th_max>=i_b_1st_max)begin //a b c
o_all_1st_max=i_t_1st_max;
o_all_2nd_max=i_t_2nd_max;
o_all_3th_max=i_t_3th_max;
o_index_1stmax=i_Num_t_1st_max;
o_index_2ndmax=i_Num_t_2nd_max;
o_index_3thmax=i_Num_t_3th_max;
end
else if(i_t_2nd_max>=i_b_1st_max)begin // a b d
o_all_1st_max=i_t_1st_max;
o_all_2nd_max=i_t_2nd_max;
o_all_3th_max=i_b_1st_max;
o_index_1stmax=i_Num_t_1st_max;
o_index_2ndmax=i_Num_t_2nd_max;
o_index_3thmax=i_Num_b_1st_max;
end
else if(i_t_1st_max>=i_b_1st_max)begin//a d b
if(i_t_2nd_max>=i_b_2nd_max)begin
o_all_1st_max=i_t_1st_max;
o_all_2nd_max=i_b_1st_max;
o_all_3th_max=i_t_2nd_max;
o_index_1stmax=i_Num_t_1st_max;
o_index_2ndmax=i_Num_b_1st_max;
o_index_3thmax=i_Num_t_2nd_max;
end
else begin//a d e
o_all_1st_max=i_t_1st_max;
o_all_2nd_max=i_b_1st_max;
o_all_3th_max=i_b_2nd_max;
o_index_1stmax=i_Num_t_1st_max;
o_index_2ndmax=i_Num_b_1st_max;
o_index_3thmax=i_Num_b_2nd_max;
end
end
else if(i_t_1st_max>=i_b_2nd_max)begin
if(i_t_2nd_max>=i_b_2nd_max)begin// d a b
o_all_1st_max=i_b_1st_max ;
o_all_2nd_max=i_t_1st_max;
o_all_3th_max=i_t_2nd_max;
o_index_1stmax=i_Num_b_1st_max;
o_index_2ndmax=i_Num_t_1st_max;
o_index_3thmax=i_Num_t_2nd_max;
end
else begin //d a e
o_all_1st_max=i_b_1st_max;
o_all_2nd_max=i_t_1st_max;
o_all_3th_max=i_b_2nd_max;
o_index_1stmax=i_Num_b_1st_max;
o_index_2ndmax=i_Num_t_1st_max;
o_index_3thmax=i_Num_b_2nd_max;
end
end
else if(i_t_1st_max>=i_b_3th_max)begin //d e a
o_all_1st_max=i_b_1st_max;
o_all_2nd_max=i_b_2nd_max;
o_all_3th_max=i_t_1st_max;
o_index_1stmax=i_Num_b_1st_max;
o_index_2ndmax=i_Num_b_2nd_max;
o_index_3thmax=i_Num_t_1st_max;
end
else begin // d e f
o_all_1st_max=i_b_1st_max;
o_all_2nd_max=i_b_2nd_max;
o_all_3th_max=i_b_3th_max;
o_index_1stmax=i_Num_b_1st_max;
o_index_2ndmax=i_Num_b_2nd_max;
o_index_3thmax=i_Num_b_3th_max;
end
end
end
endmodule
//顶层模块
module six_comparator(
input i_rst_n,
i_start,
input signed[`EXTENDWORDSIZE:0] i_t_1st_max,
i_t_2nd_max,
i_t_3th_max,
i_b_1st_max,
i_b_2nd_max,
i_b_3th_max,
input [`INDEX:0] i_Num_t_1st_max,
i_Num_t_2nd_max,
i_Num_t_3th_max,
i_Num_b_1st_max,
i_Num_b_2nd_max,
i_Num_b_3th_max,
output [`INDEX:0] o_index_1stmax,
o_index_2ndmax,
o_index_3thmax,
output signed[`EXTENDWORDSIZE:0] o_all_1st_max,
o_all_2nd_max,
o_all_3th_max
);
wire signed[`EXTENDWORDSIZE:0] tt_1st_max, tt_2nd_max,tt_3th_max, bb_1st_max, bb_2nd_max,bb_3th_max;
wire[`INDEX:0] Num_tt_1st_max, Num_tt_2nd_max, Num_tt_3th_max, Num_bb_1st_max, Num_bb_2nd_max,Num_bb_3th_max;
rank3data top(
.i_rst_n(i_rst_n),
.i_start(i_start),
.i_t_1st_max(i_t_1st_max),
.i_t_2nd_max(i_t_2nd_max),
.i_t_3th_max(i_t_3th_max),
.i_Num_t_1st_max(i_Num_t_1st_max),
.i_Num_t_2nd_max(i_Num_t_2nd_max),
.i_Num_t_3th_max(i_Num_t_3th_max),
.o_index_1stmax(Num_tt_1st_max),
.o_index_2ndmax(Num_tt_2nd_max),
.o_index_3thmax(Num_tt_3th_max),
.o_all_1st_max(tt_1st_max),
.o_all_2nd_max(tt_2nd_max),
.o_all_3th_max(tt_3th_max)
);
rank3data bottom(
.i_rst_n(i_rst_n),
.i_start(i_start),
.i_t_1st_max(i_b_1st_max),
.i_t_2nd_max(i_b_2nd_max),
.i_t_3th_max(i_b_3th_max),
.i_Num_t_1st_max(i_Num_b_1st_max),
.i_Num_t_2nd_max(i_Num_b_2nd_max),
.i_Num_t_3th_max(i_Num_b_3th_max),
.o_index_1stmax(Num_bb_1st_max),
.o_index_2ndmax(Num_bb_2nd_max),
.o_index_3thmax(Num_bb_3th_max),
.o_all_1st_max(bb_1st_max),
.o_all_2nd_max(bb_2nd_max),
.o_all_3th_max(bb_3th_max)
);
two3rank_compare findMAX(
.i_rst_n(i_rst_n),
.i_start(i_start),
.i_t_1st_max(tt_1st_max),
.i_t_2nd_max(tt_2nd_max),
.i_t_3th_max(tt_3th_max),
.i_b_1st_max(bb_1st_max),
.i_b_2nd_max(bb_2nd_max),
.i_b_3th_max(bb_3th_max),
.i_Num_t_1st_max(Num_tt_1st_max),
.i_Num_t_2nd_max(Num_tt_2nd_max),
.i_Num_t_3th_max(Num_tt_3th_max),
.i_Num_b_1st_max(Num_bb_1st_max),
.i_Num_b_2nd_max(Num_bb_2nd_max),
.i_Num_b_3th_max(Num_bb_3th_max),
.o_index_1stmax(o_index_1stmax),
.o_index_2ndmax(o_index_2ndmax),
.o_index_3thmax(o_index_3thmax),
.o_all_1st_max(o_all_1st_max),
.o_all_2nd_max(o_all_2nd_max),
.o_all_3th_max(o_all_3th_max)
);
endmodule
matlab code
function softMetricout2=Nord_Ronbin_Decoder_MST_43(pvecMetric_in,Integer_bit,Decimal_bit)
%==========================================================================================
% Project Code : DTMB Link Level Simulation
% Create Date : 2014/04/09
% Description : Nordstrom-Robinson decoder, used in the 4QAM-NR demapping of DTMB
% Version : V0.1
% Modification History:
%====================================================================
%%%------------------------parameters------------------------------%%%
format longG;
global S
pvecMetric11=pvecMetric_in(:,3)';
% Quantification :8bit wide = 1bit Symbol + 4bit Integer + 3bit Decimal
demapperout2=[];
for t=1:1:7488
demapperout2(t)=floor(pvecMetric11(1,t)*2^(Decimal_bit));
if(demapperout2(t)>(2^(Integer_bit+Decimal_bit)-1)) demapperout2(t)=2^(Integer_bit+Decimal_bit)-1;
else if(demapperout2(t)<=-(2^(Integer_bit+Decimal_bit))) demapperout2(t)=-(2^(Integer_bit+Decimal_bit));
end
end
end
pvecMetric=demapperout2;
y=zeros(1,8);
softMetric_out=zeros(1,3744);
NR_decoder_output=zeros(1,3744);
data_length=length(pvecMetric);
fp=fopen('input7488_611_mat.txt','w');
for i=1:7488
fprintf(fp,'%d ',pvecMetric(i));
end
fclose(fp);
%%%------------------------calculate the 256 states----------------------------------------%%%
% for i=1:1:256
% x=fliplr(de2bi(i-1,8));
% x7=x(1);
% x6=x(2);
% x5=x(3);
% x4=x(4);
% x3=x(5);
% x2=x(6);
% x1=x(7);
% x0=x(8);
% y(1) = mod(x7+x6+x0+x1+x3+(xor(x0,x4)&(mod(x1+x2+x3+x5,2)))+(xor(x1,x2)&(xor(x3,x5))),2);
% y(2) = mod(x7+x0+x1+x2+x4+(xor(x1,x5)&(mod(x2+x3+x4+x6,2)))+(xor(x2,x3)&(xor(x4,x6))),2);
% y(3) = mod(x7+x1+x2+x3+x5+(xor(x2,x6)&(mod(x3+x4+x5+x0,2)))+(xor(x3,x4)&(xor(x5,x0))),2);
% y(4) = mod(x7+x2+x3+x4+x6+(xor(x3,x0)&(mod(x4+x5+x6+x1,2)))+(xor(x4,x5)&(xor(x6,x1))),2);
% y(5) = mod(x7+x3+x4+x5+x0+(xor(x4,x1)&(mod(x5+x6+x0+x2,2)))+(xor(x5,x6)&(xor(x0,x2))),2);
% y(6) = mod(x7+x4+x5+x6+x1+(xor(x5,x2)&(mod(x6+x0+x1+x3,2)))+(xor(x6,x0)&(xor(x1,x3))),2);
% y(7) = mod(x7+x5+x6+x0+x2+(xor(x6,x3)&(mod(x0+x1+x2+x4,2)))+(xor(x0,x1)&(xor(x2,x4))),2);
% y(8) = mod((x0+x1+x2+x3+x4+x5+x6+x7+y(1)+y(2)+y(3)+y(4)+y(5)+y(6)+y(7)),2);
% %calculate 256 states 16bit
% for j=1:1:8
% C(i,j)=x(9-j);
% end
% for k=1:1:8
% C(i,8+k)=y(k);
% end
%
% %s=1-2c
% for m=1:1:16
% S(i,m)=1-2*C(i,m);
% end
% end
% load C;
% load S;
%%%-------------Find the maximum correlation rm, the sub-maximum correlation rsm----------------%%%
for i=0:16:(data_length-1)
% calculate tne rc
rc=zeros(1,256);
for state_num=1:1:256
for j=1:1:16
% if ((S(state_num,j)==-1)&(pvecMetric(i+j)==-128))
% rc(state_num)=rc(state_num)+ 127;
% else
% rc(state_num)=rc(state_num)+ S(state_num,j)*pvecMetric(i+j);
% en
if (S(state_num,j)==-1)
if(pvecMetric(i+j)== -128)
rc(state_num)=rc(state_num)+ 127;
else
rc(state_num)=rc(state_num)-pvecMetric(i+j);
end
else
rc(state_num)=rc(state_num)+pvecMetric(i+j);
end
end
end
% pvecMetric(1,1:16)
% save rc;
%
% fp=fopen('rc_256state.txt','w');
% for i=1:256
% fprintf(fp,'%d \n',rc(i));
% end
% fclose(fp);
%
max_rc=min(rc);%
smax_rc=min(rc);
tmax_rc=min(rc);
C_max=zeros(1,16);
C_smax=zeros(1,16);
C_tmax=zeros(1,16);
max_pos=0;
smax_pos=0;
%Find the maximum correlation max_rc and its corresponding codewords
for t=1:1:256
if(rc(t)>max_rc)
max_rc=rc(t);
C_max(1,:)=S(t,:);
max_pos=t;
end
end
%Find the sub-maximum correlation smax_rc and its corresponding codewords
for l=1:1:256
if l~=max_pos %skip the maximum piont
if(rc(l)>smax_rc)
smax_rc=rc(l);
C_smax(1,:)=S(l,:);
smax_pos=l;
end
end
end
%Find the Third Maximum correlation tmax_rc and its corresponding codewords
for ll=1:1:256
if ll~=max_pos %skip the maximum piont
if ll~=smax_pos %skip the sub-maximum piont
if(rc(ll)>tmax_rc)
tmax_rc=rc(ll);
C_tmax(1,:)=S(ll,:);
end
end
end
end
%%%-------------calculate the output LLR----------------%%%
for num=1:8
%For those information bits different between C_max C_smax
if(C_smax(num)~=C_max(num))
softMetric_out(floor(i/2)+num)= (1/2)*((C_max(num)*max_rc)-(C_max(num)*smax_rc));
else if(C_tmax(num)~=C_max(num)) %for the same bits within C_max and C_smax, but different from C_tmax
softMetric_out(floor(i/2)+num)= (1/2)*((C_max(num)*max_rc)- (C_max(num)*tmax_rc));
else % for the same bits within C_max, C_smax and C_tmax
softMetric_out(floor(i/2)+num)= (1/2)*((C_max(num)*max_rc)- 0.875*(C_max(num)*tmax_rc));
end
end
end
% softMetric_out(1:8)
end
for t=1:1:3744
if ((softMetric_out(1,t)>0))
softMetricout2(t)=floor(softMetric_out(1,t));
else
softMetricout2(t)=round(softMetric_out(1,t));
end
%
% else
% end
if(softMetricout2(t)>(2^(Integer_bit+Decimal_bit)-1)) softMetricout2(t)=2^(Integer_bit+Decimal_bit)-1;
else if(softMetricout2(t)<=-(2^(Integer_bit+Decimal_bit))) softMetricout2(t)=-(2^(Integer_bit+Decimal_bit));
end
end
end
end
