求达人给个解我来我网·5come5 Forum

我来我网·5come5 Forum » 学业有成 » 电子设计·数学建模 » 求达人给个解

本页主题: 求达人给个解	显示签名 \| 打印 \| 加为IE收藏 \| 收藏主题 \| 上一主题 \| 下一主题

huyuge512

∷ 性别:

∷ 状态:

∷ 等级: 栋梁之材

∷ 发贴: 807

∷ 威望: 0

∷ 浮云: 1129

∷ 在线等级:

∷ 注册时间: 2005-06-11

∷ 最后登陆: 2007-06-25

【复制此帖地址只看此人回复】

5come5帮你背单词 [ miss /mis/ B n. 小姐 ]

求达人给个解

n1=1.55;n2=1.67;n3=1.50;
H=1.5;
q=0;
w0=1.55;
syms x;
k0=2*pi/w0;
r1=k0*sqrt(x^2-n1^2);
r2=k0*sqrt(n2^2-x^2);
r3=k0*sqrt(x^3-n3^2);
f=r2*H-atan(n2^2*r1/(n1^2*r2))-atan(n2^2*r3/(n3^2*r2))-q*pi;

我用牛顿迭代法得一个解是1.62529822143391.
而正确解是1.6312728221.
下面是我用的牛顿迭代法算法
function y=newton(x0)
x1=x0-fc(x0)/df(x0);%fc(x)是原函数,df(x)是导函数;
n=1;
while abs(x1-x0)>=1.0e-6
x0=x1;
x1=x0-fc(x0)/df(x0);n=n+1;
end
x1
n
不知道哪里出了问题,请达人赐教