%用傅立葉變換法進行數(shù)據(jù)處理的Matlab計算程序：
clear;
clc;
clf;
global?Co
Co=299792458;???%light?velocity?constant

%INPUT?%
widt=1400;??????%sample?points
win_ct=790;?????%filter?center
win_wd=30;???????%filter?half?width

x1=1260;?????????%start?wavelength
x2=1400;?????????%end?wavelength


L=750;?????????%fiber?length?unit:m
order=3;????????%Order?of?polynomial?Fit
nFit=20;??????%fit?points
sign=-1;????????%sign?of?dispersion

%Main?Program?%
?[filename?pathname]?=?uigetfile（‘*.txt‘?‘Pick?an?Data?file‘）;
????????if?isequal（filename0）
???????????disp（‘User?selected?Cancel‘）
????????else
???????????disp（[‘User?selected‘?fullfile（pathname?filename）]）
????????end
FilePath=[pathname?filename];
fid?=?fopen（FilePath‘r‘）;
[Acount]?=?fscanf（fid‘%f‘）;
t=length（A）/2;

wavelength=p;
omega1=2*pi*Co./wavelength*10^9;
signal1=q;
l1=length（signal1）;

figure（1）
plot（wavelengthsignal）;
xlabel（‘\fontsize{12}\bfWavelength（nm）‘）;
ylabel（‘\fontsize{12}\bfIntensity（a.u）‘）;

figure（2）
plot（omegasignal）;
xlabel（‘\fontsize{12}\bfomega（nm）‘）;
ylabel（‘\fontsize{12}\bfIntensity（a.u）‘）;

figure（3）
fsi=fft（signal