%矢量

format long

clear

n1=1.45

n2=1.447

a=4.5

k=(a^2)*(n1^2-n2^2);

u1=(a^2)*(n1^2);

w1=(a^2)*(n2^2);

XiShu1=(n1^2+n2^2)/(2*n1^2);

XiShu2=(n1^2-n2^2)/(2*n1^2);

i=1;

n=n2;

for V=0:0.001:1

k02=V^2/k;

k01=sqrt(k02);

for BeiTa=n:0.000001:n1

%n

nn=n+0.000001

U2=u1*k02-(a^2)*k02*(BeiTa^2);

U=sqrt(U2)

W2=(a^2)*k02*(BeiTa^2)-w1*k02;

W=sqrt(W2)

if(U==0 || W==0)

break;

%disp('sss')

else

z1=besselj(0,U)/(U*besselj(1,U));

z2=1/(U^2)+XiShu1*(1/W^2+besselk(0,W)/(W*besselk(1,W)))-sqrt((XiShu2*(1/W^2+besselk(0,W)/(W*besselk(1,W))))^2+(BeiTa/n1)^2*(V^2/(U2*W2))^2);

if(abs(z1-z2)<0.01)

x(i)=V;

y(i)=BeiTa;

i=i+1;

z1;

z2;

%disp('nnn')

n=BeiTa;

break;

end

end

end

end

for V=1:0.01:6

k02=V^2/k;

k01=sqrt(k02);

for BeiTa=n:0.00001:n1

%n

U2=u1*k02-(a^2)*k02*(BeiTa^2);

U=sqrt(U2);

W2=(a^2)*k02*(BeiTa^2)-w1*k02;

W=sqrt(W2);

if(U==0 || W==0)

break;

else

z1=besselj(0,U)/(U*besselj(1,U));

z2=1/(U^2)+XiShu1*(1/W^2+besselk(0,W)/(W*besselk(1,W)))-sqrt((XiShu2*(1/W^2+besselk(0,W)/(W*besselk(1,W))))^2+(BeiTa/n1)^2*(V^2/(U2*W2))^2);

if(abs(z1-z2)<0.01)

x(i)=V;

y(i)=BeiTa;

i=i+1;

z1;

z2;

n=BeiTa;

break;

end

end

end

end

plot(x,y)

save SL_V.mat x

save SL_N.mat y

运行结果:

这个程序是自己学习用的,程序中有一些中间参数仅为了调试方便,这里未做删改。最后的运行结果也没有添加任何坐标说明。这是单模光纤特性分析中的第一部分,待全部做完后,写论文时再做整理。

