課程設(shè)計---信號分析與處理c語言編程

上傳人：奔*** IP屬地：河北更新時間：2024-04-03 格式：doc 頁數(shù)：24 大?。?.35MB 人氣指數(shù)：12 舉報 版權(quán)申訴

已閱讀1頁，還剩23頁未讀，繼續(xù)免費(fèi)閱讀

版權(quán)說明：本文檔由用戶提供并上傳，收益歸屬內(nèi)容提供方，若內(nèi)容存在侵權(quán)，請進(jìn)行舉報或認(rèn)領(lǐng)

文檔簡介

1、　　勘查技術(shù)課程設(shè)計：信號分析與處理基礎(chǔ)　　對于勘查技術(shù)與工程專業(yè)的學(xué)生來說，《信號分析與處理基礎(chǔ)》是一門專業(yè)基礎(chǔ)課，我是2010級的，我們是在大三第一學(xué)期上的，這門課數(shù)學(xué)與物理知識要求比較高，不過一開認(rèn)真仔細(xì)學(xué)的話，也會學(xué)的好的，起碼要比那空洞、生奧、蛋疼《彈性波動力學(xué)》好學(xué)些。　　隨著課程的結(jié)束，《信號分析與處理基礎(chǔ)》的課程設(shè)計也隨

2、之而來，我們是老師布置了4 個題目，分單號與雙號各自做2道，我是單號，做的是濾波與相關(guān)。　　這次課程設(shè)計，注意考驗(yàn)大家的編程能力，目前我們學(xué)過得就只有C語言，可以用Fortran,Matlab等等，Matlab可以現(xiàn)學(xué)現(xiàn)用，上手快。但是大家也可以挑戰(zhàn)下自己的C語言，提高下自己的編程能力，這是一次很好地機(jī)會，真正實(shí)用的時刻。我就是用C語言編的。　　其余2題

3、，我也把程序與結(jié)果圖收集到了這里，以供學(xué)弟、學(xué)妹們參考之用！我的QQ：593066480，有什么不懂的，或者好的見教，歡迎來信息交流！　　題目如下：　　濾波　　已知原始地震記錄x(t),要求：設(shè)計濾波器，消除x(t)中10Hz以下，80H

4、z以上的干擾信號。　　建議參數(shù)： A1=1，A2=0.8，A3=0.5　　f1=25Hz，f2=45Hz，f3=5Hz，f4=80　　取樣點(diǎn)數(shù)：N=200　　抽樣間隔：Δ=0.004，ß=100　　時域

5、濾波：　　由 (h(t)為濾波因子)　　建議參數(shù)：第一參數(shù)：ff1=20Hz，ff2=50Hz　　第二參數(shù)：ff1=25Hz，ff2=45Hz　　抽樣點(diǎn)數(shù)：M=60

6、抽樣間隔：Δ=0.004　　，　　，　　單獨(dú)計算： 　　要求：畫出x(t),h(t),y(t)圖形，為了分析方便，也可以畫出有效波s(t),干擾波n(t)及

7、其頻譜進(jìn)行分析，如下圖: 　　最后就是答辯，老師問問題，學(xué)生回答。主要注意幾點(diǎn)就行了　　熟悉課本濾波部分知識；　　第二參數(shù)要比第一參數(shù)濾波效果好，因?yàn)殚T第一參數(shù)開大了，進(jìn)來的干擾波也多了，從第一參數(shù)：Fy、第二參數(shù)Fy圖形上可以看出來，干擾波頻譜被壓小了。　　第二

8、參數(shù)壓制了干擾波，突顯了有效波，所有好。　　頻域?yàn)V波　　由公式：　　對x(t)進(jìn)行補(bǔ)0，28或者29總之必須是2的次方，因?yàn)橐玫紽FT公式與IFFT公式，進(jìn)行FFT變換得到X（k）；　　H（k）

9、的求法：H（k）也必須與X（k）點(diǎn)數(shù)相同　　，　　單門：20, 50　　雙門：20, 50,=N-,=N-,　　3, 在用IFFT反變換得y(t),存在實(shí)部與虛部，需要分析與處理

