-
大小: 2KB文件類型: .m金幣: 1下載: 0 次發(fā)布日期: 2021-05-17
- 語言: Matlab
- 標(biāo)簽: matlab??色散??耦合強(qiáng)度??
資源簡(jiǎn)介
matlab色散補(bǔ)償程序,光譜儀采集到的數(shù)據(jù)進(jìn)行一系列變換,畫出耦合曲線的圖
代碼片段和文件信息
%用傅立葉變換法進(jìn)行數(shù)據(jù)處理的Matlab計(jì)算程序:
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
評(píng)論
共有 條評(píng)論