矢量的模 matlab,单模光纤基模色散曲线矢量分析的Matlab实现相关推荐

  1. 异步电动机直接转矩matlab仿真,异步电机直接转矩控制原理MATLAB仿真.docx

    异步电机直接转矩控制原理MATLAB仿真.docx 异步电机直接转矩控制原理MATLAB仿真分析引言直接转矩控制变频调速技术是近年来继矢量控制变频调速技术之后发展起来的一种新型的具有高性能的交流变频调 ...

  2. matlab 课程设计循环码性能分析,matlab课程设计--循环码的性能分析.docx

    matlab课程设计--循环码的性能分析.docx 课程设计任务书学生姓名 专业班级 指导教师 工作单位 题目 循环码的性能分析 初始条件 MATLAB,速率为100Bd的矩形输入信号,AWGN信道要 ...

  3. matlab中的abc-dq相坐标变换,MATLAB中的abc-dq相坐标变换

    坐标变换总结 姓名: 日期:2011.11.4 坐标变换的总结 一. 由三项坐标系变换到两相旋转坐标系 1. 三相到两相静止坐标系的变换 首先,确定三相电压的相序: 在坐标图上表示三相到两相静止坐标系 ...

  4. matlab里面板有什么作用,MATLAB轻松享受GPU的强大功能

    MATLAB轻松享受GPU的强大功能 MATLAB的GPU支持为活跃于许多学科的大量研究人员(不一定是CUDA编程专家)提供了一种加速科学计算的新方法.考虑到MATLAB主要是用于科学计算和工程计算, ...

  5. matlab 坐标不用科学计数法,matlab不用科学计数法

    『壹』 matlab中怎么才能不是科学计数法表示结果.比如1.0e+003 * 2.7581,怎么使它显示为2758.1谢谢了,很急啊 format long (小数位14) 或 format sho ...

  6. matlab噪音的消除办法,基于MATLAB的噪声消除方法.ppt

    基于MATLAB的噪声消除方法 基于MATLAB的噪声消除方法 答 辩 人 :徐 苏 美 指导教师:杨卫平教授 论文的指导思想 21世纪的社会是信息化社会,我们生活中的每一天都离不开数字信号.随着信息 ...

  7. MATLAB 学习笔记(3)MATLAB 矩阵的进阶操作

    目录 MATLAB 矩阵标量操作 实际例子 MATLAB 矩阵的转置 实际例子 MATLAB 串联矩阵 实际例子 MATLAB 矩阵的行列式 MATLAB 逆矩阵 详细例子 MATLAB 矩阵标量操作 ...

  8. 在matlab中ATB是什么意思,Matlab—线性代数入门

    <Matlab-线性代数入门>由会员分享,可在线阅读,更多相关<Matlab-线性代数入门(48页珍藏版)>请在人人文库网上搜索. 1.线性代数Matlab入门,作者刘进生太原 ...

  9. 如何用matlab画电场强度三维图,浅析matlab模拟静电场三维图与恒定电流场模拟静电场...

    浅析matlab模拟静电场三维图与恒定电流场模拟静电场 1 浅析 Matlab 模拟静电场三维图 与恒定电流场模拟静电场 摘要:本文利用恒定电流场模拟静电场,通过寻找等势点的方法描绘出点 电荷对的电力 ...

  10. matlab制作动画并保存,[转载]Matlab中动画的实现、制作和保存(转载Matlab技术论坛)...

    Matlab中动画的实现.制作和保存(转载Matlab技术论坛) 本帖由Matlab技术论坛原创,原帖参见 http://www.matlabsky.com/thread-592-1-1.html 光 ...

最新文章

  1. 摄像头标定GML Camera Calibration
  2. 从换脸到换姿势,AI 在图像处理的道路上越走越魔幻
  3. MySQL事务隔离级别介绍
  4. 重庆两江新区将建国内最大“云计算”数据基地
  5. android 开发 命名规范
  6. c语言求成绩标准差,C程序计算标准偏差
  7. .NET Framework 4.8预览
  8. 摄像头分辨率怎么调整_网络监控摄像头怎么选择 网络监控摄像头选择方法【介绍】...
  9. c语言开发实际,21实际c语言教程-1 (转)
  10. Elasticsearch 磁盘使用率超过警戒水位线,怎么办?
  11. CC2530单片机入门学习
  12. 51单片机用三种方法实现流水灯(超详细)
  13. 安装第三方库,出现“error: Microsoft Visual C++ 14.0 is required. Get it with “Microsoft Visual……”解决方案
  14. MySQL 简单查询语句练习题(4+18)
  15. 硕士毕业论文框架怎么搭建?
  16. 家用汽车蓄电池亏电解决方案
  17. beats耳机红白交替闪烁三次_beats耳机红白灯交替闪怎么解决
  18. Exchange Server 2013 运维系列——恢复已删除邮件
  19. latex 跳转标签_LaTeX 之 \label 的运用
  20. 【ArcGIS Pro二次开发】(15):用地用海名称和代码互转

热门文章

  1. 华为交换机链路聚合使用ENSP模拟器进行实验
  2. 机械臂动力学建模(2)- Kane凯恩算法
  3. 知道焊缝长度如何确定节点板尺寸_钢结构节点板(钢结构节点板尺寸如何选取)...
  4. 几款实用Chrome插件
  5. Redis 实战案例总结
  6. postsql时间计算
  7. postSQL hash分表
  8. 中美线径对照表_美标线径与国标对照及导线介绍
  9. 如何设置视频的帧率、分辨率、码率来调节音视频即时通讯的质量
  10. 复数基2 DIT FFT程序