10、;　　要求：畫出H(k)、X(k)、Y(k)圖形，并且分析X(k)、Y(k)的區(qū)別，還有開單雙門的區(qū)別與差異，時域與頻域?yàn)V波誰好、為什么？　　分析：從單雙門實(shí)部圖形看出，與有效波是完全一樣的，但是幅值變大，這體現(xiàn)了濾波突顯有效波的特性；從圖形看出虛部對結(jié)果沒有影響；　　單雙門虛部完全不同，這是由于開單雙門效果不同引起的。單門虛部變化

11、大，而且幅值與起伏變化也大，而雙門幅值很小，小到可以忽略不計。按理說反變化IFFT后應(yīng)該只有實(shí)部，沒有虛部，至于為什么會產(chǎn)生虛部，希望讀者自己下去研究下，希望大家相互交流、交流！　　第一個注意問題：單門與雙門誰好？答案當(dāng)然是：雙門好。因?yàn)樵加行Р▁(t)是實(shí)數(shù)信號，對應(yīng)的頻譜是偶對稱的，單門時，只是讓一部分通過了；而雙門則全部通過，肯定效果要好！而且實(shí)部波形，明顯看出雙門是干擾部分起伏比單門時

12、要平緩，小得多了。　　第二個注意問題：時域與頻域哪個好？頻域要好，首先從兩者波形上可以看出，其次就是頻域?yàn)V波，只讓有效波通過，而之外的完全被濾掉，可謂是真正的理想狀態(tài)。　　第三個注意問題：頻域?yàn)V波y(t)實(shí)部圖200點(diǎn)以后，為什么有波形起伏？因?yàn)橛捎跁r域離散，必將導(dǎo)致頻域周期化，這是由于頻域周期化的結(jié)果。<

13、b>　　相關(guān)　　建議參數(shù)：1，30Hz，抽樣點(diǎn)數(shù)：M=100　　0.8，50Hz，抽樣點(diǎn)數(shù)：M=150　　抽樣間隔：Δ=0.004s　　要求：畫出，，的圖形并分析驗(yàn)證書上自互相關(guān)性質(zhì)的正確性！<b&

14、gt;　　快速褶積　　建議參數(shù)：1，25Hz，抽樣點(diǎn)數(shù)：M=250　　0.7，55Hz，抽樣點(diǎn)數(shù)：M=200　　抽樣間隔：Δ=0.004s　　地震記錄的生成和頻譜分析　　地震記錄的生成和頻譜分析、信噪比計算以及

15、補(bǔ)0對頻譜的影響，給定地震子波的數(shù)學(xué)表達(dá)式：　　和反射系數(shù)序列：0.2，0.4，0.15，0.5　　0.35， 0.1， 0.2　　產(chǎn)生一個含有隨機(jī)干擾信號的地震信號：　?。ㄆ渲衝(t)可以用-0.4—+0.4之間的隨機(jī)數(shù)代替）&

16、lt;b>　　要求：　　制作合成地震記錄x(t),并對b(t)和s(t)做頻譜分析（用FFT），其中：b(t)（N=41,Δ=0.004s，f=35Hz，α=100，A1=1）反射序列和隨機(jī)干擾N=200；　　計算x(t)信噪比，改變n(t)(用-0.8—+0.8之間的隨機(jī)數(shù)代替)的大小再計算。<p

17、>　　信噪比：　　分別對b(t)后邊、前邊、中間補(bǔ)0，計算補(bǔ)0前后頻譜的變化及補(bǔ)0多少對頻譜的影響。　　附帶程序：　　最麻煩的就是編程序了，開始的時候，是很麻煩，不過只有去啃，就一定會有收獲，比如：我開始編褶積程序的時候，整理了幾天

18、，上網(wǎng)查，翻圖書，后來突然明白，靠公式就可以編出來了。只是FFT需要些功夫,其余都很快！書上介紹的是基2FFT，這個代碼網(wǎng)上到處都有，想提升自己的編程能力，就去嘗試編基4FFT，以及考慮基nFFT吧，祝大家好運(yùn)。　　1、時域?yàn)V波程序代碼　　#include<stdio.h>　　#include<

19、math.h>　　#include"conv.cpp"　　#include"dft.cpp"　　#include"gyh.cpp"　　#define T 200//************T表示x(t)點(diǎn)數(shù)

20、，R表示n(t)點(diǎn)數(shù)　　#define R 60　　void main()　　{　　int i,j,B=100;//*********************B表示貝塔的值

21、　FILE *fp;　　double A1=1.0,A2=0.8,A3=0.5;//********產(chǎn)生x(t)　　double f1=25,f2=45,f3=5,f4=80;//******頻率取值　　double data=0.004;//******************取樣點(diǎn)數(shù)<

