| 查看: 1800 | 回复: 11 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
watertxf管理员
|
[求助]
求帮忙，关于含有bessel函数的方程的编程求解问题
|
|
|
我现在有个含有贝塞尔函数的方程，方程只有一个未知数k,我需要求方程等于0时k随l的变化。我编的程序如下，请各位高手帮忙看下哪里出了问题，非常感谢! ps：附件是要求的方程及所用到的关系式 a=0.01159; b=0.02317; e=10; l=7; syms k jleka=((pi/(2*e^0.5*k*a))^0.5)*besselj(l+0.5,e^0.5*k*a); jlka=((pi/(2*k*a))^0.5)*besselj(l+0.5,k*a); jlkb=((pi/(2*k*b))^0.5)*besselj(l+0.5,k*b); jldka=0.5*jlka+(k*a*pi/8)^0.5*(besselj(l-0.5,k*a)-besselj(l+1.5,k*a)); jldeka=0.5*jleka+(e^0.5*k*a*pi/8)^0.5*(besselj(l-0.5,e^0.5*k*a)-besselj(l+1.5,e^0.5*k*a)); ylka=((pi/(2*k*a))^0.5)*bessely(l+0.5,k*a); ylkb=((pi/(2*k*b))^0.5)*bessely(l+0.5,k*b); yldka=0.5*ylka+(k*a*pi/8)^0.5*(bessely(l-0.5,k*a)-bessely(l+1.5,k*a)); h11=jleka; h12=jlka/jlkb-ylka/ylkb; h21=jldeka; h22=jldka/jlkb-yldka/ylkb; y=h11*h22-h12*h21; k=solve(y)  |
回复此楼» 猜你喜欢
求调剂0703
已经有5人回复
-
一志愿天津大学339材料与化工求调剂
已经有4人回复
-
一志愿华南师范361分，化学求调剂
已经有3人回复
-
282求调剂
已经有6人回复
-
一志愿 南京航空航天大学 ，080500材料科学与工程学硕
已经有10人回复
-
273求调剂
已经有13人回复
-
一志愿085600中科院宁波所276分求调剂
已经有21人回复
-
求调剂
已经有6人回复
-
339求调剂
已经有5人回复
-
265求调剂
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求matlab求解一个方程组的问题
已经有6人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助完毕】用ode逆向求解微分方程问题
已经有4人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】求助PDEPE函数求解高阶偏微分方程 ：D
已经有5人回复
- 【求助】指数函数方程求解
已经有5人回复
watertxf
专家顾问
- 应助: 0 (幼儿园)
- 金币: 504.8
- 散金: 42
- 红花: 1
- 帖子: 46
- 在线: 16.5小时
- 虫号: 1289373
- 注册: 2011-05-08
- 性别: MM
- 专业: 同步辐射技术及其应用
5楼2011-07-14 14:38:35
fengsj9898
无虫
- 应助: 3 (幼儿园)
- 金币: 285.1
- 散金: 3600
- 红花: 5
- 帖子: 371
- 在线: 82.8小时
- 虫号: 425920
- 注册: 2007-07-28
- 专业: 机械结构强度学
【答案】应助回帖
★ ★
xzhdty(金币+2): 欢迎常来程序语言看看 2011-07-14 15:41:58
watertxf(金币+5): 您好！能不能麻烦您在跟我详细说一下？ 2011-07-15 14:07:49
xzhdty(金币+2): 欢迎常来程序语言看看 2011-07-14 15:41:58
watertxf(金币+5): 您好！能不能麻烦您在跟我详细说一下？ 2011-07-15 14:07:49
|
ki=0:0.1:1; nki=length(k); yi=zeros(nki,1); for i=1:nk k=ki(i); jleka=((pi/(2*e^0.5*k*a))^0.5)*besselj(l+0.5,e^0.5*k*a); jlka=((pi/(2*k*a))^0.5)*besselj(l+0.5,k*a); jlkb=((pi/(2*k*b))^0.5)*besselj(l+0.5,k*b); jldka=0.5*jlka+(k*a*pi/8)^0.5*(besselj(l-0.5,k*a)-besselj(l+1.5,k*a)); jldeka=0.5*jleka+(e^0.5*k*a*pi/8)^0.5*(besselj(l-0.5,e^0.5*k*a)-besselj(l+1.5,e^0.5*k*a)); ylka=((pi/(2*k*a))^0.5)*bessely(l+0.5,k*a); ylkb=((pi/(2*k*b))^0.5)*bessely(l+0.5,k*b); yldka=0.5*ylka+(k*a*pi/8)^0.5*(bessely(l-0.5,k*a)-bessely(l+1.5,k*a)); h11=jleka; h12=jlka/jlkb-ylka/ylkb; h21=jldeka; h22=jldka/jlkb-yldka/ylkb; y=h11*h22-h12*h21; yi(i)=y; end |
2楼2011-07-14 11:57:37
watertxf
管理员
- 应助: 0 (幼儿园)
- 金币: 504.8
- 散金: 42
- 红花: 1
- 帖子: 46
- 在线: 16.5小时
- 虫号: 1289373
- 注册: 2011-05-08
- 性别: MM
- 专业: 同步辐射技术及其应用
4楼2011-07-14 14:20:49
watertxf
管理员
- 应助: 0 (幼儿园)
- 金币: 504.8
- 散金: 42
- 红花: 1
- 帖子: 46
- 在线: 16.5小时
- 虫号: 1289373
- 注册: 2011-05-08
- 性别: MM
- 专业: 同步辐射技术及其应用
6楼2011-07-14 14:40:37
