%观察群速度色散和kell效应同时扰动作用下基态孤子的传输特性clear
clc%?注意时间域与频域的选择问题
T = 20; %时间窗口(-10:10)
tn= 2^10;  % 时间区域个数(2^M个)
dt = T/tn;   %选择的时间步长
t = ((1:tn)'-(tn+1)/2)*dt;
w = wn(T,tn);%傅里叶域频率,调用计算频率函数 m=ran(0,1)%n为0-1之间随机变化数
r=ran(-1,1)%l为-1和1之间随机变化数
%上面两个式子调用了随机数产生函数b=0.03;
c=-i*(1+m*r)*w'.^2/2+i*b*w'.^3;
%步长的选择影响计算效率
dz=0.1;  %迭代时选择的步长
z=0:dz:100;  % 位置坐标
%z=0:dz:50;  % 位置坐标  n=ran(0,1)%n为0-1之间随机变化数
l=ran(-1,1)%l为-1和1之间随机变化数
%上面两个式子调用了随机数产生函数% 输入初始脉冲
u(1,:)=sech(t);  %?注意复数的表示,傅里叶变换函数的使用
%由于输入负频,fft对应逆变换,ifft对应傅里叶变换
for k=1:1:length(z)
% 分布傅里叶变换解方程u1=fft(exp(c*dz/2).*ifft(u(k,:)));u2=u1.*exp(i*dz*(1+n*l)*abs(u(k,:)).^2);u(k+1,:)=fft(exp(c*dz/2).*ifft(u2));%二次迭代for p=1:2u2=u1.*exp(i*dz/2*(1+n*l)*(abs(u(k,:)).^2+abs(u(k+1,:)).^2));u(k+1,:)=fft(exp(c*dz/2).*ifft(u2));end
%缺少迭代仿真将出现很多问题
end subplot(1,2,1);
plot(t(1:1:length(t)),u(1,:))
xlabel('t')
ylabel('u')
title('输入基态孤子u=secht')subplot(1,2,2);
% waterfall 用于画流水线图样,mesh 用于画彩色曲面,mesh计算速度快
%waterfall(t(1:1:length(t)),z(1:(length(z))/50:length(z)),abs(u(1:(length(z))/50:length(z),:)))
mesh(t(1:1:length(t)),z(1:(length(z))/50:length(z)),abs(u(1:(length(z))/50:length(z),:)))
%注意采样点少点好画,
xlabel('t')
ylabel('z')
zlabel('|u|')
title('GVD和KELL同时扰动作用输出光脉冲')% subplot(1,2,2);
%waterfall(t(1:1:length(t)),z(1:(length(z))/50:length(z)),abs(u(1:(length(z))/50:length(z),:)))
%mesh(t(1:1:length(t)),z(1:(length(z))/50:length(z)),abs(u(1:(length(z))/50:length(z),:)))
% %注意采样点少点好画,
% view(-30,45);
% xlabel('t')
% ylabel('z')
% zlabel('|u|')
% title('输出光脉冲旋转一个角度观察两个孤子相互作用')

非线性光纤光学——光孤子相关推荐

  1. 非线性光纤光学——光孤子5

    %在色散位移光纤中研究多孤子的相互作用clear clc%?注意时间域与频域的选择问题 % T = 20; %时间窗口(-10:10) T = 40; tn= 2^10; % 时间区域个数(2^M个) ...

  2. 非线性光纤光学中分步傅里叶算法(SSFFT)的matlab代码实现

    非线性光纤光学中分步傅里叶算法(SSFFT)的matlab代码实现 SSFFT(分布傅里叶算法)函数代码如下: function [waveform, f_spectrum] = SSFFT_arra ...

  3. matlab中taufft,非线性光纤光学中分步傅里叶算法(SSFFT)的matlab代码实现

    非线性光纤光学中分步傅里叶算法(SSFFT)的matlab代码实现 SSFFT(分布傅里叶算法)函数代码如下: function [waveform, f_spectrum] = SSFFT_arra ...

  4. 【非线性光纤光学】,第四章第4题,我的解答,画出KDP晶体的角度调谐曲线

    <非线性光纤光学>第四章第4题,我的解答 作者:最甜甜甜甜的小草莓莓 题目:画出KDP晶体的角度调谐曲线 先声明,不保证就是正确答案!!!! (但是某J大佬觉得ok,小声) (1) 前辈有 ...

  5. 非线性光纤光学_一个锁模光纤激光器,居然也有这么多种类

    锁模光纤激光器原理 激光腔的边界条件决定了激光以一系列分立的谐振模式稳定存在于腔中.依据光的传播方向,这些模式可分成纵模和横模.仅考虑单个横模的情况下,如果诸多纵模间的相位关系被锁定(即锁模),则可以 ...

  6. 非线性光纤光学_北京大学刘忠范院士/刘开辉教授等《自然·纳米技术》:光纤内二维材料的均匀生长实现超高非线性...

    非线性光纤已被广泛应用于光学变频.超快激光和光通信等领域.在目前的制造技术中,非线性是通过将非线性材料注入到纤维或制造微结构纤维等途径来实现的.然而,这两种策略都存在低的光学非线性或设计灵活性差的问题 ...

  7. 非线性光纤光学_深度解析:光纤随机激光器及其应用研究进展!

    本文作者饶云江,电子科技大学信息与通信工程学院,光纤传感与通信教育部重点实验室,光纤光学研究中心,仅供交流学习之用,感谢分享! 引言 自从2010年分布反馈式光纤随机激光器的概念提出以来,光纤随机激光 ...

  8. 非线性光纤光学_1.56 m波段高能量百飞秒光纤激光器

    华东师范大学曾和平教授和上海理工大学郝强副教授带领的研究团队成功研制了参数为1.56μm/ 60nJ/ 100fs/ 10MHz的高稳定性.可靠耐用的超快激光光源. "2019中国光学十大进 ...

  9. 特种光纤所需学习知识(光纤光学)

    光纤光学学习知识点 感谢 第一章 第二章 第三章 第四章 第五章 第六章 第七章 光纤连接与耦合 第八章 简答题: 感谢 第一章 光纤通信的优点 宽带宽,低损耗,保密性好,易铺设 光纤 介质圆柱光波导 ...

  10. 计算机仿真平台的优点包括,光纤通信系统中的计算机仿真

    摘要: 随着通信系统复杂性的增加,传统的分析设计方法已经不能适应发展的需要.如果每次都直接用真实系统进行实验,不仅耗资昂贵,费工费时,有时甚至难于找到问题症结所在.因此,解决上述问题的有效方法是采用计 ...

最新文章

  1. Python发送邮件以及对其封装
  2. ajax将数据显示在class为content的标签中_python selenium:自动化爬取某鱼数据
  3. 函数式编程 - 组合compose
  4. c语言指针字符串逆序存放,C语言指针实现字符串的反转
  5. php 请求url没有扩展名,apache – 将.php重定向到无扩展名,并在没有.php的情况下创建url...
  6. 网络matlab程序_【Matlab】官网资源盘点
  7. 【C++】C++中头文件的使用
  8. 【SICP归纳】6 副作用与环境模型
  9. 什么是Python脚本?
  10. python tornado websocket_Python Tornado实现WEB服务器Socket服务器共存并实现交互的方法...
  11. 项目实施流程和规范模板(测试方向)
  12. Lync 小技巧-8-超速-更新-地址簿
  13. 基于matlab和lingo的数学实验,MATLAB和LINGO软件在数学建模竞赛中的应用-精选教育文档...
  14. html数据透视,高级Excel – 数据透视表工具
  15. Diamond简介和快速应用
  16. ppt太大如何压缩到最小
  17. linux ps 主进程,Linux ps 查看进程
  18. C语言数据结构猜数小游戏
  19. 用matlab抽奖,现在社会上销售彩票的很多。一家三口在抽奖时,常常喜欢让孩子来抽,请问这是遵循了什么决策原则(   )...
  20. bat批处理静默自动安装证书

热门文章

  1. 广义表存储结构图的绘制
  2. batchplot3.5.9如何使用_Flink 1.9 实战:使用 SQL 读取 Kafka 并写入 MySQL
  3. xp系统启动sql服务器,XP系统,sqlserver服务启动了又自己停止,服务管理器点了启动随后马上就已停止!...
  4. GTD工具 Wunderlist使用心得总结
  5. 金士顿优盘突然出现写保护,无法删除更改数据也没有办法格式化
  6. 资源搜索神器 海量精品资源教程 大大啦啦资源神器
  7. 如何设置.net控件SplitContainer平均分配
  8. logo语言是计算机语言吗,LOGO语言的编程
  9. 游戏资源提取软件介绍
  10. 【计算机组成原理】Verilog语言编写32位并行加法器的理解