| 2 | 1/1 | 返回列表 |
| 查看: 583 | 回復(fù): 1 | |||
jiasen1980金蟲 (正式寫手)
|
[交流]
【求助】那位大俠幫看下哪里出問題了啊
|
|
Po;%輸入光強(qiáng),單位W alpha;%光纖損耗值,單位dB/km gamma;%光纖非線性參數(shù) to;%初始脈沖寬度,單位秒 C;%第一次計算輸入的啁啾參數(shù) b2;%波數(shù)的倒數(shù) cputime=0; tic; ln=1; i=sqrt(1); pi=3.1415926535; alph=alpha/(4.343); Ld=(to^2)/(abs(b2)); %擴(kuò)散長度,單位是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的指數(shù)表達(dá)式 f=ifft(spectrum); %快速離散反傅立葉變換 f=f.*exp(g2);%g2為非線性算符e的指數(shù)表達(dá)式 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' |
金蟲 (正式寫手)
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +6 | 白QF 2026-03-21 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 275求調(diào)劑 +10 | jjjjjjjjjjl 2026-03-27 | 10/500 |
|
|
[有機(jī)交流]
高溫高壓反應(yīng)求助
10+4
|
chibby 2026-03-25 | 4/200 |
|
|
[考研] 化學(xué)調(diào)劑 +4 | 愛吃番茄的旭 2026-03-24 | 5/250 |
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫 2026-03-26 | 5/250 |
|
|
[考研] 一志愿211,335分,0856,求調(diào)劑院校和導(dǎo)師 +4 | 傾____蕭 2026-03-27 | 5/250 |
|
|
[考研]
|
孅華 2026-03-22 | 16/800 |
|
|
[考研] 324求調(diào)劑 +5 | hanamiko 2026-03-26 | 5/250 |
|
|
[考研]
|
鐘llll 2026-03-26 | 4/200 |
|
|
[考研]
|
平樂樂樂 2026-03-26 | 4/200 |
|
|
[考研] 機(jī)械學(xué)碩310分,數(shù)一英一,一志愿211本科雙非找調(diào)劑信息 +3 | @357 2026-03-25 | 3/150 |
|
|
[考研]
|
WWW西西弗斯 2026-03-24 | 8/400 |
|
|
[考研] 一志愿天津大學(xué)339材料與化工求調(diào)劑 +3 | 江往賣魚 2026-03-26 | 3/150 |
|
|
[考研] 07化學(xué)303求調(diào)劑 +5 | 睿08 2026-03-25 | 5/250 |
|
|
[考研] 293求調(diào)劑 +7 | 加一一九 2026-03-24 | 7/350 |
|
|
[考研] 0854電子信息求調(diào)劑 324 +4 | Promise-jyl 2026-03-23 | 4/200 |
|
|
[考研] 求調(diào)劑 +6 | 研研,接電話 2026-03-24 | 7/350 |
|
|
[考研] 材料/農(nóng)業(yè)專業(yè),07/08開頭均可,過線就行 +3 | 呵唔哦豁 2026-03-23 | 4/200 |
|
|
[考研] 276求調(diào)劑 +3 | YNRYG 2026-03-21 | 4/200 |
|
|
[考研] 一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑 +5 | vv迷 2026-03-21 | 8/400 |
|