22、;/p>　　double s[T],n[T], x[T],h[R],y[T+R-1];　　double si[T]={0},s1[T]={0},s1i[T]={0},Fs[T]={0};　　double ni[T]={0},n1[T]={0},n1i[T]={0},Fn[T]={0};

23、　　double xi[T]={0},x1[T]={0},x1i[T]={0},Fx[T]={0};　　double hi[T]={0},h1[T]={0},h1i[T]={0},Fh[T]={0};　　double y1[T+R-1]={0},yb[T+R-1]={0},yb1[T+R-1]={0},Fy[T+R-1]={0};

24、　　for(i=0;i<T;i++)　　{　　s[i]=A1*exp(-B*pow((i*data),2)) * sin(2*PI*f1*i*data) + A2*exp(-B*pow((i*data),2)) * sin(2*PI*f2*i*data);<p&g

25、t;　　n[i]=A3*(sin(2*PI*f3*i*data)+cos(2*PI*f4*i*data));　　x[i]=s[i]+n[i];　　}　　//對s(t)進(jìn)行頻譜分析，用DFT　　dft(s,si,s1,s1i,Fs,T,

26、1);　　dft(n,ni,n1,n1i,Fn,T,1);　　dft(x,ni,x1,x1i,Fx,T,1);　　//統(tǒng)一導(dǎo)出　　fp=fopen("D:\\sh\\time1\\snxFsFnFx1.txt",&

27、quot;w");　　for(i=0;i<T;i++)　　fprintf(fp,"%f\t%f\t%f\t%f\t%f\t%f\n",s[i],n[i],x[i],Fs[i],Fn[i],Fx[i]);　　fclose(fp);<

28、b>　　//產(chǎn)生h(t)　　float ff1=20,ff2=50;//*****************可以修改的h(t)參數(shù)　　double w1,w2,dataw,w0;　　w1=2*PI*ff1,w2=2*PI*ff2;　　dataw=(

29、w2-w1)/2,w0=(w2+w1)/2;　　h[0]=2*dataw/PI;　　for(j=1;j<R;j++)　　{　　h[j]=(2.0/(PI*j*data))*cos(w0*j*data)*sin(dataw*j*d

30、ata);　　}　　dft(h,hi,h1,h1i,Fh,R,1);　　gyh(Fh,R);　　fp=fopen("D:\\sh\\time1\\hFh.txt","w");

31、　　for(j=0;j<R;j++)　　fprintf(fp,"%f\t%f\n",h[j],Fh[j]);　　fclose(fp);　　//褶積濾波得到y(tǒng)(t)　　con(x,h,y,T,R);&l

32、t;p>　　//對y(t)作傅里葉變換　　dft(y,y1,yb,yb1,Fy,T+R-1,1);　　//導(dǎo)出y及頻譜Y(k)　　fp=fopen("D:\\sh\\time1\\yFy1.txt","w");　　for(i=

33、0;i<T+R-1;i++)　　fprintf(fp,"%f\t%f\n",y[i],Fy[i]);　　fclose(fp); 　　printf("\nover\n\n");}　　頻域?yàn)V波程序代碼：</

34、b>　　#include<stdio.h>　　#include<math.h>　　#include"fft.cpp"　　#include"ifft.cpp"

35、#define G 256　　void main()　　{　　int i,B=100;　　FILE *fp;//定義指針文件　　double A1=1.0,A2=0.8,A3=0.5;//產(chǎn)生

36、x(t)　　double f1=25,f2=45,f3=5,f4=80;//參數(shù)　　double data=0.004;//抽樣間隔　　double s[200],n[200];　　double x[G]={0},x1[G]={0},h[G]={0},h1[G]={0},

37、y[G]={0},y1[G]={0};　　double Fx[G]={0},Fh1[G]={0},Fh2[G]={0},Fy1[G]={0},Fy2[G]={0};　　double Y1[G]={0},Y1i[G]={0},Y2[G]={0},Y2i[G]={0};　　for(i=0;i<200;i++)&

38、lt;/p>　　{　　s[i]=A1*exp(-B*pow((i*data),2)) * sin(2*PI*f1*i*data) + A2*exp(-B*pow((i*data),2)) * sin(2*PI*f2*i*data);　　n[i]=A3*(sin(2*PI*f3*i*data

39、)+cos(2*PI*f4*i*data));　　x[i]=s[i]+n[i];　　}　　//導(dǎo)出x[t]　　fp=fopen("D:\\sh\\py\\x.txt",&quo

