亚洲精品四区麻豆文化传媒,欧美精品一区二区大全,尤物邪恶午夜女性爽爽影院

首頁射頻工程師培訓 HFSS培訓 ADS培訓 CST培訓天線設計培訓微波測量操作培訓射頻專業(yè)基礎課官方淘寶店

,圓形陣列的泰勒綜合

文章來源: 互聯(lián)網錄入: mweda.com

n=input('陣元個數(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為公式推導中和主瓣電平有關的常量
n_1=round(2*A^2+1/2);          %******    n_1必須大于2*A^2+1/2
            %**********  判斷想控制的旁瓣數(shù)是否滿足條件   ***********%
n_2=input('想控制第幾對零點：')
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('輸入錯誤','Data message');
elseif  n_2<n_1
    errordlg('數(shù)太小了','Data message');
end
gama=n_1/(sqrt(A^2+(n_1-1/2)^2));     %********  展寬因子，用于計算半功率波瓣寬度
                                      %********  根據公式求泰勒方向圖的各個根
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根的平方，利于方程的計算
                    %*********  a_n即為各個根的向量
                    %********  構造公式中(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ū)別其他幾個綜合法(藍色三角形)
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]);
                      %*****    分奇偶討論電流激勵分布     *********%
if mod(n,2)==0          %*******    陣元個數(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;
                      %*******  利用第二個利于編程的歸一化方向圖函數(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);           %******* 翻轉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          %**********    陣元個數(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]
                           % 激勵分布圖
    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)與波長之比(lamda)之比：')
q=0.4%input('陣元距離(d)與波長之比(lamda)之比：')
n=9%input('陣元個數(shù)：')
r1=-25%input('要求所有副瓣歸一化最大值(R0/DB)：')
theta0=90%input('發(fā)射方向(度數(shù))：')
D=10^((abs(r1)+3)/20);      %******   最大的方向性系數(shù)(L無限長時)
fai1=0%input('步進激勵相位(度數(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));       %******   提取多項式系數(shù)
r1=roots(t); %多項式求根
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));    %******   提取需要的最后一個參量
m2=abs(m);                 %%%以上所有求線陣激勵，即環(huán)陣的遠場模
m2=m2/max(m2);              %%%%%%%歸一化
M=(n-1)/2;                   %%%%%%%%%%%%求近場模
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個波長)'},'FontSize',18)
grid on

hold on                   %%%%在同一圖中畫出遠場模（虛線）和近場模（實線）
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('模式標號','color','k')
ylabel('幅度','color','k')
%title('使用9個相模的圓形陣列的遠場模（虛線）和激勵模（實線）的歸一化幅值')
grid on

fi=linspace(0,2*pi,1000);           %%%%%%%%%%求環(huán)形天線的激勵并畫圖
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ù)激勵分布')
grid on

theta=linspace(0,2*pi,1000);            %%%%%%%%%%%%%實際方向圖求解
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({'由離散激勵構成的環(huán)形陣列方向圖（虛線）和理想方向圖（實線）';'(10個輻射單元，間距為0.4個波長)'},'FontSize',18);
grid on

%%%%%%%%%%%%%%%%%%%改變陣元數(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({'由離散激勵構成的環(huán)形陣列方向圖（虛線）和理想方向圖（實線）';'(11個輻射單元，間距為0.4個波長)'},'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)在想要圓陣的泰勒綜合,聽說只要用第一個直線泰勒程序換第二個程序里開頭的切比雪夫的內容
求大神幫忙搞一下,做出圓陣的泰勒綜合

網友回復:

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

申明：網友回復良莠不齊，僅供參考。如需專業(yè)解答，請咨詢本站專家，或者學習本站天線設計視頻培訓課程。

上一篇文章：小型化圓極化微帶天線

下一篇文章：環(huán)形天線的匹配問題。

天線設計培訓教程推薦

5門天線設計視頻課程和一本圖書，將天線設計理論和實踐相結合，幫助您全面、系統(tǒng)、深入地學習天線設計，讓天線設計不再難...【詳細介紹】

1天線設計培訓課程套裝
2HFSS天線設計培訓套裝
313.56MHz線圈天線設計培訓
4HFSS天線設計入門 — 中文視頻教程
5CST天線設計詳解 — 中文視頻教程
6PCB天線設計和HFSS仿真分析實例
7HFSS陣列天線設計實例視頻
8HFSS-IE仿真器應用詳解和天線分析實例
6天線技術原理 — 中文視頻課程
9HFSS天線設計 — 圖書