Question

MATLAB仿真图出问题了？

此图为我想要得到的图，但以程序没有实线
clc
clear
P=5.2786e-007;
L1=4e-3;
h1=12e-3;
v=1;
neff=1.468155;
bolongD=2*neff*P;
bolong=[1549.5:0.001:1550.5]*1e-9;
detan=4e-5;%n调制
for j=1:1001
Bolong=1549.5*1e-9+(j-1)*0.001*1e-9;
dbeta=2*pi*neff*(1/Bolong-1/bolongD);%失谐量
beta=dbeta+2*pi*detan/bolongD;%总的直流自耦合系数
k=pi*detan*v/bolongD;
q=sqrt(k^2-beta^2);
r1(j)=(sinh(L1*q))^2/((cosh(L1*q))^2-(beta^2)/(k^2));
f_p1(j)=2*r1(j)*(1+cos(4*pi*neff*(L1+h1)/Bolong));
end
h=figure
plot(bolong,f_p1,'r-');
%axis ([1549.5*1e-9,1550.5*1e-9,0,1]);
set(gca,'XTickLabel',{'1549.5' '1550.0' '1550.5' });
set(gca,'yTickLabel',{'0' '0.1' '0.2' '0.3' '0.4' '0.5' '0.6' '0.7' '0.8' '0.9' '1.0'});
legend(['L1=',num2str(1e003*L1),'mm']);
xlabel(' \lambda/nm ')
ylabel('P/Reflectivity')
grid on
P=5.2786e-007;
L2=5e-3;
h2=10e-3;
v=1;
neff=1.468095;
bolongD=2*neff*P;
bolong=[1549.5:0.001:1550.5]*1e-9;
detan=1e-4;%n调制
for j=1:1001
Bolong=1549.5*1e-9+(j-1)*0.001*1e-9;
dbeta=2*pi*neff*(1/Bolong-1/bolongD);%失谐量
beta=dbeta+2*pi*detan/bolongD;%总的直流自耦合系数
k=pi*detan*v/bolongD;
q=sqrt(k^2-beta^2);
r2(j)=(sinh(L2*q))^2/((cosh(L2*q))^2-(beta^2)/(k^2));
f_p2(j)=2*r2(j)*(1+cos(4*pi*neff*(L2+h2)/Bolong));
end
figure(h);hold on;
plot(bolong,f_p2,'b--');
%axis ([1549.5*1e-9,1550.5*1e-9,0,1]);
legend(['L1=',num2str(1e003*L1),'mm'],['L2=',num2str(1e003*L2),'mm']);
xlabel(' \lambda/nm ')
ylabel('P/Reflectivity')
grid on
最后的图片是我想要的，我没得出来，上传的我变的程序，麻烦大侠帮解决，事后还有感谢？

Answer 1

这个好像是的啊。。

追问

幅值不一样 告诉我你ＱＱ　我详细跟你说

Answer 2

你没得出来，图哪里来的？？我帮你调试一下

追问

这个图是论文中的，是我想要仿真成那个样子，我的程序却得不出来。

追答

又是个学术不端的，哎。。。我帮你看看吧

追问

这倒不是，我是刚学，试着练习仿真，谢谢了！

Answer 3

z这是天线阵列吧