24小时热门版块排行榜

返回列表

媛小圈爱吃鱼

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 9.5
帖子: 3
在线: 3小时
虫号: 4338281
注册: 2016-01-06
专业: 电磁场与波

[求助] 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

回复此楼