40、t;w");　　for(i=0;i<G;i++)　　fprintf(fp,"%f\n",x[i]);　　fclose(fp);　　//x[t]頻譜計算　　fft(x,x1,Fx,G,1);<

41、;/p>　　//導(dǎo)出X(k)　　fp=fopen("D:\\sh\\py\\Fx.txt","w");　　for(i=0;i<G;i++)　　fprintf(fp,"%f\n",Fx

42、[i]);　　fclose(fp);　　//處理h[t]　　double dataf;　　double m1,m2,m3,m4;　　dataf=1.0/(data*G);

43、　　m1=20.0/dataf,m2=50.0/dataf;　　m3=G-m2,m4=G-m1;　　//計算,導(dǎo)出單門　　for(i=0;i<G;i++)　　{</

44、p>　　if(i>=int(m1)&&i<=int(m2))　　Fh1[i]=1;　　else Fh1[i]=0;　　}

45、　//計算,導(dǎo)出雙門　　for(i=0;i<G;i++)　　{　　if((i>=int(m1)&&i<=int(m2))||(i>=int(m3)&&i<=int(m4)))

46、　　Fh2[i]=1;　　else Fh2[i]=0;　　}　　fp=fopen("D:\\sh\\py\\Fh12.txt","w");<p&

47、gt;　　for(i=0;i<256;i++)　　fprintf(fp,"%f\t%f\n",Fh1[i],Fh2[i]);　　fclose(fp);　　for(i=0;i<G;i++)//單門濾波　　{

48、;　　Y1[i]=x[i]*Fh1[i];　　Y1i[i]=x1[i]*Fh1[i];　　Fy1[i]=sqrt(pow(Y1[i],2)+pow(Y1i[i],2));　　}　　for(i=

49、0;i<G;i++)//開雙門　　{　　Y2[i]=x[i]*Fh2[i];　　Y2i[i]=x1[i]*Fh2[i];　　Fy2[i]=sqrt(pow(Y2[i],2)+pow(Y2i[i],2));}

50、　　fp=fopen("D:\\sh\\py\\Fy12.txt","w");//導(dǎo)出單雙門　　for(i=0;i<G;i++)　　fprintf(fp,"%f\t%f\n",Fy1[i],Fy2[i]);　　fclose(fp

51、);　　//Y[k]作反變換IFFT,并且導(dǎo)出實(shí)部，虛步　　ifft(Y1,Y1i,G,-1);　　ifft(Y2,Y2i,G,-1);　　fp=fopen("D:\\sh\\py\\yt12.txt","w");

52、;　　for(i=0;i<G;i++)　　fprintf(fp,"%f\t%f\t%f\t%f\n",Y1[i],Y1i[i],Y2[i],Y2i[i]);　　fclose(fp);　　printf("\nover\n\n");}</p&g

53、t;　　相關(guān)程序代碼　　#include<stdio.h>　　#include<math.h>　　#include"correl.cpp"　　#define PI 3.

54、14159265　　void main()　　{　　int i,j;　　double A1=1.0,A2=0.8;　　double f1=30,f2=50

55、;　　double data=0.004,p=100;　　double x[100]={0},y[150]={0},xy[249]={0},yx[249]={0};　　double xx[199],yy[299],xyft[249];　　for(i=0;i<100;

56、i++)　　x[i]=A1*sin(2*PI*f1*i*data);　　for(j=0;j<150;j++)//計算y(t)　　y[j]=A2 * exp(-p*j*data*j*data)*sin(2*PI*f2*j*data);　　correl(x,y,100,

57、150,xy);　　correl(x,x,100,100,xx);　　correl(y,y,150,150,yy);　　correl(y,x,150,100,yx);　　FILE *fp,*fp1,*fp2,*fp3,*fp4,*fp5,*fp6;//導(dǎo)出數(shù)據(jù)</p

58、>　　int i1,j1,k,l,m,n;　　fp1=fopen("D:\\sh\\data\\x.txt","w");　　fp2=fopen("D:\\sh\\data\\y.txt","w");　　fp

