频谱分析仪模拟仿真

数字信号处理课程设计题目：频谱分析仪模拟仿真专业：电子信息工程班级：姓名：日期：2015-0 目录设计目的------------------------------------------------------------------------1程序流程------------------------------------------------------------------------2软件设计------------------------------------------------------------------------3心得体会------------------------------------------------------------------------4参考文献------------------------------------------------------------------------5N=32;for(i=0;i<N/2;i++){a=1;printf("%d",a);}for(i=N/2;i<N;i++){A=-1printf("%d",a);}system("pause");}软件截图2、DFT程序#include<stdio.h>#include<math.h>voidmain(){floatt,re=0,im=0,sum,mod,P=3.1415926;intn,k,N=32,*x;for(n=0;n<N/2;n++){x[n]=1;}for(n=N/2;n<N;n++){x[n]=-1;}for(k=0;k<N;k++){ for(n=0;n<N;n++) { t=(-2*P/N)*n*k; re+=x[n]*cos(t); } for(n=0;n<N;n++) { t=(-2*P/N)*n*k; im+=x[n]*sin(t); }sum=re*re+im*im;mod=sqrt(sum);printf("x[%d]=%f,\n",k,mod);re=0,im=0;sum=0;mod=0;}}软件截图3、倒位序程序#include<stdio.h>#include<math.h>#defineN32voidmain(){intnv2,nm1,i,k,j=0,x[100],t;nv2=N/2;nm1=N-1;printf("请输入自然序：\n");for(i=0;i<N;i++)scanf("%d",&x[i]);for(i=0;i<nm1;i++){ if(i<j) { t=x[j]; x[j]=x[i]; x[i]=t; } k=nv2; while(k<=j) { j=j-k; k=k/2; } j=j+k; }printf("输出倒位序：\n");for(j=0;j<N;j++)printf("%d",x[j]);printf("\n");}软件截图4、FFT程序#include<stdio.h>#include<math.h>structcmplex{ floatr; floati;};structcmplexss(structcmplexa1,structcmplexa2){ structcmplexc; c.r=a1.r*a2.r-a1.i*a2.i; c.i=a1.r*a2.i+a1.i*a2.r; return(c);}voidmain(){ inti,k,N=32,nv,nm,j=0,m,le,le0,ip; floatpai=3.1415926,tmp,sum,*a; structcmplexu,w,t,x[32]; for(i=0;i<N;i++) { x[i].r=1; x[i].i=0; } for(i=N/2;i<N;i++) { x[i].r=-1; x[i].i=0; } nv=N/2; nm=N-1; for(i=0;i<nm;i++) { if(i<j) { t=x[j]; x[j]=x[i]; x[i]=t; } k=nv; while(k<=j) { j-=k; k=k/2; } j+=k; } for(m=1;m<=5;m++) { le=pow(2,m); le0=le/2; u.r=1; u.i=0; tmp=pai/le0;w.r=cos(-tmp); w.i=sin(-tmp); for(j=0;j<le0;j++) { for(i=j;i<N;i+=le) { ip=i+le0; t=ss(x[ip],u); x[ip].r=x[i].r-t.r;x[ip].i=x[i].i-t.i; x[i].r+=t.r;x[i].i+=t.i; } u=ss(u,w); }} printf("\n"); for(i=0;i<N;i++) { sum=x[i].r*x[i].r+x[i].i*x[i].i; a[i]=sqrt(sum); printf("a[%d]=%f\n",i,a[i]); }}软件截图三、心得体会该实验巩固了我们上学期所学的数字信号处理和信号与系统知识，通过C语言上机操作仿真信号，不仅能够加深我们对各种信号抽样和变换的认识，也能让我们真真切切体会到计算机仿真系统带给我们学生和技术人员的巨大便利。总而言之