| 查看: 756 | 回复: 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. |
回复此楼» 猜你喜欢
filecode变化情况
已经有15人回复
-
关于豆爷回答的JTJC与%2F数量
已经有6人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有31人回复
-
filecode
已经有19人回复
-
基金中了
已经有11人回复
-
听说今天filecode变了
已经有44人回复
-
固定端突然变了，今天
已经有8人回复
-
filecode
已经有11人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有6人回复
-
一个有机合成实验室都需要哪些设备？
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求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
50