文件名:f1.m
内容：
function F=f1(k,w)
F=tanh(pi/0.65*w*sqrt(k^2-1))-1/(17.36-0.715*i)*sqrt(k^2+17.36-0.715*i)/sqrt(k^2-1);
---------------------------------------------------------
文件名:*.m
内容：
clear all;
w=0:0.1:100;%这里设了步长0.1，可自行调节
for n=1:size(w,2)
k(n)=fsolve(@(k)f1(k,w(n)),6,optimset('Display','off','TolFun',1e-6));%用了fsolve，实际上是当作优化问题处理了
end
kr=real(k);
ki=imag(k);
subplot(1,2,1);%左边的画实部
plot(w,kr);
subplot(1,2,2);%右边的画虚部
plot(w,ki);
---------------------------------------------------------
ps：获得的图形很奇特，不连续，如果w代表频率，想不出k代表什么。
LZ可以用命令
tanh(pi/0.65*w.*sqrt(k.^2-1))-1./(17.36-0.715*i).*sqrt(k.^2+17.36-0.715*i)./sqrt(k.^2-1);
验证计算结果
|