| 查看: 757 | 回复: 1 | ||
我不懂我不懂金虫 (初入文坛)
山大王
|
[求助]
tayltap工具箱所用算法
|
|
Matlab的工具箱tayltap，里面用到的是什么算法呀？？看不懂代码，代码如下： function [g,h] = tayltap(SLL,nbar,N) % ***************************************************************************** % * DBT, A Matlab Toolbox for Radar Signal Processing * % (c) FOA 1994-2000. See the file dbtright.m for copyright notice. % % function [g,h] = tayltap(SLL,nbar,N) % % Computes a Taylor taper for a specified sidelobe level % SLL>0 and no of elements = N % The function has been verified by comparing the results % with tables and plots in the following literature. % [R.J.Mailloux, "Phased Array Antenna Handbook", 1994, fig. 3.6] % [D.K.Barton, H.R.Ward,"Handbook of Radar Measurement,1984, % fig. A.14] % [M. Skolnik, "Radar Handbook", 2nd ed.,1990,table 10.9] % % % Start : 9xxxxx Fredrik Athley (freath). % Latest change: $Date: 2000/10/16 15:22:14 $ $Author: svabj $. % $Revision: 1.7 $ % ***************************************************************************** A = 1/pi*acosh(10^(SLL/20)); sig2 = (nbar^2)/(A^2+(nbar-.5)^2); F = zeros(1,nbar-1); nbar_min = ceil(2*A^2+0.5); if nbar<nbar_min nbar = nbar_min; sig2 = (nbar^2)/(A^2+(nbar-.5)^2); disp('The specified nbar is too small.') disp(['Using nbar = ',num2str(nbar)]) end for m=1:nbar-1 denprod=1; numprod=1; for n=1:nbar-1 numprod=numprod*(1-(m^2)/(sig2*(A^2+(n-0.5)^2))); if m~=n denprod=denprod*(1-(m/n)^2); end end F(m) = (-1)^(m+1)*numprod/(2*denprod); end % The Taylor distribution function must be sampled at points one-half % element spacing from the end of the distribution function x = 1/N*(-(N-1)/2  N-1)/2);g = zeros(1,N); num2 = 0; for m=1:nbar-1 num2 = num2 + F(m)*cos(2*pi*m*x); end g = (1+2*num2)/(1+2*sum(F)); % h=abs(fft(g,2^10)); % h = fftshift(h); % h = 20*log10(h/max(h)); % plot(h);grid % axis([0 1024 -SLL-20 5]) |
回复此楼
我不懂我不懂
金虫 (初入文坛)
山大王
- 应助: 1 (幼儿园)
- 金币: 1340
- 散金: 10
- 帖子: 47
- 在线: 26.9小时
- 虫号: 2029124
- 注册: 2012-09-25
- 性别: GG
- 专业: 通信理论与系统
2楼2015-09-22 21:19:36
10