59、3=fopen("D:\\sh\\data\\xy.txt","w");　　fp4=fopen("D:\\sh\\data\\xx.txt","w");　　fp5=fopen("D:\\sh\\data\\yx.txt","w");&l

60、t;/p>　　fp6=fopen("D:\\sh\\data\\yy.txt","w");　　for(i=0;i<249;i++)　　xyft[i]=yx[i];　　fp=fopen("D:\\sh\\data\\xyt.txt&

61、quot;,"w");　　for(i=0;i<249;i++)　　fprintf(fp,"%f\n",xyft[i]);　　fclose(fp);　　for(i1=0;i1<100;i1++)&l

62、t;p>　　fprintf(fp1,"%f\n",x[i1]);　　fclose(fp1);　　for(j1=0;j1<150;j1++)　　fprintf(fp2,"%f\n",y[j1]);　　fclose(fp2)

63、;　　for(k=0;k<249;k++)　　fprintf(fp3,"%f\n",xy[k]);　　fclose(fp3);　　for(l=0;l<199;l++)　　fprintf(fp4,&q

64、uot;%f\n",xx[l]);　　fclose(fp4);　　for(m=0;m<249;m++)　　fprintf(fp5,"%f\n",yx[m]);　　fclose(fp5);　　fo

65、r(n=0;n<299;n++)　　fprintf(fp6,"%f\n",yy[n]);　　fclose(fp6);}　　4、地震子波及頻譜分析　　#include<stdio.h>　　#in

66、clude"conv.cpp"　　#include"fft.cpp"　　#include"uni.cpp"　　#define V 64//*********宏定義的參數(shù)，可以修改的點(diǎn)數(shù)　　#define W 25

67、6　　void main()　　{　　int i;　　FILE *fp;　　double A1=1,f=35,a=1

68、00;//********************a代表阿爾法α　　double data=0.004;//************************* 抽樣間隔　　double b[V]={0},bi[V]={0},Fb[V]={0},ks[200]={0};　　double sz[W]={0},s[20

69、0]={0},szi[W]={0},Fsz[W]={0};　　double x[W]={0},xi[W]={0},Fx[W]={0},n[200]={0};　　double xwp[V];　　for(i=0;i<41;i++)　　b[i]=A1*exp(-a*i*d

70、ata)*sin(2*PI*f*i*data);　　//*****產(chǎn)生ξ(t)　　ks[29]=0.2,ks[64]=0.4,ks[80]=0.15,ks[102]=0.5,ks[114]=0.35,ks[145]=0.1,ks[156]=0.2;　　fp=fopen("D:\\da\\b.txt&qu

71、ot;,"w");　　for(i=0;i<V;i++)　　fprintf(fp,"%f\n",b[i]);　　fclose(fp);　　fp=fopen("D:\\da\\ks.txt","w&q

72、uot;);　　for(i=0;i<200;i++)　　fprintf(fp,"%f\n",ks[i]);　　fclose(fp);　　con(b,ks,sz,41,200);　　fp=fopen(&qu

73、ot;D:\\da\\s.txt","w");　　for(i=0;i<W;i++)　　fprintf(fp,"%f\n",sz[i]);　　fclose(fp);　　for(i=0;i<200;i++)//**

74、*****************甩掉后40個樣點(diǎn)　　s[i]=sz[i];　　uni(-0.4,0.4,200,n);//*****************產(chǎn)生隨機(jī)數(shù)　　fp=fopen("D:\\da\\n.txt","w");

75、　　for(i=0;i<200;i++)　　fprintf(fp,"%f\n",n[i]);　　fclose(fp);　　for(i=0;i<200;i++)　　x[i]=s[i]+n[i];　　f

76、p=fopen("D:\\da\\x.txt","w");　　for(i=0;i<W;i++)　　fprintf(fp,"%f\n",x[i]);　　fclose(fp);　　fft(sz,szi,Fsz,

77、W,1);//**********　　fft(b,bi,Fb,V,1);//************************用fft對b[t],x[t]頻譜分析　　fft(x,xi,Fx,W,1);　　fp=fopen("D:\\da\\Fs.txt","w");</

78、p>　　for(i=0;i<W;i++)　　fprintf(fp,"%f\n",Fsz[i]);　　fclose(fp);　　fp=fopen("D:\\da\\Fb.txt","w");<

79、p>　　for(i=0;i<V;i++)　　fprintf(fp,"%f\n",Fb[i]);　　fclose(fp);　　fp=fopen("D:\\da\\Fx.txt","w");　　for(i

80、=0;i<W;i++)　　fprintf(fp,"%f\n",Fx[i]);　　fclose(fp);　　for(i=0;i<V;i++)　　xwp[i]=atan(bi[i]/b[i]);　　fp

81、=fopen("D:\\da\\xw.txt","w");　　for(i=0;i<V;i++)　　fprintf(fp,"%f\n",xwp[i]);　　fclose(fp);}　　5、信

