资源简介
用有限差分法求右图所示金属矩形波导的基模及第一个高阶模的色散特性曲线。
代码片段和文件信息
clear;clc;close all;
sz_X=41;sz_Y=23;
sz_M=sz_X-2;sz_N=sz_Y-2;
h=0.254;
% sz_X=81;sz_Y=45;
% sz_M=sz_X-2;sz_N=sz_Y-2;
% h=0.127
% sz_X=21;sz_Y=12;
% sz_M=sz_X-2;sz_N=sz_Y-2;
% h=0.508;
%TM波计算
DM=diag(4*ones(sz_M1)0)+diag(-ones(sz_M-11)1)+diag(-ones(sz_M-11)-1);
for i=1:sz_N
D1(::i)=DM;
end
sz_D1 = size(D1);
c = mat2cell(reshape(D1sz_D1(1)[])sz_D1(1)repmat(sz_D1(2)1sz_D1(3)));
K1=feval(@blkdiagc{:});
KM=K1+diag(-ones(1sz_M*sz_N-sz_M)sz_M)+diag(-ones(1sz_M*sz_N-sz_M)-sz_M);
[v_TMd_TM]=eig(KM);
[dx_TMdy_TM]=find(d_TM>0);
l_dxTM=length(dx_TM);
for i=1:l_dxTM
d2_TM(i)=d_TM(dx_TM(i)dy_TM(i));
end
for i=1:l_dxTM-1
minIndex=i;
for j=i+1:l_dxTM
if(d2_TM(j) minIndex=j;
end
end
temp=d2_TM(i);
d2_TM(i)=d2_TM(minIndex);
d2_TM(minIndex)=temp;
end
kc_TM=sqrt(d2_TM)/h;
%TE波计算
DE=diag(4*ones(sz_X1)0)+diag(-ones(sz_X-11)1)+diag(-ones(sz_X-11)-1);
DE(12)=-2;DE(sz_Xsz_X-1)=-2;
for i=1:sz_Y
D2(::i)=DE;
end
sz_D2 = size(D2);
c = mat2cell(reshape(D2sz_D2(1)[])sz_D2(1)repmat(sz_D2(2)1sz_D2(3)));
K1=feval(@blkdiagc{:});
KE=K1+diag(-1*ones(1sz_X*sz_Y-sz_X)sz_X)+diag(-1*ones(1sz_X*sz_Y-sz_X)-sz_X);
for i=1:sz_X
KE(ii+sz_X)=-2;
KE(i+sz_X*(sz_Y-1)i+sz_X*(sz_Y-2))=-2; %KE矩阵
end
[v_TEd_TE]=eig(KE);
[dx_TEdy_TE]=find(d_TE>1e-6);
l_dxTE=length(dx_TE);
for i=1:l_dxTE
d2_TE(i)=d_TE(dx_TE(i)dy_TE(i));
end
for i=1:l_dxTE-1
minIndex=i;
for j=i+1:l_dxTE
if(d2_TE(j)
- 上一篇:并网仿真模型matlab
- 下一篇:多通道滤波器组
相关资源
- Matlab路面裂缝识别69319
- 高灵敏度GPS接收机MATLAB仿真,附捕获
- 基于MATLAB的质点弹道计算与外弹道优
- 阵列天线的matlab仿真
- MATLAB 经典程序源代码大全
- MATLAB小波软阈值去噪代码33473
- 天线阵的波束形成在MATLAB仿真程序及
- 非线性SVM算法-matlab实现
- 《MATLAB 智能算法超级学习手册》-程序
- 组合导航matlab程序
- 读取txt文件内容matlab代码实现
- Matlab实现基于相关的模板匹配程序
- matlab优化工具箱讲解
- 基于MATLAB的快速傅里叶变换
- 光纤传输中的分布傅立叶算法matlab实
- 基于matlab的图像处理源程序
- matlab 椭圆拟合程序
- 算术编码解码matlab源代码
- optical_flow 光流法 matlab 实现程序
- 引导图像滤波器 Matlab实现
- 分形几何中一些经典图形的Matlab画法
- OFDM系统MATLAB仿真代码
- SVM工具箱(matlab中运行)
- 图像小波变换MatLab源代码
- LU分解的MATLAB实现
- 冈萨雷斯数字图像处理matlab版(第三
- 替代数据法的matlab程序
- 用matlab实现的多站定位系统性能仿真
- 通过不同方法进行粗糙集属性约简m
- k近邻算法matlab实现
评论
共有 条评论