要用DSP做一个频响函数，网上找到了一个IIR系统的频率响应算法实例，有点看不太懂我是不是能使用

2829

问答对人有帮助，内容完整，我也想知道答案 0 我做的项目是用激振锤锤击悬臂梁，DSP采集激振锤返回的力信号的电压值和安装在悬臂梁上的加速度传感器返回的加速度信号的电压值，将这两路信号做频响函数，但不知道这个程序能不能应用，各位大神有空的帮忙看看哈，万分感谢。我采样频率设置的是6K，采样点数是512个点。程序如下： #include #include #include #include "msp.h" void miirres(float a[],float b[],int lb,int la,complex h[],int n) { /--------------------------------------------------------------------- Routinue miirres:to obtain the frequence responce (Amplitude and phase) H(exp(jw)) of system function H(z)=B(z)/A(z); b(0)+b(1)z^(-1)+ ... + b(lb)z^(-lb) H(z)=------------------------------------ 1.0 +a(1)z^(-1)+ ... + a(la)z^(-la) input parameters: a :la+1 dimensioned real array. a(0)=1.0, a(n) (n=1 to la) is the coef. of z^(-n). b :lb+1 dimensioned real array. b(n) (n=0 to la) is the coef. of z^(-n). n :data point number of frequency response from 0.0 to pi.the frequency resolution is pi/n. output parameters: H :n dimensioned complex array. frequency response is stored in H(0) to H(n-1), which is normalized from 0.0 to pi. Note: If H(k)=0/0,then let H(K)=H(k+1) (This is for the case of that H(exp(jw)) is a sinc function) : If la=0: for FIR system. in Chapter 2 ---------------------------------------------------------------------/ complex z,bsum,asum; int k,i; int isign; float pole,peak,pi,s,sf,sa,**,temp; pole=1.e-10; peak=1.e+10; pi=4.atan(1.); s=pi/(float)(n); isign=0; for(k=0;k { sf=sk; bsum.real=b[0];bsum.imag=0.0; if(lb>0) for(i=1;i<=lb;i++) { *=sfi; z.real=cos(*);z.imag=-sin(); bsum.real+=bz.real; bsum.imag+=bz.imag; } asum.real=1.; asum.imag=0.0; if(la!=0) for(i=1;i<=la;i++) { sa=sfi; z.real=cos(sa);z.imag=-sin(sa); asum.real+=az.real; asum.imag+=az.imag; } if(mabs(asum) { if(mabs(bsum) isign=1; else { h[k].real=peak; h[k].imag=0.0; } } else { temp=pow(mabs(asum),2.0); h[k].real=(bsum.realasum.real+bsum.imagasum.imag)/temp; h[k].imag=(bsum.imagasum.real-bsum.realasum.imag)/temp; if(isign==1) { h[k-1]=h[k]; isign=0; } } } return; } 0 已退回10积分
2016-6-3 10:21:04　　评论淘帖0 邀请回答您可以邀请以下用户，快速回答问题 × 柑橘与柠檬该类别下有 19 个回答。邀请回答 jinyi7016 该类别下有 11 个回答，其中被选为最佳答案 2 次。邀请回答 ruanjianming 该类别下有 11 个回答。邀请回答翟荣会该类别下有 10 个回答。邀请回答 24不可说该类别下有 9 个回答，其中被选为最佳答案 9 次。邀请回答 zuqiuxiaozi_srx 该类别下有 8 个回答。邀请回答陈该类别下有 8 个回答。邀请回答 dsphahaha 该类别下有 8 个回答。邀请回答 zercola 该类别下有 8 个回答。邀请回答武-之-魂该类别下有 8 个回答。邀请回答大蜜该类别下有 8 个回答。邀请回答云中月1 该类别下有 8 个回答。邀请回答 xxdjx0206 该类别下有 7 个回答。邀请回答 sunyu817 该类别下有 7 个回答。邀请回答北纬该类别下有 7 个回答。邀请回答 lz_1995 该类别下有 7 个回答。邀请回答 pannybei888 该类别下有 7 个回答。邀请回答王立峰该类别下有 7 个回答。邀请回答 forlovefy 该类别下有 6 个回答。邀请回答小_艺人该类别下有 6 个回答。邀请回答举报王相关推荐 • 绘制信号处理系统的频率响应曲线时，为何DAC和ADC响应下降？ 1348 • 关于运算放大器的频率响应异常的疑问：OPA847 4228 • 频率响应分析仪的使用问题 3737 • [求助]关于氢火焰的一些问题 2700 • 求助一下！一个电路容阻抗的问题~ 2843 • LSM9DS1的频率响应 1722 • 不稳定系统如何稳定？ 4798 • 如何画出这个传递函数的波特图或者频率响应特性 7222 • 哪些参数会影响麦克风频率响应？哪种评估频率响应的最佳方法？ 3351 • 在MCFW框架下，有OSD（SWOD）的图像叠加的算法在DSP端执行，请问这样的算法是不是应该符合Ti的XDM标准？ 2165