%
Jm=[2.40482555769577,5.52007811028631,8.65372791291101,11.7915344390143,14.9309177084878,18.0710639679109,21.2116366298793,24.3524715307493,27.4934791320403,30.6346064684320;
3.83170597020751,7.01558666981562,10.1734681350627,13.3236919363142,16.4706300508776,19.6158585104682,22.7600843805928,25.9036720876184,29.0468285349169,32.1896799109744;
5.13562230184068,8.41724414039986,11.6198411721491,14.7959517823513,17.9598194949878,21.1169970530218,24.2701123135731,27.4205735499846,30.5692044955164,33.7165195092227;
6.38016189592398,9.76102312998167,13.0152007216984,16.2234661603188,19.4094152264350,22.5827295931044,25.7481666992950,28.9083507809218,32.0648524070977,35.2186707386101;
7.58834243450380,11.0647094885012,14.3725366716176,17.6159660498048,20.8269329569624,24.0190195247711,27.1990877659813,30.3710076671172,33.5371377118192,36.6990011287447;
8.77148381595995,12.3386041974669,15.7001740797117,18.9801338751799,22.2177998965613,25.4303411542227,28.6266183072911,31.8117167240478,34.9887812945593,38.1598685619671;
9.93610952421768,13.5892901705412,17.0038196678160,20.3207892135665,23.5860844355814,26.8201519834114,30.0337223865705,33.2330417628471,36.4220196682585,39.6032394160754;
11.0863700192451,14.8212687270132,18.2875828324817,21.6415410198484,24.9349278876730,28.1911884594832,31.4227941922656,34.6370893520693,37.8387173828536,41.0307736915855;
12.2250922640047,16.0377741908877,19.5545364309971,22.9451731318746,26.2668146411766,29.5456596709986,32.7958000373415,36.0256150638696,39.2404479951781,42.4438877432736;
13.3543004774353,17.2412203824891,20.8070477892641,24.2338852577506,27.5837489635730,30.8853789676967,34.1543779238551,37.4000999771566,40.6285537189645,43.8438014203374;
14.4755006865545,18.4334636669666,22.0469853646978,25.5094505541828,28.8873750635305,32.2118561997127,35.4999092053739,38.7618070178817,42.0041902366718,45.2315741035350];
%Jm为bessel函数的根。 Jm是11*10。竖是从0阶到10阶,横的是从第1个根到第10个根。
%
Jmp=[3.83170597020751,7.01558666981562,10.1734681350627,13.3236919363142,16.4706300508776,19.6158585104682,22.7600843805928,25.9036720876184,29.0468285349169,32.1896799109744;
1.84118378134066,5.33144277352503,8.53631636634629,11.7060049025921,14.8635886339090,18.0155278626818,21.1643698591888,24.3113268572108,27.4570505710592,30.6019229726691;
3.05423692822714,6.70613319415846,9.96946782308759,13.1703708560161,16.3475223183218,19.5129127824882,22.6715817724774,25.8260371417853,28.9776727729937,32.1273270204435;
4.20118894121053,8.01523659837595,11.3459243107430,14.5858482861670,17.7887478660665,20.9724769365377,24.1448974329093,27.3100579302043,30.4702688062904,33.6269491827967;
5.31755312608399,9.28239628524161,12.6819084426389,15.9641070377316,19.1960288000489,22.4010322676890,25.5897596813867,28.7678362176665,31.9385393409728,35.1039166773468;
6.41561637570024,10.5198608737723,13.9871886301403,17.3128424878846,20.5755145213869,23.8035814765939,27.0103078977777,30.2028490789817,33.3854439010101,36.5607776868804;
7.50126614468415,11.7349359530427,15.2681814610979,18.6374430096662,21.9317150178022,25.1839255994996,28.4097763625101,31.6178757161050,34.8133929842974,37.9996408977153;
8.57783648971407,12.9323862370896,16.5293658843669,19.9418533665273,23.2680529264576,26.5450320618236,29.7907485831966,33.0151786413751,36.2243805487872,39.4222745789393;
9.64742165199722,14.1155189078946,17.7740123669153,21.2290626228531,24.5871974863177,27.8892694279551,31.1553265561883,34.3966285542722,37.6200780441971,40.8301786818220;
10.7114339706999,15.2867376673330,19.0045935379461,22.5013987267773,25.8912772768391,29.2185634999361,32.5052473523755,35.7637929288088,39.0019028115142,42.2246384307533;
11.7708766749556,16.4478527484865,20.2230314126817,23.7607158603274,27.1820215271905,30.5345047540071,33.8419657751357,37.1180004236656,40.3710689053339,43.6067649013795];
%Jmp为bessel函数导数的根。 Jmp是11*10。竖是从0阶到10阶,横的是从第1个根到第10个根。
%
%输入相关常数
e0=8.85e-12;
u0=4.0*pi*1.0e-7;
er=1.0;
ur=1.0;
epslion=e0*er;
mur=u0*ur;
a=0.03; %波导半径
z=0; %显示z=0平面的场值
%
%选择传输模式类型
disp('TM=-1 or TE=1,m and n');
wavetype=input('wavetype is:');
m=input('input m: ');
n=input('input n: ');
%
%计算截止频率
if wavetype==1 %TE,导数根Jmp
Jmpo=Jmp(m+1,n);
wc=Jmpo/(a*sqrt(epslion*mur));
elseif wavetype==-1 %TM,函数根Jm
Jmo=Jm(m+1,n);
wc=Jmo/(a*sqrt(epslion*mur));
end
fc=wc/(2*pi);
fc=fc/1e9;
disp('cutoff fc is (GHz):');
%输入工作频率
disp('input f');
f=input('f is (GHz):');
w=2*pi*f*1e9;
%
%画图
r=linspace(0,a,300)';
phi=linspace(0,2*pi,300);
x=r*cos(phi);
y=r*sin(phi);
%
if wavetype==1 %TE
kz=sqrt(w^2*epslion*mur-(Jmpo/a)^2); %计算z方向传播常数
krou=sqrt(w^2*epslion*mur-kz^2);
%Jmd=(besselj((m-1),r*krou)-besselj(m,r*krou)*m./(r*krou));
Hz=besselj(m,r*krou)*cos(phi*m)*cos(z*kz);
surf(x,y,Hz),
colorbar
%
elseif wavetype==-1 %TM
kz=sqrt(w^2*epslion*mur-(Jmo/a)^2); %计算z方向传播常数
krou=sqrt(w^2*epslion*mur-kz^2);
Ez=besselj(m,r*krou)*cos(phi*m)*cos(z*kz);
surf(x,y,Ez),
colorbar
end
shading flat
shading interp
%
%
评论0