| 查看: 1712 | 回复: 11 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
求帮忙，关于含有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)  |
回复此楼» 猜你喜欢
青椒八年已不青，大家都被折磨成啥样了？
已经有5人回复
-
青年基金C终止
已经有3人回复
-
26申博求博导推荐-遥感图像处理方向
已经有4人回复
-
限项规定
已经有7人回复
-
救命帖
已经有8人回复
-
招博士
已经有4人回复
-
西南交通大学国家级人才团队2026年博士研究生招生（考核制）—机械、材料、力学方向
已经有3人回复
-
英文综述是否需要润色及查重
已经有5人回复
-
为什么nbs上溴 没有产物点出现呢
已经有9人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求matlab求解一个方程组的问题
已经有6人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助完毕】用ode逆向求解微分方程问题
已经有4人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】求助PDEPE函数求解高阶偏微分方程 ：D
已经有5人回复
- 【求助】指数函数方程求解
已经有5人回复
6楼2011-07-14 14:40:37
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
4楼2011-07-14 14:20:49
5楼2011-07-14 14:38:35