| 查看: 1016 | 回复: 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 |
回复此楼» 猜你喜欢
本人考085602 化学工程 专硕
已经有23人回复
-
求调剂院校信息
已经有4人回复
-
085600材料与化工306
已经有4人回复
-
286求调剂
已经有10人回复
-
328求调剂，英语六级551，有科研经历
已经有9人回复
-
一志愿北京化工大学070300 学硕336求调剂
已经有4人回复
-
286分人工智能专业请求调剂愿意跨考！
已经有8人回复
-
资源与环境 调剂申请(333分)
已经有5人回复
-
280求调剂
已经有12人回复
-
269专硕求调剂
已经有5人回复
chendequan
铁虫 (小有名气)
- 应助: 152 (高中生)
- 金币: 918.6
- 红花: 17
- 帖子: 262
- 在线: 208.6小时
- 虫号: 994474
- 注册: 2010-04-11
- 专业: 结构工程
2楼2017-11-05 09:27:56
