MATLAB实现带通滤波:从时域到频域的信号处理

需积分: 19 28 下载量 130 浏览量 更新于2024-09-08 收藏 4KB TXT 举报
该资源是关于使用MATLAB进行带通滤波器设计和信号处理的实践案例。通过创建和分析不同类型的信号,如白噪声和调频信号,并应用带通滤波器,来展示如何在时域和频域中理解和评估滤波效果。 在MATLAB中,带通滤波器通常用于从信号中选取特定频率范围内的成分,去除不需要的高频或低频噪声。在这个例子中,首先生成一个长度为1024的白噪声序列,然后调整其均值和标准差以规范化信号。白噪声在时域中的表示通过`plot(noise)`展示,并通过快速傅里叶变换(FFT)转换到频域,以`plot(freq22, abs(g))`显示频谱。 接着,生成一个调频信号`x1`,然后使用调频-modulation函数`modulate`将其调制到中心频率`fc`,采样率`fs`。在时域和频域中分别展示调频信号,以便观察调制的效果。 为了设计带通滤波器,选择3000Hz至4000Hz的频率范围,并利用`fir1`函数创建一个汉明窗函数的有限冲击响应(FIR)滤波器。滤波器的截止频率`wlp`和`whp`设置为3000*2*pi和4000*2*pi,以对应于3000Hz和4000Hz。滤波器的系数通过`fir1`计算得到,然后通过FFT转换为频域表示。 在最后的步骤中,应用这个滤波器到调频信号`x2`上,然后再次查看时域和频域的表现,以评估滤波器的性能。在时域中,可以观察到滤波器如何平滑信号并减少噪声;在频域中,可以看到滤波器有效地选择了特定频率范围内的信号分量,而抑制了其他频率。 总结来说,这个示例涵盖了以下几个关键知识点: 1. MATLAB信号处理基础:包括生成随机信号、时域和频域分析。 2. 快速傅里叶变换(FFT):用于将信号从时域转换到频域,以便分析其频率成分。 3. 调频信号的生成:通过`modulate`函数实现信号的调频。 4. 带通滤波器设计:使用`fir1`函数创建FIR滤波器,指定汉明窗函数和截止频率。 5. 滤波器的应用:将设计的滤波器应用于信号,以提取特定频率范围内的信息。 6. 信号滤波效果的评估:通过比较滤波前后的时域和频域图像,评估滤波器的性能。 通过这个案例,学习者可以深入理解带通滤波器的工作原理及其在实际信号处理中的应用。
2019-04-16 上传
float DigFil(invar, setic) float invar; int setic; /******************************************************************************/ /* Filter Solutions Version 2009 Nuhertz Technologies, L.L.C. */ /* www.nuhertz.com */ /* +1 602-279-2448 */ /* 3rd Order Band Pass Butterworth */ /* Bilinear Transformation with Prewarping */ /* Sample Frequency = 5.000 KHz */ /* Standard Form */ /* Arithmetic Precision = 4 Digits */ /* */ /* Center Frequency = 300.0 Rad/Sec */ /* Pass Band Width = 20.00 Rad/Sec */ /* */ /******************************************************************************/ /* */ /* Input Variable Definitions: */ /* Inputs: */ /* invar float The input to the filter */ /* setic int 1 to initialize the filter to zero */ /* */ /* Option Selections: */ /* Standard C; Initializable; Internal States; Not Optimized; */ /* */ /* There is no requirement to ever initialize the filter. */ /* The default initialization is zero when the filter is first called */ /* */ /******************************************************************************/ /* */ /* This software is automatically generated by Filter Solutions */ /* no restrictions from Nuhertz Technologies, L.L.C. regarding the use and */ /* distributions of this software. */ /* */ /******************************************************************************/ { float sumnum=0.0, sumden=0.0; int i=0; static float states[6] = {0.0,0.0,0.0,0.0,0.0,0.0}; static float znum[7] = { -7.968e-09, 0.0, 2.39e-08, 0.0, -2.39e-08, 0.0, 7.968e-09 }; static float zden[6] = { .992, -5.949, 14.88, -19.86, 14.92, -5.981 }; if (setic==1){ for (i=0;i<6;i++) states[i] = [i] = [i]*invar; return 0.0; } else{ sumnum = sumden = 0.0; for (i=0;i<6;i++){ sumden += states[i]*zden[i]; sumnum += states[i]*znum[i]; if (i<5) states[i] = states[i+1]; } states[5] = invar-sumden; sumnum += states[5]*znum[6]; return sumnum; } }