%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](https://img-mid.csdnimg.cn/release/static/image/mid/ask/652232610136158.png?x-oss-process=image/auto-orient,1/resize,w_320,m_lfit)