24小时热门版块排行榜    

查看: 453  |  回复: 1

fenglisha

铜虫 (初入文坛)

[求助] 求:牛顿迭代法matlab程序及输出结果.求指点

x=xlsread('d:\x');
y=xlsread('d:\y');
r=xlsread('d:\curvature');
m(1)=0;n(1)=0;

u(1)=-826.9388,v(1)=7517.6;
for  i=2:150
[m(i) n(i)]=fsolve(‘(m(i)-u(i-1))^2+(n(i)-u(i-1))^2=r(i-1)^2','(m(i)-m(i-1))^2+(n(i)-n(i-1))^2=(x(i)-x(i-1))^2+(y(i)-y(i-1))^2’,’m(i)',’n(i)’)
k(i)=(m(i)-u(i-1),n(i)-v(i-1));
u(i)=m(i)-r(i)*abs(m(i)-u(i-1))/norm(k(i));
v(i)=m(i)-r(i)*abs(m(i)-v(i-1))/norm(k(i));

end
总体就是已知点(0,0)和(-826.9388,7517.6),循环求解后面的149个相应的点。能不能帮我看看程序要怎么写。
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : curvature.xls
  • 2013-03-28 11:38:52, 23.5 K
  • 附件 2 : x.xls
  • 2013-03-28 11:38:54, 23.5 K
  • 附件 3 : y.xls
  • 2013-03-28 11:38:56, 23.5 K

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fenglisha

铜虫 (初入文坛)

不要沉啊
2楼2013-03-28 17:17:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fenglisha 的主题更新
信息提示
请填处理意见