西邮MATLAB光波场的时域频谱分析实验报告.docxVIP

. 光学仿真课程设计 实验报告 课程名称： 光学仿真课程设计 姓 名： 学 院： 电子工程学院 系 部： 光电子技术系 专 业： 年 级： 学 号： 指导教师： 李晓莉 职 称： 讲师 时 间： 2013-11-18 至 2013-11-29 . . 光波场的时域频谱分析 一、实验目的 对一些常见光波进行傅里叶变换计算并绘出频谱图， 总结影响频谱宽窄的因素。 二、实验原理 无限长时间的等幅振荡 其表达式为 ( ) = E0 e -i2 πυt - ∞ ?? ∞ 0 E t 式中 ，E 、υ为常数，且 E 可以取复数值。由 0 0 0 ∞ -i2 πυt ( υ) = [ ( )] = ∫ ( ) E F E t E t e dt 式，它的频谱为 - ∞ ∞ -i2 πυ t i2 ( υ) ∫ πυt = E0 e 0 e dt - ∞ = E ∞ -i2 π(υ - υ) t dt ∫ e 0 0 - ∞ = E δυ υ 0 ( - 0 ) 该式说明，等幅振荡光场对应的频谱只含有一个频率成分 2 为理想单色振动。其功率谱为 |E(υ)|，如图所示。  υ，我们称其 0 持续有限时间的等幅振荡其表达式为 ( 设振幅等于 1) -i2 πυt T T ( ) e 0 , - 2 ≤t ≤ 2 E t = { 0, 其它 这时 T 2 -i2 πυt 2πυt E(υ) = ∫T e 0 e dt - 2 . . sin πT(υ- υ) 0 = T πT(υ- υ0) 或表示成 E( υ) = Tsinc[T ( υ- υ )] 0 其相应功率谱为 |E( υ)| 2 = T2 sinc 2 [T ( υ- υ0) ] 如图所示。 可见，这种光场频谱的主要部分集中在从 υ到υ 的频率范围之内，主峰 1 2 中心位于 υ处， υ是振荡的表现频率，或称为中心频率。 0 0 为表征频谱分布特性，定义最靠近 υ的两个强度为零的点所对应的频率 0 υ和υ 之差的一半为这个有限正弦波的频谱宽度 △υ。 2 1 | ( υ)| 2 2 2 υ υ 式，当 υ υ 时， | ( υ 2 2 = T sinc [T( 0 )| = T ; - 0 )] = 0 当υ υ 1 | ( υ)| = 0 ，所以有 △υ 1 。因此，振荡持续的时间越 = ± 时， E = 0 T T 长，频谱宽度愈窄。 衰减振荡 其表达式可写为 ( ) - βt -i2 πυ t , t ≥ 0 3) = { e e 0 E t 0, ?? 0 ( ) 为 ( ) ∞ - βt -i2 πυ0t i2 πυt = ∫ 相应的 E v E v e e e dt - ∞ ∞ [ ( ) ] = ∫ e dt 0 0 = 功率谱为  i 2π(υ- υ0) + i β | ( )| 2 ( υ) ? (υ) E v = E E = 1 2 υ υ) 2 2 4 π ( + β - 0 如图所示。 . . 可见，这个衰减振荡也可视为无限多个振幅不同、频率连续变化的简谐振荡的叠加， υ为其中心频率。这时，把最大强度一半所对应的两个频 0 率υ和 υ之差 △υ，定义为这个衰减振荡的频谱宽度。 2 1 由于 υ= υ ( 或 υ) ，|E( υ)| 2 = |E(υ0 )| 2 ，即 2 1 2 2 = 1 1 化简后 2 1 2 2 +β 2 4π (υ - υ) 2 β 0 得( υ υ) β所以 △υ= υ - υ = ( υ υ ) + ( υ - υ) β 2 -0 = 2π 2 1 2 - 0 0 1 = π . . 三、实验流程及程序 流程图： 开始 定义变量： t1,w1,f,Eo,B,tao 无限长时间的等幅 持续有限时间等幅 衰 减 震 荡 震荡 震荡 yt3=Eo*exp(-B*t3)*e yt1=Eo*exp(-2i* yt2=Eo*exp(-2i*pi xp(-2i*pi*f.*t3)*(h pi*f*t1) *f*t2)*(heaviside eaviside(t3)) (t2+tao)-heavisid Yw1=fft(yt1) e(t2-tao)) Yw3=fourier(yt3 Yw2=fourier(yt2) 显示函数 显示函数 显示函数 plot(t1,yt1) ezplot(t2,yt2,[ ezplot(t3,yt3,[- -2,2]) 1,10]) plot(f,abs(Yw1) ezplot(w3,abs(Yw ) ezplot(w2,abs(Y 3),[-10,-2]) w2),[-22,10]) 结束 . . 程序： clear all; close all; clc;