基于CCS实现带通滤波器的DSP实现

带通数字滤波器的DSP实现

设计内容及要求:

已知x1(n)=sin(2pif11nT);

x2(n)=0.5sin(2pif12n*T)

x3(n)=0.8sin(2pif13n*T);

x(n) = x1(n)+ x2(n)+ x3(n)

其中:f11=100; f12=400; f13=1000; fs=10000

要求:采用c语言编程,分别设计FIR带通滤波器,把f11和f13滤掉,保留f12。

MATLAB部分
1)设置好FIR滤波系数
在这里插入图片描述
2)将系数导出为MATLAB变量
在这里插入图片描述
覆盖原数据
在这里插入图片描述
导出Num数据
在这里插入图片描述
CCS部分
1)将FIR滤波系数粘贴到CCS中在这里插入图片描述
代码部分

#include "math.h"
#define Fs 10000         /* 采样频率 */
#define T  1/Fs         /* 采样时间 */
#define f1 100          /* 正弦信号1频率 */
#define f2 400         /* 正弦信号2频率 */
#define f3 1000         /* 正弦信号3频率 */
#define PI 3.1415926
#define w1 (2*PI*f1*T)  /* 正弦信号1数字频率=2*pi*f1/Fs =2*pi/100*/
#define w2 (2*PI*f2*T)  /* 正弦信号2数字频率=2*pi*f2/Fs =2*pi*4/10*/
#define w3 (2*PI*f3*T)  /* 正弦信号3数字频率=2*pi*f3/Fs =2*pi/10*/
#define a1 1        /* 正弦信号1幅度 */
#define a2 0.5         /* 正弦信号2幅度 */
#define a3 0.8         /* 正弦信号3幅度 */
#define FIRNUMBER 100//滤波阶数

#define FIRNUMBER 100//滤波阶数
#define FIRNUMBER 100//滤波阶数
extern int fir(int *,int *,unsigned int,int );
/* Low-pass FIR fi
  • 16
    点赞
  • 146
    收藏
    觉得还不错? 一键收藏
  • 8
    评论
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; } }

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 8
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值