,圓形陣列的泰勒綜合

n=input('陣元個(gè)數(shù)：') 
f=3e8; 
lamda=2*pi/f; 
d=lamda/2; 
L=n*d; 
theta=linspace(0,pi,500); 
x=L/lamda*cos(theta); 
r1=input('要求所有副瓣歸一化最大值R0(DB)：') 
R0=10^(abs(r1)/20); 
A=acosh(R0)/pi;                %******    A為公式推導(dǎo)中和主瓣電平有關(guān)的常量 
n_1=round(2*A^2+1/2);          %******    n_1必須大于2*A^2+1/2 
            %**********  判斷想控制的旁瓣數(shù)是否滿足條件   ***********% 
n_2=input('想控制第幾對(duì)零點(diǎn)：') 
if n_2>=n_1&abs(n_2-round(n_2))==0 %判斷輸入是否正確 
     n_1=n_2 
elseif  abs(n_2-round(n_2))<1&abs(n_2-round(n_2))>0 
    errordlg('輸入錯(cuò)誤','Data message'); 
elseif  n_2<n_1 
    errordlg('數(shù)太小了','Data message'); 
end 
gama=n_1/(sqrt(A^2+(n_1-1/2)^2));     %********  展寬因子，用于計(jì)算半功率波瓣寬度 
                                      %********  根據(jù)公式求泰勒方向圖的各個(gè)根 
t=A^2+(n_1-1/2)^2;                
for i=1:n_1-1;%正的各根n<=n_1-1 
      a_n1(i)=n_1*sqrt((A^2+(i-1/2)^2)/(A^2+(n_1-1/2)^2)); 
end 
a_n=a_n1.^2;        %*********  Xn根的平方，利于方程的計(jì)算 
                    %*********  a_n即為各個(gè)根的向量 
                    %********  構(gòu)造公式中(1-X^2/Xn) 
for   j=1:n_1-1; 
       i=1:length(theta)     
        t1(j,i)=(1-x(i).^2./a_n(j))./(1-x(i).^2/j^2); 
end 
t1=prod(t1);  %%%%  prod()為連乘(1*500)的向量 
                    %*********  歸一化方向圖函數(shù)  *********% 
for i=1:length(theta) 
    G(i)=cosh(pi*A)*sin(pi.*x(i))./x(i)/pi.*t1(i); 
end 
g=G/max(G); 
                    %********        ********% 
g=20*log10(abs(g)); 
hold on 
figure(1) 
                       %*********  擬合,區(qū)別其他幾個(gè)綜合法(藍(lán)色三角形) 
theta=theta/pi*180; 
theta2=linspace(0,pi,40); 
theta2=theta2/pi*180' 
g2=interp1(theta,g,theta2); 
plot(theta,g,theta2,g2,'b^')  
grid on 
axis([0 180 -55 0]);                                       
                      %*****    分奇偶討論電流激勵(lì)分布     *********% 
if mod(n,2)==0          %*******    陣元個(gè)數(shù)為偶數(shù)個(gè)   *******% 
     for   i=1:n_1-1; 
       j=1:n_1-1; 
    s1(j,i)=1-(i^2./a_n(j)); 
     end 
    s2=prod(s1); 
      for i=1:n_1-1; 
                      %*******  利用第二個(gè)利于編程的歸一化方向圖函數(shù) *****% 
    s3(i)=(factorial(n_1-1))^2/(factorial(n_1-1+i))/(factorial(n_1-1-i)).*s2(i); 
      end 
       for i=1:n_1-1; 
     j=0:n/2-1; 
       In(i,j+1)=2*s3(i).*cos(i*pi*(2*j+1)/n); 
       end 
       In2=1.+sum(In); 
       In3=rot90(In2,2);           %******* 翻轉(zhuǎn)In2 但是要去掉公共的元素 
        for i=2:length(In2) 
         In5(i-1)=In2(i);     %*******     去掉i=0的元素 
        end 
        In4=[In3 In5] 
        hold on 
        figure(2) 
        plot(In4/max(In2),'-bs'); 
        grid on  
