基于MATLAB的拍频仿真

在本文中，我们模拟了光学拍频现象，利用窄带激光和带有边带的激光进行拍频。

具体参数如下:

信号x，为频率100M的激光。

信号x1，为中心频率为100M，且具有三阶边带的激光，FSR（自由谱范围）为1MHz。

信号x2，为中心频率为100M，且具有FSR为20M的一阶边带激光。

拍频可以得到

原理如下

在测量过程中，两束激光的电场复振幅 E1(t)和 E2(t) 同时进人光电探测器时，由于光电探测器响应的是光强，是电场复振幅平方和的时间积分【1】，即

用低带宽的滤波器只能检测到频差信号，即w1 - w2。即实验中常用到的拍频技术，来检测两个频率相差不多的激光的频差。

clc
clear

c =  3e8;
fs = 400e6;                             %采样率
f0 = 100e6;                             %中心频率
fm = 1e6;                               %X1的FSR
fn = 20e6;                              %X2的FSR
N  = 6;                                 %阶数*2
f1 = f0-fm*N/2:fm:f0+fm*N/2;            %X1信号频率
f2 = [f0-fn f0 f0+fn];                  %X2信号频率

T = 1e-3;                                         %时宽1s
n = round(T*fs);                                  %采样点个数
t = linspace(0,T,n);                              %时域横坐标

x = 0.5*sin(2*pi*f2(3)*t);                            %无边带信号 

x1 = 0;
for i = 1:length(f1)
    x1 = x1 + 0.5^abs(i-N/2-1)*sin(2*pi*f1(i)*t);             %x1的3阶边带信号
end

x1 = x1 - 0.2*sin(2*pi*f2(2)*t);

x2  = 0.2*sin(2*pi*f2(1)*t) + 0.6*sin(2*pi*f2(2)*t) + 0.25*sin(2*pi*f2(3)*t);

y  = (x + x1).^2;
y2 = (x2 + x1).^2;
y3 = x2.^2;

X = fft(x./(n),n); %用fft得出离散傅里叶变换
X = X(1:floor(n/2+1));%单边取一半
X1 = fft(x1./(n),n); %用fft得出离散傅里叶变换
X1 = X1(1:floor(n/2+1));%单边取一半
X2 = fft(x2./(n),n); %用fft得出离散傅里叶变换
X2 = X2(1:floor(n/2+1));%单边取一半
Y =  fft(y./(n),n); %用fft得出离散傅里叶变换
Y = Y(1:floor(n/2+1));%单边取一半
Y2 =  fft(y2./(n)，n); %用fft得出离散傅里叶变换
Y2 = Y2(1:floor(n/2+1));%单边取一半
Y3 =  fft(y3./(n)，n); %用fft得出离散傅里叶变换
Y3 = Y3(1:floor(n/2+1));%单边取一半

f = fs*(0:(n/2))/n;%频率刻度
figure(1)
plot(f,abs(X1));%画单侧频谱幅度图
title('x1')

figure(2)
plot(f,abs(X2));%画单侧频谱幅度图
title('x2')


figure(3);
plot(f,abs(Y2),'-b')
xlabel("f/Hz")
ylabel("幅度");

参考文献：

【1】曹翔科, 何耀, 张蓉竹. 两种不同类型独立激光器的拍频实验[J]. Chinese Journal of Lasers, 2009, 36(2): 285-289.