24小时热门版块排行榜    

查看: 1468  |  回复: 12
本帖产生 1 个 程序强帖 ,点击这里进行查看
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

mafuwu

银虫 (小有名气)

[交流] 【求助】求助编程,急 已有7人参与

有一函数2/(2 + 1/4 (a + Sqrt[a^2 + 8])^2),画出该函数的图像,随后将函数里面的所有a替换为a (a + Sqrt[a^2 + 8])^2/16,以此进行,迭代5次。sqrt是对[a^2 + 8]开根号。用For循环,while循环,do循环都可。但要在mathematica的程序下。谢谢了

[ Last edited by senlia on 2010-4-13 at 13:44 ]
回复此楼

» 猜你喜欢

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

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

resonant

木虫 (正式写手)

独自为政

mafuwu(金币+8):谢谢,但是我运行不出来,得不到你说的这个结果。我用的是mathematica 7.0 2010-04-20 13:46
楼主这个东西的解答应该为
For[i=1;expr=2/(2+1/4*(a+Sqrt[a^2+8])^2),i<5,i++,expr=expr/.(a->a*(a+Sqrt[a^2+8])^2/16)]
这样得到的结果就正确了——

\!\(2/\((2 +
      1\/4\ \((\(\(1\/65536\)\((a\ \((a + \@\(8 + a\^2\))\)\^2\ \((1\/16\ a\ \
\((a + \@\(8 + a\^2\))\)\^2 + \@\(8 + 1\/256\ a\^2\ \((a + \@\(8 + \
a\^2\))\)\^4\))\)\^2\ \((1\/256\ a\ \((a + \@\(8 + a\^2\))\)\^2\ \((1\/16\ a\ \
\((a + \@\(8 + a\^2\))\)\^2 + \@\(8 + 1\/256\ a\^2\ \((a + \@\(8 + \
a\^2\))\)\^4\))\)\^2 + \@\(8 + \(a\^2\ \((a + \@\(8 + a\^2\))\)\^4\ \((1\/16\ \
a\ \((a + \@\(8 + a\^2\))\)\^2 + \@\(8 + 1\/256\ a\^2\ \((a + \@\(8 + a\^2\))\
\)\^4\))\)\^4\)\/65536\))\)\^2\ \((\(\(1\/4096\)\((a\ \((a + \@\(8 + \
a\^2\))\)\^2\ \((1\/16\ a\ \((a + \@\(8 + a\^2\))\)\^2 + \@\(8 + 1\/256\ a\^2\
\ \((a + \@\(8 + a\^2\))\)\^4\))\)\^2\ \((1\/256\ a\ \((a + \@\(8 + a\^2\))\)\
\^2\ \((1\/16\ a\ \((a + \@\(8 + a\^2\))\)\^2 + \@\(8 + 1\/256\ a\^2\ \((a + \
\@\(8 + a\^2\))\)\^4\))\)\^2 + √\((8 + \(\(1\/65536\)\((a\^2\ \((a + \@\(8 + \
a\^2\))\)\^4\ \((1\/16\ a\ \((a + \@\(8 + a\^2\))\)\^2 + \@\(8 + 1\/256\ a\^2\
\ \((a + \@\(8 + a\^2\))\)\^4\))\)\^4)\)\))\))\)\^2)\)\) + √\((8 + \
\(\(1\/16777216\)\((a\^2\ \((a + \@\(8 + a\^2\))\)\^4\ \((1\/16\ a\ \((a + \@\
\(8 + a\^2\))\)\^2 + \@\(8 + 1\/256\ a\^2\ \((a + \@\(8 + \
a\^2\))\)\^4\))\)\^4\ \((1\/256\ a\ \((a + \@\(8 + a\^2\))\)\^2\ \((1\/16\ a\ \
\((a + \@\(8 + a\^2\))\)\^2 + \@\(8 + 1\/256\ a\^2\ \((a + \@\(8 + \
a\^2\))\)\^4\))\)\^2 + √\((8 + \(\(1\/65536\)\((a\^2\ \((a + \@\(8 + \
a\^2\))\)\^4\ \((1\/16\ a\ \((a + \@\(8 + a\^2\))\)\^2 + \@\(8 + 1\/256\ a\^2\
\ \((a + \@\(8 + a\^2\))\)\^4\))\)\^4)\)\))\))\)\^4)\)\))\))\)\^2)\)\) + \
√\((8 + \(\(1\/4294967296\)\((a\^2\ \((a + \@\(8 + a\^2\))\)\^4\ \((1\/16\ a\ \
\((a + \@\(8 + a\^2\))\)\^2 + \@\(8 + 1\/256\ a\^2\ \((a + \@\(8 + \
a\^2\))\)\^4\))\)\^4\ \((1\/256\ a\ \((a + \@\(8 + a\^2\))\)\^2\ \((1\/16\ a\ \
\((a + \@\(8 + a\^2\))\)\^2 + \@\(8 + 1\/256\ a\^2\ \((a + \@\(8 + \
a\^2\))\)\^4\))\)\^2 + √\((8 + \(a\^2\ \((a + \@\(8 + a\^2\))\)\^4\ \((1\/16\ \
a\ \((a + \@\(8 + a\^2\))\)\^2 + \@\(8 + 1\/256\ a\^2\ \((a + \@\(8 + a\^2\))\
\)\^4\))\)\^4\)\/65536)\))\)\^4\ \((\(\(1\/4096\)\((a\ \((a + \@\(8 + a\^2\))\
\)\^2\ \((1\/16\ a\ \((a + \@\(8 + a\^2\))\)\^2 + \@\(8 + 1\/256\ a\^2\ \((a \
+ \@\(8 + a\^2\))\)\^4\))\)\^2\ \((1\/256\ a\ \((a + \@\(8 + a\^2\))\)\^2\ \
\((1\/16\ a\ \((a + \@\(8 + a\^2\))\)\^2 + \@\(8 + 1\/256\ a\^2\ \((a + \@\(8 \
+ a\^2\))\)\^4\))\)\^2 + √\((8 + \(\(1\/65536\)\((a\^2\ \((a + \@\(8 + \
a\^2\))\)\^4\ \((1\/16\ a\ \((a + \@\(8 + a\^2\))\)\^2 + \@\(8 + 1\/256\ a\^2\
\ \((a + \@\(8 + a\^2\))\)\^4\))\)\^4)\)\))\))\)\^2)\)\) + √\((8 + \
\(\(1\/16777216\)\((a\^2\ \((a + \@\(8 + a\^2\))\)\^4\ \((1\/16\ a\ \((a + \@\
\(8 + a\^2\))\)\^2 + \@\(8 + 1\/256\ a\^2\ \((a + \@\(8 + \
a\^2\))\)\^4\))\)\^4\ \((1\/256\ a\ \((a + \@\(8 + a\^2\))\)\^2\ \((1\/16\ a\ \
\((a + \@\(8 + a\^2\))\)\^2 + \@\(8 + 1\/256\ a\^2\ \((a + \@\(8 + \
a\^2\))\)\^4\))\)\^2 + √\((8 + \(\(1\/65536\)\((a\^2\ \((a + \@\(8 + \
a\^2\))\)\^4\ \((1\/16\ a\ \((a + \@\(8 + a\^2\))\)\^2 + \@\(8 + 1\/256\ a\^2\
\ \((a + \@\(8 + \
a\^2\))\)\^4\))\)\^4)\)\))\))\)\^4)\)\))\))\)\^4)\)\))\))\)\^2)\)\)
交朋识友,猛搞科研。两手抓,两手都要硬。
5楼2010-04-20 13:10:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答

