Last report for Lab of Physics
本文最后更新于:1 年前
同轴光子晶体实验报告
实验仪器
信号发生器、示波器、测试电缆(约 1m
实验原理
同轴电缆电学参数
对于两种特征阻抗的同轴电览交替连接构成的同轴光子晶体,当两种电缆的特征阻抗之比
Z0H/Z0L
常用同轴电缆的特征阻抗为 75Ω
本实验只使用一种特征阻抗的同轴电缆。通过并联两根电缆的方法获得另一特征阻抗以及更高的阻抗比。
同轴电缆的等效电路如图 5
透射与反射
对于半无限长的同轴电缆,信号从同轴电缆始端输入后,只有沿 +
两根特征阻抗分别为 Z0
测量同轴光子晶体特性
由群速度公式 vg=cn+ωdndω
可测得某一频率下示波器 Ch2
电学参数理论计算(A 部分)
衰减系数和幅度
由电压幅度值每米减少 0.571%≪1
距离/m | 0 | 10 | 20 | 30 |
---|---|---|---|---|
幅度/V | 1.000 | 0.944 | 0.892 | 0.842 |
电缆并联
由图 1 可知,记每一小段 Δx
高低阻抗交替连接
由上问可知,高特征阻抗使用单根电缆,低特征阻抗使用并联电缆,且有
Z0HZ0L=\textcolorblue2
频率计算
由电磁波公式 f=vλ=c√ϵrμrλ
间接测量色散关系(B 部分)
正常光子晶体
测量数据
以 1MHz
测量数据整理如下
f (MHz) | V1pp (V) | V2pp (V) | φ (°) |
---|---|---|---|
1 | 10.250 | 9.220 | -72.10 |
2 | 10.260 | 9.380 | -139.80 |
3 | 10.270 | 8.750 | 147.70 |
4 | 10.280 | 8.730 | 80.80 |
5 | 10.270 | 8.250 | 2.70 |
6 | 10.320 | 7.590 | -66.40 |
7 | 10.480 | 5.890 | -166.90 |
8 | 10.510 | 2.320 | 138.10 |
9 | 10.570 | 1.450 | 112.60 |
10 | 10.640 | 1.360 | 92.10 |
11 | 10.750 | 1.850 | 68.20 |
12 | 10.860 | 4.160 | 27.00 |
13 | 10.980 | 8.000 | -63.80 |
14 | 11.110 | 8.160 | -145.60 |
15 | 11.280 | 8.770 | 140.40 |
16 | 11.450 | 8.700 | 71.10 |
17 | 11.630 | 9.190 | 2.40 |
18 | 11.810 | 8.670 | -67.30 |
19 | 12.060 | 9.340 | -131.90 |
20 | 12.310 | 9.030 | 155.50 |
传输效率和等效折射率
同轴晶体未接入、阻抗匹配状态下 50Ω
f (MHz) | φ' (°) | η | n |
---|---|---|---|
1 | -72.10 | 80.91% | 1.71548 |
2 | -139.80 | 83.58% | 1.66313 |
3 | -212.30 | 72.59% | 1.68375 |
4 | -279.20 | 72.12% | 1.66076 |
5 | -357.30 | 64.53% | 1.70025 |
6 | -426.40 | 54.09% | 1.69089 |
7 | -526.90 | 31.59% | 1.79094 |
8 | -581.90 | 4.87% | 1.73065 |
9 | -607.40 | 1.88% | 1.60577 |
10 | -627.90 | 1.63% | 1.49397 |
11 | -651.80 | 2.96% | 1.40985 |
12 | -693.00 | 14.67% | 1.37405 |
13 | -783.80 | 53.09% | 1.43454 |
14 | -865.60 | 53.95% | 1.47109 |
15 | -939.60 | 60.45% | 1.49040 |
16 | -1008.90 | 57.73% | 1.50030 |
17 | -1077.60 | 62.44% | 1.50820 |
18 | -1147.30 | 53.89% | 1.51654 |
19 | -1211.90 | 59.98% | 1.51762 |
20 | -1284.50 | 53.81% | 1.52811 |
其效率和频率关系如下
等效折射率与频率的变化如下
由图表可得 f=\textcolorblue1,2,3,4,5,13,14,15,16,17,18,19,20MHz
群速度
代入公式 (9) ,其中离散化处理 dndω≈nk+1−nkωk+1−ωk,k=1,⋯,19
f (MHz) | ωdn/dω | vg/c |
---|---|---|
1 | -0.05234 | 0.60127 |
2 | 0.04124 | 0.58673 |
3 | -0.06900 | 0.61929 |
4 | 0.15799 | 0.54983 |
5 | -0.04679 | 0.60479 |
6 | 0.60026 | 0.43646 |
7 | -0.42203 | 0.73051 |
8 | -0.99904 | 1.36686 |
9 | -1.00621 | 1.66789 |
10 | -0.84119 | 1.53193 |
11 | -0.39378 | 0.98418 |
12 | 0.72587 | 0.47621 |
13 | 0.47518 | 0.52364 |
14 | 0.27029 | 0.57426 |
15 | 0.14856 | 0.61015 |
16 | 0.12638 | 0.61475 |
17 | 0.14183 | 0.60605 |
18 | 0.01941 | 0.65106 |
19 | 0.19927 | 0.58245 |
20 | N/A | N/A |
将群速度和频率作图如下
可以在处于反常色散频率时群速度可能超光速(对应 f=8,9,10MHz
缺陷光子晶体
测量数据
在同轴光子晶体中增加一根同轴电缆,形成有缺陷的同轴电缆,按照 4.1 节内容测量数据结果如下
f (MHz) | V1pp (V) | V2pp (V) | φ (°) |
---|---|---|---|
1 | 10.250 | 9.280 | -81.60 |
2 | 10.260 | 9.210 | -163.30 |
3 | 10.270 | 7.860 | 119.80 |
4 | 10.280 | 8.570 | 44.90 |
5 | 10.260 | 8.540 | -45.00 |
6 | 10.310 | 6.110 | -148.40 |
7 | 10.380 | 2.660 | 150.00 |
8 | 10.430 | 1.980 | 118.80 |
9 | 10.490 | 2.950 | 84.30 |
10 | 10.560 | 4.990 | -34.40 |
11 | 10.670 | 2.340 | -90.10 |
12 | 10.790 | 2.340 | -122.70 |
13 | 10.910 | 4.550 | -172.50 |
14 | 11.030 | 7.790 | 96.50 |
15 | 11.190 | 8.470 | 4.30 |
16 | 11.370 | 7.760 | -76.20 |
17 | 11.530 | 8.570 | -151.80 |
18 | 11.720 | 8.490 | 126.90 |
19 | 11.960 | 8.930 | 51.80 |
20 | 12.220 | 8.670 | -27.30 |
传输效率和等效折射率
同理代入公式 (16)(12) 列表如下
f (MHz) | φ' (°) | η | n |
---|---|---|---|
1 | -81.60 | 81.97% | 1.94151 |
2 | -163.30 | 80.58% | 1.94270 |
3 | -240.20 | 58.57% | 1.90503 |
4 | -315.10 | 69.50% | 1.87430 |
5 | -405.00 | 69.28% | 1.92724 |
6 | -508.40 | 35.12% | 2.01606 |
7 | -570.00 | 6.57% | 1.93743 |
8 | -601.20 | 3.60% | 1.78805 |
9 | -635.70 | 7.91% | 1.68058 |
10 | -754.40 | 22.33% | 1.79495 |
11 | -810.10 | 4.81% | 1.75225 |
12 | -842.70 | 4.70% | 1.67087 |
13 | -892.50 | 17.39% | 1.63348 |
14 | -983.50 | 49.88% | 1.67146 |
15 | -1075.70 | 57.29% | 1.70628 |
16 | -1156.20 | 46.58% | 1.71935 |
17 | -1231.80 | 55.25% | 1.72402 |
18 | -1313.10 | 52.48% | 1.73570 |
19 | -1388.20 | 55.75% | 1.73840 |
20 | -1467.30 | 50.34% | 1.74558 |
其效率和频率关系如下,其中透射峰使用 \textcolorred↙
计算表 6和表 3计算的传输效率之比,列表如下
f (MHz) | η1 | η2 | η2/η1 |
---|---|---|---|
1 | 80.91% | 81.97% | 1.013 |
2 | 83.58% | 80.58% | 0.964 |
3 | 72.59% | 58.57% | 0.807 |
4 | 72.12% | 69.50% | 0.964 |
5 | 64.53% | 69.28% | 1.074 |
6 | 54.09% | 35.12% | 0.649 |
7 | 31.59% | 6.57% | 0.208 |
8 | 4.87% | 3.60% | 0.740 |
9 | 1.88% | 7.91% | 4.202 |
10 | 1.63% | 22.33% | 13.667 |
11 | 2.96% | 4.81% | 1.624 |
12 | 14.67% | 4.70% | 0.321 |
13 | 53.09% | 17.39% | 0.328 |
14 | 53.95% | 49.88% | 0.925 |
15 | 60.45% | 57.29% | 0.948 |
16 | 57.73% | 46.58% | 0.807 |
17 | 62.44% | 55.25% | 0.885 |
18 | 53.89% | 52.48% | 0.974 |
19 | 59.98% | 55.75% | 0.929 |
20 | 53.81% | 50.34% | 0.935 |
将两者的传输效率之比绘图如下
可以看出在 f=9,10,11MHz
群速度
同理离散化代入公式 (9) 计算群速度,列表如下
f (MHz) | ωdn/dω | vg/c |
---|---|---|
1 | 0.00119 | 0.51475 |
2 | -0.07534 | 0.53552 |
3 | -0.09220 | 0.55162 |
4 | 0.21176 | 0.47937 |
5 | 0.44414 | 0.42170 |
6 | -0.47178 | 0.64755 |
7 | -1.04570 | 1.12142 |
8 | -0.85972 | 1.07721 |
9 | 1.02929 | 0.36902 |
10 | -0.42698 | 0.73101 |
11 | -0.89521 | 1.16681 |
12 | -0.44859 | 0.81815 |
13 | 0.49371 | 0.47010 |
14 | 0.48744 | 0.46320 |
15 | 0.19600 | 0.52569 |
16 | 0.07474 | 0.55739 |
17 | 0.19867 | 0.52011 |
18 | 0.04846 | 0.56049 |
19 | 0.13645 | 0.53338 |
20 | N/A | N/A |
其群速度和频率的曲线绘制如下
由图像可知,在 f=9MHz
色散关系直接测量(C 部分)
特征频率
在 50 kHz 至 150 kHz
范围内改变信号发生器频率,观察脉冲信号的波形、幅度和延时变化情况。从上到下分别为
114、130、160 kHz
群速度计算
对图片截图计算中心波腹的延迟像素与示波器背景单位格点像素比例,例如对第一张图两者像素为
333 和 543,故其延迟时间为 t1=333 pix543 pix⋅500 ns=306.6~ns
f (MHz) | 基准像素 | 波腹像素 | Δt/ns |
---|---|---|---|
11.4 | 543 | 333 | 306.6 |
13.0 | 493 | 99 | 100.4 |
16.0 | 498 | 164 | 164.7 |
由群速度代表波包的运动,vg=lΔt
f (MHz) | vg/c |
---|---|
11.4 | 0.3807 |
13.0 | 1.1628 |
16.0 | 0.7090 |
将该结果标注在图5中如下,其中 \textcolorred∙
可以发现,除中间点外,测量结果基本符合预期。
传输特性理论计算(D.1−D.4 部分)
叠加电压求解
由公式 (6) 的两个方程反解得 {V+0=12(v2+Z0i2)eγlV−0=12(v2−Z0i2)e−γl
传递矩阵计算
取 x=0
传递矩阵性质
将两个同轴电缆串联连接如下
由上节传递矩阵的性质 [v1i1]=A1[v2i2],[v2i2]=A2[v3i3]
负载电压表达式
令总传递矩阵 AT=[ABCD]
编程计算(D.5 部分)
使用 matlab 对给定数据进行模拟,已知单根同轴电缆的特征阻抗 Z0=50Ω
由连接图可知,令 A1
程序设计
将转移矩阵对应的元素设计在函数中,例如第一个元素获取函数
A1()
如下
1 |
|
主函数针对传输效率、等效折射率、群速度计算,纵坐标标注不同结果
1 |
|
其中使用 unwrap()
函数对相位差有跃变的情况加以处理,自动增加 nπ
正常光子晶体
取定 \textcolorbluef∈[0,8]~MHz
传输效率
从左往右依次为通带、带隙、通带交替产生,并且随着频率地增加,多次反射和透射产生的损耗增加,从而引起效率极值点下降。
等效折射率
可以看出在频率逐渐增大的过程中存在多个反常色散区域,这与反常色散曲线吻合。
群速度
可以看出在反常色散区域群速度超过光速 c
缺陷光子晶体
假定增加的电缆位于左侧,其总传递矩阵更改为 A′T=A2A2A2A1A2A1A2
传输效率
仅仅更改矩阵,其余代码不变,传输效率关系图如下
与图
6对比,在每个通带和带隙之间增加了若干个透射峰(确切来说有 4
等效折射率
群速度
有缺陷的情形均与实际情况吻合,同时符合能带理论。