|
回复 8 # sunhu2003 的帖子
谢谢,我现在直接应用的资料上的结果,即把他的行列式直接写出,然后令他为0.得到的应该是一个关于w和cp相速度的高阶方程。请问下面应该这么画出w-cp的曲线呢?我得到的行列式为disperse。他的表达式为下面的形式:比较复杂
isperse=-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 |
|