| 查看: 1932 | 回复: 10 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
利用matlab画出级联OPO波长调谐曲线
|
||
|
求高手指点！如题，模拟PPLN晶体级联OPO+OPA温度调谐曲线，但是不知道该怎么模拟！要得到最后面的那个图。 下面是我的程序，希望哪位高人能帮忙调一下！十分感激！ clear clc n=1; m=1; a1=5.756; a2=0.0983; a3=0.2020; a4=189.32; a5=12.52; a6=1.32e-2; b1=2.860e-6; b2=4.700e-8; b3=6.113e-8; b4=1.516e-4; A=31.5; lamdap=1.064; T = [25:0.1:250]; lamdas = [1.065:0.01:3.2]; lamdass = [2.5:0.01:4.5]; for i_T = 1:length(T) f=(T(i_T)-24.5)*(T(i_T)+570.82); ni = zeros(size(lamdas)); np = ni; ns = ni; nss = ni; nii = ni; for i_lamdas = 1:length(lamdas) lamdai(i_lamdas)=1./(1./lamdap-1./lamdas(i_lamdas)); ni(i_lamdas)=(a1+b1.*f+(a2+b2.*f)./(lamdai(i_lamdas).^2-(a3+b3.*f)^2)+(a4+b4.*f)./(lamdai(i_lamdas).^2-a5^2)-a6.*lamdai(i_lamdas).^2).^(1/2); np(i_lamdas)=(a1+b1.*f+(a2+b2.*f)./(lamdap.^2-(a3+b3*f)^2)+(a4+b4.*f)./(lamdap.^2-a5^2)-a6.*lamdap.^2).^(1/2); ns(i_lamdas)=(a1+b1.*f+(a2+b2.*f)./(lamdas(i_lamdas).^2-(a3+b3.*f)^2)+(a4+b4.*f)./(lamdas(i_lamdas).^2-a5^2)-a6.*lamdas(i_lamdas).^2).^(1/2); end for i_lamdass = 1:length(lamdass) lamdaii(i_lamdass)=1/(1/lamdas-1/lamdass(i_lamdass)); nii(i_lamdass)=(a1+b1.*f+(a2+b2.*f)./(lamdai2(i_lamdass).^2-(a3+b3.*f)^2)+(a4+b4.*f)./(lamdai2(i_lamdass).^2-a5^2)-a6.*lamdai2(i_lamdass).^2).^(1/2); ns(i_lamdass)=(a1+b1.*f+(a2+b2.*f)./(lamdas(i_lamdass).^2-(a3+b3.*f)^2)+(a4+b4.*f)./(lamdas(i_lamdass).^2-a5^2)-a6.*lamdas(i_lamdass).^2).^(1/2); nss(i_lamdass)=(a1+b1.*f+(a2+b2.*f)./(lamdass(i_lamdass).^2-(a3+b3.*f)^2)+(a4+b4.*f)./(lamdass(i_lamdass).^2-a5^2)-a6.*lamdass(i_lamdass).^2).^(1/2); end temp1 = abs(np./lamdap-ns./lamdas-ni./lamdai-1/A); [x1,i_x1] = min(temp1); lamdai_opt(i_T) = lamdai(i_x1); lamdas_opt(i_T) = lamdas(i_x1); temp2 = abs(ns./lamdas-nss./lamdass-nii./lamdaii-1/A); [x2,i_x2] = min(temp2); lamdaii_opt(i_T) = lamdaii(i_x2); lamdass_opt(i_T) = lamdass(i_x2); end plot(T,lamdai_opt,'b.') hold on plot(T,lamdas_opt,'b.') hold on plot(T,lamdaii_opt,'b.') hold on plot(T,lamdass_opt,'b.')  级联OPO+OPA温度调谐曲线.jpg |
回复此楼» 猜你喜欢
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有177人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有0人回复
-
北工大宋晓艳教授团队招收2026年9月入学博士【学博和专博名额均有】！
已经有1人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收10人
已经有0人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收调剂，10人
已经有0人回复
-
计算与实验科学工程国际会议(8月7-11日,香港) 二维材料分会S5与声学换能器分会S8
已经有0人回复
8楼2013-01-06 14:41:30
2楼2013-01-05 10:44:09
3楼2013-01-05 15:55:28
4楼2013-01-05 16:49:51
