关于C代码实现差分方程滤波filter
来源:互联网 发布:小区网络装iptv服务器 编辑:程序博客网 时间:2024/06/06 00:59
很久之前就已经把代码编译好了,没有时间整理:今天整理一下,已备后查;
整个差分方式的实现过程如下:公式不好编辑:
matlab实现了以后,用C照着写不是什么难事,找到2种实现方法,一种是照着filter差分的方式去写,误差说是在10^-6,另一种说是没有误差,实现过程没有去仔细研究,觉得小数点后面6位的误差也是可以的了。
代码可见:点击打开链接
核心 过程不过是两个for循环:
for (i=0;i<DataLen;i++)
{
double sum_b_x=0;
double sum_a_y=0;
k=0;
for( j=i;j<=i+ord;j++)
{
sum_b_x=sum_b_x+bz1[k]*up[j];
/*printf("bz1[%d]*up[%d]=\n",k,j);*/
k++;
}
/*printf("%f\n",sum_b_x);*/
//double sum_a_y=0;
k=0;
for( j=i;j<=i+ord-1;j++)
{
sum_a_y=sum_a_y+aa[k]*yp[j];
/*printf("aa[%d]*yp[%d]=\n",k,j);*/
k++;
}
/*printf("%f\n",sum_a_y);*/
yp[i+ord]=sum_b_x-sum_a_y;
/*printf("---%d\n,%e\n",i+ord,yp[i+ord]);*/
}
0 0
- 关于C代码实现差分方程滤波filter
- 差分方程建模
- 微分方程差分方程
- 热传导方程的差分格式原理与matlab实现
- 等差数列和的差分方程
- 等比数列和的差分方程
- 对差分方程应用积分方法
- 对差分方程应用微分方法
- 特殊形式的差分方程
- 特殊形式的差分方程(改变)
- 差分方程与exp(At)
- 椭圆形方程的差分解法
- 双边滤波(bilateral filter)灰度图 matlab实现代码
- 双边滤波(bilateral filter)彩色图 matlab实现代码
- 双边滤波(bilateral filter)快速算法 matlab实现代码
- 关于差分信号
- 关于差分信号
- 关于差分
- JZOJ3997. 树
- 文章标题Fire Net
- Intellij IDEA常用快捷键整理
- Fragments
- vim
- 关于C代码实现差分方程滤波filter
- css HACK
- 导入可执行文件
- build.gradle常见配置方式
- 基于Spring框架的简单多数据源切换解决办法
- HTML之DocType的几种类型
- Unity学习之Tanks官网教程
- 17 - 03 - 02 图解HTTP(06)
- Iterable 迭代器