写一段基于matlab的仿真光学频率梳的代码

代码如下：

% 仿真光学频率梳

lambda = 1.55; % 波长(μm) delta = 0.1; % 频率间隔(THz)

% 生成频率点 f = (-2:delta:2)*1e12;

% 计算振幅 A = exp(-0.5*(f/1e12).^2).(2pi*f/lambda).^2;

% 画出振幅图 plot(f,A); xlabel('频率(THz)'); ylabel('振幅');