| 1 | 1/1 | 返回列表 |
| 查看: 2052 | 回復(fù): 0 | ||
Albert Feng新蟲 (初入文壇)
|
[求助]
matlab非線性方程組數(shù)值求解(負(fù)雙軸晶體非共線相位匹配角計(jì)算)
|
|
本人要計(jì)算BiBO晶體的非共線的相位匹配角,物理模型為一束泵浦光產(chǎn)生滿足相位條件會(huì)產(chǎn)生兩束子光束,其有一定的夾角。求泵浦光在晶體中的相位匹配角度。 本人主方程 clear all clc lambdap=0.39 %um lambdai=0.78 lambdas=0.78 nxp=f(lambdap) nxi=f(lambdai) nxs=f(lambdai) nyp=g(lambdap) nyi=g(lambdai) nys=g(lambdas) nzp=h(lambdap) nzi=h(lambdai) nzs=h(lambdas) alphaphi=0 %3/180*pi syms alpha beta thetap phip thetas phis thetai phii thetapie np=k(thetap,phip,lambdap); ns=l(thetas,phis,lambdas); ni=l(thetai,phii,lambdai); eq0=thetap-pi/2-thetapie eq1=phip-pi/2 eq2=2*np-cos(alpha)*ns-cos(beta)*ni eq3=sin(alpha)*ns-sin(beta)*ni eq4=ns*sin(alpha)-sin(alphaphi) eq5=cos(thetas)+sin(thetapie)*cos(alpha) eq6=cos(thetai)+sin(thetapie)*cos(beta) eq7=tan(phii)*sin(beta)-cos(thetapie)*cos(beta) eq8=tan(phis)*sin(alpha)+cos(alpha)*cos(thetapie) [alpha,beta,thetap,phip,thetas,phis,thetai,phii,thetapie]=solve(eq0,eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8,alpha,beta,thetap,phip,thetas,phis,thetai,phii,thetapie) 調(diào)用函數(shù): function f=f(x) %BiBO nx f=sqrt(3.0740+0.0323/(x^2-0.0316)-0.01337*x^2); end %nx function g=g(x) %BIBO ny g=sqrt(3.1685+0.0373/(x^2-0.0346)-0.01750*x^2); end function h=h(x) %BIBO nz h=sqrt(3.6545+0.0511/(x^2-0.0371)-0.0226*x^2); end function k=k(theta,phi,lambda) kx=sin(theta)*cos(phi) ; ky=sin(theta)*sin(phi); kz=cos(theta); nx=f(lambda); ny=g(lambda); nz=h(lambda); a=nx^-2; b=ny^-2; c=nz^-2; B=-(sin(theta)*cos(phi))^2*(b+c)-(sin(theta)*sin(phi))^2*(a+c)-(cos(theta))^2*(a+b); C=(sin(theta)*cos(phi))^2*b*c+(sin(theta)*sin(phi))^2*a*c+(cos(theta))^2*a*b; k=sqrt(2)/sqrt(-B-sqrt(B^2-4*C)) %慢光 function l=l(theta,phi,lambda) kx=sin(theta)*cos(phi) ; ky=sin(theta)*sin(phi); kz=cos(theta); nx=f(lambda); ny=g(lambda); nz=h(lambda); a=nx^-2; b=ny^-2; c=nz^-2; B=-(sin(theta)*cos(phi))^2*(b+c)-(sin(theta)*sin(phi))^2*(a+c)-(cos(theta))^2*(a+b); C=(sin(theta)*cos(phi))^2*b*c+(sin(theta)*sin(phi))^2*a*c+(cos(theta))^2*a*b; l=sqrt(2)/sqrt(-B+sqrt(B^2-4*C)) %e2-快光 解除來無解 可能不能用solve |

| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求b區(qū)院校調(diào)劑 +3 | 周56 2026-03-24 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 材料調(diào)劑 +3 | iwinso 2026-03-23 | 3/150 |
|
|
[考研] 278求調(diào)劑 +4 | 我可以上岸的對(duì)?/a> 2026-03-19 | 4/200 |
|
|
[考研] 292求調(diào)劑 +4 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 4/200 |
|
|
[考研] 07化學(xué)280分求調(diào)劑 +7 | 722865 2026-03-23 | 7/350 |
|
|
[考研] 化學(xué)調(diào)劑 +6 | yzysaa 2026-03-21 | 6/300 |
|
|
[考研] 300分,材料,求調(diào)劑,英一數(shù)二 +5 | 超贊的 2026-03-24 | 5/250 |
|
|
[有機(jī)交流]
20+3
|
FENGSHUJEI 2026-03-23 | 5/250 |
|
|
[考研] 材料專碩找調(diào)劑 +5 | 哈哈哈吼吼吼哈 2026-03-23 | 5/250 |
|
|
[考研]
|
pk3725069 2026-03-19 | 17/850 |
|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +10 | 紙魚ly 2026-03-21 | 11/550 |
|
|
[考研] 277分求調(diào)劑,跨調(diào)材料 +3 | 考研調(diào)劑lxh 2026-03-24 | 3/150 |
|
|
[考研] 一志愿吉大化學(xué)322求調(diào)劑 +4 | 17501029541 2026-03-23 | 6/300 |
|
|
[考博] 26申博自薦 +3 | whh869393 2026-03-24 | 3/150 |
|
|
[考研] 327求調(diào)劑 +5 | prayer13 2026-03-23 | 5/250 |
|
|
[考研] 361求調(diào)劑 +3 | Glack 2026-03-22 | 3/150 |
|
|
[考研] 336求調(diào)劑 +4 | 收到VS 2026-03-20 | 4/200 |
|
|
[考研] 321求調(diào)劑 +9 | 何潤(rùn)采123 2026-03-18 | 11/550 |
|
|
[考研]
|
然11 2026-03-19 | 4/200 |
|
|
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|