else          %**********    陣元個(gè)數(shù)為奇數(shù)  % 
     for i=1:n_1-1; 
          j=1:n_1-1; 
           s1(j,i)=1-(i^2./a_n(j)); 
     end 
         s2=prod(s1); 
     for i=1:n_1-1; 
         s3(i)=(factorial(n_1-1))^2/(factorial(n_1-1+i))/(factorial(n_1-1-i)).*s2(i); 
     end 
     for i=1:n_1-1; 
         j=0:(n-1)/2; 
           In(i,j+1)=2*s3(i).*cos(i*pi*2*j/n); 
     end 
    In2=1.+sum(In); 
    In3=rot90(In2,2); 
    for i=2:length(In2) 
       In5(i-1)=In2(i);        % 
     end 
       In4=[In3,In5] 
                           % 激勵(lì)分布圖 
    hold on 
    figure(2) 
    plot(In4/max(In2),'-bs'); 
    grid on 
end 
Bw=2*asin(lamda*gama/pi/L*((acosh(R0))^2-(acosh(R0/sqrt(2)))^2)^(1/2)); 
disp(['半功率寬度=(°)',num2str(Bw/pi*180)]); 
 
這是直線的泰勒綜合 
 
 
%%%%%%%%%%%% 均勻分布圓形陣列方向圖的切比雪夫綜合   %%%%%%%%%%%% 
clear;clc; 
R=1%input('環(huán)陣半徑(R)與波長(zhǎng)之比(lamda)之比：') 
q=0.4%input('陣元距離(d)與波長(zhǎng)之比(lamda)之比：') 
n=9%input('陣元個(gè)數(shù)：') 
r1=-25%input('要求所有副瓣歸一化最大值(R0/DB)：') 
theta0=90%input('發(fā)射方向(度數(shù))：') 
D=10^((abs(r1)+3)/20);      %******   最大的方向性系數(shù)(L無限長(zhǎng)時(shí)) 
fai1=0%input('步進(jìn)激勵(lì)相位(度數(shù))：') 
syms w tp s_w 
x=sym('x'); 
f=3e8; 
lamda=2*pi/f; 
L=n*q*lamda; 
k=2*pi/lamda; 
theta=linspace(0,pi,500); 
fai=fai1/180*pi; 
R0=10^(abs(r1)/20); 
for kk=1:n-2              %******   泰勒P階方程遞推 
  tp(1)=1; 
  tp(2)=x; 
  tp(kk+2)=2*x*tp(kk+1)-tp(kk) 
