| 查看: 333 | 回复: 0 | |||
zdzhang木虫 (正式写手)
|
[求助]
我用这种方法解了这个方程，但是算出来的lambta和预想得到的差别有点大
|
|
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了，都解不了。这个问题我弄了好长时间了，这段时间把握折腾的头晕脑胀了，希望大侠们最好能写一下具体的编程过程，感激不尽！  v=(k_0^2-k_z^2)^0.5,这是我求解的方程 |
回复此楼» 猜你喜欢
存款400万可以在学校里躺平吗
已经有36人回复
-
英文综述是否需要润色及查重
已经有5人回复
-
救命帖
已经有5人回复
-
限项规定
已经有5人回复
-
为什么nbs上溴 没有产物点出现呢
已经有9人回复
-
招博士
已经有3人回复
-
最失望的一年
已经有18人回复
-
求推荐英文EI期刊
已经有5人回复
找到一些相关的精华帖子，希望有用哦~
我的这种情况能否报考执业药师？
已经有9人回复
- 为什么有的样品批不出来东西，DNA检测时挺好的啊
已经有26人回复
- RT-PCR时GAPDH能跑出来，但目的条带跑不出来~~
已经有7人回复
- 高温固相法合成材料高低温前后差别之原因
已经有7人回复
- 能请教一下金相图像与SEM图像的差别
已经有13人回复
- 药学学硕和专硕差别在哪啊
已经有39人回复
- 腐蚀的深浅，组织差别很大的原因？？
已经有17人回复
- matlab解微分方程
已经有10人回复
- 青基申请书字体的问题——win7和xp有差别！
已经有7人回复
- 求绳子的曲线方程
已经有6人回复
- 求matlab 解方程的问题
已经有5人回复
- 正反向引物长度差别较大可用吗？
已经有8人回复
- 用DPM模型时，solution里面会不会有 discrece model方程
已经有4人回复
- Matlab求解偏微分方程组
已经有22人回复
- 分子量300~500的 化学物质，如何透析出来？请点帖详看，得到您的帮助，十分感谢！
已经有1人回复
- fslove解非线性方程组解不出来
已经有3人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人