*A+c0*c00*B+C+D*c0);%反射率Р r1=r*conj(r);Р s=s+1;Р kr(1,s)=r1;Р endРd=1.25:0.001:2.35;Рhold onРaxis([1.25 2.35 0 1])Рset(gca,'XTick',[1.25:0.1:2.35])Р switch dbР case 0.35Р subplot(3,1,1)Р plot(d,kr)Р xlabel('入射波长(um)'); ylabel('反射系数')Р legend('db=0.35')Р case 0.36Р subplot(3,1,2)Р plot(d,kr,'r:')Р xlabel('入射波长(um)');ylabel('反射系数')Р legend('db=0.36')Р case 0.37Р subplot(3,1,3)Р plot(d,kr,'g-.')Р xlabel('入射波长(um)');ylabel('反射系数')Р legend('db=0.37')РendРendР三、程序结果与分析Р该程序在matlab6.5中运行结果如下图所示Р 图4 反射率与入射波长的关系图谱Р观察结果发现,当保持A层厚度不变,而增加B层厚度发现带隙红移(光谱向长波段(红)段的位移);反观结果,随着膜层变薄,带隙会出现蓝移(光谱向短波长波段(紫)段的位移)。在物理上可以理解为增加(或减小)结构单元厚度,导致光在周期性结构中传输时光程的相应增加(或减小),必然导致反射加强的部分向波长长(或短)的位置移动。Р四、心得体会Р此次设计还有很多不完善的地方,还有待提高。比如程序的细想,我想还有很大的提升空间。Р程序设计过程中,有许多难题,上网参考了诸多资料,在此感谢浩瀚的互联网资源。其次,也感谢系主任黄仙山老师,给我们提供了这样一个动手、实际操作的机会,在此表示感谢,3Q!
全文预览