82、噪比　　#include<stdio.h>　　#include<math.h>　　#include"conv.cpp"　　#include"uni.cpp"<

83、p>　　#define PI 3.14159265　　#define V 64　　#define W 256　　void main()　　{　　int i;</b

84、>　　double A1=1,f=35,a=100;　　double data=0.004,Sr,tps=0,tpn=0;　　double b[V]={0},sz[240]={0},s[200]={0},ks[200]={0};　　double x[W]={0},n

85、[200]={0};　　for(i=0;i<41;i++)　　b[i]=A1*exp(-a*i*data)*sin(2*PI*f*i*data);　　ks[29]=0.2,ks[64]=0.4,ks[80]=0.15;　　ks[102]=0.5,ks[114]=0.

86、35,ks[145]=0.1,ks[156]=0.2;　　con(b,ks,sz,41,200);　　for(i=0;i<200;i++)//甩掉后40個樣點(diǎn)　　s[i]=sz[i];　　uni(-0.4,0.4,200,n);//隨機(jī)數(shù)&l

87、t;p>　　for(i=0;i<200;i++)　　x[i]=s[i]+n[i];　　for(i=0;i<200;i++)//信噪比計算　　{　　tps+=pow(s[i],2);<p&

88、gt;　　tpn+=pow(n[i],2);　　}　　Sr=tps/tpn;　　printf("SNR1=%f\n",Sr);}　　補(bǔ)0影響<p&

89、gt;　　#include<stdio.h>　　#include"fft.cpp"　　#define V 64　　void main()　　{

90、　　int i;　　FILE *fp;　　double A1=1,f=35,a=100;　　double data=0.004;　　double b[V]={0},bi[V]={0},Fb[V]={0};</

91、p>　　for(i=0;i<41;i++)　　b[i]=A1*exp(-a*i*data)*sin(2*PI*f*i*data);　　fp=fopen("D:\\daa\\b1.txt","w");　　for(i=0;i<V;i++)

92、　　fprintf(fp,"%f\n",b[i]);　　fclose(fp);　　fft(b,bi,Fb,V,1);//**************用fft對b[t]頻譜分析　　fp=fopen("D:\\daa\\Fb1.txt"

93、,"w");　　for(i=0;i<V;i++)　　fprintf(fp,"%f\n",Fb[i]);　　fclose(fp);　　}　　《最后》

94、所用到的子程序　　相關(guān)：void correl(double x[],double y[],int N,int M,double z[])　　{　　int i,j,l,ts=0;//ts代表相關(guān)序列的第一個數(shù)的序列號　　double

95、t=0;　　l=1-M;　　for(i=-M+1;i<N;i++)　　{　　for(j=0;j<N;j++)　　if(j-i>=0 &a

96、mp;& (j-i)<M)　　t+=x[j]*y[j-i];　　z[i+M-1]=t;　　t=0;　　}　　printf("\

97、n**********\n");　　printf("ts=%d\n",l);　　printf("**********\n");}　　褶積：//傳入數(shù)組1,2以及存儲數(shù)組3；1,2的長度　　void con(double a

98、[],double b[],double c[],int M,int L)　　{　　int i,j,N;　　N=M+L-1;　　for(i=0;i<N;i++)<

99、;p>　　{　　double tp=0.0;　　for(j=0;j<M;j++)　　{　　if((i-j)>=0&&(i-j)<L)

100、　　tp+=a[j]*b[i-j];　　}　　c[i]=tp,tp=0.0;　　}　　}　　DFT

101、：//dft子函數(shù)，傳遞實(shí)部，虛部，存儲sx數(shù)組，振幅數(shù)組。點(diǎn)數(shù)，正反　　#include<math.h>　　#define PI 3.14159265　　void dft(double x[],double y[],double x1[],double y1[],double Fn[],int N,

102、int sign)　　{　　int i,j;　　double Q,D,C,c,s;　　Q=2*PI/N;　　for(

103、i=0;i<N;i++)　　{　　for(j=0;j<N;j++)　　{　　D=Q*i*j;<b&

104、gt;　　c=cos(D);　　s=sin(D)*sign;　　x1[i]+=c*x[j] + s*y[j];　　y1[i]+=c*y[j] - s*x[j];　　}<

105、;b>　　}　　if(sign==-1)　　{　　C=1.0/N;　　for(i=0;i<N;i++)

106、　{　　x1[i]=C*x1[i],y1[i]=C*y1[i];　　}}　　for(i=0;i<N;i++)　　Fn[i]=sqrt(pow(x1[i],2) + pow(y1[i],2));</p&g

107、t;　　}　　FFT://fft子函數(shù)，傳遞實(shí)部，虛部，存儲振幅數(shù)組，點(diǎn)數(shù)，正反　　#include<math.h>　　#define PI 3.14159265　　void fft(double x

