| 查看: 1115 | 回复: 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 |
回复此楼» 猜你喜欢
40+中年人失业后，还能找到好工作吗？
已经有13人回复
-
用NBS溴化时溶剂除了CCl4外还有什么
已经有6人回复
-
申请科研项目纯靠实力能成功吗？
已经有11人回复
-
极限拉扯，如何抉择
已经有15人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
chendequan
铁虫 (小有名气)
- 应助: 152 (高中生)
- 金币: 918.6
- 红花: 17
- 帖子: 262
- 在线: 209小时
- 虫号: 994474
- 注册: 2010-04-11
- 专业: 结构工程
2楼2017-11-05 09:27:56
200