| 查看: 1073 | 回复: 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 |
回复此楼» 猜你喜欢
关于filecode，很负责任的告诉大家
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
-
好奇怪的filecode
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
chendequan
铁虫 (小有名气)
- 应助: 152 (高中生)
- 金币: 918.6
- 红花: 17
- 帖子: 262
- 在线: 208.7小时
- 虫号: 994474
- 注册: 2010-04-11
- 专业: 结构工程
2楼2017-11-05 09:27:56
