| 查看: 828 | 回复: 8 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
yongbinz金虫 (著名写手)
|
[交流]
【讨论】一个小程序，过来讨论讨论！
|
||
|
我需要一个小程序，求解一个方程！ X/L=Sin(θ/2)÷( θ/2) 公式中的θ为弧度，已知X与L，输入X与L后，要能求出θ的大小。 这个程序怎么写呢，大家能否给点意见？ |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有8人回复
-
面上意见出来了
已经有12人回复
-
为什么到现在没收到通知？
已经有4人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有10人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有6人回复
-
有没有仍没收到信息的
已经有5人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有9人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有10人回复
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
suton(金币+1,VIP+0):谢谢！ 6-7 09:10
小木虫(金币+0.5):给个红包，谢谢回帖交流
suton(金币+1,VIP+0):谢谢！ 6-7 09:10
|
MATLAB程序如下： function y=f6(x) y=Sin(x)/( x)-X/L; end ---------------------------------------------------- a=-1;b=1.5; y1=f6(a);y2=f6(b); while abs(b-a)>0.001 x0=0.5*(a+b);y0=f6(x0); if abs(y0)< 0.0001 x0 elseif y0*y1<0 b=x0;y2=y0; else a=x0;y1=y0; end end x0 |
6楼2009-06-06 22:28:45
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
幻影无痕(金币+0,VIP+0):呵呵，亲自上阵！鼓励一下！ 6-3 11:04
小木虫(金币+0.5):给个红包，谢谢回帖交流
幻影无痕(金币+0,VIP+0):呵呵，亲自上阵！鼓励一下！ 6-3 11:04
|
这个是一个超越方程求解的问题，可以看一些《计算方法》或者《数值分析》的书，比如牛顿迭代什么的方法都可以找到数值解的。如果是Matlab或者Mathematic等软件直接有求解方程的函数，相对比较简单。 就是不知道你要用什么语言去写。给你一个简单的例子吧 用牛顿迭代法求方程X3+4X2-10=0 的根 #include main() {float x,x0,f,f1; x=1.5; do {x0=x; f=x0*x0*(x0+4)-10; f1=x0*(3*x0+8); /* f的导数，即它的切线 */ x=x0-f/f1; /* 切线与x轴的交点 */ } while(fabs(x-x0)>=1e-5); printf("The root is %8.7f\n",x); } [ Last edited by gjliu on 2009-6-3 at 11:21 ] |
2楼2009-06-03 11:00:49
yongbinz
金虫 (著名写手)
- 应助: 6 (幼儿园)
- 金币: 1180.2
- 散金: 1415
- 红花: 3
- 帖子: 1995
- 在线: 892.5小时
- 虫号: 428019
- 注册: 2007-08-01
- 专业: 文学理论
3楼2009-06-03 11:40:43
5楼2009-06-06 00:14:43
