| 查看: 704 | 回复: 2 | ||
jiasen1980金虫 (正式写手)
|
[求助]
急求MATLAB解色散方程精度问题
|
|
运行一个有缘色散的超越方程求数值解， 下边是我的程序，用MATLAB7.1 运行只能结出前四个解，for循环到w=5时显示 没有精确解，不知问题出在了那里，求大虫指点迷津 clc; clear; syms b; neff=zeros(1,31); digits(200); for w=1:31 d=vpa(w*5*1e-9); k0=vpa(2*pi/0.8e-6); e1=1; e2=13.425; w z=vpa(solve(exp(sqrt(b-e1*k0^2).*d)-(e2*sqrt(b-e1*k0^2)-e1*sqrt(b-e2*k0^2))/(e2*sqrt(b-e1*k0^2)+e1*sqrt(b-e2*k0^2)),b)); neff(1,w)=vpa(sqrt(z)/k0) end 用算结果显示如下 neff = Columns 1 through 4 3.3286 + 1.6971i 2.2151 + 1.4575i 1.8060 + 1.2174i 1.5881 + 1.0458i Columns 5 through 8 0 0 0 0 Columns 9 through 12 0 0 0 0 Columns 13 through 16 0 0 0 0 Columns 17 through 20 0 0 0 0 Columns 21 through 24 0 0 0 0 Columns 25 through 28 0 0 0 0 Columns 29 through 31 0 0 0 w = 5 Warning: Explicit solution could not be found. > In solve at 140 In sym.solve at 49 ??? Subscripted assignment dimension mismatch. |
回复此楼» 猜你喜欢
【考研调剂】化学专业 281分，一志愿四川大学，诚心求调剂
已经有8人回复
-
一志愿华南师大 070300（化学）304分求调剂
已经有3人回复
-
一志愿 西北大学 ，070300化学学硕，总分287，双非一本，求调剂。
已经有3人回复
-
一志愿重庆大学085700资源与环境专硕，总分308求调剂
已经有3人回复
-
354求调剂
已经有8人回复
-
304求调剂
已经有6人回复
-
一志愿西南交大，求调剂
已经有5人回复
-
299求调剂
已经有5人回复
-
材料专业求调剂
已经有6人回复
-
311求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求matlab求解一个方程组的问题
已经有6人回复
- 求各位前辈帮忙！matlab解复数方程，出现图中的警告，怎来修改呢？谢谢了啊
已经有5人回复
- 求助Matlab解方程问题
已经有6人回复
- MATLAB解方程出现??? Subscripted assignment dimension mismatch.
已经有3人回复
- 急切求助用matlab曲线拟合色散方程
已经有10人回复
- 【求助】MATLAB导入数据精度问题。。。
已经有4人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】matlab 解矩阵方程 重金悬赏解决方案
已经有10人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
xianpingwang 
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 361.1
- 散金: 217
- 红花: 1
- 帖子: 72
- 在线: 94.7小时
- 虫号: 905711
- 注册: 2009-11-17
- 专业: 光学
★
jjdg(金币+1): 感谢参与 2011-06-23 23:03:02
jjdg(金币+1): 感谢参与 2011-06-23 23:03:02
| 最好把方程写出来，要不大家不知道编的是什么程序 |
» 本帖已获得的红花（最新10朵）
jiasen19802楼2011-06-23 12:12:34
jiasen1980
金虫 (正式写手)
- 应助: 12 (小学生)
- 金币: 1391.2
- 散金: 1600
- 红花: 7
- 帖子: 459
- 在线: 494.7小时
- 虫号: 580859
- 注册: 2008-07-18
- 专业: 光学
3楼2011-06-23 17:02:10
10