?? lcxeight.m
字號(hào):
clear;
a=9e-3; %內(nèi)導(dǎo)體半徑
b=22.8e-3; %外導(dǎo)體半徑
epsir=1.25; %介質(zhì)的介電常數(shù)
sita=pi/4; %開(kāi)槽的傾角
l0=20e-3; %開(kāi)槽長(zhǎng)度的一半
w=5e-3; %開(kāi)槽寬度
V0=0.01; %激勵(lì)電壓
f=24e8; %工作頻率900MHz
womiga=2*pi*f; %角頻率
c=3e8; %自由空間電磁波傳播速度
lamda=c/f; %波長(zhǎng)
k0=2*pi*f/c; %自由空間波數(shù)
mu0=4e-7*pi;
beta=k0*sqrt(epsir); %電纜內(nèi)電波傳播常數(shù)
beta1=beta*cos(sita); %開(kāi)縫方向傳播常數(shù)分量
P=0.2; %開(kāi)縫周期必須大于半波長(zhǎng)
syms l;
Fi=double(int(exp(-j*k0*sqrt(l^2+(w/2)^2))/(2*pi*sqrt(l^2+(w/2)^2)),l,-20e-3,20e-3))
H0=V0*sqrt(epsir)/(b*log(b/a)); %外導(dǎo)體內(nèi)側(cè)的周向磁場(chǎng)
A=-0.5j*womiga*mu0*H0*sin(sita)/[(k0^2-beta1^2)*Fi*sin(2*k0*l0)]; %磁流的系數(shù)
Nh=floor(k0*(1-sqrt(epsir))*P/(2*pi))
Nl=ceil(-k0*(1+sqrt(epsir))*P/(2*pi))
for n=Nl:Nh;
syms z1;
J(n-Nl+1)=double(1/P*int(2*A*cos(sita)*(sin((beta1-k0)*l0)*exp(j*k0*z1/cos(sita))-sin((beta1+k0)*l0)*exp(-j*k0*z1/cos(sita))...
+sin(2*k0*l0)*exp(-j*beta1*z1/cos(sita)))*exp(-j*n*2*pi/P*z1),z1,-l0*cos(sita),l0*cos(sita)));
end; %計(jì)算軸向磁流的傅里葉系數(shù)
s=1;
for z=0:0.1:200;
Z(s)=z;
Efi=0;
for n=Nl:Nh;
kesai(n-Nl+1)=sqrt(k0^2-(beta+2*pi*n/P)^2)
besselh(1,2,kesai(n-Nl+1)*1)
Efi=Efi+j/4*J(n-Nl+1)*kesai(n-Nl+1)*besselh(1,2,kesai(n-Nl+1)*2)*exp(-j*(beta+2*pi*n/P)*z);
end;
Pr(s)=abs(Efi);
s=s+1;
end;
Z
Pr
plot(Z,Pr);hold on;
?? 快捷鍵說(shuō)明
復(fù)制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號(hào)
Ctrl + =
減小字號(hào)
Ctrl + -