| 查看: 708 | 回复: 8 | |||
| 当前主题已经存档。 | |||
yongbinz金虫 (著名写手)
|
[交流]
【讨论】一个小程序，过来讨论讨论！
|
||
|
我需要一个小程序，求解一个方程！ X/L=Sin(θ/2)÷( θ/2) 公式中的θ为弧度，已知X与L，输入X与L后，要能求出θ的大小。 这个程序怎么写呢，大家能否给点意见？ |
回复此楼» 猜你喜欢
271求调剂
已经有37人回复
-
一志愿中科大材料与化工，353分还有调剂学校吗
已经有8人回复
-
申博自荐
已经有3人回复
-
300求调剂
已经有10人回复
-
26年电池方向博士申请
已经有5人回复
-
急需调剂
已经有6人回复
-
RY：中国产出的科学垃圾论文，绝对数量和比例都世界第一
已经有18人回复
-
291求调剂
已经有8人回复
-
0854调剂
已经有9人回复
-
322求调剂
已经有15人回复
★
小木虫(金币+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
4楼2009-06-03 12:31:30
5楼2009-06-06 00:14:43
★ ★
小木虫(金币+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
yongbinz
金虫 (著名写手)
- 应助: 6 (幼儿园)
- 金币: 1180.2
- 散金: 1415
- 红花: 3
- 帖子: 1995
- 在线: 892.5小时
- 虫号: 428019
- 注册: 2007-08-01
- 专业: 文学理论
7楼2009-06-06 23:42:20
8楼2009-06-07 08:40:46
sitonyl
铁杆木虫 (文坛精英)
- 应助: 33 (小学生)
- 金币: 12935.9
- 红花: 3
- 帖子: 10448
- 在线: 338.8小时
- 虫号: 722261
- 注册: 2009-03-14
- 性别: GG
- 专业: 导航、制导与传感技术
9楼2009-06-09 09:12:21
5