基于FPGA的数字信号处理系统开发笔记（二）FIR滤波器

原创于 2025-04-11 22:34:39 发布 · 1.8k 阅读 ·

本内容遵循CC 4.0 BY-SA版权协议

GEO检测

收录于

基于FPGA的数字信号处理

代码可运行

想深耕嵌入式？这个专辑值得收藏

MCU、FPGA、工控、传感器一站式学习，实战项目直接抄

一键订阅

一，matlab产生初始信号并混合

主文件代码如下，引用的fft_plot函数为自己定义的函数，在第二段代码

特别注意：考虑到上次结尾遗留的图像突变的问题，本次采样点数（数据深度）N改为1000

这样就可以保证上一个周期最后一个点，和下一个周期第一个点都为0，此时图像会连续不跳变。

clear;
    Fs=10000;%采样频率
    N=1000; %采样点数
    Ts=1/Fs; %采样周期

%%输入信号参数
    A = 3;              %幅值
    t =(0:N-1)*Ts;   %采样时间长度
    ft_0 = 100;         %信号1频率
    vi_0 = zeros(1,N);  %创建一个1*N的全零数组
    ft_1 = 2000;        %信号1频率
    vi_1 = zeros(1,N);  %创建一个1*N的全零数组
    vi = zeros(1,N);    %创建一个1*N的全零数组
    
    for n = 1:N
        vi_0(n) = A * sin(2*pi*ft_0*t(n)); %输入信号1
        vi_1(n) = A * sin(2*pi*ft_1*t(n)); %输入信号2
        vi(n) = vi_0(n) + vi_1(n);         %混合信号1和信号2
    end
    
    %%作图
    figure(1)
    subplot(3,1,1);% 3行1列第一个
    plot(t,vi_0,'k');
    subplot(3,1,2);% 3行1列第二个
    plot(t,vi_1,'k');
    subplot(3,1,3);% 3行1列第三个
    plot(t,vi,'k');
    
    %%快速傅里叶变换
    [Fre0,Amp0] = fft_plot(vi_0,Fs,N); %%fft_plot是自己定义的函数
    [Fre1,Amp1] = fft_plot(vi_1,Fs,N);
    [Fre,Amp] = fft_plot(vi,Fs,N);

    %%作图
    figure(2)
    subplot(3,1,1);% 3行1列第一个
    plot(Fre0,Amp0,'k');
    subplot(3,1,2);% 3行1列第二个
    plot(Fre1,Amp1,'k');
    subplot(3,1,3);% 3行1列第三个
    plot(Fre,Amp,'k');

fft_plot代码如下

function [Fre, Amp] = fft_plot(signal, fs, num_samples)
    % 计算FFT
    f_fft = fft(signal);
    N = length(f_fft); % FFT之后的序列长度
    
    % 提取频谱正半部分
    Y = f_fft(1:N/2+1);
    Y(2:end-1) = 2*Y(2:end-1); % 除了直流和Nyquist频率，其他频率成分乘以2
    
    % 计算幅值
    Y_amp = abs(Y)/N;
    
    % 计算频率向量
    df = fs/N;  %%fs采样频率 N为FFT的点数
    Fre = (0:1:N/2)*df;
    
    % 输出幅值
    Amp = Y_amp;
    
    % 如果需要在此函数中绘制图形，可以取消注释以下代码
    % 但通常最好在主脚本中处理绘图
    % figure;
    % plot(Fre, Amp);
    % title('Frequency Spectrum');
    % xlabel('Frequency (Hz)');
    % ylabel('Amplitude');
end

运行结果如下