主题:奇怪的图片?程序有错?请高手帮忙!
原始公式:
y=-2arctan[tan(x)*(1-r)/(1+r)]
x=2π/λ*d*cos(θ)
其中d=1.55mm,r1=0.1,r2=0.3
程序为:
>> r1=0.1;r2=0.3;
>> x=-2*pi:2*pi;
>> R1=(1-r1)/(1+r1);
>> s=cos(x);
>> M1=R1*tan(s);
>> y1=-2*atan(M1);
>> R2=(1-r2)/(1+r2);
>> M2=R2*tan(s);
>> y2=-2*atan(M2);
>> plot(x,y1,'r',x,y2,'g',x,y1,'+',x,y2,'*')
>> xlabel('Δ')
>> ylabel('Θ(Δ)')
>> gtext('r1=0.1'),gtext('r2=0.3')
上传图片中,图1为我用上述程序做出的图,图2为要求的图,不知道是哪里出错了,会差这么多,请高手帮忙看看,谢谢了哈!
y=-2arctan[tan(x)*(1-r)/(1+r)]
x=2π/λ*d*cos(θ)
其中d=1.55mm,r1=0.1,r2=0.3
程序为:
>> r1=0.1;r2=0.3;
>> x=-2*pi:2*pi;
>> R1=(1-r1)/(1+r1);
>> s=cos(x);
>> M1=R1*tan(s);
>> y1=-2*atan(M1);
>> R2=(1-r2)/(1+r2);
>> M2=R2*tan(s);
>> y2=-2*atan(M2);
>> plot(x,y1,'r',x,y2,'g',x,y1,'+',x,y2,'*')
>> xlabel('Δ')
>> ylabel('Θ(Δ)')
>> gtext('r1=0.1'),gtext('r2=0.3')
上传图片中,图1为我用上述程序做出的图,图2为要求的图,不知道是哪里出错了,会差这么多,请高手帮忙看看,谢谢了哈!