| 查看: 1023 | 回复: 1 | ||
[求助]
matlab求解色散方程，非线性方程 已有1人参与
|
|
求教如何计算色散方程，自己写的代码如下，但是求解不出来，不知道什么原因，希望有大神来帮忙，这里的y是一个数组，求解出的x也应该是一个数组，但是程序报错： %%%输入参数%%%% f=1e+12; q=1.602e-19; hb=6.626e-34/(2*pi); T=300; kb=1.38e-23; tao=1e-12; uo=1.257e-6; eo=8.854e-12; %%自由空间介电常数 ed=3.9; %%二氧化硅相对介电常数 ea=1; %%空气相对介电常数 Vf=10e+6; Vdc=20; %%%衬底厚度%%%% H=0.1e-6:0.1e-6:20e-6; mu=hb.*Vf.*((pi.*eo.*ed.*Vdc./(q.*H)).^(0.5)); %%%求解化学势 ko=2.*pi.*f/3e+8; %%%自由空间波数 y=-1i*q.*q.*kb.*T./(pi.*hb.*hb.*(2.*pi.*f-1i/tao)).*(mu./(kb.*T)+2.*log(exp(-mu./(kb.*T))+1)); %%面电导率 x0=zeros(1,200); x=fsolve(@(x)func(x,ea,eo,ed,ko,y,f),x0); function fb=func(x,ea,eo,ed,ko,y,f) ya=1i*ko*y./2*pi*f; p1=ea*eo/((x.^2-ea)^0.5); p2=ed*eo/((x.^2-ed)^0.5); yaa=p1+p2; fb=yaa-ya; end |
回复此楼
chendequan
铁虫 (小有名气)
- 应助: 152 (高中生)
- 金币: 918.6
- 红花: 17
- 帖子: 262
- 在线: 208.6小时
- 虫号: 994474
- 注册: 2010-04-11
- 专业: 结构工程
2楼2017-11-05 09:27:56
20