应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 【求助】求助编程，急
131/2返回列表12下一页
查看: 1194  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 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 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
奋斗！！！！！！
1楼 2010-04-13 09:25:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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 ]
赞一下(1人) 回复此楼
交朋识友，猛搞科研。两手抓，两手都要硬。
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 ]
赞一下(2人) 回复此楼
交朋识友，猛搞科研。两手抓，两手都要硬。
3楼2010-04-13 21:23:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1377429145

★
mafuwu(金币+2):谢谢参与
好像也不是很复杂啊
赞一下(1人) 回复此楼
4楼2010-04-13 21:56:41
已阅   回复此楼   关注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的回帖

resonant

木虫 (正式写手)

独自为政
mafuwu(金币+1):解这个迭代程序，不管用什么方法都可 2010-04-20 13:44
楼主的要求很奇怪，又要求使用For等循环，又不能是个迭代程序...
赞一下 回复此楼
交朋识友，猛搞科研。两手抓，两手都要硬。
6楼2010-04-20 13:26:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

resonant

木虫 (正式写手)

独自为政
mafuwu(金币+4):恩，谢谢 2010-04-20 21:44
hehe,楼主真是抱歉，我用5.2
原来装过7，结果由于乱改Preference，导致出现了意想不到的错误，卸掉以后只能等下次系统重装才能再安装这个了。

其实你要不用多核心的cpu，一般功能的话估计5.2和6.0也差不多。哈哈，只要能实现功能不妨用个低版本，占系统资源还少点。
赞一下 回复此楼
交朋识友，猛搞科研。两手抓，两手都要硬。
7楼2010-04-20 13:50:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

resonant

木虫 (正式写手)

独自为政
mafuwu(金币+5):我用5.2再试试，谢谢了，优秀 2010-04-20 21:47
wangen994(程序强帖+1):辛苦了 2010-05-09 21:02:30
其实一般还是要有点版本兼容的。
要不楼主使用这个：

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)]
采用延迟赋值方式呢？
在5.2种这两个方式得到的结果是一样的，但是不知道7里面是不是也一样。按理是应该用延迟的这种的。

[ Last edited by resonant on 2010-4-20 at 14:04 ]
赞一下(1人) 回复此楼
交朋识友，猛搞科研。两手抓，两手都要硬。
8楼2010-04-20 14:00:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mshwangg

至尊木虫 (正式写手)
★
mafuwu(金币+2):谢谢参与
引用回帖:
Originally posted by resonant at 2010-04-20 13:26:35:
楼主的要求很奇怪，又要求使用For等循环，又不能是个迭代程序...

版主很强大！稀饭
回复此楼
9楼2010-04-30 15:07:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

旦复旦

金虫 (正式写手)
★ ★
mafuwu(金币+2):谢谢参与
resonant(金币+1):这里的Mathematica程序都比较简单，使用的朋友还比较少，我也是急用现学。欢迎讨论交流：-）提出您的问题，共同解决，共同提高。 2010-04-30 20:04
这个我看了一下，好像也不是很复杂啊。
赞一下(2人) 回复此楼
10楼2010-04-30 15:10:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 mafuwu 的主题更新
131/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