>> n=1.5;
>> [x fval]=fsolve(@(x) tan(x)-  n/(1+x^2*n/((n-1)^2)  ),0)

2楼: Originally posted by 月只蓝 at 2015-08-12 11:25:56
程序中以x表示k*L，由于你没有给出N的数值，我就以n=1.5为例，代码如下：
>> n=1.5;
>> =fsolve(@(x) tan(x)-  n/(1+x^2*n/((n-1)^2)  ),0)

计算结果：
解x =

    0.5198


方程残差fval = ...

3楼: Originally posted by Fanny675 at 2015-08-12 12:14:44
你把方程弄错了，丢掉了（k*L)^2...

2楼: Originally posted by 月只蓝 at 2015-08-12 11:25:56
程序中以x表示k*L，由于你没有给出N的数值，我就以n=1.5为例，代码如下：
>> n=1.5;
>> =fsolve(@(x) tan(x)-  n/(1+x^2*n/((n-1)^2)  ),0)

计算结果：
解x =

    0.5198


方程残差fval = ...

4楼: Originally posted by 月只蓝 at 2015-08-12 12:20:07
x表示 k*L，x^2即（k*L)^2。
你若不信，把 k*L=0.5198代回方程两边，看两边相等不。

此外，你给出的方程 “tan（k*L）=（k*L)/[1+N/(N-1)^2*(k*L)^2 ” 和你图片中给出的不一致，这个需要你核对。

我是按照 ...

5楼: Originally posted by Fanny675 at 2015-08-12 12:27:09
我有一个例子，其中N=2.89,求得k*L=3.455。我用楼主的方法试了下不对...

7楼: Originally posted by 月只蓝 at 2015-08-12 12:33:10
你要求解的方程，到底是  “tan（k*L）=（k*L)/[1+N/(N-1)^2*(k*L)^2 ” 还是你给出的图片中的 “tan（k*L）=N/[1+N/(N-1)^2*(k*L)^2 ”？
我给出的代码是针对图片中方程的。...

8楼: Originally posted by Fanny675 at 2015-08-12 12:36:23
对不住楼主了，我仔细看了下，我要求的应该是tan（k*L=（k*L)/[1+N/(N-1)^2*(k*L)^2，图片中的分子部分打错了...

n=2.89;
[x fval]=fsolve(@(x) tan(x)-  x/(1+x^2*n/((n-1)^2)  ),1)