??眾所周知，Matlab 中的 Filter Designer 可以直接生成 FIR 濾波器的 verilog 代碼，可以方便地生成指定階數(shù)、指定濾波器參數(shù)的高通、低通、帶通濾波器，生成的 verilog 代碼也可以指定輸入輸出信號的類型和位寬。然而其生成的代碼實(shí)在算不上美觀，復(fù)用性也很差，要實(shí)現(xiàn)不同濾波特性的切換就要生成多個濾波器的代碼。
??出于以上考慮，自己設(shè)計(jì)實(shí)現(xiàn)了 FIR 濾波器的通用 verilog 代碼，其濾波器參數(shù)通過接口輸入，從而可以通過輸入不同的參數(shù)獲得相應(yīng)的濾波結(jié)果。verilog 代碼如下：
/*
* file : FIR_filter.v
* author : 今朝無言
* date : 2023-07-03
* version : v1.0
* description : FIR 濾波器
*/
module FIR_filter(
input clk,
input rst_n,
input [16*N-1:0] filter_params,
input signed [15:0] data_in,
output reg signed [15:0] data_out
);
parameter N = 32; //濾波器參數(shù)個數(shù)
parameter div_N = 16; //sum結(jié)果除 2^div_N，作為 filter 的輸出
//FIR 濾波器參數(shù)
reg signed [15:0] b[0:N-1];
integer m;
always @(*) begin
for(m=0; m<N; m=m+1) begin
b[m] <= filter_params[(m << 4) +: 16];
end
end
reg signed [15:0] shift_reg[0:N-1];
integer i;
always @(posedge clk) begin
if(~rst_n) begin
for(i=N-1; i>=0; i=i-1) begin
shift_reg[i] <= 16'd0;
end
end
else begin
for(i=N-1; i>0; i=i-1) begin
shift_reg[i] <= shift_reg[i-1];
end
shift_reg[0] <= data_in;
end
end
reg signed [31:0] multi[0:N-1];
integer j;
always @(*) begin
for(j=0; j<N; j=j+1) begin
multi[j] <= shift_reg[j] * b[j];
//這里可以考慮使用multiplier IP核，使用LUT搭建（而這里直接乘使用的是DSP資源，一般的FPGA芯片只有幾百個）
end
end
reg signed [47:0] sum;
integer k;
always @(*) begin
sum = 0;
for(k=0; k<N; k=k+1) begin
sum = sum + multi[k];
end
end
always @(posedge clk) begin
data_out <= sum[47-div_N : 32-div_N];
end
endmodule
Lowpass Filter示例
??當(dāng)濾波器階數(shù)較高時，濾波器參數(shù)如何給出無疑是個麻煩事，因此又編寫了 matlab 代碼，可以一鍵生成所需的 .v 文件以實(shí)現(xiàn)參數(shù)的配置：
%-----------FIR濾波器參數(shù)（生成.v）-----------------
clc,clear,close all
fs=1e6;
N=20;
Wn=0.1;
b = fir1(N, Wn); % 默認(rèn)Hamming窗
freqz(b,1,512)
%% pramas
B=floor(b*65536);
B=B';
%% test
t=0:1/fs:1e-3;
s=(mod(t,1e-4)<5e-5)*1.0;
%s_filt=filter(B,1,s)/65536;
for i=1:size(s,2)-N-1
s_filt(i)=s(i:i+N)*double(B)/65536;
end
figure
subplot(2,1,1)
plot(t,s)
subplot(2,1,2)
plot(t(1:end-N-1),s_filt)
%% 生成.v
filename='FIR_params';
fid=fopen(['./v/',filename,'.v'],'w');
fprintf(fid,['/* ','\n',...
' * file\t\t\t: ',filename,'.v','\n',...
' * author\t\t: 今朝無言','\n',...
' * date\t\t\t: 2023-07-04','\n',...
' * version\t\t: v1.0','\n',...
' * description\t: FIR 濾波器','\n',...
' */','\n']);
fprintf(fid,['module ',filename,'(','\n',...
'output\t[',num2str(size(B,1)*16-1),':0]\tparams\n',...
');\n\n']);
for i=1:size(B,1)
if(B(i)>=0)
hex=dec2hex(B(i),4);
else
hex=dec2hex(65536+B(i),4);
end
fprintf(fid,['assign\t','params[',...
num2str(i*16-1),':',num2str((i-1)*16),...
']\t= 16','''','h',hex,';\n']);
end
fprintf(fid,'\nendmodule\n');
fclose(fid);
??testbench與測試結(jié)果如下
`timescale 1ns/100ps
module FIR_filter_tb();
reg clk_100M = 1'b1;
always #5 begin
clk_100M <= ~clk_100M;
end
localparam N = 20; //FIR濾波器階數(shù)
wire [16*(N+1)-1:0] filter_params;
FIR_params_0d1 FIR_params_inst(
.params (filter_params)
);
reg [15:0] data_in;
wire signed [15:0] data_out;
FIR_filter #(.N(N+1))
FIR_filter_inst2(
.clk (clk_100M),
.filter_params (filter_params), //濾波器參數(shù)
.data_in (data_in),
.data_out (data_out)
);
reg [7:0] cnt = 8'd0;
always @(posedge clk_100M) begin
cnt <= cnt + 1'b1;
if(cnt<100) begin
data_in <= -10000;
end
else if(cnt<200) begin
data_in <= 10000;
end
else begin
data_in <= 0;
end
end
initial begin
#10000;
$stop;
end
endmodule
Hilbert 示例
??使用以上 FIR 濾波器代碼，還可以實(shí)現(xiàn)許多其他濾波功能，比如常用的 90 度相移，可以使用 Hilbert 變換實(shí)現(xiàn)，Hilbert 濾波器參數(shù)的 matlab 生成代碼如下
%-----------------Hilbert----------------------
clc,clear,close all
%% Hilbert
N=200;
% method 1 這種直接通過 h(n) 表達(dá)式生成的更為精確，推薦
n=(1:floor((N-1)/2));
b1=(1-(-1).^n)./(pi.*n);
if mod(N,2)==0
b1=[0,b1,0,-b1(end:-1:1)]';
else
b1=[0,b1,-b1(end:-1:1)]';
end
% method 2 構(gòu)造 Hilbert 的頻域特性，經(jīng) IFFT 獲得
H=[-1j*ones(1,floor((N+1)/2)),1j*ones(1,floor(N/2))];
b2=ifft(H);
b2=real(b2)';
b=b1;
freqz(b,1,100)
%% Filter
fs=1e3;
t=0:1/fs:1;
s=5*sin(2*pi*10*t);
% f >= fs/N 時，可以由很好的90度移相
s2=filter(b,1,s);
figure
hold on
plot(t,s,'r-')
plot(t,s2,'b--')
hold off
%% 量化
B=floor(b*32768);
s3=filter(B,1,s)/32768;
figure
hold on
plot(t,s,'r-')
plot(t,s3,'b--')
hold off
%% 生成params.v
filename='Hilbert_params';
fid=fopen(['./v/',filename,'.v'],'w');
fprintf(fid,['/* ','\n',...
' * file\t\t\t: ',filename,'.v','\n',...
' * author\t\t: 今朝無言','\n',...
' * date\t\t\t: 2023-08-04','\n',...
' * version\t\t: v1.0','\n',...
' * description\t: FIR濾波器參數(shù)（Hilbert）',...
' N=',num2str(N),'\n',...
' */','\n']);
fprintf(fid,['module ',filename,'(','\n',...
'output\t[',num2str(size(B,1)*16-1),':0]\tparams\n',...
');\n\n']);
for i=1:size(B,1)
if(B(i)>=0)
hex=dec2hex(B(i),4);
else
hex=dec2hex(65536+B(i),4);
end
fprintf(fid,['assign\t','params[',...
num2str(i*16-1),':',num2str((i-1)*16),...
']\t= 16','''','h',hex,';\n']);
end
fprintf(fid,'\nendmodule\n');
fclose(fid);
??仿真結(jié)果如下文章來源：http://www.zghlxwxcb.cn/news/detail-638141.html
文章來源地址http://www.zghlxwxcb.cn/news/detail-638141.html
到了這里，關(guān)于通用FIR濾波器的verilog實(shí)現(xiàn)（內(nèi)有Lowpass、Hilbert參數(shù)生成示例）的文章就介紹完了。如果您還想了解更多內(nèi)容，請?jiān)谟疑辖撬阉鱐OY模板網(wǎng)以前的文章或繼續(xù)瀏覽下面的相關(guān)文章，希望大家以后多多支持TOY模板網(wǎng)！
