| 查看: 1940 | 回复: 3 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
xiaojixiong金虫 (小有名气)
|
[求助]
求助：如何求含贝塞尔函数的超越方程的复数解
|
|
|
如题，小弟最近碰到一含贝塞尔函数的超越方程。程序如下： i=1; for x=0:0.001:5 % x=kc/Wp for y=0:0.001:5 % y=w/Wp if (y-x)<=4.2&&0<(y-x) t=0.6*1000*sqrt(y.*y-x.*x); ta=0.6*sqrt(y.*y-x.*x); tc=0.6*2* sqrt(y.*y-x.*x); t1=(y.*y-1); t2=x.*x*(3.4596-y.*y); t3=y.*y*(y*y-4.4596); T=0.6*1000*sqrt((t1.*(t2+t3))/t3); Ta=0.6*sqrt((t1.*(t2+t3))/t3); pTa=(Ta/4).*(Ta/4); pta=(ta/4).*(ta/4); ptc=(tc/4).*(tc/4); % ta有关的贝塞尔函数（4个）； cta0=besselj(0,ta(：)); cta1=besselj(1,ta(：)); dta0=bessely(0,ta(：)); dta1=bessely(1,ta(：)); % tc 有关的贝塞尔函数(2个) ctc0=besselj(0,tc(：)); dtc0=bessely(0,tc(：)); A=ctc0.*dta0-cta0.*dtc0; B=-ctc0.*dta1+cta1.*dtc0; cTa0=besselj(0,Ta(：)); cTa1=besselj(1,Ta(：)); f=((y.*y-1)/T).*cTa1.*A+(y.*y/t)*cTa0 *B; if abs(f)<=0.00000001 a0(i)=x;b0(i)=y; i=i+1; disp([num2str(x) ' ' num2str(y)]); end end end end 我通过限制f绝对值的精度，去求出了满足f=0的实数解，但是我若想求出满足f=0，x和对应的y的复数解。这个程序应该怎么改啊，谢谢大家了。 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意！
已经有13人回复
-
同一个分子，一条来自文献，一条来自AI——不告诉你答案，你会选哪条？
已经有7人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 高数关于不显含x和y的微分方程求解？
已经有8人回复
- 如何求解含复数系数的贝塞尔微分方程
已经有8人回复
xiaojixiong
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 3482.5
- 散金: 3
- 帖子: 239
- 在线: 67.6小时
- 虫号: 425899
- 注册: 2007-07-28
- 性别: GG
- 专业: 等离子体物理
3楼2014-09-25 08:35:26
0605102034
铜虫 (著名写手)
- 应助: 2 (幼儿园)
- 金币: 4272.9
- 散金: 89
- 红花: 4
- 帖子: 2445
- 在线: 670小时
- 虫号: 1514410
- 注册: 2011-11-28
- 性别: GG
- 专业: 碳素材料与超硬材料
2楼2014-09-24 22:44:15
xiaojixiong
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 3482.5
- 散金: 3
- 帖子: 239
- 在线: 67.6小时
- 虫号: 425899
- 注册: 2007-07-28
- 性别: GG
- 专业: 等离子体物理
4楼2014-09-25 08:43:09
