24小时热门版块排行榜    

查看: 744  |  回复: 4

baobiao007

木虫 (职业作家)

中国特色

[交流] 【分享】线性卷积c程序 已有3人参与

#include
/*  y(n)=x(n)*h(n)
     m--length of x(n);
     n--length of h(n);
     l=m+n-1  length of y(n)
*/
void conv(float x[],int m,float h[],int n,float y[],int l)
{
        int i,j,k;
        for(i=0; i         {
                y=0.0;
                for(j=0; j                 {   k=i-j;
                    if(k>=0 && k                         y += x[j]*h[k];
                }
        }
}
void main()
{              //参数m,n,a1,a2,b1,b2,x[],h[]可任意改动
        const int m=3;//length of x(n)
        const int n=3;//length of h(n)
        const int l=m+n-1;//length of y(n)
        const int a1=0;//x(n)左下标,表示x(n)={x(0),x(1),x(2)}
        const int a2=2;//x(n)右下标
        const int b1=4;//h(n)左下标,h(n)={h(4),h(5),h(6)}
        const int b2=6;//h(n)右下标
        int c1;//卷积结果序列的左下标
        int c2;//卷积结果序列的右下标
        int i;
        float x[m]={1,1,1};
        float h[n]={1,1,1};
        float y[l];
        c1=a1+b1;
        c2=a2+b2;
        conv(x,m,h,n,y,l);
                printf("result:\n";
        for(i=0; i              printf("y(%d)=%f\n",c1+i,y);
}
回复此楼
我同意叔本华的观点,人们投身艺术和科学领域的强烈愿望之一就是逃离痛苦、残酷和枯燥无味的现实生活,逃离自己飘忽不定的七情六欲的桎梏。--爱因斯坦
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

baobiao007

木虫 (职业作家)

中国特色

主要用于合成地震记录,搞地震数据处理必须掌握的第一个程序
我同意叔本华的观点,人们投身艺术和科学领域的强烈愿望之一就是逃离痛苦、残酷和枯燥无味的现实生活,逃离自己飘忽不定的七情六欲的桎梏。--爱因斯坦
2楼2011-02-12 19:52:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xuliforever

木虫 (正式写手)

木虫村村委书记


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by baobiao007 at 2011-02-12 19:52:06:
主要用于合成地震记录,搞地震数据处理必须掌握的第一个程序

可以用到其他地方吗?
比如说
用于重金属的污染评价中
只要不下流,我们就主流!
3楼2011-02-13 17:29:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

baobiao007

木虫 (职业作家)

中国特色

不行吧  相差太远了
我同意叔本华的观点,人们投身艺术和科学领域的强烈愿望之一就是逃离痛苦、残酷和枯燥无味的现实生活,逃离自己飘忽不定的七情六欲的桎梏。--爱因斯坦
4楼2011-02-13 17:49:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

04lrcg

新虫 (初入文坛)

还不错
5楼2011-06-25 12:36:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 baobiao007 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见