|
本帖最后由 Nexus 于 2012-11-10 20:01 编辑
用MatLab的fdatool设计了滤波器,然后导出头文件,参数什么的先不管,有以下疑问:
1.NL[MWSPT_NSEC]这个数组是干嘛用的?
2.const real32_T NUM[MWSPT_NSEC][3] = {....这个为什么要二维数组呢?如果单纯是系数的话,一维的就够了。
3.const real32_T NUM[MWSPT_NSEC][3] 和 const real32_T DEN[MWSPT_NSEC][3] 这两个数组如何配合使用?分别什么意思啊·???
4.这个图片里Order是22,那这个section是什么意思?
求指教大神,数字信号处理学过,都忘记了。。。
最后贴一下代码:
/*
* Filter Coefficients (C Source) generated by the Filter Design and Analysis Tool
*
* Generated by MATLAB(R) 7.12 and the Signal Processing Toolbox 6.15.
*
* Generated on: 10-Nov-2012 19:57:29
*
*/
/*
* Discrete-Time IIR Filter (real)
* -------------------------------
* Filter Structure : Direct-Form II, Second-Order Sections
* Number of Sections : 11
* Stable : Yes
* Linear Phase : No
*/
/* General type conversion for MATLAB generated C-code */
#include "tmwtypes.h"
/*
* Expected path to tmwtypes.h
* C:\Program Files\MATLAB\R2011a\extern\include\tmwtypes.h
*/
/*
* Warning - Filter coefficients were truncated to fit specified data type.
* The resulting response may not match generated theoretical response.
* Use the Filter Design & Analysis Tool to design accurate
* single-precision filter coefficients.
*/
#define MWSPT_NSEC 23
const int NL[MWSPT_NSEC] = { 1,3,1,3,1,3,1,3,1,3,1,3,1,3,1,3,1,3,1,3,1,3,1 };
const real32_T NUM[MWSPT_NSEC][3] = {
{
0.09686300904, 0, 0
},
{
1, 2, 1
},
{
0.08955428004, 0, 0
},
{
1, 2, 1
},
{
0.08345051855, 0, 0
},
{
1, 2, 1
},
{
0.07838561386, 0, 0
},
{
1, 2, 1
},
{
0.07421958447, 0, 0
},
{
1, 2, 1
},
{
0.07083738595, 0, 0
},
{
1, 2, 1
},
{
0.06814649701, 0, 0
},
{
1, 2, 1
},
{
0.06607432663, 0, 0
},
{
1, 2, 1
},
{
0.06456578523, 0, 0
},
{
1, 2, 1
},
{
0.06358129531, 0, 0
},
{
1, 2, 1
},
{
0.06309531629, 0, 0
},
{
1, 2, 1
},
{
1, 0, 0
}
};
const int DL[MWSPT_NSEC] = { 1,3,1,3,1,3,1,3,1,3,1,3,1,3,1,3,1,3,1,3,1,3,1 };
const real32_T DEN[MWSPT_NSEC][3] = {
{
1, 0, 0
},
{
1, -1.530096531, 0.9175485373
},
{
1, 0, 0
},
{
1, -1.414644241, 0.7728613615
},
{
1, 0, 0
},
{
1, -1.318226099, 0.652028203
},
{
1, 0, 0
},
{
1, -1.238218307, 0.5517607927
},
{
1, 0, 0
},
{
1, -1.172409654, 0.4692879915
},
{
1, 0, 0
},
{
1, -1.118982673, 0.4023322761
},
{
1, 0, 0
},
{
1, -1.076476097, 0.349062115
},
{
1, 0, 0
},
{
1, -1.043743134, 0.3080404103
},
{
1, 0, 0
},
{
1, -1.019913435, 0.2781765163
},
{
1, 0, 0
},
{
1, -1.004361987, 0.2586871684
},
{
1, 0, 0
},
{
1, -0.9966852069, 0.249066487
},
{
1, 0, 0
}
};
|
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有帐号?注册
x
阿莫论坛20周年了!感谢大家的支持与爱护!!
如果天空是黑暗的,那就摸黑生存;
如果发出声音是危险的,那就保持沉默;
如果自觉无力发光,那就蜷伏于牆角。
但是,不要习惯了黑暗就为黑暗辩护;
也不要为自己的苟且而得意;
不要嘲讽那些比自己更勇敢的人。
我们可以卑微如尘土,但不可扭曲如蛆虫。
|