24小时热门版块排行榜    

CyRhmU.jpeg
查看: 608  |  回复: 4

hedan1215

银虫 (小有名气)

[交流] 【求助】求解方程程序出错已有2人参与

方程为  8*J1(X)*J2(X)/(X^2)=0 求解x
其中J1(X)为一阶贝塞尔函数,J2(X)为二阶贝塞尔函数

程序如下:
syms x y1 y2 ;
y1=besselj(1,x);
y2=besselj(2,x);
x=solve('8*(y1)*(y2)/(x^2)=0')

运行结果:x =

[ empty sym ]             请帮忙改一下程序!多谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
robert2020(金币+2):多谢提供解决方案,欢迎常来仿真模拟版交流! 2010-06-04 22:27:35
>> syms x
x=solve('8*besselj(1,x)*besselj(2,x)/(x^2)=0')

x =

matrix([[-226.97841709642947178848018795482]])

>> eval(x)

ans =

   -2.269784170964295e+002


>> x=fsolve('8*besselj(1,x)*besselj(2,x)/(x^2)',3)


x =

   3.831697904464683

>> x=fsolve('8*besselj(1,x)*besselj(2,x)/(x^2)',-3)

x =

-3.831697904100898 + 0.000000000000689i

[ Last edited by change0618 on 2010-6-4 at 22:11 ]
2楼2010-06-04 22:06:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
robert2020(金币+1):多谢提供解决方案,欢迎常来仿真模拟版交流! 2010-06-04 22:28:04
先画一下图看看它是什么样的图形,然后再决定怎么解,能不能接,需要哪个范围的的解。
>> x=-300:300;
>> y=8*besselj(1,x).*besselj(2,x)./(x.^2);
>> plot(x,y)
3楼2010-06-04 22:13:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hedan1215

银虫 (小有名气)

引用回帖:
Originally posted by change0618 at 2010-06-04 22:06:48:
>> syms x
x=solve('8*besselj(1,x)*besselj(2,x)/(x^2)=0')

x =

matrix([[-226.97841709642947178848018795482]])

>> eval(x)

ans =

   -2.269784170964295e+002


>>  ...

请问楼主:
syms x
x=solve('8*besselj(1,x)*besselj(2,x)/(x^2)=0')
运行之后显示x =

0.
0.
另外,
“x =

matrix([[-226.97841709642947178848018795482]])”这句是什么意思呢?多谢!
4楼2010-06-05 14:35:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师


小木虫(金币+0.5):给个红包,谢谢回帖交流
我不懂,希望你搞懂
5楼2010-06-05 15:02:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hedan1215 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见