resonant

木虫 (正式写手)

独自为政


mafuwu(金币+2):谢谢参与
画个函数还是很简单的吧?
fa[a_]=2/(2 + 1/4 (a + Sqrt[a^2 + 8])^2);
Plot[fa[a], {a, -10, 10}, DisplayFunction -> $DisplayFunction];

[ Last edited by resonant on 2010-4-14 at 17:22 ]
交朋识友,猛搞科研。两手抓,两手都要硬。
2楼2010-04-13 21:14:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

resonant

木虫 (正式写手)

独自为政

★ ★ ★ ★ ★
wangen994(金币+2):感谢参与应助 2010-04-14 08:08
wangen994(金币+3):活动期间额外奖励 2010-04-14 08:08
mafuwu(金币+6):谢谢。但不对,我想让他们出现在一个程序里面,这是一个迭代程序 2010-04-14 15:16
后一半也不难

For[i = 0, i < 5, fa[a_] = 2/(2 + 1/
  4 (a + Sqrt[a^2 + 8])^2) /. a -> a (a + Sqrt[a^2 + 8])^2/16, i++]
?fa
楼主自己运行一下看看。我用的是5.2

my result:
\!\(\*
  InterpretationBox[GridBox[{
        {GridBox[{
              {\(fa[a_] = 2\/\(2 +
        1\/4\ \((1\/16\
            a\ \((a + \@\(8 +
              a\^2\))\)\^2 + \@\(8 +
                  1\/256\ a\^2\ \((a + \@\(8 + a\^2\))\)\^4\))\)\^2\)\)}
              },
            GridBaseline->{Baseline, {1, 1}},
            ColumnWidths->0.999,
            ColumnAlignments->{Left}]}
        },
      GridBaseline->{Baseline, {1, 1}},
      ColumnAlignments->{Left}],
    Definition[ "fa"],
    Editable->False]\)

