《数字信号实验报告.docxVIP

数字信号实验报告姓名：学号：专业： 2012·12·27实验一:数字离散信号实验原理：采用一定的时间间隔，对连续信号进行抽样，得到离散信号，即离散序列。根据离散序列运算规则，对离散信号进行反褶及求和求差运算。实验任务：生成离散信号并计算其振幅谱；数字离散序列的运算。 实现程序段：生成离散信号程序段：void signal(float *h, float f, float dt, float m, int N){int i;float t,a;for(i=0;iN;i++){h[i]=0;t=i*dt;a=(float)(2*f*f*log(m));h[i]=(float)((exp(-(a)*t*t))*sin(2*3.14159*f*t));}}计算序列程序段：float h1[64], h2[64],h3[64], h4[64];{for(i=0;iN*2;i++){h1[i]=0;h2[i]=0;h3[i]=0;h4[i]=0;}for(i=0;iN;i++){h[i+32]=h[i]; h[32-i]=h[i]; }for(i=0;iN*2;i++){h3[i]=h1[i]+h2[i]; h4[i]=h1[i]-h2[i];}}实验结果： 采样间隔 = 0.004秒频谱图 采样间隔 = 0.009秒频谱图 采样间隔 = 0.011秒频谱图四种采样间隔的信号图 信号h(n)图 信号h(-n)图 信号h(n)+h(-n)图 信号h(n)-h(-n)图实验分析：采样信号越小，信号恢复得越好，恢复后的信号失真越小。采样间隔过大，导致信号失真，恢复出来在波峰波谷处表现最明显。实验二：卷积及相关运算实验原理：根据卷积及相关运算的定义，用计算机编程实现离散序列的卷积及相关运算。实验任务：计算序列h(n)和x(n)的卷积；计算序列y1(n)和y2(n)的相关函数。实现程序段：计算卷积的程序段：void convolution(float *y, float *h, int N, float *x, int L){int i,j;for(i=0;iL;i++){ y[i]=0; for(j=0;jN;j++) { if( (i-j)=0) { y[i]+=h[j]*x[i-j]; }}}}计算相关函数的程序段：void relation(float *r, int M, float *y1, float *y2, int L){int i,j;for(i=0;iL;i++){r[i]=0;for(j=0;jM;j++){if(i-j=0)r[i]+=y1[j]*y2[j-i];}}}实现结果： 卷积结果图 相关结果图 实验分析： 相关与卷积的运算是类似的，基本运算过程都包括平移、相乘、求和。二者的差别是相关不对系统函数进行反褶，卷积则需要对系统函数进行反褶。实验三：离散信号振幅谱和相位谱的计算实验原理：通过傅里叶变换，把时间域的信号变换成幅频域的信号，研究信号在复频域的特征，即振幅谱及相位谱。实验任务：编写计算DFT的程序；编写计算振幅谱和相位谱。实现程序段：DFT计算程序段：void DFT(float *X, float *Xr, float *Xi, int N, float s){int i,j;for(i=0;iN;i++){ Xr[i]=0; Xi[i]=0;}float Pi=3.14159;if(s==1)for(i=0;iN;i++){for(j=0;jN;j++){Xr[i]+=X[j]*cos(j*i*2*Pi/N);Xi[i]+=X[j]*sin(j*i*2*Pi/N);}}if(s==-1){for(i=0;iN;i++){for(j=0;jN;j++){X[i]+=Xr[j]*cos(i*j*2*Pi/float(N))-Xi[j]*sin(i*j*2*Pi/float(N));}}}}振幅谱和相位谱计算程序段：int i;for(i=0;iN;i++) {Amp[i]=sqrt(Xr[i]*Xr[i]+Xi[i]*Xi[i]);}for(i=0;iN;i++){ phi[i]=atan(Xi[i]/Xr[i]); }实验结果： 信号图 振幅谱图 相位谱图实验分析：快速傅里叶变换（FFT）相比于离散傅里叶变换（DFT）,在数据量小的时候，基本没有多大差别。但到数据量很大的时候，就能明显地感觉到FFT的速度远快于DFT，并且这种优势随着数理量的增大而越发突出。实验四：FIR滤波器设计和滤波处理实验原理：根据所要设计的滤波器的系统特性，即所满足的系统因子，设计出相应的滤波器，并用编程实现。实验任务：编写滤波处理的程序片断;设计高通及低通滤波器。实现程序段：频率域的滤波处理程序段：for(k=0;