一、⛄获取代码方式
获取代码方式1:
完整代码已上传我的资源:【波导仿真】基于matlab模式匹配法分析波导H面不连续性【含Matlab源码 3557期】
点击上面蓝色字体,直接付费下载,即可。
获取代码方式2:
付费专栏Matlab物理应用(初级版)
备注:
点击上面蓝色字体付费专栏Matlab物理应用(初级版),扫描上面二维码,付费29.9元订阅海神之光博客付费专栏Matlab物理应用(初级版),凭支付凭证,私信博主,可免费获得1份本博客上传CSDN资源代码(有效期为订阅日起,三天内有效);
点击CSDN资源下载链接:1份本博客上传CSDN资源代码
⛄二、部分源代码
clear all;
close all;
clc;
M=20;
a=7.12/1000;b=3.56/1000;%单位m
a1=1e-3*[3.9,2.85,2.65,2.65,2.85,3.9];
L=1e-3*[2,4.2,4.83,4.9,4.83,4.2,2];%前后对称,主波导的长度;
D=1e-31.2;%每个阶梯波导的长度;
mu=4pi1e-7;%磁导率
epsilon=1/36/pi1e-9;%介电常数
f=301e9:0.0251e9:401e9;%频率范围
for m=1:length(f) %做频率点数。
ST11=zeros(M);ST12=eye(M);ST21=ST12;ST22=ST11;%初始化广义散射矩阵S;
[SL11,SL12,SL21,SL22] = Transline( L(1),a,f(m),M );%计算第一个波导传输线的S参数矩阵;
[S11,S12,S21,S22]= Transline( L(1),a,f(m),M );%这句代码和下句代码都可以;
% [S11,S12,S21,S22] = Cascade(ST11,ST12,ST21,ST22,SL11,SL12,SL21,SL22);%第一次级联,可以不要,因为刚开始ST全部为零;
% [S11,S12,S21,S22] = HPlaneStepGSM( a,a1(1),b,f(m),D,M ) ;
for n=1:length(L)-1
[SL11,SL12,SL21,SL22] = HPlaneStepGSM( a,a1(n),b,f(m),D,M ) ;%计算6次不连续面S矩阵
[SR11,SR12,SR21,SR22] = Transline( L(n+1),a,f(m),M );%计算6次传输波导S矩阵;
[ST11,ST12,ST21,ST22] = Cascade(S11,S12,S21,S22,SL11,SL12,SL21,SL22);%依次级联
[S11,S12,S21,S22] = Cascade(ST11,ST12,ST21,ST22,SR11,SR12,SR21,SR22);%依次级联
end
%最关心的两个参数!!!
S11_result(m)=S11(1,1);
S21_result(m)=S21(1,1);
%S11=ST22;S12=ST21;S21=ST12;S22=ST11;
%[S11,S12,S21,S22] = Cascade(ST11,ST12,ST21,ST22,S11,S12,S21,S22);
%输出S1和S2和HFSS仿真进行对比。
S1= 20log10(abs(S11_result’));
S2= 20log10(abs(S21_result’));
end
plot(1e-9f,20log10(abs(S11_result)),1e-9f,20*log10(abs(S21_result)));
title(‘H面双边不连续性分析实例—Ka波段波导H面膜片滤波器’);
xlabel(‘f/GHz’);ylabel(‘S21/S11衰减曲线’);
grid on;
⛄三、运行结果
⛄四、matlab版本及参考文献
1 matlab版本
2014a
2 参考文献
[1]徐光宪.电磁场与电磁波[M].高等教育出版社
3 备注
简介此部分摘自互联网,仅供参考,若侵权,联系删除