| 2 | 1/1 | 返回列表 |
| 查看: 582 | 回復(fù): 1 | |||
jiasen1980金蟲 (正式寫手)
|
[交流]
【求助】那位大俠幫看下哪里出問題了啊
|
|
Po;%輸入光強(qiáng),單位W alpha;%光纖損耗值,單位dB/km gamma;%光纖非線性參數(shù) to;%初始脈沖寬度,單位秒 C;%第一次計(jì)算輸入的啁啾參數(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 %不同的光纖長度不同,這個(gè)量可變 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); %零延遲對(duì)中的譜 w=fftshift(w); %零延遲對(duì)中的譜 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);%保存在所有間隔點(diǎn)上的輸出脈沖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ā)表 | |
|---|---|---|---|---|
|
[考研] 085602 化工專碩 338分 求調(diào)劑 +6 | 路癡小琪 2026-03-27 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 0856調(diào)劑 +5 | 求求讓我有書讀?/a> 2026-03-26 | 6/300 |
|
|
[考研] 考研調(diào)劑 +10 | 呼呼?~+123456 2026-03-24 | 10/500 |
|
|
[考研] 333求調(diào)劑 +3 | question挽風(fēng) 2026-03-23 | 3/150 |
|
|
[考研] 考研292求調(diào)劑 +3 | YIXIAO789 2026-03-26 | 3/150 |
|
|
[考研] 一志愿吉大071010,316分求調(diào)劑 +3 | xgbiknn 2026-03-27 | 3/150 |
|
|
[考研] 一志愿陜師大生物學(xué)071000,298分,求調(diào)劑 +5 | SYA! 2026-03-23 | 5/250 |
|
|
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +4 | 柳溪邊 2026-03-26 | 4/200 |
|
|
[考研]
|
平樂樂樂 2026-03-26 | 4/200 |
|
|
[考研] 一志愿北京化工大學(xué)材料與化工(085600)296求調(diào)劑 +9 | 稻妻小編 2026-03-26 | 9/450 |
|
|
[考研] 尋找調(diào)劑 +5 | 倔強(qiáng)芒? 2026-03-21 | 8/400 |
|
|
[考研] 一志愿國科過程所081700,274求調(diào)劑 +3 | 三水研0水立方 2026-03-23 | 3/150 |
|
|
[考研] 336求調(diào)劑 +4 | 收到VS 2026-03-20 | 4/200 |
|
|
[論文投稿] 急發(fā)核心期刊論文 +3 | 賢達(dá)問津 2026-03-23 | 5/250 |
|
|
[考研] 280分求調(diào)劑 一志愿085802 +4 | PUMPT 2026-03-22 | 7/350 |
|
|
[考研] 石河子大學(xué)(211、雙一流)碩博研究生長期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
[考研] 求調(diào)劑 +5 | Zhangbod 2026-03-21 | 7/350 |
|
|
[考研] 一志愿華中科技大學(xué)071000,求調(diào)劑 +4 | 沿岸有貝殼6 2026-03-21 | 4/200 |
|
|
[考研] 一志愿深大,0703化學(xué),總分302,求調(diào)劑 +4 | 七月-七七 2026-03-21 | 4/200 |
|