经再次检验,如果单次手动计算的话,似乎结果要复杂的多。可能我错了。期待高手上场。

[ Last edited by resonant on 2010-4-13 at 22:03 ]
交朋识友,猛搞科研。两手抓,两手都要硬。
3楼2010-04-13 21:23:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1377429145


mafuwu(金币+2):谢谢参与
好像也不是很复杂啊
4楼2010-04-13 21:56:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:50 by 5BDX0d0WFp7t
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-15 05:34 by 5BDX0d0WFp7t
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 05:22 by 5BDX0d0WFp7t
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 6/300 2026-09-15 05:13 by 5BDX0d0WFp7t
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 QUjhNVAcOSff 2026-09-13 8/400 2026-09-15 05:01 by mibUvS8DDCwf
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 00:21 by mibUvS8DDCwf
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 00:12 by vZfe6xYu34yj
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-14 23:33 by vZfe6xYu34yj
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 s3fFTmArrBt6 2026-09-13 3/150 2026-09-14 23:09 by vZfe6xYu34yj
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 LwdutQ8HoqWP 2026-09-13 3/150 2026-09-14 22:24 by vZfe6xYu34yj
[找工作] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-14 22:13 by vZfe6xYu34yj
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 21:37 by vZfe6xYu34yj
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:36 by vZfe6xYu34yj
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +4 两三岁ss 2026-09-14 4/200 2026-09-14 20:29 by 班主任的同桌
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-14 18:01 by Aj1rhIDL5ixY
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 15:09 by Aj1rhIDL5ixY
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[找工作] 浙江师范大学是怎么坑我的 +6 recruit123 2026-09-09 6/300 2026-09-14 12:00 by 萨仁其其格
[考博] 西北工业大学材料学院电化学传感与催化课题组招收2027年推免研究生、博士生、博士后 +3 马方园 2026-09-10 4/200 2026-09-14 09:29 by Chem张zz
[教师之家] 浙江师范大学是怎么坑我的 +10 recruit123 2026-09-09 15/750 2026-09-13 10:45 by seaskyy
信息提示
请填处理意见