取样光纤光栅传输函数源程序

《取样光纤光栅传输函数源程序》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、源程序：#include#include#include#definepi3.14159265#defineN18#defineN28typedefstd::complexcomplex;/*----二维基2FFT算法----*/voidfft2(complexf[N1][N2]){inti,j,k,m,n,l,r,M;intla,lb,lc;complext;/*----逐行进行一维FFT变换----*/for(i=0;i

2、*//*----计算分解的级数M=log2(N2)----*/for(m=N2,M=1;(m=m/2)!=1;M++);/*----按照倒位序重新排列原信号----*/for(j=1,k=N2/2;j<=N2-2;j++){if(j

3、示每个蝶形单元上下节点之间的距离/*----蝶形运算----*/for(l=1;l<=lb;l++){r=(l-1)*pow(2,M-m);for(n=l-1;n

4、像矩阵第j列数据进行FFT变换----*//*----计算分解的级数M=log2(N1)----*/for(m=N1,M=1;(m=m/2)!=1;M++);/*----按照倒位序重新排列原信号----*/for(i=1,k=N1/2;i<=N1-2;i++){if(i

5、级每个分组所含蝶形单元数//同时它也表示每个蝶形单元上下节点之间的距离/*----蝶形运算----*/for(l=1;l<=lb;l++){r=(l-1)*pow(2,M-m);for(n=l-1;n

6、[N2]){inti,j;for(i=0;i

7、*/voidmain(){inti,j;/*----原始图像矩阵数据x[N1][N2]----*/doublex[N1][N2]={{0,0,0,0,0,0,0,0},{0,0,0,1,1,0,0,0},{0,0,2,4,4,2,0,0},{0,1,4,8,8,4,1,0},{0,1,4,8,8,4,1,0},{0,0,2,4,4,2,0,0},{0,0,0,1,1,0,0,0},{0,0,0,0