| 查看: 1101 | 回复: 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 |
回复此楼» 猜你喜欢
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
面上没中，邀请各位路过的虫友分析一下分数
已经有13人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有8人回复
-
两块石头
已经有8人回复
-
初秋的晨风
已经有6人回复
-
学科评审组评审是指会评吗？
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有14人回复
-
广西大学-广州大学招聘博士后 欢迎广大优秀人才！！！
已经有7人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有17人回复
chendequan
铁虫 (小有名气)
- 应助: 152 (高中生)
- 金币: 918.6
- 红花: 17
- 帖子: 262
- 在线: 208.7小时
- 虫号: 994474
- 注册: 2010-04-11
- 专业: 结构工程
2楼2017-11-05 09:27:56