【代码更新】标准差

 

【代码更新】标准差

每513个点中，对前512个取其标准差，最后一个点不在计算范围内；

具体的过程是，将512点的平方和均值，减去，512个点的均值的平方

方差=用平方和的均值减去均值的平方得到方差

//用平方和的均值减去均值的平方得到方差

module downSampling(
    input  wire                 clk,
    input  wire               rst_n,
    
    input   wire           din_valid,
    input   wire    [7:0]        din,    
    
    // output reg [7:0] data_in,
    
    output  wire   down_sampling_dout_valid,
    output  wire   [7:0]  down_sampling_dout    
);
//test
// reg [31:0] cnt0_test;
// wire      add_cnt0_test;
// wire      end_cnt0_test;
// always@(posedge clk or negedge rst_n)begin
    // if(rst_n==1'b0)begin
        // cnt0_test<=0;
    // end
    // else if(add_cnt0_test)begin
        // if(end_cnt0_test)
            // cnt0_test<=0;
        // else
            // cnt0_test<=cnt0_test+1;
    // end 
// end 
// assign add_cnt0_test = din_valid;
// assign end_cnt0_test = add_cnt0_test && cnt0_test == 100000-1;

// reg [7:0] data_in;
// always@(posedge clk or negedge rst_n)begin
    // if(rst_n==1'b0)begin
        // data_in<=0;
    // end
    // else if(add_cnt0_test && cnt0_test < 'd20000)begin
        // data_in<='d1;
    // end 
    // else if(add_cnt0_test && cnt0_test == 'd20000 )begin
        // data_in<='d35;
    // end 
    // else if(add_cnt0_test && cnt0_test > 'd20000 && cnt0_test < 'd20010 )begin
        // data_in<='d1;
    // end    
    // else if(add_cnt0_test && cnt0_test == 'd20010 )begin
        // data_in<='d35;
    // end    
    // else if(add_cnt0_test && cnt0_test > 'd20010 && cnt0_test < 'd60000 )begin
        // data_in<='d1;
    // end        
// end 
// wire    [7:0]        din;
// assign din = data_in;
//-----------------------------------------------------------------
//8位
reg     [16:0]  sum;
reg     [24:0]  square_sum;
wire    [15:0]  square_mean;
wire    [15:0]  mean_square;
wire    [33:0]  sum_square; 
wire    [15:0]  variance;
reg             valid_variance; 
//
reg [9:0] cnt0;
wire      add_cnt0;
wire      end_cnt0;
always@(posedge clk or negedge rst_n)begin
    if(rst_n==1'b0)begin
        cnt0<=0;
    end
    else if(add_cnt0)begin
        if(end_cnt0)
            cnt0<=0;
        else
            cnt0<=cnt0+1;
    end 
end 
assign add_cnt0 = din_valid;//1'b1 din_valid
assign end_cnt0 = add_cnt0 && cnt0 == 513-1;   //----512 ---9
//
always@(posedge clk or negedge rst_n)begin
    if(rst_n==1'b0)begin
        sum<='d0;
    end
    else if(add_cnt0)begin
        if(end_cnt0)
            sum<='d0;
        else
            sum<=sum + din;        
    end         
end 
//
always@(posedge clk or negedge rst_n)begin
    if(rst_n==1'b0)begin
       square_sum <= 'd0; 
    end
    else if(add_cnt0)begin
        if(end_cnt0)
            square_sum <= 'd0;
        else
            square_sum <= square_sum + din*din;            
    end 
    else begin
        square_sum <= 'd0;    
    end 
end 
//
always@(posedge clk or negedge rst_n)begin
    if(rst_n==1'b0)begin
       valid_variance<='b0; 
    end
    else if(add_cnt0 && cnt0 == 'd511)begin    
        valid_variance<='b1;
    end 
    else begin
        valid_variance<='b0;    
    end 
end
//
assign sum_square       = (add_cnt0&& end_cnt0)?(sum*sum):'d0; 
//8位
assign mean_square      = sum_square[33-:16]; //----------[33:18]
assign square_mean      = square_sum[24-:16];//-----------[24: 9]
assign variance         = (add_cnt0&&end_cnt0&&(square_mean > mean_square))?(square_mean - mean_square):'d0;
//
wire       m_axis_dout_tvalid;
wire [15:0] m_axis_dout_tdata;
cordic_0 cordic_0_uut (
  .aclk(clk),                                         // input wire aclk
  .s_axis_cartesian_tvalid(valid_variance),           // input wire s_axis_cartesian_tvalid
  .s_axis_cartesian_tdata(variance),                  // input wire [15 : 0] s_axis_cartesian_tdata
  .m_axis_dout_tvalid(m_axis_dout_tvalid),            // output wire m_axis_dout_tvalid
  .m_axis_dout_tdata(m_axis_dout_tdata)               // output wire [15 : 0] m_axis_dout_tdata
);
assign down_sampling_dout_valid=     m_axis_dout_tvalid;
assign      down_sampling_dout = m_axis_dout_tdata[7:0];
//8位
// ila_3 ila_3_uut (
    // .clk(clk), // input wire clk

    // .probe0(din_valid          ), // input wire [0:0]  probe0  
    // .probe1(din                ), // input wire [7:0]  probe1  //din_test    
    // .probe2(m_axis_dout_tvalid ), // input wire [0:0]  probe2 
    // .probe3(m_axis_dout_tdata  ), // input wire [15:0]  probe3 
    // .probe4(cnt0               ), // input wire [9:0]  probe4 
    // .probe5(sum                ), // input wire [16:0]  probe5 
    // .probe6(square_sum         ), // input wire [24:0]  probe6 
    // .probe7(variance           ), // input wire [15:0]  probe7 
    // .probe8(valid_variance     )  // input wire [0:0]  probe8
// );
endmodule
//test
// wire  [7:0] din_test;
// assign      din_test=(add_cnt0 && cnt0 == 'd127 && cnt0 < 'd150)?'d0:din; 
//
//16位
// reg     [24:0]  sum;
// reg     [40:0]  square_sum;
// wire    [31:0]  square_mean;
// wire    [31:0]  mean_square;
// wire    [49:0]  sum_square; 
// wire    [31:0]  variance;
// reg             valid_variance; 
//16位
// assign mean_square      = sum_square[49-:32]; //----------[49:18]
// assign square_mean      = square_sum[40-:32];//-----------[40: 9]
//16位
// ila_3 ila_3_uut (
    // .clk(clk), // input wire clk

    // .probe0(din_valid          ), // input wire [0:0]  probe0  
    // .probe1(din                ), // input wire [15:0]  probe1 
    // .probe2(m_axis_dout_tvalid ), // input wire [0:0]  probe2 
    // .probe3(m_axis_dout_tdata  ), // input wire [16:0]  probe3 
    // .probe4(cnt0               ), // input wire [9:0]  probe4 
    // .probe5(sum                ), // input wire [24:0]  probe5 
    // .probe6(square_sum         ), // input wire [40:0]  probe6 
    // .probe7(variance           ), // input wire [31:0]  probe7 
    // .probe8(valid_variance     )  // input wire [0:0]  probe8
// );