Index exceeds the number of array elements. Index must not exceed 1. Error in shiv7 (line 7) da = D(j);

Question

shiv gaur 2022-1-7

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1624030-index-exceeds-the-number-of-array-elements-index-must-not-exceed-1-error-in-shiv7-line-7-da

评论： shiv gaur 2022-1-7

function shiv7

%K0= (2*pi/1530)*1e6:(2*pi/10)*1e6:(2*pi/1630)*1e6;

%K0= linspace((2*pi/1530)*1e6,(2*pi/1630)*1e6,5);

D=linspace(0,300,10);

for j=1:length(D)

da = D(j);

p0 = 1;

p1 = 1.5;

p2 = 2;

TOL = 10^-8;

N0 = 100; format long

h1 = p1 - p0;

h2 = p2 - p1;

DELTA1 = (f(p1,da) - f(p0,da ))/h1;

DELTA2 = (f(p2,da ) - f(p1,da ))/h2;

d = (DELTA2 - DELTA1)/(h2 + h1);

i=3;

while i <= N0

b = DELTA2 + h2*d;

D = (b^2 - 4*f(p2,da )*d)^(1/2);

if abs(b-D) < abs(b+D)

E = b + D;

else

E = b - D;

end

h = -2*f(p2,da)/E;

p = p2 + h;

if abs(h) < TOL

disp(p)

break

end

p0 = p1;

p1 = p2;

p2 = p;

h1 = p1 - p0;

h2 = p2 - p1;

DELTA1 = (f(p1,da ) - f(p0,da ))/h1;

DELTA2 = (f(p2,da ) - f(p1,da ))/h2;

d = (DELTA2 - DELTA1)/(h2 + h1);

i=i+1;

end

P(j)=real(p);

%R(j)=real(r);

end

plot(D,P)

end

function y=f(kx,da )

l=6330*(10)^(-10);

k0=2*pi/l;

nc=1.33;

ec=nc^2;

na=1.43;

ea=na^2;

nf=1.59;

ef=nf^2;

nm=0.064+1i*4;

em=nm^2;

ns=1.495;

es=ns^2;

df=0.5e-6;

dm=0.03e-6;

%da=0.06e-6;

%kx=(2*pi/l)*np*sin(a);

kc=sqrt(kx^2-ec*k0^2);

ka=sqrt(ea*k0^2-kx^2);

kf=sqrt(ef*k0^2-kx^2);

km=sqrt(em*k0^2-kx^2);

ks=sqrt(kx^2-es*k0^2);

rfm=(kf-km)/(kf+km);

rms=(km-ks)/(ks+km);

rfa=(kf-ka)/(kf+ka);

rac=(ka-kc)/(ka+kc);

a=(1-rfm)/(1+rfm);

b=(1-rms*exp(2*1i*dm*km))/(1+rms*exp(2*1i*dm*km));

c=(1-rfa)/(1+rfa);

f=(1-rac*exp(2*1i*da*ka))/(1+rac*exp(2*1i*da*ka));

pfms=2*atan(1i*a*b);

pfac=2*atan(1i*c*f);

y=2*kf*df-pfms-pfac;

end

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

shiv gaur 2022-1-7

this is the problem ie showing

请先登录，再进行评论。

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

Answer 1

Walter Roberson 2022-1-7

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1624030-index-exceeds-the-number-of-array-elements-index-must-not-exceed-1-error-in-shiv7-line-7-da#answer_869870

在 MATLAB Online 中打开

D = (b^2 - 4*f(p2,da )*d)^(1/2);

That line is overwriting your vector D

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

shiv gaur 2022-1-7

function shiv7

%K0= (2*pi/1530)*1e6:(2*pi/10)*1e6:(2*pi/1630)*1e6;

%K0= linspace((2*pi/1530)*1e6,(2*pi/1630)*1e6,5);

D1=linspace(1,300,10);

for j=1:numel(D1)

da = D1(j);

p0 = 1;

p1 = 1.5;

p2 = 2;

TOL = 10^-8;

N0 = 100; format long

h1 = p1 - p0;

h2 = p2 - p1;

DELTA1 = (f(p1,da) - f(p0,da ))/h1;

DELTA2 = (f(p2,da ) - f(p1,da ))/h2;

d = (DELTA2 - DELTA1)/(h2 + h1);

i=3;

while i <= N0

b = DELTA2 + h2*d;

D = (b^2 - 4*f(p2,da )*d)^(1/2);

if abs(b-D) < abs(b+D)

E = b + D;

else

E = b - D;

end

h = -2*f(p2,da)/E;

p = p2 + h;

if abs(h) < TOL

disp(p)

break

end

p0 = p1;

p1 = p2;

p2 = p;

h1 = p1 - p0;

h2 = p2 - p1;

DELTA1 = (f(p1,da ) - f(p0,da ))/h1;

DELTA2 = (f(p2,da ) - f(p1,da ))/h2;

d = (DELTA2 - DELTA1)/(h2 + h1);

i=i+1;

end

P(j)=real(p);

%R(j)=real(r);

end

plot(D1,P)

end

function y=f(kx,da )

l=6330*(10)^(-10);

k0=2*pi/l;

nc=1.33;

ec=nc^2;

na=1.43;

ea=na^2;

nf=1.59;

ef=nf^2;

nm=0.064+1i*4;

em=nm^2;

ns=1.495;

es=ns^2;

df=0.5e-6;

dm=0.03e-6;

%da=0.06e-6;

%kx=(2*pi/l)*np*sin(a);

kc=sqrt(kx^2-ec*k0^2);

ka=sqrt(ea*k0^2-kx^2);

kf=sqrt(ef*k0^2-kx^2);

km=sqrt(em*k0^2-kx^2);

ks=sqrt(kx^2-es*k0^2);

rfm=(kf-km)/(kf+km);

rms=(km-ks)/(ks+km);

rfa=(kf-ka)/(kf+ka);

rac=(ka-kc)/(ka+kc);

a=(1-rfm)/(1+rfm);

b=(1-rms*exp(2*1i*dm*km))/(1+rms*exp(2*1i*dm*km));

c=(1-rfa)/(1+rfa);

f=(1-rac*exp(2*1i*da*ka))/(1+rac*exp(2*1i*da*ka));

pfms=2*atan(1i*a*b);

pfac=2*atan(1i*c*f);

y=2*kf*df-pfms-pfac;

end

shiv gaur 2022-1-7

but plot between da and kx is not showing

shiv gaur 2022-1-7

please help to plot the change in kx/change in na vs da

请先登录，再进行评论。

Index exceeds the number of array elements. Index must not exceed 1. Error in shiv7 (line 7) da = D(j);

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

回答（1 个）

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

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

Index exceeds the number of array elements. Index must not exceed 1. Error in shiv7 (line 7) da = D(j);

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

回答（1 个）

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

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

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

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