why not able to plot size of both same

Question

shiv gaur 2022-2-16

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1651375-why-not-able-to-plot-size-of-both-same

回答： Image Analyst 2022-2-16

k0=(2*pi/0.6328)*1e6;
t2=1.5e-6;
n1=1.512;n2=1.521;n3=4.1-1i*0.211;
n4=1;
m=0;
x=2;
t3=linspace(1e-9,1e-6,20)
    
    k1=k0*sqrt(n1^2-x.^2);k2=k0*sqrt(n2^2-x.^2);
    k3=k0*sqrt(n3^2-x^2);
    k4=k0*sqrt(n4^2-x.^2);
    tol = 1e-12; 
    n = 1;
    while (n<=200)
        x_new =-(k2).*t2+atan(k1./1i*k2)+atan((k3./k2).*tan(atan(k4./1i*k2)-k3.*t3))+m*pi;
        if abs(x_new-x) < tol, break; end
        x =x_new;
        k1=k0*sqrt(n1^2-x.^2);
        k2=k0*sqrt(n2^2-x.^2);
        k3=k0*sqrt(n3^2-x.^2);
        k4=k0*sqrt(n4^2-x.^2);
        n=n+1;
    end
    
disp(x_new)
plot(t3,real(x_new))

size is same but not able to plot.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

shiv gaur 2022-2-16

k0=(2*pi/0.6328)*1e6;

t2=1.5e-6;

n1=1.512;n2=1.521;n3=4.1-1i*0.211;

n4=1;

m=0;

x=2;

t3=linspace(1e-9,1e-6,20)

k1=k0*sqrt(n1^2-x.^2);k2=k0*sqrt(n2^2-x.^2);

k3=k0*sqrt(n3^2-x^2);

k4=k0*sqrt(n4^2-x.^2);

tol = 1e-12;

n = 1;

while (n<=200)

x_new =-(k2).*t2+atan(k1./1i*k2)+atan((k3./k2).*tan(atan(k4./1i*k2)-k3.*t3))+m*pi;

if abs(x_new-x) < tol, break; end

x =x_new;

k1=k0*sqrt(n1^2-x.^2);

k2=k0*sqrt(n2^2-x.^2);

k3=k0*sqrt(n3^2-x.^2);

k4=k0*sqrt(n4^2-x.^2);

n=n+1;

end

disp(x_new)

plot(t3,real(x_new))

edit program

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Image Analyst 2022-2-16

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1651375-why-not-able-to-plot-size-of-both-same#answer_897700

在 MATLAB Online 中打开

Break this up into simpler terms before combining them.

x_new =-(k2).*t2+atan(k1./1i*k2)+atan((k3./k2).*tan(atan(k4./1i*k2)-k3.*t3))+m*pi;

It will help you figure out what the problem is.

k0=(2*pi/0.6328)*1e6;
t2=1.5e-6;
n1=1.512;n2=1.521;n3=4.1-1i*0.211;
n4=1;
m=0;
x=2;
t3=linspace(1e-9,1e-6,20)
k1=k0*sqrt(n1^2-x.^2);
k2=k0*sqrt(n2^2-x.^2);
k3=k0*sqrt(n3^2-x^2);
k4=k0*sqrt(n4^2-x.^2);
tol = 1e-12;
n = 1;
while (n<=200)
    term1 =-(k2) .* t2 + atan(k1 ./ 1i*k2);
    term2 = +atan(k3./k2)
    term3 = atan(k4./1i*k2) -k3 .* t3;
    term4 = tan(term3) + m*pi;
    x_new = term1 + term2 .* term4;
    %         x_new =-(k2).*t2+atan(k1./1i*k2)+atan((k3./k2).*tan(atan(k4./1i*k2)-k3.*t3))+m*pi;
    if abs(x_new-x) < tol, break; end
    x =x_new;
    k1=k0*sqrt(n1^2-x.^2);
    k2=k0*sqrt(n2^2-x.^2);
    k3=k0*sqrt(n3^2-x.^2);
    k4=k0*sqrt(n4^2-x.^2);
    n=n+1;
end
disp(x_new)
plot(t3,real(x_new))