24小时热门版块排行榜    

查看: 2194  |  回复: 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
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

惜梦寻草

新虫 (初入文坛)

引用回帖:
6楼: Originally posted by csgt0 at 2013-01-06 11:12:53
能算,就是不知道图对不对,不对就看你的算法错没错

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; ...

非常感谢!
不过我还有两个问题:
第一,为什么我的程序里lamdas = [1.2:0.01:3.2];lamdass = [2.5:0.01:4.5];修改以后还是出不了图呢?
第二,为什么你写的程序理论上没有问题,但是输出的图形和我要的不一样呢?
7楼2013-01-06 14:39:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
惜梦寻草: 金币+3, 有帮助 2013-01-05 15:56:02
你这问题好多,比如
lamdai2是多少
lamdas和lamdass长度不一样怎么一对一的计算
  lamdai(i_lamdas)=1./(1./lamdap-1./lamdas(i_lamdas));到底是除元素还是除矩阵?我估计除元素
showmethemoney
2楼2013-01-05 10:44:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

惜梦寻草

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by csgt0 at 2013-01-05 10:44:09
你这问题好多,比如
lamdai2是多少
lamdas和lamdass长度不一样怎么一对一的计算
  lamdai(i_lamdas)=1./(1./lamdap-1./lamdas(i_lamdas));到底是除元素还是除矩阵?我估计除元素

lamdai2就是lamdaii,我在编程序的时候疏忽了
lamdas和lamdass需要长度一样才能算吗?这个我不太清楚
lamdai(i_lamdas)=1./(1./lamdap-1./lamdas(i_lamdas));正如你想的,是要除元素。
我知道这个程序的错误很多,我不知道应该怎么算这个,你能帮我写一下能得到下面那张图的程序吗?谢谢啦!
3楼2013-01-05 15:55:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

引用回帖:
3楼: Originally posted by 惜梦寻草 at 2013-01-05 15:55:28
lamdai2就是lamdaii,我在编程序的时候疏忽了
lamdas和lamdass需要长度一样才能算吗?这个我不太清楚
lamdai(i_lamdas)=1./(1./lamdap-1./lamdas(i_lamdas));正如你想的,是要除元素。
我知道这个程序的错误很多 ...

因为你的temp2 = abs(ns./lamdas-nss./lamdass-nii./lamdaii-1/A);
在一个公式里面必须保证lamdas,lamdass,lamdaii都是一样长得
showmethemoney
4楼2013-01-05 16:49:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +27 医学老男孩 2026-08-13 58/2900 2026-08-15 15:57 by GZUZM
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 k0dTPqJtl0jt 2026-08-14 5/250 2026-08-15 15:10 by QXtWNz7PI7wZ
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 7lpolszZVXgi 2026-08-14 3/150 2026-08-15 14:53 by QXtWNz7PI7wZ
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 k0dTPqJtl0jt 2026-08-14 3/150 2026-08-15 12:41 by KxMI1BYBxWX1
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 HFw0lei2R37i 2026-08-14 9/450 2026-08-15 11:33 by KxMI1BYBxWX1
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 k0dTPqJtl0jt 2026-08-14 5/250 2026-08-15 07:09 by 4wMiSEwB6436
[基金申请] filecode +7 cratir 2026-08-14 11/550 2026-08-15 06:19 by 学员8dgXkO
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 7lpolszZVXgi 2026-08-14 5/250 2026-08-15 04:41 by 4wMiSEwB6436
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
信息提示
请填处理意见