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

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

已知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