屏蔽室的固有谐振频率计算及程序代码
摘要：介绍了屏蔽室的固有谐振频率计算公式，并提供了计算的程序代码，*后给出了一个计算案例。
为了避免屏蔽室谐振引起的测量误差，应通过理论计算和实际测量来获得屏蔽室的主要谐振频率点，并记录下来，以便在今后的测试中，避开这些谐振频率。
谐振频率与屏蔽室的几何尺寸有关，平行六面体结构的屏蔽室谐振频率计算按下面公式计算：
式中：
f0―谐振频率，单位为兆赫兹（MHz)；
m,n,k―为0,1,2,3⋯,m,n,k只能有一个取0；
l,w,h―为屏蔽室的边长，单位为米（m)，且l≥w≥h。
可以通过Excel中的VBA编程功能，在其中录制一个宏，将下面代码拷贝到其中，即可计算谐振频率。
下面是VBA的代码，表中数据是l=9，w=6，h=4条件下的部分数据，图是Excel的计算结果截图。
m | n | k | f/MHz | ||||||
1 | 1 | 0 | 30.05 | 2 | 52.70 | 3 | 76.83 | ||
2 | 41.67 | 60.09 | 82.07 | ||||||
3 | 55.90 | 70.71 | 90.14 | 41.04 | 113.73 | 50.17 | 117.33 | 62.50 | 123.11 |
0 | 45.07 | 79.06 | 115.24 | 83.85 | 106.07 | 135.21 |
Sub f() '计算长方体屏蔽室的谐振频率
Dim m, n, k, i As Integer
Dim l, w, h, f As Double
Dim No As Integer
l = CDbl(Range("B1").Value)
w = CDbl(Range("B2").Value)
h = CDbl(Range("B3").Value)
No = CInt(Range("B4").Value)
If l < w Then
MsgBox ("l需大于等于w,请修改。")
End If
If w < h Then
MsgBox ("w需大于等于h,请修改。")
i = 6 '从第六行开始输入数值
k = 0
For m = 1 To No
For n = 1 To No
f = 150 * ((m / l) ^ 2 + (n / w) ^ 2+ (k / h) ^ 2) ^ 0.5
Range("A" & i).Value = m
Range("B" & i).Value = n
Range("C" & i).Value = k
Range("D" & i).Value = f
i = i + 1
Next
Next
n = 0
For k = 1 To No
m = 0
For n = 1 To No
For k = 1 To No
Range("A" & i).Value =m
Range("B" & i).Value =n
Range("C" & i).Value =k
Range("D" & i).Value =f
i = i + 1
Next
MsgBox ("计算结束。")
End Sub