在《自適應天線與相控陣》這門課中，我了解到了關於理想低副瓣陣列設計的一些方法，其中切比雪夫等副瓣陣列設計方法是一種基礎的方法，故將其設計流程寫成maltab程序供以后學習使用。在此分享一下。
此方法全稱為道爾夫-切比雪夫綜合法，簡稱為切比雪夫綜合法，是一種工程實際中常用的可控制副瓣電平的陣列天線綜合方法。切比雪夫陣列的特點是：
（1）等副瓣電平；
（2）在相同副瓣電平和相同陣列長度下主瓣最窄，為最佳陣列；
（3）單元數過多時，陣列兩端單元激勵幅度跳變大，使饋電困難。
一般在雷達系統中，為了使其具有較高的抗干擾、抗反輻射導彈的能力，往往要求雷達天線的副瓣盡量低，而采用道爾夫-切比雪夫綜合法以及進一步的泰勒綜合法等設計的陣列天線就可以實現低副瓣。
最早，道爾夫（C.L.Dolph）利用切比雪夫函數來逼近天線陣列的陣因子函數，得到了這種嚴謹規范的綜合方法。
而且，經過前人研究，當天線單元N≤13時，切比雪夫陣列從中間到兩端的激勵分布是單調減小的；而當N>13時，陣列兩端單元的激勵開始出現跳變。所以對於大型陣列來說一般不宜采用切比雪夫方法綜合陣列。所以下面的Matlab程序正常工作在天線單元數N為3到13這個范圍內。
關於如何采用切比雪夫多項式去設計陣因子的具體技術步驟，另一位博主較為詳細地介紹了，此處不再贅述，只放該文鏈接：切比雪夫低副瓣陣列設計 MATLAB
下面是可以綜合設計天線單元從3到13單元的切比雪夫綜合法的Matlab程序：
%% --------------------------------------------------------------------------
% 切比雪夫低副瓣陣列綜合
% 設計一個間距為d,單元數為N，主副瓣電平比為RdB，掃描角度為theta0的切比雪夫陣列。
% 2019.11.10
%--------------------------------------------------------------------------
%% 初始數據賦值
clear
clc
N = 13; %單元數N（3<N<=13,N取整數）
if rem(N,2)==0 %求和項數M(奇偶不同)
M = N/2;
else
M = (N-1)/2+1;
end
RdB = 26; % 主副瓣比(dB值)
lamuda = 10; % 波長
d = 0.6*lamuda; % 單元間距
theta0 = 80/180*pi; % 掃描角度，相對於陣列排布方向的夾角
A = [1,0,0,0,0,0,0,0,0,0,0,0,0,0; % chebyshev多項式Tn(x) = cos(nu)= f(x)系數矩陣A
0,1,0,0,0,0,0,0,0,0,0,0,0,0; % 系數矩陣A每一行表示n，從n = 0開始
-1,0,2,0,0,0,0,0,0,0,0,0,0,0; % 列表示x的冪次方，從0次方開始
0,-3,0,4,0,0,0,0,0,0,0,0,0,0;
1,0,-8,0,8,0,0,0,0,0,0,0,0,0;
0,5,0,-20,0,16,0,0,0,0,0,0,0,0;
-1,0,18,0,-48,0,32,0,0,0,0,0,0,0;
0,-7,0,56,0,-112,0,64,0,0,0,0,0,0;
1,0,-32,0,160,0,-256,0,128,0,0,0,0,0;
0,9,0,-120,0,432,0,-576,0,256,0,0,0,0;
-1,0,50,0,-400,0,1120,0,-1280,0,512,0,0,0;
0,-11,0,220,0,-1232,0,2816,0,-2816,0,1024,0,0;
1,0,-72,0,840,0,-3584,0,6912,0,-6144,0,2048,0;
0,13,0,-364,0,2912,0,-9984,0,16640,0,-13312,0,4096];
% 初始矩陣賦值
I = zeros(1,M); % 電流幅度矩陣
S = zeros(M,M); % 陣因子系數矩陣
S_compare = zeros(1,M); % 系數比對矩陣
R = 10^(RdB/20); % 非dB 值的主副瓣比
x0 = 1/2*( (R+sqrt(R^2-1))^(1/(N-1))+...% 變量代換值x0
(R-sqrt(R^2-1))^(1/(N-1)) );
%% 求S、S_compare和I
% 從系數矩陣中擇選出M個求和項對應的系數S(奇偶分開討論)
for i = 1:M
if rem(N,2)==0 % 偶數情況
for j = 1:M % 第i行表示x的i次方，
S(i,j) = A(2*j,2*i); % 第j列表示第j個求和項系數（未除x0)
end
S_compare(i) = A(N,2*i); % 比對矩陣，即下標為N-1的chebyshev多項式的系數
else % 奇數情況
for j = 1:M
S(i,j) = A(2*j-1,2*i-1);
end
S_compare(i) = A(N,2*i-1);
end
end
% 通過S和S_compare系數比對求出電流幅度
for k = 1:M
i = M-k+1;
if rem(N,2)==0 % 偶數
I(i) = (S_compare(i)*x0^(2*i-1) -...
I*S(i,:)')/S(i,i);
else % 奇數
I(i) = (S_compare(i)*x0^(2*(i-1)) -...
I*S(i,:)')/S(i,i);
end
end
I = I/max(I); % 對I歸一化
if rem(N,2)==0
I_final = [fliplr(I),I]; % 最終的單元排列（左右對稱）
else
I_final = [fliplr(I),I(2:end)];
end
sprintf('天線單元歸一化電流幅度：')
sprintf('%.3f ',I_final)
%% 獲得最終陣列方向圖S_P
theta_rad = 0:0.01:pi;
theta = theta_rad*180/pi;
u = pi*d/lamuda*( cos(theta_rad)- cos(theta0));
S_P = zeros(1,length(theta_rad)); % 最終方向圖
for k = 1:M
if rem(N,2)==0
S_P = S_P + I(k)*cos((2*k-1)*u);% 偶數
else
S_P = S_P + I(k)*cos(2*(k-1)*u);% 奇數
end
end
S_P_abs = abs(S_P); % 對S_P取絕對值
S_PdB = 20*log10(S_P_abs/max(S_P_abs)); % 對S_P取dB值
%% 繪圖
H = -ones(1,length(S_P_abs))*26; % 根據預先設置的主副瓣比得到的參考曲線
% 直角坐標系
figure('NumberTitle', 'off', 'Name', 'S Parameter (abs)-Plot');
plot(theta,S_P_abs,'b','LineWidth',1.5)
xlabel('theta(°)')
ylabel('|S| ')
title('chebyshev低副瓣陣列直角坐標圖')
figure('NumberTitle', 'off', 'Name', 'S Parameter (dB)-Plot');
plot(theta,H,'r--','LineWidth',1.5)
hold on
plot(theta,S_PdB,'b','LineWidth',1.5)
xlabel('theta(°)')
ylabel('|S| dB')
title('chebyshev低副瓣陣列直角坐標圖')
legend('預設副瓣參考曲線','方向圖')
% 極坐標系
figure('NumberTitle', 'off', 'Name', 'S Parameter (dB)-Polar');
polarplot(theta_rad,H,'r--','LineWidth',1.5)
hold on
polarplot(theta_rad,S_PdB,'b','LineWidth',1.5)
thetalim([0 180]);
rmin = S_PdB(1,1);
rmax = max(S_PdB);
rlim([-50 rmax]);
title('chebyshev低副瓣陣列極坐標圖')
legend('預設副瓣參考曲線RdB','方向圖(dB)')
下面即為一個示例：單元間距d=0.6λ、單元數13、主副瓣電平比26dB、掃描角度80度（相對於單元排布方向）的切比雪夫陣列設計。
歸一化單元電流幅度比為：0.406 0.432 0.604 0.770 0.908 1.000 0.516 1.000 0.908 0.770 0.604 0.432 0.406
文末小詩：
離殤
——紅葉
絲絲秋雨思秋語，秋雨無情。秋語無情，落盡殘花，揮手贈余情。
戚戚南風別南楓，南風忽停。南楓忽停，終是離殤，望月感零丁。