| 查看: 2568 | 回复: 12 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
Illusionist银虫 (正式写手)
|
[求助]
matlab求解色散方程，求助
|
|
|
m=1; h=1.5E-6; n0=1.504; n1=2.2899; n2=1; K0=5.927E+6; pi=3.141592653; K=((K0×n1)^2-B^2)^(1/2) P0=(B^2-(K0×n0)^2)^(1/2) P2=(B^2-(K0×n2)^2)^(1/2) hK=m×pi+artan(P0/K)+artan(P2/K) 精确求解B 我写的代码，但是出现一些问题 clear clc m=1; h=1.5E-6; n0=1.504; n1=2.2899; n2=1; K0=5.927E+6; pi=3.141592653; syms B K=((K0*n1)^2-B.^2).^0.5; P0=(B.^2-(K0*n0)^2).^0.5; P2=(B.^2-(K0*n2)^2).^0.5; ya=m*pi+atan(P0./K)+atan(P2./K); f=fsolve(@(B)ya-K.*h,30); 提示错误信息： ??? Undefined function or method 'full' for input arguments of type 'sym'. Error in ==> trustnleqn at 28 Fvec = full(Fvec); Error in ==> fsolve at 378 [x,FVAL,JACOB,EXITFLAG,OUTPUT,msgData]=... Error in ==> Untitled at 16 f=fsolve(@(B)ya-K.*h,30); |
回复此楼» 猜你喜欢
三甲基碘化亚砜的氧化反应
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
孩子确诊有中度注意力缺陷
已经有12人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
-
论文投稿，期刊推荐
已经有4人回复
-
硕士和导师闹得不愉快
已经有13人回复
Illusionist
银虫 (正式写手)
- 应助: 20 (小学生)
- 金币: 232.8
- 散金: 1001
- 红花: 15
- 帖子: 429
- 在线: 328.1小时
- 虫号: 1123329
- 注册: 2010-10-15
- 专业: 半导体晶体与薄膜材料
13楼2012-05-30 08:00:44
xcj3630385
木虫 (小有名气)
长江11号
- 应助: 6 (幼儿园)
- 金币: 5104.3
- 帖子: 205
- 在线: 217.1小时
- 虫号: 1359358
- 注册: 2011-08-02
- 性别: GG
- 专业: 金属非晶态、准晶和纳米晶
2楼2012-05-28 23:30:36
小刘8885
木虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 10041.3
- 散金: 200
- 帖子: 468
- 在线: 160.5小时
- 虫号: 735175
- 注册: 2009-03-30
- 性别: MM
- 专业: 材料物理化学
3楼2012-05-29 10:23:19
4楼2012-05-29 10:56:13