近场天线测量作业.docx

作业一 ：分别采用直接求和与快速Fourier变换（FFT）两种方法计算出F(?)，并与理论计算结果比较，并比较两种方法所用时间。

已知f(x)?e?x

求 F(?)?

??e?xej2??xdx

??

直接积分：F(?)? 2 (1-1)

1?(2??)2

当?很大时，F(?)?0

取??100时，F(?)?10?5?0

故近似认为当???时，F(?)?0，即可以近似认为f（x）是一个谱宽有限得函数，带限为2?，取?x??2???0.005，则由取样定理有

F(?)?

??

m???

e?m?x

?ej2??m?x??x

令??

nN?x，

N?1

F(?)? ?2

e?m?xej

2?mn

N ?x

令m?N?k,则有

2

m??N

2

n ?N?1 ?(k?N)?x

2?n N

j (k? )

F( )?

N?x

e

k?0

2 e N 2?x

?(?1)n

?N?1

k?0

N?1

?xe

?(k?N)?x2

2?kn

2?kn

je N

j

?(?1)n? fej

k

k?0

N (n?0,1, ,N?1) FFT形式 (1-2)

? N

??xe(k? )?x

?

近场天线测量大作业 2

k?0,1,?,N?1

2

??其中： f 2

?

?

(1-3)

k ??xe(N?k)?x

k?N

,?,N?1

?? 2 2

取N=2048，则

1024*0.005?5，

N?x

e2

e

??1

Matlab程序清单如下：clc

clearallWp=100;

dx=1/(2*Wp);N=2048;n=0:N-1;

w1=n./(N*dx);

%%%理论值w=0:.001:Wp

ticFP=2./(1+(2*pi*w).^2);toc

%%%直接求和tic

forn1=0:N-1

FS1=0;

fork=0:N-1;

FS2=dx.*exp(-abs((k-N/2).*dx)).*exp(i*2*pi*k*n1/N);FS1=FS2+FS1;

end

FS(1,n1+1)=(-1).^n1.*FS1;

endFS=abs(FS);

toc

%%%FFTk=0:N-1;

w_=n./(N*dx);

Fk=dx.*exp(-abs((k-N/2).*dx));tic

FF=fft(Fk,N);

FF=abs(FF);

toc

近场天线测量大作业 3

figure(1)subplot(1,2,1)

plot(w,FP,--,w_,FS,s);gridon;

title(直接求和计算结果与理论结果)legend(理论值,直接求和)

axis([0302]);subplot(1,2,2)

plot(w,FP,--,w_,FF,s);gridon;title(FFT计算结果与理论结果)legend(理论值,FFT计算)axis([0302])

figure(2)

plot(w,FP,--,w_,FS,s,w_,FF,o);gridon;legend(理论值,直接求和,FFT计算)axis([0302]);

近场天线测量大作业 4

Elapsedtimeis0.006355seconds. （理论值计算见式（1-1）Elapsedtimeis10.422329seconds.（直接求和计算见式（1-2）Elapsedtimeis0.001044seconds. （FFT见式（1-3）

4.结果与讨论

由计算结果图可以看出：用直接求和计算和FFT算法得到的结果均与理论结果吻合很好，几乎重合；由计算所用时间可以得出：FFT算法比直接求和法具有明显的优势，当N=1024

时，直接计算需要N2=1048576次乘法，然而FFT算法只需要1024(log21024) 5120次

2

乘法，算法次数减小自然能节约系统资源缩短计算时间，从而比直接求和法更实用。

近场天线测量大作业PAGE5

近场天线测量大作业

PAGE5

n作业二： 利用一维驻相法推导天线的远场方向函数与柱面波谱an(h)，b(h)的关系式。

n

由远区场的表达式：

E ??? ???a(h)M(4)(r)?b(h)N(4)(r)?dh

?n

n?????

n,h

n n,h ?