| 查看: 871 | 回复: 2 | ||
[求助]
关于MATLAB程序设计的一些问题 已有1人参与
|
|
我需要的程序是：计算一元二维ABAB型啁啾结构的光子晶体，在A层为SiO2折射率为1.4,B层为Te折射率为4.6,,d1=d2=350nm,其中A层厚度不变，B层厚度以5nm为步长减少时，其波长-透射率图。帮我看一下我错在哪了，要怎么改，我用的是传输矩阵法 clc clear close all n0=1;%n0为空气折射率 a0=input('请输入入射角a0:'); N=input('请输入需要计算的周期数N:'); C=0.05；%啁啾系数（步长） for i = 1:N if(i>N) db=0.35-c; i=i+1; else end for d=1.25:0.001:2.35;%波长范围 na=1.4; nb=4.6; a1=asin(n0*sin(a0)/na);%第一分界面折射角 a2=asin(na*sin(a1)/nb);%第二分界面折射角 da=0.35; p1=2*pi*na*da*cos(a1)/d;%相位变化 p2=2*pi*nb*db*cos(a2)/d; u0=4*pi*1e-7;%真空磁导率 e0=1e-9/(36*pi);%真空介电常数 b=sqrt(e0/u0); c1=b*cos(a1)*na;%波阻抗 c2=b*cos(a2)*nb; A1=cos(p1);B1=-i*sin(p1)/c1;C1=-i*c1*sin(p1);D1=cos(p1); M1=[A1 B1;C1 D1];%介质1传输矩阵 A2=cos(p2);B2=-i*sin(p2)/c2;C2=-i*c2*sin(p2);D2=cos(p2); M2=[A2 B2;C2 D2];%介质2传输矩阵 M=(M1*M2)^N;%N个周期总的传输矩阵 A=M(1,1);B=M(1,2);C=M(2,1);D=M(2,2); a3=asin(nb*sin(a2)/n0);%射出介质时的出射角 c0=b*cos(a0)*n0;c00=b*cos(a3)*n0; r=(A*c0+B*c0*c00-C-D*c00)/(c0*A+c0*c00*B+C+D*c0);%反射率 r1=r*conj(r); tr=1-r1; end plot(d,tr,'k')% 画图 title('(a) 0');% 为当前坐标系添加标题 xlabel('入射波长(um)');ylabel('透射系数') |
回复此楼» 猜你喜欢
314求调剂
已经有5人回复
-
314求调剂
已经有11人回复
-
还有化工二轮调剂的学校吗
已经有28人回复
-
材料化工总分334求调剂
已经有12人回复
-
296求调剂
已经有5人回复
-
297求调剂
已经有14人回复
-
考研二轮调剂
已经有5人回复
-
求机械专硕297第二批调剂
已经有5人回复
-
本科郑州大学，一志愿华东师范大学282求调剂
已经有25人回复
-
283求调剂
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 《MATLAB编程与工程应用（第二版）》---Stormy Attaway
已经有284人回复
- 求解相算法的相关matlab程序块
已经有3人回复
- 优化工具箱GUI编程求助
已经有3人回复
- matlab中BP神经网络问题
已经有6人回复
- Matlab心得及学习方法（不断更新）
已经有26人回复
- matlab编程约束条件问题
已经有10人回复
- 我用matlab编写了一个风力机叶片的优化设计程序，有点问题
已经有13人回复
- matlab的m文件关联问题
已经有17人回复
- 求助：资源优化问题的遗传算法matlab编程
已经有5人回复
- 请matlab高手看一下这个问题怎样编程才好，附程序代码，谢谢！
已经有9人回复
- 学习matlab算法编程，推荐一些有用的资料或者书籍吧
已经有21人回复
- 哪一本matlab程序设计的书比较好？？
已经有12人回复
- matlab编程中遇到的问题
已经有4人回复
- 求助matlab编程
已经有13人回复
- 关于应用MATLAB工具箱BP网络编程问题
已经有1人回复
- 【求助】用MATLAB编程求级数的和
已经有16人回复
我来蹭经验
铁虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 26.8
- 帖子: 14
- 在线: 7.4小时
- 虫号: 3014244
- 注册: 2014-03-03
- 性别: GG
- 专业: 数理逻辑和与计算机相关的
2楼2014-05-07 20:57:06
3楼2014-05-30 08:25:01
