www.pudn.com > NearFieldMeasurement.rar > RecHorn.m, change:2016-01-18,size:3247b

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 此文件为角锥喇叭理论E面、H面方向图
%
% E面方向图：做出考虑波形因子和不考虑情况下的对比图
% 扫描角度： -60到+60度
% 参考文献：魏文元，宫德明，天线原理[M]. 国防工业出版社, 1985.
% 17/1/2016
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clc
%clear
%% %%%%%%%%%%%%    定义角锥喇叭结构参数      %%%%%%%%%%%%%%%%%%%%%%%
a = 23*10.^(-3);  %波导输入端宽度 m（H面）
b = 10*10.^(-3);  %波导输入端宽度 m（E面）
D1 = 238*10.^(-3);%角锥口径宽度 m（H面）
D2 = 176*10.^(-3);%角锥口径宽度 m（E面）
h = 465*10.^(-3); %角锥喇叭长度

%% %%%%%%%%%%%%%%      定义工作参数  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
f = 9.375*10.^9;         %工作频率 Hz
lamd = 3*10.^8/f;        %工作波长 m
R1 = h/(1-a/D1);         %H面喇叭虚顶点到口径中心距离
R2 = h/(1-b/D2);         %E面喇叭虚顶点到口径中心距离
theta1 = -60:0.2:60;     %观测范围
theta  = theta1.*pi/180; %角度转化为弧度

%% %%%%%%%%%%%%%%   H面方向图作图     %%%%%%%%%%%%%%%%%%%%%%%%%%%%%

M = 0.5*sqrt(lamd*R1/2)*exp( 1i*(pi/4)*lamd*R1*(1/D1 + 2*sin(theta)/lamd).^2 );
N = 0.5*sqrt(lamd*R1/2)*exp( 1i*(pi/4)*lamd*R1*(1/D1 - 2*sin(theta)/lamd).^2 );
v1 = 0.707*(sqrt(lamd*R1)*(1/D1 + 2*sin(theta)/lamd) + D1/sqrt(lamd*R1));
v2 = 0.707*(sqrt(lamd*R1)*(1/D1 + 2*sin(theta)/lamd) - D1/sqrt(lamd*R1));
v3 = 0.707*(sqrt(lamd*R1)*(1/D1 - 2*sin(theta)/lamd) + D1/sqrt(lamd*R1));
v4 = 0.707*(sqrt(lamd*R1)*(1/D1 - 2*sin(theta)/lamd) - D1/sqrt(lamd*R1));
%%一个符号错误就导致完全错误！认真检查

FH = (1+cos(theta)).*(M.*Fresnel(v1,v2) + N.*Fresnel(v3,v4));  %方向图函数

FH_M=abs(FH);         %取模值
FH_1=FH_M./max(FH_M); %归一化
FHdB=20*log10(FH_1);  %采用分贝形式

figure(1)
hold on
plot(theta1,FHdB,'b','LineWidth',1.5);
%plot(t,FHA,'--r','LineWidth',1.6);  % %需要比较实测去掉注释 并已经在变量空间加载实测值
%legend('理论H方向图','实测H方向图');
grid on
title('角锥喇叭H面方向图(理论实测比较)','FontName','华文隶书',...
'FontWeight','Bold','FontSize',16)
xlabel('\bf角度θ   单位：度','FontSize',15);
ylabel('\bff(θ)    单位：dB','FontSize',15);

%% %%%%%%%%%%%%%%   E面方向图作图     %%%%%%%%%%%%%%%%%%%%%%%%%%%%%

G =sqrt(1-((lamd/(2*a))^2)); %波形因子
lamdg=lamd/G;

for n=1:2   %循环做出波导波长近似为空气波长 和无近似的方向图比较
if n==2
lamdg=lamd;
G=1;
end

v5 = 0.707*(sqrt(lamdg*R2).*(2*sin(theta)/lamd) + D2/sqrt(lamdg*R2));
v6 = 0.707*(sqrt(lamdg*R2).*(2*sin(theta)/lamd) - D2/sqrt(lamdg*R2));

FE=(1+G*cos(theta)).*Fresnel(v5,v6); %方向图函数
FEM=abs(FE);        %取模值
FE1=FEM./max(FEM);  %归一化
FEdB=20*log10(FE1); %采用分贝形式

figure(2)
hold on
if n==1
plot(theta1,FEdB,'--g','LineWidth',1.5);
else
plot(theta1,FEdB,'b','LineWidth',1.5);
end
%plot(tt,FEA,'--r','LineWidth',2);  %需要比较实测去掉注释
%legend('理论E方向图','理论E方向图(G=1)','实测E方向图');
grid on
title('角锥喇叭E面方向图（理论实测比较）','FontName','华文隶书',...
'FontWeight','Bold','FontSize',16)
xlabel('\bf角度θ   单位：度','FontSize',15);
ylabel('\bff(θ)    单位：dB','FontSize',15);
end