²é¿´: 578  |  »Ø¸´: 1

jiasen1980

½ð³æ (ÕýʽдÊÖ)

[½»Á÷] ¡¾ÇóÖú¡¿ÄÇλ´óÏÀ°ï¿´ÏÂÄÄÀï³öÎÊÌâÁ˰¡

Po£»%ÊäÈë¹âÇ¿£¬µ¥Î»W
alpha£»%¹âÏËËðºÄÖµ£¬µ¥Î»dB/km
gamma£»%¹âÏË·ÇÏßÐÔ²ÎÊý
to£»%³õʼÂö³å¿í¶È£¬µ¥Î»Ãë
C£»%µÚÒ»´Î¼ÆËãÊäÈëµÄßúà±²ÎÊý
b2£»%²¨ÊýµÄµ¹Êý
cputime=0;
tic;
ln=1;
i=sqrt(1);
pi=3.1415926535;
alph=alpha/(4.343);
Ld=(to^2)/(abs(b2)); %À©É¢³¤¶È£¬µ¥Î»ÊÇm
Ao=sqrt(Po); %¹âÕñ·ù
tau =4096e12:
1e12:
4095e12;
dt=1e12;
h=1000;%²½³¤
for ii=0.1:0.1:1.5 %²»Í¬µÄ¹âÏ˳¤¶È²»Í¬£¬Õâ¸öÁ¿¿É±ä
z=ii*Ld;
u=Ao*exp(((
1+i*(C))/
2)*(tau/to).^2);
figure(1)
plot(abs(u),'r');
title('Input Pulse'); xlabel('Time'); ylabel('Amplitude');
grid on;
hold on;
l=max(size(u));
fwhm1=find(abs(u)>abs(max(u)/2));
fwhm1=length(fwhm1);
dw=1/l/dt*2*pi;
w=(1*
l/2:1:l/21)*
dw;
u=fftshift(u); %ÁãÑÓ³Ù¶ÔÖÐµÄÆ×
w=fftshift(w); %ÁãÑÓ³Ù¶ÔÖÐµÄÆ×
spectrum=fft(fftshift(u)); %¿ìËÙÀëÉ¢¸µÁ¢Ò¶±ä»»
for jj=h:h:z
spectrum=spectrum.*exp(g1) ; %g1ΪÏßÐÔËã·ûeµÄÖ¸Êý±í´ïʽ
f=ifft(spectrum); %¿ìËÙÀëÉ¢·´¸µÁ¢Ò¶±ä»»
f=f.*exp(g2);%g2Ϊ·ÇÏßÐÔËã·ûeµÄÖ¸Êý±í´ïʽ
spectrum=fft(f); %¿ìËÙÀëÉ¢¸µÁ¢Ò¶±ä»»
spectrum=spectrum.*exp(g1) ;
end
f=ifft(spectrum); %¿ìËÙÀëÉ¢·´¸µÁ¢Ò¶±ä»»
op_pulse(ln,=abs(f);%±£´æÔÚËùÓмä¸ôµãÉϵÄÊä³öÂö³å
fwhm=find(abs(f)>abs(max(f)/2));
fwhm=length(fwhm);
ratio=fwhm/fwhm1;
pbratio(ln)=ratio;
dd=atand((abs(imag(f)))/(abs(real(f))));
phadisp(ln)=dd;%±£´æÂö³åÏàλ
ln=ln+1;
end
toc;
cputime=toc;
figure(2);
mesh(op_pulse(1:1:ln1,
title('Pulse Evolution');
xlabel('Time'); ylabel('distance'); zlabel('amplitude');
figure(3)
plot(pbratio(1:1:ln1),'
k');
xlabel('Number of steps');
ylabel('Pulse broadening ratio');
grid on;
hold on;
figure(5)
plot(phadisp(1:1:ln1),'
k');
xlabel('distance travelled');
ylabel('phase change');
grid on;
hold on;
disp('CPU time:'), disp(cputime);
??? Undefined function or variable 'Po'
»Ø¸´´ËÂ¥
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

jiasen1980

½ð³æ (ÕýʽдÊÖ)

лл ´ó¼Ò  ÎÊÌâÒѾ­½â¾ö, ÎÒÊÇ´ÓÒ»¸öÂÛÎÄÉÏ¿´µ½µÄÕâ¸ö´úÂë, ºóÀ´·¢ÏÖÕâ¸öÂÛÎÄÒ²Êdz­±ðÈ˵Ä, ÎÒÕÒµ½ÁËԭʼµÄ´úÂë, ÊÇ¿ÉÒÔÔËÐеÄ.
Çë°ßÖñËøÌû°É
2Â¥2010-09-09 14:16:39
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ jiasen1980 µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 070300»¯Ñ§319Çóµ÷¼Á +5 ½õÀð0909 2026-03-17 5/250 2026-03-18 10:21 by macy2011
[¿¼ÑÐ] 298-Ò»Ö¾Ô¸Öйúũҵ´óѧ-Çóµ÷¼Á +4 ÊÖ»úÓû§ 2026-03-17 4/200 2026-03-18 09:30 by b²»Ò×
[¿¼ÑÐ] 334Çóµ÷¼Á +3 Ö¾´æ¸ßÔ¶ÒâÔÚ»úÐ 2026-03-16 3/150 2026-03-18 08:34 by lm4875102
[¿¼ÑÐ] 268Çóµ÷¼Á +7 ºÃÔËÁ¬Ã಻¾ø 2026-03-12 8/400 2026-03-17 20:28 by xilongliang
[¿¼ÑÐ] 328Çóµ÷¼Á£¬Ó¢ÓïÁù¼¶551£¬ÓпÆÑо­Àú +3 ÉúÎ﹤³Ìµ÷¼Á 2026-03-16 8/400 2026-03-17 19:03 by Wangjingyue
[¿¼ÑÐ] 326Çóµ÷¼Á +5 Éϰ¶µÄСÆÏ 2026-03-15 6/300 2026-03-17 17:26 by ruiyingmiao
[¿¼ÑÐ] 085601Çóµ÷¼Á +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[ÂÛÎÄͶ¸å] ÓÐûÓдóÀз¢Ð¡ÂÛÎÄÄÜ´øÎÒ¸ö¶þ×÷ +3 ÔöÈñ©ÈË 2026-03-17 4/200 2026-03-17 09:26 by xs74101122
[¿¼ÑÐ] 0854¿ØÖƹ¤³Ì 359Çóµ÷¼Á ¿É¿çרҵ +3 626776879 2026-03-14 9/450 2026-03-16 17:42 by 626776879
[¿¼ÑÐ] 321Çóµ÷¼Á +5 ´óÃ×·¹£¡ 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[¿¼ÑÐ] 0703 ÎïÀí»¯Ñ§µ÷¼Á +3 ÎÒ¿ÉÒÔÉϰ¶µÄ¶Ô 2026-03-13 5/250 2026-03-16 10:50 by ÎÒ¿ÉÒÔÉϰ¶µÄ¶ÔÂ
[¿¼ÑÐ] 085601²ÄÁϹ¤³Ì315·ÖÇóµ÷¼Á +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[¿¼ÑÐ] ²ÄÁϹ¤³Ì327Çóµ÷¼Á +3 xiaohe12w 2026-03-11 3/150 2026-03-14 20:20 by ms629
[¿¼ÑÐ] Öпƴó²ÄÁÏÓ뻯¹¤319Çóµ÷¼Á +3 ÃÏöβÄÁÏ 2026-03-14 3/150 2026-03-14 20:10 by ms629
[¿¼ÑÐ] 265Çóµ÷¼Á +4 Íþ»¯±ý07 2026-03-12 4/200 2026-03-14 17:23 by userper
[¿¼ÑÐ] 330Çóµ÷¼Á +3 ?½´¸øµ÷¼Á¹òÁË 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[¿¼ÑÐ] ¹¤¿Æ278·ÖÇóµ÷¼Á +5 ÖÜÂýÈȰ¡ 2026-03-12 7/350 2026-03-13 15:49 by JourneyLucky
[¿¼ÑÐ] 274Çóµ÷¼Á +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
[¿¼ÑÐ] 0856»¯Ñ§¹¤³Ì280·ÖÇóµ÷¼Á +4 shenzxsn 2026-03-11 4/200 2026-03-13 11:55 by ymwdoctor
[¿¼ÑÐ] 289Çóµ÷¼Á +3 ÀîÕþÓ¨ 2026-03-12 3/150 2026-03-13 11:02 by Çóµ÷¼Ázz
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û