end 
t=sym2poly(tp(n));       %******   提取多項(xiàng)式系數(shù) 
r1=roots(t); %多項(xiàng)式求根 
x0=cosh(acosh(R0)/(n-1)); 
u=2*acos(r1'/x0); 
w1=exp(i*u);             %******   謝坤諾夫圓根 
for kk=1:n-2              %******   切比雪夫遞歸求解方程 
  s_w(1)=w-w1(1); 
  s_w(kk+1)=s_w(kk)*(w-w1(kk+1)); 
end 
m=sym2poly(s_w(n-1));    %******   提取需要的最后一個(gè)參量 
m2=abs(m);                 %%%以上所有求線陣激勵(lì)，即環(huán)陣的遠(yuǎn)場(chǎng)模 
m2=m2/max(m2);              %%%%%%%歸一化 
M=(n-1)/2;                   %%%%%%%%%%%%求近場(chǎng)模 
z=4 ;  
for kk=1:1:n 
    p=kk-M-1; 
    J(kk)=besselj(p,z); 
    C(kk)=m2(kk)/(J(kk)*(i^p)); 
    D(kk)=abs(C(kk)); 
end 
 
 
E=0;                           %%%%%%環(huán)陣的方向圖 
theta=linspace(-pi,pi,1000) 
for kk=1:n 
    E=E+m2(kk)*exp(i*kk*theta); 
end 
E=abs(E); 
hold on 
figure(1) 
E=E/max(E); 
E=20*log10(E); 
theta=theta/pi*180; 
plot(theta,E); 
axis([-180 180 -55 0]); 
xlabel('角度','color','k') 
ylabel('dB','color','k') 
%title({'均勻分布圓形陣列的方向圖';'(相模數(shù)為9，副瓣電平為-25dB，間距為0.4個(gè)波長(zhǎng))'},'FontSize',18) 
grid on 
 
hold on                   %%%%在同一圖中畫出遠(yuǎn)場(chǎng)模（虛線）和近場(chǎng)模（實(shí)線） 
figure(2) 
D=D/max(D); 
plot(m2,':bo'); 
hold on 
D=D/max(D); 
plot(D,'-bo'); 
set(gca,'XTickLabel',{'-4','-3','-2','-1','0','1','2','3','4'}); 
xlabel('模式標(biāo)號(hào)','color','k') 
ylabel('幅度','color','k') 
%title('使用9個(gè)相模的圓形陣列的遠(yuǎn)場(chǎng)模（虛線）和激勵(lì)模（實(shí)線）的歸一化幅值') 
grid on 
 
fi=linspace(0,2*pi,1000);           %%%%%%%%%%求環(huán)形天線的激勵(lì)并畫圖 
syms V  
V=0; 
for kk=1:n  
    p=kk-M-1; 
    V=V+C(kk)*exp(i*p*fi); 
end 
V2=V; 
V=abs(V); 
V=V/max(V); 
V=20*log10(V); 
fi=fi/pi*180; 
hold on 
figure(3) 
plot(fi,V); 
axis([0 360 -20 0]); 
xlabel('單元位置','color','k') 
ylabel('dB','color','k') 
%title('圓形陣列的連續(xù)激勵(lì)分布') 
grid on 
 
theta=linspace(0,2*pi,1000);            %%%%%%%%%%%%%實(shí)際方向圖求解 
R2=z/k; 
for kk=1:10                           
    kkk=kk-1; 
    Y(kk)=V2(50+kkk*100); 
end 
E2=0;                           
theta=linspace(-pi,pi,1000); 
for kk=1:10 
    kkk=kk-1; 
    E2=E2+Y(kk).*exp(i*z*(cos(0.1*pi+0.2*kkk*pi-theta))); 
end 
E2=abs(E2); 
E2=E2/max(E2); 
E2=20*log10(E2); 
theta=theta/pi*180; 
hold on                    
figure(4) 
plot(theta,E2,'b:'); 
hold on   
plot(theta,E) 
axis([-180 180 -55 0]); 
xlabel('角度','color','k') 
ylabel('dB','color','k') 
%title({'由離散激勵(lì)構(gòu)成的環(huán)形陣列方向圖（虛線）和理想方向圖（實(shí)線）';'(10個(gè)輻射單元，間距為0.4個(gè)波長(zhǎng))'},'FontSize',18); 
grid on 
 
%%%%%%%%%%%%%%%%%%%改變陣元數(shù)目后的實(shí)際方向圖求解%%%%%%%%%%%%%% 
fi=linspace(0,2*pi,1100);            
syms V  
V=0; 
for kk=1:n  
    p=kk-M-1; 
    V=V+C(kk)*exp(i*p*fi); 
end 
V3=V; 
theta=linspace(0,2*pi,1100); 
R2=z/k; 
for kk=1:11                          
    kkk=kk-1; 
    Y3(kk)=V3(50+kkk*100); 
end 
E3=0;                           
theta=linspace(-pi,pi,1100); 
for kk=1:11 
    kkk=kk-1; 
    E3=E3+Y3(kk).*exp(i*z*(cos(1/11*pi+2/11*kkk*pi-theta))); 
end 
E3=abs(E3); 
E3=E3/max(E3); 
E3=20*log10(E3); 
theta=theta/pi*180; 
hold on                    
figure(5) 
plot(theta,E3,'b:'); 
axis([-180 180 -55 0]); 
xlabel('角度','color','k') 
ylabel('dB','color','k') 
%title({'由離散激勵(lì)構(gòu)成的環(huán)形陣列方向圖（虛線）和理想方向圖（實(shí)線）';'(11個(gè)輻射單元，間距為0.4個(gè)波長(zhǎng))'},'FontSize',18); 
grid on 
E=0;                            
theta=linspace(-pi,pi,1100) 
for kk=1:n 
    E=E+m2(kk)*exp(i*kk*theta); 
end 
E=abs(E); 
hold on 
E=E/max(E); 
E=20*log10(E); 
theta=theta/pi*180; 
plot(theta,E); 
 
這是圓陣的切比雪夫綜合 
現(xiàn)在想要圓陣的泰勒綜合,聽說只要用第一個(gè)直線泰勒程序換第二個(gè)程序里開頭的切比雪夫的內(nèi)容 
求大神幫忙搞一下,做出圓陣的泰勒綜合 
 

網(wǎng)友回復(fù):

求大神幫個(gè)忙啊，弄出來了我的積分都送出以示感謝

申明：網(wǎng)友回復(fù)良莠不齊，僅供參考。如需專業(yè)解答，請(qǐng)咨詢本站專家，或者學(xué)習(xí)本站天線設(shè)計(jì)視頻培訓(xùn)課程。