求助：如何用matlab画含bessel函数的频散曲线？

yuandianqi

已知disperse是含有w和cp的表达式，若disperse =0，则应可以画出w关于cp的关系曲线，只是不知道这么复杂的式子如何来做？
disperse=-w^6*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2)*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2)*(.669529e-9*besselk(1.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*besseli(0.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*cp^4*besseli(1.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))*besselk(0.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))-.1798160e-5*besselk(1.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*besseli(0.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*cp^2*besseli(1.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))*besselk(0.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))+.669529e-9*besselk(1.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*besseli(0.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*cp^4*besselk(1.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))*besseli(0.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))-.1798160e-5*besselk(1.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*besseli(0.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*cp^2*besselk(1.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))*besseli(0.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))+.975975e-3*besselk(1.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*besseli(0.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*besseli(1.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))*besselk(0.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))+.975975e-3*besselk(1.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*besseli(0.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*besselk(1.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))*besseli(0.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))+.669529e-9*besseli(1.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*besseli(1.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))*cp^4*besselk(0.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*besselk(0.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))+.669529e-9*besseli(1.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*besseli(0.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))*besselk(1.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))*besselk(0.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*cp^4-.1798160e-5*besseli(1.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*besseli(0.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))*besselk(1.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))*besselk(0.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*cp^2+.975975e-3*besseli(1.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*besseli(0.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))*besselk(1.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))*besselk(0.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))-.1798160e-5*besseli(1.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*besseli(1.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))*besselk(0.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*cp^2*besselk(0.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))+.975975e-3*besseli(1.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*besseli(1.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2))*besselk(0.,.1050e-1*(-w^2*(1.*cp^2-.3546e8)/cp^2)^(1/2))*besselk(0.,.1964e-1*(-w^2*(1.*cp^2-.1013e8)/cp^2)^(1/2)))/cp^4

ChaChing

水平/时间有限但好奇! 这麼复杂, 不能先设法简化下!?

yuandianqi

回复 2 # ChaChing 的帖子

比如呢？怎么化简呢？一直发愁哦

ChaChing

個人水平/时间有限, 看到这麼复杂的東西, 真的懶得...
怎简化, LZ應該比別人容易, 畢竟是自己東西, 所以...
或許試試simple...

yuandianqi

谢谢您的解答。请问：在做自由板的频散曲线的时候，是一个模态一个模态的分别求取吗？怎么确定是哪个模态啊？如果是通过方程直接得到的，那最后你又怎么能判断众多的点中哪一个是哪个模态呢？您知道吗

ChaChing

yuandianqi 发表于 2012-5-29 23:49

                               
登录/注册后可看大图

谢谢您的解答。请问：在做自由板的频散曲线的时候，是一个模态一个模态的分别求取吗？怎么确定是哪个模态啊 ...

LZ这问题都快一年, 还没能解决吗? 水平有限专业不同, 建议LZ先这帖!
建议提问的网友分清 编程问题 和 专业问题 http://www.chinavib.com/thread-36746-1-1.html

还有这些是否应该在声学版块?