| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 418 人关注过本帖
标题:求助编程高手:设计滤波器软件时遇到一个问题
只看楼主 加入收藏
guquan
Rank: 1
等 级:新手上路
帖 子:2
专家分:0
注 册:2012-3-31
收藏
 问题点数:0 回复次数:1 
求助编程高手:设计滤波器软件时遇到一个问题
求助编程高手:设计滤波器软件时遇到一个问题,对高手来说应该不难,可我暂时难以解决,求各位高手专家、仁人志士指教,不胜感激…
设计的是FIR滤波器软件, 问题是窗长度(程序中变量m)只能从小到大选取,且不能返回,否则就会出错(如下图示)。

确定窗长度m即可求得数组coef中的前m项,也是实际有用m个数,数组中其他元素为0。
GetCtrlVal (panelHandle, PANEL_ORDER, &m);   // 获取窗长度值
……
……
Ksr_LPF(samplefrequency,cutofffrequency,m,coef,beta) ; // coef即单位脉冲响应,由此函数产生。
数组coef产生无误,并始终显示正确。

以下为对coef进行快速傅里叶变换(FFT)的程序:
double *img_coef,*mag_coef,*phase_coef;
img_coef=malloc (n*sizeof (double));
mag_coef=malloc (n*sizeof (double));
phase_coef=malloc (n*sizeof (double));
    for(a=0;a<N;A++)
    {
     img_coef[a]=0;
    }
Copy1D(coef,n,coef_1);  // 复制数组元素
FFT (coef_1,img_coef,n);  // 进行FFT变换
ToPolar1D(coef_1,img_coef,n,mag_coef,phase_coef);  // 极坐标变换
DeleteGraphPlot (panelHandle, PANEL_GRAPH9, -1, VAL_IMMEDIATE_DRAW); //清除graph控件中图形
PlotY(panelHandle,PANEL_GRAPH9,mag_coef,400,VAL_DOUBLE,VAL_THIN_LINE,VAL_EMPTY_SQUARE, VAL_SOLID, 1, VAL_BLACK);  //在graph控件上绘出滤波器频率响应曲线
    free(img_coef);
    free(mag_coef);
    free(phase_coef);

 图传不上去,详见附件,谢谢各位!!
 
搜索更多相关主题的帖子: 滤波器 软件 编程高手 不胜感激 
2012-03-31 20:32
guquan
Rank: 1
等 级:新手上路
帖 子:2
专家分:0
注 册:2012-3-31
收藏
得分:0 
2012-03-31 20:33
快速回复:求助编程高手:设计滤波器软件时遇到一个问题
数据加载中...
 
   



关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.017163 second(s), 10 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved