24小时热门版块排行榜

返回列表

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

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

zhaoyanqiao

铁虫 (正式写手)

应助: 10 (幼儿园)
金币: 1578.5
散金: 187
红花: 5
帖子: 473
在线: 101.2小时
虫号: 1763251
注册: 2012-04-18
性别: GG
专业: 光学

[求助] 我想请问一个matlab编程的问题，sin（v/2）=3*v/14/pi

我想求解“sin（v/2）=3*v/14/pi；”，不知道matlab程序如何编写，恳请高人解答，并麻烦附上程序代码，谢谢了

回复此楼

» 猜你喜欢

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

向大家学习

1楼2012-07-06 10:32:38

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

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