24小时热门版块排行榜    

查看: 829  |  回复: 8
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

yongbinz

金虫 (著名写手)

[交流] 【讨论】一个小程序,过来讨论讨论!

我需要一个小程序,求解一个方程!
X/L=Sin(θ/2)÷( θ/2) 公式中的θ为弧度,已知X与L,输入X与L后,要能求出θ的大小。
这个程序怎么写呢,大家能否给点意见?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

duzhuqi

★ ★
小木虫(金币+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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

gjliu

银虫 (小有名气)


小木虫(金币+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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yongbinz

金虫 (著名写手)

谢谢楼上,不过,本人就是不会C语言呀!谢谢支持!
3楼2009-06-03 11:40:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhuzhs

金虫 (正式写手)

maple,mathmatics
当所有一切都已看平淡 是否有一种精致还留在心田
5楼2009-06-06 00:14:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-31 01:38 by hZiFeudZyoGR
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 gy1nBQXYQJqL 2026-08-29 8/400 2026-08-31 01:26 by hZiFeudZyoGR
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-31 01:07 by hZiFeudZyoGR
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +7 ASdOkHsho7FD 2026-08-28 10/500 2026-08-31 00:42 by hZiFeudZyoGR
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-30 22:24 by xiangy672
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 4FFAWE8HcgUD 2026-08-29 4/200 2026-08-30 20:40 by hZiFeudZyoGR
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 jCd0dEvKHShX 2026-08-29 6/300 2026-08-30 19:13 by hZiFeudZyoGR
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 G6APbkg8SA6w 2026-08-29 6/300 2026-08-30 19:02 by hZiFeudZyoGR
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 G6APbkg8SA6w 2026-08-29 5/250 2026-08-30 18:40 by hZiFeudZyoGR
[博后之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 18:28 by hZiFeudZyoGR
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 gy1nBQXYQJqL 2026-08-29 9/450 2026-08-30 17:56 by hZiFeudZyoGR
[论文投稿] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 17:03 by hZiFeudZyoGR
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
请填处理意见