108、[],double y[],double Fn[],int N,int sign)　　{　　int i,j,r,m1,m2,m3,m4,k1,k2,l,k;　　double u,v; //暫時存儲中間實(shí)部與虛部　　double

109、 c,s,tr,ti,C;　　l=(int)(log10(N)/log10(2));//計算l的值　　for(j=0,i=0;i<N-1;i++)//變址運(yùn)算，倒序，采用雷德SF　　{　　if(i<j)&

110、lt;/b>　　{　　tr=x[j];ti=y[j];x[j]=x[i];y[j]=y[i];x[i]=tr;y[i]=ti;　　}　　k=N/2;&l

111、t;/p>　　while(k<(j+1))　　{　　j=j-k;　　k=k/2;　　}

112、　　j=j+k;　　}　　for(i=1;i<=l;i++)//第一重循環(huán)，遞推迭代次數(shù)　　{　　m1=int

113、(pow(2,i-1));　　m2=2*m1;　　m3=N/m2;　　for(j=0;j<m3;j++)　　{

114、　　m4=j*m2;　　for(r=0;r<m1;r++)　　{　　k1=m4+r;　　k2=k1+m1;<

115、/p>　　c=cos(2*PI*r/m2);　　s=sign*sin(2*PI*r/m2);　　u=x[k2]*c + y[k2]*s;　　v=y[k2]*c - x[k2]*s;　　x[k2]=x[k1]-u;y[k2]=y[k1]-v;&l

116、t;/p>　　x[k1]=x[k1]+u;y[k1]=y[k1]+v;　　}　　}　　}　　if(sign==-1)</p

117、>　　{　　C=1.0/N;　　for(i=0;i<N;i++)　　{　　x[i]=C*x[i];&

118、lt;p>　　y[i]=C*y[i];　　}　　}　　for(i=0;i<N;i++)　　{　　Fn[i

119、]=sqrt(pow(x[i],2) + pow(y[i],2));　　} 　　}　　void gyh(double x[],int N)　　{</p

120、>　　int i,j;　　double t=0;　　for(i=0;i<N;i++)　　if(t<x[i])　　t=x[i];<

121、;p>　　for(j=0;j<N;j++)　　{x[j]=x[j]/t;}}　　隨機(jī)數(shù)：　　//隨機(jī)數(shù)子函數(shù),需要傳遞區(qū)間，產(chǎn)生個數(shù)，以及存儲數(shù)組　　#include<stdlib.h>

122、　　#include<time.h>　　void uni(double x,double y,int Z,double a[])　　{　　int i;　　unsig

123、ned seed=time(NULL);　　printf("seed=%d\n",seed);　　printf("\n");　　srand(seed);　　for(i=0;i<Z;i++)<p&g

眾賞文庫> 全部分類> 畢業(yè)設(shè)計

溫馨提示

1. 本站所有資源如無特殊說明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
3. 本站RAR壓縮包中若帶圖紙，網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽，若沒有圖紙預(yù)覽就沒有圖紙。
4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
5. 眾賞文庫僅提供信息存儲空間，僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理，對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯，并不能對任何下載內(nèi)容負(fù)責(zé)。
6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容，請與我們聯(lián)系，我們立即糾正。
7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

課程設(shè)計---信號分析與處理c語言編程

文檔簡介

溫馨提示

最新文檔

評論

課程設(shè)計---信號分析與處理c語言編程

文檔簡介

溫馨提示

最新文檔

評論

免費(fèi)下載