| 查看: 1613 | 回复: 9 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
MATLAB程序调试
|
|||
|
i=1; %设定从1开始计数 lamda=0.0019; LengthCrystal=14; RefractiveIndex=1.44; l1=5; h= LengthCrystal/2/ RefractiveIndex; l2=10; l31=20; l32=40; l33=60; d1=l1+h; d2=l2+h; Radius=300; FundamentalModeRadius11=zeros(1,10000); %设定第一基模有10000个矩阵元，初值都为0 FundamentalModeRadius12=zeros(1,10000); FundamentalModeRadius13=zeros(1,10000); FundamentalModeRadius21=zeros(1,10000); %设定第二基模有10000个矩阵元，初值都为0 FundamentalModeRadius22=zeros(1,10000); FundamentalModeRadius23=zeros(1,10000); y=zeros(1,10000); for f=1:0.1:1000 M1=[1,d1;0,1]*[1,0;0,1]*[1,d1;0,1]*[1-h/f,2*h-h*h/f;-1/f,1-h/f]*[1,d2;0,1]*... [1,0;0,1]*[1,l31;0,1]*[1,0;-2/ Radius,1]*[1,l31;0,1]*[1,0;0,1]*[1,d2;0,1]*... [1-h/f,2*h-h*h/f;-1/f,1-h/f]; A1=M1(1,1); %第一行第一列矩阵元 B1=M1(1,2); %第一行第二列矩阵元 D1=M1(2,2); %第二行第二列矩阵元 M2=[1,d1;0,1]*[1,0;0,1]*[1,d1;0,1]*[1-h/f,2*h-h*h/f;-1/f,1-h/f]*[1,d2;0,1]*... [1,0;0,1]*[1,l32;0,1]*[1,0;-2/ Radius,1]*[1,l32;0,1]*[1,0;0,1]*[1,d2;0,1]*... [1-h/f,2*h-h*h/f;-1/f,1-h/f]; %第一基模的ABCD矩阵 A2=M2(1,1); %第一行第一列矩阵元 B2=M2(1,2); %第一行第二列矩阵元 D2=M2(2,2); %第二行第二列矩阵元 M3=[1,d1;0,1]*[1,0;0,1]*[1,d1;0,1]*[1-h/f,2*h-h*h/f;-1/f,1-h/f]*[1,d2;0,1]*... [1,0;0,1]*[1,l33;0,1]*[1,0;-2/ Radius,1]*[1,l33;0,1]*[1,0;0,1]*[1,d2;0,1]*... [1-h/f,2*h-h*h/f;-1/f,1-h/f]; %第一基模的ABCD矩阵 A3=M3(1,1); %第一行第一列矩阵元 B3=M3(1,2); %第一行第二列矩阵元 D3=M3(2,2); %第二行第二列矩阵元 N1=[1,d2;0,1]*[1,0;0,1]*[1,l31;0,1]*[1,0;-2/ Radius,1]*[1,l31;0,1]*[1,0;0,1]*... [1,d2;0,1]*[1-h/f,2*h-h*h/f;-1/f,1-h/f]*[1,d1;0,1]*[1,0;0,1]*... [1,d1;0,1]*[1-h/f,2*h-h*h/f;-1/f,1-h/f]; a1=N1(1,1); %第一行第一列矩阵元 b1=N1(1,2); %第一行第二列矩阵元 d1=N1(2,2); %第二行第二列矩阵元 N2=[1,d2;0,1]*[1,0;0,1]*[1,l32;0,1]*[1,0;-2/ Radius,1]*[1,l32;0,1]*[1,0;0,1]*... [1,d2;0,1]*[1-h/f,2*h-h*h/f;-1/f,1-h/f]*[1,d1;0,1]*[1,0;0,1]*... [1,d1;0,1]*[1-h/f,2*h-h*h/f;-1/f,1-h/f]; a2=N2(1,1); %第一行第一列矩阵元 b2=N2(1,2); %第一行第二列矩阵元 d2=N2(2,2); %第二行第二列矩阵元 N3=[1,d2;0,1]*[1,0;0,1]*[1,l33;0,1]*[1,0;-2/ Radius,1]*[1,l33;0,1]*[1,0;0,1]*... [1,d2;0,1]*[1-h/f,2*h-h*h/f;-1/f,1-h/f]*[1,d1;0,1]*[1,0;0,1]*... [1,d1;0,1]*[1-h/f,2*h-h*h/f;-1/f,1-h/f]; a3=N3(1,1); %第一行第一列矩阵元 b3=N3(1,2); %第一行第二列矩阵元 d3=N3(2,2); %第二行第二列矩阵元 if abs((A1+D1)/2)<=1&abs((A2+D2)/2)<=1&abs((A3+D3)/2)<=1 FundamentalModeRadius11(i)=abs((2*lamda*B1)/(pi*sqrt(4-(A1+D1)^2)))^(1/2); FundamentalModeRadius21(i)=abs((2*lamda*b1)/(pi*sqrt(4-(a1+d1)^2)))^(1/2); FundamentalModeRadius12(i)=abs((2*lamda*B2)/(pi*sqrt(4-(A2+D2)^2)))^(1/2); FundamentalModeRadius22(i)=abs((2*lamda*b2)/(pi*sqrt(4-(a2+d2)^2)))^(1/2); FundamentalModeRadius13(i)=abs((2*lamda*B3)/(pi*sqrt(4-(A3+D3)^2)))^(1/2); FundamentalModeRadius23(i)=abs((2*lamda*b3)/(pi*sqrt(4-(a3+d3)^2)))^(1/2); y(i)=f; i=i+1; end end FundamentalModeRadius11n= FundamentalModeRadius11(1:i-1); FundamentalModeRadius21n= FundamentalModeRadius21(1:i-1); FundamentalModeRadius12n= FundamentalModeRadius12(1:i-1); FundamentalModeRadius22n= FundamentalModeRadius22(1:i-1); FundamentalModeRadius13n= FundamentalModeRadius13(1:i-1); FundamentalModeRadius23n= FundamentalModeRadius23(1:i-1); yn=y(1:i-1); wp=0.4*ones(1,i-1) plot(yn, FundamentalModeRadius11n,yn, FundamentalModeRadius21n,yn,FundamentalModeRadius12n,':',yn,FundamentalModeRadius22n,':',yn,FundamentalModeRadius13n,'-.',yn, FundamentalModeRadius23n,'-.',yn,wp) 运行结果什么也没有 求助高手调试 |
回复此楼» 收录本帖的淘帖专辑推荐
 Better One |
» 猜你喜欢
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有287人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有120人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有70人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有5人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 各位matlab编程时，一般怎么调试的？
已经有8人回复
- matlab调试问题
已经有3人回复
- 求帮忙啊~ matlab一多项式结果太长，怎么写到txt里？
已经有14人回复
- matlab调试问题
已经有3人回复
- while语句死循环了，求教
已经有6人回复
- 上海高等研究院（筹）信息科学与技术研究部招聘启事
已经有8人回复
- 求助matlab编程
已经有13人回复
- 【求助】VS2005数值计算程序调试【急求！！！！加金币】
已经有9人回复
9楼2012-11-28 11:25:04
2楼2012-11-27 10:34:16
3楼2012-11-27 11:01:23
4楼2012-11-27 11:08:32
