solution of transcedental equation

1 次查看（过去 30 天）
shiv gaur
shiv gaur 2021-10-14
评论： Matt J 2021-10-14
function beta1 = trial0(beta, eps1, eps2, epsm,k0,h)
S1 = sqrt(power(beta,2)-(eps1*power(k0,2)));
S2 = sqrt(power(beta,2)-(epsm*power(k0,2)));
S3 = sqrt(power(beta,2)-(eps2*power(k0,2)));
beta1 = tanh(S2*h)*((eps1*eps2*power(S2,2))+(power(epsm,2)*S1*S3)) + (S2*((eps1*S3)+(eps2*S1))*epsm);
end
eps1 = 1.5471;
eps2 = 1.5431;
eps1m = -9.894 ;
epsm2 = 1.0458;
epsm = eps1m + i*eps2m;
lambda = 633;
k0 = (2*pi)/lambda;
h = 50;
find the value of beta

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

采纳的回答

Matt J
Matt J 2021-10-14
编辑：Matt J 2021-10-14
Ran in:
eps1 = 1.5471;
eps2 = 1.5431;
eps1m = -9.894 ;
eps2m = 1.0458;
epsm = eps1m + i*eps2m;
lambda = 633;
k0 = (2*pi)/lambda;
h = 50;
[b,fval]=fminsearch(@(b) abs(trial0(complex(b(1),b(2)), eps1, eps2, epsm,k0,h)) ,[1,1]);
beta=complex(b(1),b(2)), fval
beta = -0.0145 - 0.0004i
fval = 2.7735e-05
function beta1 = trial0(beta, eps1, eps2, epsm,k0,h)
S1 = sqrt(power(beta,2)-(eps1*power(k0,2)));
S2 = sqrt(power(beta,2)-(epsm*power(k0,2)));
S3 = sqrt(power(beta,2)-(eps2*power(k0,2)));
beta1 = tanh(S2*h)*((eps1*eps2*power(S2,2))+(power(epsm,2)*S1*S3)) + (S2*((eps1*S3)+(eps2*S1))*epsm);
end
  7 个评论
显示 5更早的评论
shiv gaur
shiv gaur 2021-10-14
relates with this type of equation
Matt J
Matt J 2021-10-14
Well, if we've answered your first question, please Accept-click it and post the related question in its own thread.

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Mathematics 的更多信息

标签

产品


版本

R2021a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by