 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 1128 | 回復(fù): 1 | |||
hywdx510木蟲 (小有名氣)
|
[交流]
【討論】光纖中光信號(hào)非線性效應(yīng)的傳輸方程仿真的有關(guān)討論 已有1人參與
|
|
Agrawal在其所著Nonlinear Fiber Optics（4ed）Appendix B中給出了該方程一個(gè)利用split-step解法Matlab程序。程序中按理用快速傅里葉變換fft處都用了反變換ifft，而應(yīng)該用反變換ifft處都用了正變換fft（在程序中用黑體標(biāo)注）。 討論問(wèn)題： 1. 這樣反著用對(duì)結(jié)果有無(wú)影響？ 2. 在理論層面分析為什么這樣用？% This code solves the NLS equation with the split-step method % idu/dz - sgn(beta2)/2 d^2u/d(tau)^2 + N^2*|u|^2*u = 0 % Written by Govind P. Agrawal in March 2005 for the NLFO book %---Specify input parameters clear all; % distance = input(’Enter fiber length (in units of L_D) =’); % beta2 = input(’dispersion: 1 for normal, -1 for anomalous’);% N = input(’Nonlinear parameter N = ’); % soliton order mshape = input(’m = 0 for sech, m > 0 for super-Gaussian = ’); chirp0 = 0; % input pulse chirp (default value) %---set simulation parameters nt = 1024; Tmax = 32; % FFT points and window size step_num = round(20*distance*N^2); % No. of z steps to deltaz = distance/step_num; % step size in z dtau = (2*Tmax)/nt; % step size in tau %---tau and omega arrays tau = (-nt/2:nt/2-1)*dtau; % temporal grid omega = (pi/Tmax)*[(0:nt/2-1) (-nt/2:-1)]; % frequency grid %---Input Field profile if mshape==0 uu = sech(tau).*exp(-0.5i*chirp0*tau.^2); % soliton else % super-Gaussian uu = exp(-0.5*(1+1i*chirp0).*tau.^(2*mshape)); end %---Plot input pulse shape and spectrum temp = fftshift(ifft(uu)).*(nt*dtau)/sqrt(2*pi); % spectrum figure; subplot(2,1,1); plot(tau, abs(uu).^2, ’--k’); hold on; axis([-5 5 0 inf]); xlabel(’Normalized Time’); ylabel(’Normalized Power’); title(’Input and Output Pulse Shape and Spectrum’); subplot(2,1,2); plot(fftshift(omega)/(2*pi), abs(temp).^2, ’--k’); hold on; axis([-.5 .5 0 inf]); xlabel(’Normalized Frequency’); ylabel(’Spectral Power’); %---store dispersive phase shifts to speedup code dispersion = exp(i*0.5*beta2*omega.^2*deltaz); % phase factor hhz = 1i*N^2*deltaz; % nonlinear phase factor %*********[ Beginning of MAIN Loop]*********** % scheme: 1/2N -> D -> 1/2N; first half step nonlinear temp = uu.*exp(abs(uu).^2.*hhz/2); % note hhz/2 for n=1:step_num f_temp = ifft(temp).*dispersion; uu = fft(f_temp); temp = uu.*exp(abs(uu).^2.*hhz); end uu = temp.*exp(-abs(uu).^2.*hhz/2); % Final field temp = fftshift(ifft(uu)).*(nt*dtau)/sqrt(2*pi); %Final spectrum %***************[ End of MAIN Loop ]************** %----Plot output pulse shape and spectrum subplot(2,1,1) plot(tau, abs(uu).^2, ’-k’) subplot(2,1,2) plot(fftshift(omega)/(2*pi), abs(temp).^2, ’-k’) 歡迎大家討論！ 我試過(guò)，在循環(huán)迭代過(guò)程中，交換FFT和IFFT對(duì)結(jié)果沒(méi)任何影響，說(shuō)明FFT和IFFT只是成對(duì)出現(xiàn)時(shí)隨便那個(gè)都可以定義為FFT，另一個(gè)就成了IFFT |
回復(fù)此樓» 猜你喜歡
南京大學(xué)化學(xué)調(diào)劑
已經(jīng)有11人回復(fù)
-
材料專碩調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿華東師范大學(xué)有機(jī)化學(xué)專業(yè)，初試351分，復(fù)試被刷求調(diào)劑!
已經(jīng)有5人回復(fù)
-
283求調(diào)劑（080500）
已經(jīng)有7人回復(fù)
-
311求調(diào)劑
已經(jīng)有6人回復(fù)
-
275求調(diào)劑
已經(jīng)有20人回復(fù)
-
各位老師好，我的一志愿為北京科技大學(xué)085601材料專碩
已經(jīng)有9人回復(fù)
-
294分080500材料科學(xué)與工程求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿北京工業(yè)大學(xué)，324分求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿武漢理工，總分321，英一數(shù)二，求老師收留。
已經(jīng)有11人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 石墨烯復(fù)合材料種的協(xié)同效應(yīng)（synergistic effect）怎么理解？
已經(jīng)有8人回復(fù)
- 有關(guān)磁流體方程計(jì)算的詢問(wèn)
已經(jīng)有8人回復(fù)
- 管理應(yīng)該知道的效應(yīng)
已經(jīng)有42人回復(fù)
- 關(guān)于流體沖刷壁面的模擬，及壁面磨損的后處理
已經(jīng)有9人回復(fù)
- matlab擬合方程參數(shù)時(shí)初值的選擇
已經(jīng)有15人回復(fù)
- 求助這個(gè)非線性方程改怎么解��？？
已經(jīng)有4人回復(fù)
- matlab求解非線性方程組
已經(jīng)有16人回復(fù)
- matlab的fsove 命令求解非線性方程組
已經(jīng)有6人回復(fù)
- 【資源】數(shù)值解 光纖中 非線性薛定諤方程 軟件
已經(jīng)有27人回復(fù)
hywdx510
木蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 4515.5
- 帖子: 92
- 在線: 69.3小時(shí)
- 蟲號(hào): 865943
- 注冊(cè): 2009-10-08
- 專業(yè): 光子與光電子器件
2樓2010-12-23 17:17:17
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 hywdx510 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 【求調(diào)劑】085601材料工程專碩 | 總分272 | +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 材料與化工272求調(diào)劑 +14 | 阿斯蒂芬2004 2026-03-28 | 14/700 |
|
|
|
[考研] 085601材料工程找調(diào)劑 +13 | oatmealR 2026-03-29 | 14/700 |
|
|
|
[考研] 299求調(diào)劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
|
[考研] 調(diào)劑 +3 | 好好讀書。 2026-03-28 | 3/150 |
|
|
|
[考研] 材料277求調(diào)劑 +7 | min3 2026-03-24 | 7/350 |
|
|
|
[考研] 291求調(diào)劑 +15 | hhhhxn.. 2026-03-23 | 21/1050 |
|
|
|
[考研] 328求調(diào)劑 +7 | 嗯滴的基本都 2026-03-27 | 7/350 |
|
|
|
[考研] 一志愿上海理工能源動(dòng)力（085800）310分求調(diào)劑 +3 | zhangmingc 2026-03-27 | 4/200 |
|
|
|
[考研] 266求調(diào)劑 +11 | 陽(yáng)陽(yáng)哇塞 2026-03-27 | 12/600 |
|
|
|
[考研] 322求調(diào)劑 +4 | 我真的很想學(xué)習(xí) 2026-03-23 | 4/200 |
|
|
|
[考研] 0703化學(xué)一志愿南京師范大學(xué)303求調(diào)劑 +3 | zzffylgg 2026-03-24 | 3/150 |
|
|
|
[考研] 一志愿吉大071010，316分求調(diào)劑 +3 | xgbiknn 2026-03-27 | 3/150 |
|
|
|
[考研] 材料調(diào)劑 +8 | 匹克i 2026-03-23 | 8/400 |
|
|
|
[考研] 333求調(diào)劑 +6 | wfh030413@ 2026-03-23 | 6/300 |
|
|
|
[考研] 327求調(diào)劑 +7 | prayer13 2026-03-23 | 7/350 |
|
|
|
[考研] 一志愿天津大學(xué)339材料與化工求調(diào)劑 +3 | 江往賣魚 2026-03-26 | 3/150 |
|
|
|
[考研] 機(jī)械學(xué)碩總分317求調(diào)劑！�。。� +4 | Acaciad 2026-03-25 | 4/200 |
|
|
|
[考研] 生物技術(shù)與工程 +3 | 1294608413 2026-03-25 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +6 | 研研，接電話 2026-03-24 | 7/350 |
|