求助：如何用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
还是没有回复哦
页:
[1]