数字信号处理三级项目案例.pptxVIP

;目录;CONTENTS;谱分析过程;利用Tubor C绘制时域采样序列图及频谱图;信号的频谱分布特性;高密度谱;;CONTENTS;;01;常用的窗函数的幅频特性;窗函数设计FIR低通滤波器的频率特性;不同窗函数对滤波性能的影响;CONTENTS;之;各滤波器截止频率选择;矩形窗——低通高通;矩形窗——带阻带通;其他窗函数的低通滤波器;其他窗函数的高通滤波器;其他窗函数的带阻滤波器;其他窗函数的带通滤波器;之;滤波前的时域波形;低通滤波器;带通滤波器;带阻滤波器;高通滤波器;CONTENTS;用DFT对信号进行频谱分析时，f1=201hz,f2=208hz, f3=214hz,采样频率fs=1000hz,观测时间为Tp=0.1s 时的时域和频域波形：;;当采样点N=32,并进行32点DFT，经过计算k=5,27的谱线代表f=40Hz； k=10,22的谱线代表f=80Hz；但对于f=60Hz，对应的k=7.5,是否该频率不能显示出来呢？;;?;利用快速卷积对信号进行滤波需注意：;FIR数字滤波器具有线性相位的充要条件是 偶对称 即h(n)=h(N-1-n)或 奇对称 即h(n)=-h(N-1-n);采样点为20点（整周期）;;参考文献;THANKS!;附录——程序代码;;#includecomplex //#includedrawfunction.h using namespace std; #define pi 3.1415926 double wc=0.408*pi;//204HZ=0.408，211HZ=0.422 int N=475; void Convolution(float *Tab_A,float *Tab_B,float *cout,int la,int lb) { int ii=0,jj=0; int k=0; float out[10000]={0}; for( int p=0;p=la+lb-2;p++) //卷积后的序列号 { for(int i=0;i=la-1;i++) //序列1的序号 for( int j=0;j=lb-1;j++) //序列2的序号 { if((i+j)==p) out[k]=out[k]+Tab_A[i+abs(ii)]*Tab_B[j+abs(jj)]; } //printf(%d\n,out[k]); //输出卷积后的序列 cout[k]=out[k]; k=k+1; } } //////////////*******************************DFT函数**************************//////////////////// //////////////****************a为输入序列x(n)输入可为复数序列，c记录得到的|X(K)|,c输入必须为0***********//// /////////////****************re记录X(K)的实部，im记录X(K)的虚部，N为周期的大小************//// void DFT(float *a,float *b, float *c,float *re,float *im, int N) { complexfloat po[10000]; float op[10000]; int n; float s = float(-(2 * pi) /(float)N); complexfloat tm; for(n=0;nN;n++) { po[n]=complexfloat(a[n], b[n]);//使得输入序列为复数 } for(n=0;nN;n++) op[n]=1; for (int k = 0; k N; k++) { tm = complexfloat(c[k], 0); for ( n = 0; n N; n++) { tm += po[n]*polar(op[n],s*k*n); } re[k]=real(tm); im[k]=imag(tm);