| 查看: 2836 | 回复: 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); |
回复此楼» 猜你喜欢
材料求调剂 一志愿哈工大总分298分，前三科223分
已经有4人回复
-
086502化学工程342求调剂
已经有3人回复
-
求调剂
已经有9人回复
-
352分 化工与材料
已经有5人回复
-
085602 化工专硕 338分 求调剂
已经有10人回复
-
0856材料化工调剂 总分330
已经有10人回复
-
330一志愿中国海洋大学 化学工程 085602 有读博意愿 求调剂
已经有4人回复
-
一志愿哈尔滨工业大学材料与化工方向336分
已经有6人回复
-
材料求调剂一志愿哈工大324
已经有6人回复
-
085600 286分 材料求调剂
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助Ta2O5的色散方程
已经有4人回复
- matlab求解非线性方程组
已经有16人回复
- 急切求助用matlab曲线拟合色散方程
已经有10人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】matlab求解非线性方程组，并画图处理。要求y,z是实数解！
已经有18人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
- 【求助】用MATLAB求解偏微分方程的数值解 悬赏20
已经有8人回复
- 【求助】复超越方程的解法
已经有19人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
小刘8885
木虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 10042.8
- 散金: 200
- 帖子: 468
- 在线: 160.5小时
- 虫号: 735175
- 注册: 2009-03-30
- 性别: MM
- 专业: 材料物理化学
3楼2012-05-29 10:23:19
xcj3630385
木虫 (小有名气)
长江11号
- 应助: 6 (幼儿园)
- 金币: 5104.3
- 帖子: 205
- 在线: 217.1小时
- 虫号: 1359358
- 注册: 2011-08-02
- 性别: GG
- 专业: 金属非晶态、准晶和纳米晶
2楼2012-05-28 23:30:36
4楼2012-05-29 10:56:13
Illusionist
银虫 (正式写手)
- 应助: 20 (小学生)
- 金币: 232.8
- 散金: 1001
- 红花: 15
- 帖子: 431
- 在线: 329.1小时
- 虫号: 1123329
- 注册: 2010-10-15
- 专业: 半导体晶体与薄膜材料
|
还是有提示错误，初值改为1,5,30都不行··· Error in ==> @(B)func(B,K,k0,n0,n1,n2,m,h) Error in ==> fsolve at 254 fuser = feval(funfcn{3},x,varargin{:}); Error in ==> Untitled2 at 9 f=fsolve(@(B)func(B,K,k0,n0,n1,n2,m,h),1); Caused by: Failure in initial user-supplied objective function evaluation. FSOLVE cannot continue. |
5楼2012-05-29 11:34:49
