请问相干光通信色散补偿补偿算法怎么写？


 %D   ----->  色散系数
%lm  ----->  中心波长
%z  ----->   传输距离
%T  ----->    采样周期
%c  ----->   光速
%n  ----->   抽头变量
%N  ----->   抽头个数
%wN ----->   Nyquist 频率
Fs=80e9;
z=1000e3;
D=16e-6;
lm=1550.57/1e9;
c=3e8;
T=1/Fs;
wN=pi/T;
N=2*floor(D*lm^2*z/(2*c*T^2))+1;
n=-floor(N/2):1:floor(N/2);
bk=exp(-j*pi*c*T^2.*n.*n/(D*lm^2*z))*sqrt(j*c*T^2/(D*lm^2*z));
%画滤波器频谱响应图
freqz(bk,1);
fvtool(bk,1); 
%[H,w]=freqz(ak,1,1000000);
%Hf=abs(H);
%Hx=angle(H);
%clf
%figure(1)
%plot(w,Hf)
%title('离散系统幅频特性曲线')
%figure(2)
%plot(w,Hx)
%title('离散系统相频特性曲线')
end

写的滤波器，其不能补偿色散，求做过的分析下

img