|
|
[求助]
我用這種方法解了這個方程,但是算出來的lambta和預(yù)想得到的差別有點(diǎn)大
clear;
clc;
epsl_met=-18.36+0.46i;
%mu_met=1;
%epsl_air=1;%8.854187817e-12;
%mu_air=1;%4*pi*10^(-7)
r=60*10^(-9);
wavelength=632.8*10^(-9);
k_0=2*pi/wavelength;
k_1=k_0;
N=1000;
k_z(1)=(0.9+0.0001i)*k_0; % try value
for j=1:N
u=r*sqrt(epsl_met*k_0^2-k_z(j)^2);%we use '+'now
v=-r*sqrt(k_0^2-k_z(j)^2);%we use '+'now
temp=((1/epsl_met)*(besselk(1,v)*besselj(0,u))/(besselk(0,v)*besselj(1,u)))^2;
k_z(j+1)=sqrt(((temp*epsl_met-1)*k_0^2)/(temp-1));
re_k_z(j+1)=real(k_z(j+1));
im_k_z(j+1)=imag(k_z(j+1));
end
%a_met=1-i;
%a_air=a_met*(besselk(0,r*sqrt(epsl_met*k_0^2-k_z(j+1)^2))/besselj(0,r*sqrt(epsl_met*k_0^2-k_z(j+1)^2)))
lmbta=2*pi/real(k_z(j+1))
plot(re_k_z)
哪位能幫我看一下我用這種方法解這樣的方程是否可行?有沒有更好的解這個方程的方法呢?我試過fzero和solver了,都解不了。這個問題我弄了好長時間了,這段時間把握折騰的頭暈?zāi)X脹了,希望大俠們最好能寫一下具體的編程過程,感激不盡!
![]()
v=(k_0^2-k_z^2)^0.5,這是我求解的方程 |
|