线阵分析-matlab
05-08







%%
%linear array
% number of elements to mainlobe and sidelobe beamwidth, directivity
% different amplitude taper
% 1. W(n) = 1 ;
% 2. W(n) = abs((2*n-1)/N-1);
% 3. W(n) = abs(((2*n-1)/N)^2-1);
% 4. W(n) = power(cos(((2*n-1)/N-1)*pi/2),1)
% 5. W(n) = power(cos(((2*n-1)/N-1)*pi/2),2)
% 6. W(n) = power(cos(((2*n-1)/N-1)*pi/2),3)
% 7. W(n) = power(cos(((2*n-1)/N-1)*pi/2),4)
clc;
clear all;
close all;
N = 9; %
N = power(2,2:1:N);%elements number array
Ka =7 ; %antenna element amplitude taper conf
M = 2000; % theta sample points
AF = zeros(length(N),M);% array factor buffer
E0 = zeros(length(N),M);%
H0 = zeros(length(N),M);%
U0= zeros(length(N),M);%
f0 = 1e9;
lambda = 3e8/f0;
d = lambda/2;% elements spacing
k =2*pi/lambda;% wave constant
imp = 120*pi;%wave impdance
I0= 1;
r = 100;
l = 0.01;
BeamWidth = zeros(Ka,length(N));
D = zeros(Ka,length(N));
tau = zeros(Ka,length(N));
for m = 1:1:Ka
figure;
for i = 1:1:length(N)
theta = linspace(0,pi,M)+pi/10e10;
phi = linspace(0,2*pi,M);
af(i,:) = zeros(1,M);
a = 0;
b = 0;
for j = 1:1:N(i)
switch(m)
case 1
a= a+ 1;
b = b+1;
af(i,:) = af(i,:)+ exp(1i*(j-1)*k*d*cos(theta));
case 2
a= a+power(cos(((2*j-1)/N(i)-1)*pi/2),1);
b = b + power(cos(((2*j-1)/N(i)-1)*pi/2),2);
af(i,:) = af(i,:)+ exp(1i*(j-1)*k*d*cos(theta))*power(cos(((2*j-1)/N(i)-1)*pi/2),1);
case 3
a= a+power(cos(((2*j-1)/N(i)-1)*pi/2),2);
b = b + power(cos(((2*j-1)/N(i)-1)*pi/2),4);
af(i,:) = af(i,:)+ exp(1i*(j-1)*k*d*cos(theta))*power(cos(((2*j-1)/N(i)-1)*pi/2),2);
case 4
a= a+power(cos(((2*j-1)/N(i)-1)*pi/2),3);
b = b + power(cos(((2*j-1)/N(i)-1)*pi/2),6);
af(i,:) = af(i,:)+ exp(1i*(j-1)*k*d*cos(theta))*power(cos(((2*j-1)/N(i)-1)*pi/2),3);
case 5
a= a+power(cos(((2*j-1)/N(i)-1)*pi/2),4);
b = b + power(cos(((2*j-1)/N(i)-1)*pi/2),8);
af(i,:) = af(i,:)+ exp(1i*(j-1)*k*d*cos(theta))*power(cos(((2*j-1)/N(i)-1)*pi/2),4);
case 6
a= a+power(0.5+0.5*cos(((2*j-1)/N(i)-1)*pi/2),1);
b = b + power(0.5+0.5*cos(((2*j-1)/N(i)-1)*pi/2),2);
af(i,:) = af(i,:)+ exp(1i*(j-1)*k*d*cos(theta))*power(0.5+0.5*cos(((2*j-1)/N(i)-1)*pi/2),1);
case 7
a= a+power(0.33+0.67*cos(((2*j-1)/N(i)-1)*pi/2),1);
b = b + power(0.33+0.67*cos(((2*j-1)/N(i)-1)*pi/2),2);
af(i,:) = af(i,:)+ exp(1i*(j-1)*k*d*cos(theta))*power(0.33+0.67*cos(((2*j-1)/N(i)-1)*pi/2),1);
end
end
af(i,:) = abs(af(i,:)/N(i));
E0(i,:) = af(i,:) *imp;
Wav = real(E0(i,:) .* E0(i,:) )/2/imp;
Prad = imp*pi*3*(I0*l/lambda)^2;
theta = theta/pi*180;
af(i,:)= 10*log10(af(i,:)/max(af(i,:)));
U(i,:)= real(E0(i,:) .*E0(i,:))*r^2/2/imp;
U(i,:)= 10*log10(U(i,:)/max(U(i,:)));
plot(theta,U(i,:));
hold on;
axis([0 180 -120 0]);
error = 0.01;
af(i,:)=U(i,:);
pos1_3dB = [];
pos_max = find(max(af(i,:))==af(i,:));
while(isempty(pos1_3dB))
pos1_3dB = find(abs(((af(i,1:pos_max)-af(i,pos_max)))+3) < error);
error = error + 0.005;
end
error = 0.01;
pos2_3dB = [];
while(isempty(pos2_3dB))
pos2_3dB = find(abs(((af(i,pos_max:end)-af(i,pos_max)))+3) < error);
error = error + 0.005;
end
BeamWidth(m,i)= (theta(pos2_3dB(1)+pos_max)-theta(pos1_3dB(end)));
D(m,i) = 10*log10(a^2/b);
tau(m,i) = a^2/b/N(i);
end
switch(m)
case 1
% title('\[{W_n} = 1\]1');
text(80,5,'$$W_{n}= 1$$','interpreter','latex','fontsize',10);
% text(0.5,0.5,'$$f(x)=\frac{sin(x)}{x}$$','interpreter','latex','fontsize',20)
case 2
text(50,5,'$$W_{n}=cos\left [ \left ( \frac{2n-1}{N} -1\right )\frac{\pi}{2} \right ]$$','interpreter','latex','fontsize',10)
case 3
text(50,5,'$$W_{n}=cos^{2}\left [ \left ( \frac{2n-1}{N} -1\right )\frac{\pi}{2} \right ]$$','interpreter','latex','fontsize',10)
case 4
text(50,5,'$$W_{n}=cos^{3}\left [ \left ( \frac{2n-1}{N} -1\right )\frac{\pi}{2} \right ]$$','interpreter','latex','fontsize',10)
case 5
text(50,5,'$$W_{n}=cos^{4}\left [ \left ( \frac{2n-1}{N} -1\right )\frac{\pi}{2} \right ]$$','interpreter','latex','fontsize',10)
case 6
text(50,5,'$$W_{n}=0.5+0.5cos\left [ \left ( \frac{2n-1}{N} -1\right )\frac{\pi}{2} \right ]$$','interpreter','latex','fontsize',10)
case 7
text(50,5,'$$W_{n}=0.33+0.67cos\left [ \left ( \frac{2n-1}{N} -1\right )\frac{\pi}{2} \right ]$$','interpreter','latex','fontsize',10)
end
end










说好的可读性呢?
相关文章:
- 高速移动数据业务技术实现分析:HSPA对比EV-DO(05-08)
- 分析宽频带微带天线的辐射原理(05-08)
- 网络分析仪E5071C 租赁(05-08)
- 求助关于怎样使用矢量分析仪测天线的时域响应的问题(05-08)
- 矩量法分析线天线(05-08)
- 网络分析仪8753E 租赁/销售(05-08)
天线设计培训教程推荐