【正文】 %外電場強度no=sqrt(+(+**T.^2)./(lambda.^2(+*10.^(8)*T.^2).^2)*10.^(2)*lambda.^2)。ne=sqrt(+*10.^(7)*T.^2+(+**T.^2)./(lambda.^2(+*10.^(8)*T.^2).^2)*10.^(2)*lambda.^2)。k=2*pi*(neno)/lambda*+4*pi/(*)。%相位失配量subplot(1,2,1)。plot(T,no),xlabel(39。溫度T(K)39。),ylabel(39。no39。)subplot(1,2,2)。plot(T,ne),xlabel(39。溫度T(K)39。),ylabel(39。ne39。)%(i=1,2,3)的影響：reffi—Tlambda=。T=[270:1:350]。r232=28*。r222=*。r332=**。%LiNbO3 晶體的非零電光系數(shù)ram=4*pi/(*)。%PPLN的極化周期（準(zhǔn)相位匹配階數(shù)m=2）r=*。%晶體長度r=100*ramG0=。G2=i/pi。G20=i/pi。 k0=2*pi/lambda*。%入射o光波矢 E0=*。%外電場強度no=sqrt(+(+**T.^2)./(lambda.^2(+*10.^(8)*T.^2).^2)*10.^(2)*lambda.^2)。ne=sqrt(+*10.^(7)*T.^2+(+**T.^2)./(lambda.^2(+*10.^(8)*T.^2).^2)*10.^(2)*lambda.^2)。k=2*pi*(neno)/lambda*+4*pi/(*)。%相位失配量reff1=ne.^2.*no.^2.*r232。reff2=2*no.^4.*r222。reff3=2*ne.^4.*r332。%有效電光系數(shù)subplot(1,3,1)。plot(T,reff1),xlabel(39。溫度T(K)39。),ylabel(39。reff139。)subplot(1,3,2)。plot(T,reff2),xlabel(39。溫度T(K)39。),ylabel(39。reff239。)subplot(1,3,3)。plot(T,reff3),xlabel(39。溫度T(K)39。),ylabel(39。reff339。)%極化周期數(shù)量n=50時T改變對轉(zhuǎn)換效率η的影響：η—Tlambda=。T=[270::350]。r232=28*。r222=*。r332=**。%LiNbO3 晶體的非零電光系數(shù)ram=4*pi/(*)。%PPLN的極化周期（準(zhǔn)相位匹配階數(shù)m=2）n=50。 %極化周期數(shù)量r=**n。%晶體長度r=100*ramG0=。G2=i/pi。G20=i/pi。k0=2*pi/lambda*。%入射o光波矢 E0=*。%外電場強度no=sqrt(+(+**T.^2)./(lambda.^2(+*10.^(8)*T.^2).^2)*10.^(2)*lambda.^2)。ne=sqrt(+*10.^(7)*T.^2+(+**T.^2)./(lambda.^2(+*10.^(8)*T.^2).^2)*10.^(2)*lambda.^2)。k=2*pi*(neno)/lambda*+4*pi/(*)。%相位失配量reff1=ne.^2.*no.^2.*r232。reff2=2*no.^4.*r222。reff3=2*ne.^4.*r332。%有效電光系數(shù)v1q=reff2./(2*no)*k0*G0*E0。v2q=reff3./(2*ne)*k0*G0*E0。kq=reff1./(2*sqrt(no.*ne))*k0*G2*E0。kq0=reff1./(2*sqrt(no.*ne))*k0*G20*E0。u=*sqrt((k+v1qv2q).^2+4*kq.*kq0)。gama=*(v2qv1qk)。pro1=sqrt((cos(u*r)).^2+(gama./u).^2.*(sin(u*r).^2))。pro2=sqrt((kq0./u).^2.*(sin(u*r)).^2)。xiaolv=abs(pro2./pro1)。plot(T,xiaolv),xlabel(39。溫度T(K)39。),ylabel(39。轉(zhuǎn)換效率η39。)%極化周期數(shù)量n=150時轉(zhuǎn)換效率η—T的關(guān)系lambda=。T=[270:1:350]。r232=28*。r222=*。r332=**。%LiNbO3 晶體的非零電光系數(shù)ram=4*pi/(*)。%PPLN的極化周期（準(zhǔn)相位匹配階數(shù)m=2）n=150。 %極化周期數(shù)量r=**n。%晶體長度r=n*ramG0=。G2=i/pi。G20=i/pi。k0=2*pi/lambda*。%入射o光波矢 E0=*。%外電場強度no=sqrt(+(+**T.^2)./(lambda.^2(+*10.^(8)*T.^2).^2)*10.^(2)*lambda.^2)。ne=sqrt(+*10.^(7)*T.^2+(+**T.^2)./(lambda.^2(+*10.^(8)*T.^2).^2)*10.^(2)*lambda.^2)。k=2*pi*(neno)/lambda*+4*pi/(*)。%相位失配量reff1=ne.^2.*no.^2.*r232。reff2=2*no.^4.*r222。reff3=2*ne.^4.*r332。%有效電光系數(shù)v1q=reff2./(2*no)*k0*G0*E0。v2q=reff3./(2*ne)*k0*G0*E0。kq=reff1./(2*sqrt(no.*ne))*k0*G2*E0。kq0=reff1./(2*sqrt(no.*ne))*k0*G20*E0。u=*sqrt((k+v1qv2q).^2+4*kq.*kq0)。gama=*(v2qv1qk)。pro1=sqrt((cos(u*r)).^2+(gama./u).^2.*(sin(u*r).^2))。pro2=sqrt((kq0./u).^2.*(sin(u*r)).^2)。xiaolv=abs(pro2./pro1)。plot(T,xiaolv),xlabel(39。溫度T(K)39。),ylabel(39。轉(zhuǎn)換效率η39。)%波長λ改變時對轉(zhuǎn)換效率η的影響 lambda1=[1500:1:1650]。lambda=lambda1*。T=300。r232=28*。r222=*。r332=**。ram=4*pi/(*)。r=*。G0=。G2=i/pi。G20=i/pi。k0=2*pi./lambda*。E0=*。no=sqrt(+(+**T.^2)./(lambda.^2(+*10.^(8)*T.^2).^2)*10.^(2)*lambda.^2)。ne=sqrt(+*10.^(7)*T.^2+(+**T.^2)./(lambda.^2(+*10.^(8)*T.^2).^2)*10.^(2)*lambda.^2)。k=2*pi*(neno)./lambda*+4*pi/(*)。reff1=ne.^2.*no.^2.*r232。reff2=2*no.^4.*r222。reff3=2*ne.^4.*r332。v1q=(reff2./(2*no)).*k0*G0*E0。v2q=(reff3./(2*ne)).*k0*G0*E0。kq=(reff1./(2*sqrt(no.*ne))).*k0*G2*E0。kq0=(reff1./(2*sqrt(no.*ne))).*k0*G20*E0。u=*sqrt((k+v1qv2q).^2+4*kq.*kq0)。gama=*(v2qv1qk)。pro1=sqrt((cos(u*r)).^2+(gama./u).^2.*(sin(u*r).^2))。pro2=sqrt((kq0./u).^2.*(sin(u*r)).^2)。xiaolv=abs(pro2./pro1)。plot(lambda,xiaolv),xlabel(39。波長λ39。),ylabel(39。轉(zhuǎn)換效率η39。)%相位失配量Δk改變時對轉(zhuǎn)換效率η的影響 ：η—Δklambda1=1550。lambda=lambda1*。T=300。r232=28*。r222=*。r332=**。ram=4*pi/(*)。r=*。G0=。G2=i/pi。G20=i/pi。k0=2*pi./lambda*。E0=*。no=sqrt(+(+**T.^2)./(lambda.^2(+*10.^(8)*T.^2).^2)*10.^(2)*lambda.^2)。ne=sqrt(+*10.^(7)*T.^2+(+**T.^2)./(lambda.^2(+*10.^(8)*T.^2).^2)*10.^(2)*lambda.^2)。k=[5000:10:5000]。reff1=ne.^2.*no.^2.*r232。reff2=2*no.^4.*r222。reff3=2*ne.^4.*r332。v1q=(reff2./(2*no)).*k0*G0*E0。v2q=(reff3./(2*ne)).*k0*G0*E0。kq=(reff1./(2*sqrt(no.*ne))).*k0*G2*E0。kq0=(reff1./(2*sqrt(no.*ne))).*k0*G20*E0。u=*sqrt((k+v1qv2q).^2+4*kq.*kq0)。gama=*(v2qv1qk)。pro1=sqrt((cos(u*r)).^2+(gama./u).^2.*(sin(u*r).^2))。pro2=sqrt((kq0./u).^2.*(sin(u*r)).^2)。xiaolv=abs(pro2./pro1)。plot(k,xiaolv),xlabel(39。相位失配量Δk39。),ylabel(39。轉(zhuǎn)換效率η39。)%no,ne大小與波長λ的關(guān)系 lambda1=[1500:1:1650]。lambda=lambda1*。T=300。r232=28*。r222=*。r332=**。ram=4*pi/(*)。r=*。G0=。G2=i/pi。G20=i/pi。k0=2*pi./lambda*。E0=*。no=sqrt(+(+**T.^2)./(lambda.^2(+*10.^(8)*T.^2).^2)*10.^(2)*lambda.^2)。ne=sqrt(+*10.^(7)*T.^2+(+**T.^2)./(lambda.^2(+*10.^(8)*T.^2).^2)*10.^(2)*lambda.^2)。k=2*pi*(neno)./lambda*+4*pi/(*)。reff1=ne.^2.*no.^2.*r232。reff2=2*no.^4.*r222。reff3=2*ne.^4.*r332。v1q=(reff2./(2*no)).*k0*G0*E0。v2q=(reff3./(2*ne)).*k0*G0*E0。kq=(reff1./(2*sqrt(no.*ne))).*k0*G2*E0。kq0=(reff1./(2*sqrt(no.*ne))).*k0*G20*E0。u=*sqrt((k+v1qv2q).^2+4*kq.*kq0)。gama=*(v2qv1qk)。pro1=sqrt((cos(u*r)).^2+(gama./u).^2.*(sin(u*r).^2))。pro2=sqrt((kq0./u).^2.*(sin(u*r)).^2)。subplot(1,2,1)。plot(lambda,no),xlabel(39。波長λ39。),ylabel(39。no39。)subplot(1,2,2)。plot(lambda,ne),xlabel(39。波長λ39。),ylabel(39。ne39。)%電場E改變時對轉(zhuǎn)換效率η的影響 ：η—Elambda=。T=300。r232=28*。r222=*。r332=**。ram=4*pi/(*)。r=*。G0=。G2=i/pi。G20=i/pi。k0=2*pi./lambda*。t=[0:1:2000]。E0=**sin(*pi.*t+*pi)。no=sqrt(+(+**T.^2)./(lambda.^2(+*10.^(8)*T.^2).^2)*10.^(2)*lambda.^2)。ne=sqrt(+*10.^(7)*T.^2+(+**T.^2)./(lambda.^2(+*10.^(8)*T.^2).^2)*10.^(2)*lambda.^2)。k=2*pi*(neno)./lambda*+4*pi/(*)。reff1=ne.^2.*no.^2.*r232。reff2=2*no.^4.*r222。reff3=2*ne.^4.*r332。v1q=(reff2./(2*no)).*k0*G0*E0。v2q=(reff3./(2*ne)).*k0*G0*E0。kq=(reff1./(2*sqrt(no.*ne))).*k0*G2*E0。kq0=(reff1./(2*sqrt(no.*ne))).*k0*G20*E0。u=*sqrt((k+v1qv2q).^2+4*kq.*kq0)。gama=*(v2qv1qk)。pro1=sqr