24小时热门版块排行榜

返回列表

本帖产生 1 个博学EPI ，点击这里进行查看

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

luomia

木虫 (著名写手)

博学EPI: 1
应助: 20 (小学生)
金币: 1779
散金: 4298
红花: 14
沙发: 1
帖子: 1417
在线: 371.2小时
虫号: 1864668
注册: 2012-06-19
性别: GG
专业: 机械摩擦学与表面技术

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
zhaoyanqiao: 金币+10, 博学EPI+1, ★★★★★最佳答案, 非常感谢，我怎么就没想到迭代呢 2012-07-07 13:43:11

引用回帖:

7楼: Originally posted by zhaoyanqiao at 2012-07-07 09:40:24
v=0:0.1:10;
y1=sin(v/2);
y2=3*v/14/pi;
plot(v,y1)
hold on
plot(v,y2)
grid on

可以看到，在v5.5左右时有一个交点，也就是我所需要的值...

迭代法程序如下：
function v=diedai(tol,a,n)
v=a;
f=0;
for i=1:n
f=14*pi*sin(v/2)/3;
if abs(f-v)<=tol
      break;
else
      v=f;
end
end
end
不过迭代次数很惊人
diedai(0.0001,5,20000)

ans =

-0.3535

>> diedai(0.0001,5,200000)

ans =

  -13.6383

>> diedai(0.0001,5,2000000)

ans =

-5.5123

可以看到一共应有3个解，0，+/-5.5123