| 1 clear;clc;
2 %参数输入
3 Z0=50;ZL=100;%Z0：输入的参考阻抗，ZL：负载阻抗
4 lambda=300/3;%自定义工作频率下限的波长：mm
5 L=50;%渐变传输线的长度：mm
6 L0=linspace(0,L,251);tot_num=length(L0);
7 dL=L0(2)-L0(1);%离散化的传输线，每节的长度
8 %指数渐变传输线
9 a=1/L*log(ZL/Z0);
10 Z=Z0*exp(a*L0);
11 plot(L0,Z);%获取m文件当前路径
13p1=mfilename('fullpath');
14 i=findstr(p1,'\');
15 p1=p1(1:i(end));
16 pp=strcat(p1,'hfssapi');
17 %与m文件同级的文件夹中的hfssapi加入到路径中
18 addpath([pp '/3dmodeler']);
19 addpath([pp '/general']);
20 addpath([pp '/analysis']);
21 addpath([pp '/boundary/']);
22 %vbs脚本路径
23 tmpScriptFile = 'C:\Users\Dell\Desktop\exp_transmissionLine.vbs';
24 %计算微带线的线宽：mm
25 Er=2.2;H=0.508;Wcal=zeros(1,tot_num);
26 for num=1:tot_num
27 zcal=Z(num);
28 Wcal(num)=microstrip_Zcal(Er,H,zcal);
29 End
30
31 %vbs脚本建模：此处引用hfssapi为自己修改过后的，参数传入有差异；
32 %不可直接复制运行，当作伪代码，仅供参考
33 %建立HFSS模型文件.
34 fid = fopen(tmpScriptFile, 'wt');
35 % 建立一个新的HFSS project和design.
36 hfssNewProject(fid);
37 hfssInsertDesign(fid, 'gradually_changed_tri_transmissionLine');
38 % 创建新的Polyline
39 Points=zeros(2*tot_num,3);
40 for num=1:tot_num
41 Points(num,1:3)=[-Wcal(num)/2,(num-1)*dL,0];
42 End
43 Points((tot_num+1):end,1)=abs(flipud(Points(1:tot_num,1)));
44 Points((tot_num+1):end,2)=flipud(Points(1:tot_num,2));
45 Points((tot_num+1):end,3)=flipud(Points(1:tot_num,3));
46 hfssPolyline(fid, 'exp_transmissionLine', Points, 'mm', 1, 'Line', [255, 255, 0],0)
47 hfssSetColor(fid, 'exp_transmissionLine', [255,255,0]);
48 hfssSetTransparency(fid, {'exp_transmissionLine'}, 0);
49 hfssAssignPE(fid, 'Line', {'exp_transmissionLine'});
50
51 % Create the substrate
52 hfssBox(fid, 'substrate',[-40/2,0,0],[40,L,-H],'mm',1);
53 hfssAssignMaterial(fid, 'substrate', 'Rogers RT/duroid 5880 (tm)');
54 hfssSetTransparency(fid, {'substrate'}, 0.2);
55
56 % Create the ground
57 hfssRectangle(fid, 'GroundPlane', 'Z', [-40/2,0,-H],40,L, 'mm',1);
58 hfssSetColor(fid, 'GroundPlane', [192, 192, 192]);
59 hfssSetTransparency(fid, {'GroundPlane'}, 0.4);
60 hfssAssignPE(fid, 'PerfE_Ground', {'GroundPlane'});
61
62 %------------空气盒------------%
63 hfssBox(fid, 'AirBox', [-40/2-lambda/4,0,-H-lambda/4],[40+lambda/2, L, H+lambda/2], 'mm',1);
64 hfssSetTransparency(fid, {'AirBox'}, 0.9);
65 hfssAssignRadiation(fid, 'Rad1', 'AirBox');
66 %-----------设置端口---------------%
67 hfssRectangle(fid, 'waveport1', 'Y', [-Wcal(1)*6/2,0,-H],5*H,Wcal(1)*6, 'mm',1);
68 hfssRectangle(fid, 'waveport2', 'Y', [-Wcal(end)*6/2,L,-H],5*H,Wcal(end)*6, 'mm',1);
69 hfssSetColor(fid, 'waveport1', [255, 0, 0]); hfssSetColor(fid, 'waveport2', [255, 0, 0]);
70 hfssAssignWavePort(fid, 'Port1', 'waveport1', 1, false, [0,0,-H], [0,L,-H], 'mm');
71 hfssAssignWavePort(fid, 'Port2', 'waveport2', 1, false, [0,0,0], [0,L,0], 'mm');
72 %-----------改变积分终点--------------%
73 hfssInsertSolution(fid, 'Setup3000MHz', 3, 0.02, 15);
74 hfssInterpolatingSweep(fid, 'Sweep3to10GHz', 'Setup3000MHz',3, 10,351);
75
76 % Close the HFSS Script File.
77 fclose(fid);
78
79 function W=microstrip_Zcal(Er,H,Z0)
80 %Unit：mm，GHz，ohm
81 %计算中间变量
82 A=Z0/60*sqrt((Er+1)/2)+(Er-1)/(Er+1)*(0.23+0.11/Er);
83 B=377*pi/(2*Z0*sqrt(Er));
84 W_d_A=8*exp(A)/(exp(2*A)-2);
85 W_d_B=2/pi*(B-1-log(2*B-1)+(Er-1)/(2*Er)*(log(B-1)+0.39-0.61/Er));
86 if(W_d_A<2)
87 W_d=W_d_A;
88 else
89 W_d=W_d_B;
90 end
91 %计算微带线宽度
92 W=W_d*H;93 end |
