信号处理模块.docxVIP

实用标准文案 在今天我已经做信号处理模块差不多了， 感觉东西还是很乱， 唉，也不知道怎么搞的。我一定要好好理一理了， 我的界面设计也得马上开始了， 已经完全乱了计划，真想打自己，呜呜。 一．在信号部分 , 我得查查我做了什么。 sin 信号 double sinwn(double a,double f,double ph,double fs,int k) { double pi,result; pi=4.0*atan(1.0); f=2*pi*f/fs; ph=ph*pi/180; result=a*sin(k*f+ph); return result; } int main() { int i; double a ,f,fs,ph; double x[n]; a=2.0; f=5.0; ph=45.0; fs=150.0; for(i=0;i32;i++) { x[i]=sinwn(a,f,ph,fs,i); printf(%10.7f,x[i]); if(i%4==3) printf(\n); } return 0; } LFM信号（范围是 -T/2,T/2 ） double lfm(double a,double f0,double T,double F,double fs,int k) { double pi,y,f1,f2,f3; pi=4.0*atan(1.0); f1=2*pi*f0/fs; f2=F/T; f3=pi*f2*(1/fs)*(1/fs); y=a*cos(f1*k+f3*k*k); return y; 精彩文档 实用标准文案 } int main() { int i; double a,f0,T,F,fs; n=50; double x[n]; a=1.5; f0=10.0; T=50.0; F=200.0; fs=150.0; FILE *fp; for(i=0;i=n;i++) { x[i]=lfm(a,f0,T,F,fs,i); printf(%10.4f,x[i]); printf(\n); } fp=fopen(lfm1.txt,w+); for(i=0;i=n;i++) { fprintf(fp,%f ,x[i]); } fclose(fp); return 0; } 方波 ： 白噪声（高斯分布） ： 均匀分布随机信号 ： 双曲调频 （对数相位信号或者线性周期信号） ：（Hyperbolic FM signal 还没做，赶紧做啊） CW脉冲（相当于取一小段 0~T内的 sin 函数或者 cos 函数） 二．信号处理 1. 首先是重头戏 FFT: #include stdio.h #include math.h double fft(double x[],double y[],int n,int sign) //x[] 先 x(i) 放 fft 的实部 ,FFT 后放 X(k) 的实部；y[] 先 y(i) 放 fft 的实部 ,FFT 后放 Y(k) 的实部，n 为 2 的整数次方，sign 为 1 时做 FFT,sign 为-1 时，做 IFFT ； 精彩文档 实用标准文案 { int i,j,k,l,m,n1,n2; double c,c1,e,s,s1,t,tr,ti; j=1; for(i=1;i16;i++) // 做 n 点的 FFT，并且 N2的 16 次方 { m=i; j=j*2; if(j==n)break; } n1=n-1; j=0; for(i=0;in1;i++) // 进行码位倒置 { if(ij) { tr=x[j]; ti=y[j]; x[j]=x[i]; y[j]=y[i]; x[i]=tr; y[i]=ti; } k=n/2; while(kj+1) { j=j-k; k=k/2; } j=j+k; } n1=1; for(l=1;l=m;l++) //FFT { n1=2*n1; n2=n1/2; e=3.14159265359/n2; c=1.0; s=0.0; c1=cos(e); s1=-sign*sin(e); for(j=0;jn2;j++) { for(i=j;in;i+=n1) 精彩文档 实用标准文案 { k=i+n2; tr=c*x[k]-s*y[k]; ti=c*y[k]+s*x[k]; x[k]=x[i]-tr; y[k]=y[i]-ti; x[i]=x[i]+tr; y[i]=y[i]+ti; } t=c; c=c*c1-s*s1; s=t*s1+s*c1; } } if(sign==-1) //IFFT 的公式前面有 1/n { for(i=0;in;i++) { x[i]/=n; y[i]/=n; } }