matlab求解1隐函数
小辈求解一隐函数, 其中x 和y为自变量和因变量。别的参数值已经给出,具体程序如下:syms x y;k0=4.0516e+06;em=-104+3.68i;ed=2.152089;f=em*k0*sqrt(y^2-ed)*tanh(sqrt(y^2-ed)*k0*x/2/(10^9))+ed*k0*sqrt(y^2-em);ezplot(f,[0,500,1.5,4]);结果中,y为复数。但 根据以上程序,只得到x和y的实部 曲线图。欲求x 和y的虚部 曲线图,该怎么编程。望各位高手指教!!我来帮他解答满意回答2010-12-07 19:31%参数赋值k0=4.0516e+06;em=-127.8+3.21i;ed=2.25;?%先设定d的值,再对方程进行求解d=1:1000;z=d;??%把每一个d带入,对应求出一个xfor?i=1:500f?=?@(x)(em*k0*sqrt(x(1)^2-ed)*tanh(sqrt(x(1)^2-ed)*k0*d(i)/2/(10^9))+ed*k0*sqrt(x(1)^2-em));x=fsolve(f,1.5+0.0003i,optimset('Display','off'));z(i)=x(1);end???%至此所有的d对应的x全部求出,也就是z,下面开始画图subplot(2,2,[1?2]);plot(d,real(z),d,1.5133,'--');xlabel('d(nm)');ylabel('Real?part?of?N_{eff}');axis([0,500,1.4,4]);gtext('1.5133');legend('MIM','SP');subplot(2,2,[3?4]);plot(d,imag(z),d,0.0003,'--');xlabel('d(nm)');ylabel('Imaginary?part?of?N_{eff}');axis([0,500,-0.001,0.03]);gtext('0.0003');legend('MIM','SP');