[转帖]mathematica若干使用技巧

suffer · 发表于 2006-2-23 15:57

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

mathematica的取实部
 问：取一个复数函数的实部时，怎么定义函数其中的某些量是实数呢？ 如f=exp(i*a),a是实数，然后求re(f),但是结果总是把a当复数处理的，怎么办呢？ 
通过一个例子说明问题： AllTwo = {Exp[I a], Conjugate[Exp[I a]]}; RealPart = (ComplexExpand[AllTwo[[1]]] + ComplexExpand[AllTwo[[2]]])/2; ImaginPart = (ComplexExpand[AllTwo[[1]]] - ComplexExpand[AllTwo[[2]]])/2/I; Print["RealPart=", RealPart, "; ImaginPart=", ImaginPart]
Out[33]= RealPart="Cos[a]; ImaginPart=Sin[a]
[问题2]变量替换
如果X是x的函数，且存在如下关系：x=a+b ks，那么dX/dx如何表示成为dX/dks? In[2]:= x = a + b ks; X = fg[x]; D[X, ks]
Out[2]=b fg'[a+b ks] 
同理，可把dX/dks表示成为dX/dx的命令行如下。 In[3]:= Clear[fg, x, ks] ks = (x - a)/b; X = fg[ks]; D[X, x]
Out[3]= fg'[(-a+x)/b]/b
例1：X是ks的函数，且X=ks Z, ks=a+b x;则d2X/dx2+X表示成关于d2Z/dks2, dZ/dks，Z（Z对ks的导数等）的式子。 In[1]:= ks = a + b x; Z = fg[ks]; X = ks Z; sh = D[X, {x, 2}] + c X Out[1]=c(a+b x)fg[a+b x]+2b^2 fg'[a+b x]+b^2 (a+b x)fg''[a+b x] 其中求导为对ks的导数。这可通过如下命令验证。 D[Z, x] b fg’[a+bx] 数字填图问题
(万福永《现代数学教育》科学出版社)
数字填图问题是数学问题的一种趣味形式。早在19世纪后半期，一些数学家就在报刊中大量使用数字填图游戏和字谜游戏等，目的是使业余爱好者也能通过简单的形式去认识、理解和琢磨深奥的数学问题，这些问题中甚至包括困惑了世间智者350多年、于1994 年才刚刚被证明了的“费马大定理”。100多年来，数字填图问题对数学界所起的作用是不言而喻的。 大家都知道，数学问题一般都经过严格的逻辑证明才得以解决。而逻辑证明是指从一些公理出发，经过逻辑推理来证明问题。但随着20世纪40年代以来计算机的诞生和发展，计算机改变了整个世界，计算机已在各个领域发挥作用，并取得了许多重大进展。于是, 能否用计算机来证明数学问题便成了大家关心的话题。所谓计算机证明是指充分发挥计算机计算速度快和会“推理”的特点，用计算机程序模拟解题或进行穷举检验，最后得到问题的解。 我们将通过生活中几个常见的数字填图问题的探究，谈谈这类问题的逻辑推理解法和计算机解法。 例1：在图1的几个加法等式中，每个□表示一个非零数字，任意两个数字都不相同，问有多少个解？ □＋□＝□，□＋□＝□，□＋□＝□0 图1
计算机解法——为验证此结果，用 Mathematica编程进行模拟解题，充分利用计算机运算速度快的特点进行穷举法检验。实践表明本问题解的情况恰如上所述。我们用Timing[ ]函数对两个程序进行计时，可以看出程序2是程序1的提速版。 程序1： Timing[ (*a+b=c,d+e=f,f+h=i0*) Clear[n,a,b,c,d,e,f,g,h,i];n=0; For[a=1,a<=9,a++, For[b=1,b<=9,b++,If[b!=a, For[c=1,c<=9,c++,If[c!=a&&c!=b, For[d=1,d<=9,d++,If[d!=a&&d!=b&&d!=c, For[e=1,e<=9,e++,If[e!=a&&e!=b&&e!=c&&e!=d, For[f=1,f<=9,f++,If[f!=a&&f!=b&&f!=c&&f!=d&&f!=e, For[g=1,g<=9,g++,If[g!=a&&g!=b&&g!=c&&g!=d&&g!=e&&g!=f, For[h=1,h<=9,h++,If[h!=a&&h!=b&&h!=c&&h!=d&&h!=e&&h!=f&&h!=g, For[i=1,i<=9,i++,If[i!=a&&i!=b&&i!=c&&i!=d&&i!=e&&i!=f&&i!=g&&i!=h &&a+b==c&&d+e==f&&g+h==10*i&&a<b&&d<e&&a<d&&g<h, Print[ ++n,": ",a,"+",b,"=",c,",",d,"+",e,"=",f,",",g,"+",h,"=",i,"0"] ]]]]]]]]]]]]]]]]] (*total have 17 right ")"s *) ] 程序2： Timing[ Clear[n,a,b,c,d,e,f,g,h,i];n=0; For[a=1,a<=9,a++, For[b=a+1,b<=9,b++, For[c=b+1,c<=9,c++, For[d=a+1,d<=9,d++,If[d!=b&&d!=c, For[e=d+1,e<=9,e++,If[e!=b&&e!=c, For[f=e+1,f<=9,f++,If[f!=c, For[g=1,g<=9,g++,If[g!=a&&g!=b&&g!=c&&g!=d&&g!=e&&g!=f, For[h=g+1,h<=9,h++,If[h!=a&&h!=b&&h!=c&&h!=d&&h!=e&&h!=f, For[i=1,i<=9,i++,If[i!=a&&i!=b&&i!=c&&i!=d&&i!=e&&i!=f&&i!=g&&i!=h &&a+b==c&&d+e==f&&g+h==10*i,
 Print[ ++n,": ",a,"+",b,"=",c,",",d,"+",e,"=",f,",",g,"+",h,"=",i,"0"] ]]]]]]]]]]]]]]] (*total have 15 right ")"s *) ]
 函数说明： 1. For[...]：循环函数； 2. While[...]：循环函数； For循环函数和While循环函数适用于生成非结构循环，While循环函数只有当条件满足时才执行循环体，而For循环函数则按一定步长递增执行循环体。循环功能基本上和C语言的一样。格式：For[start, test, incr, body] 和 While[test, body]。 例如：Clear[i,t,x,n] For[i=1; t=x, i^2<10, i++, t=t^2+i; Print[t]]n=20; While[(n=Floor[n/2])!=0, Print[n]] 3. If[...]：Mathematica提供的设置单个条件的函数，一般用If函数控制分支。 格式: If[condition, t, f, u]，条件满足时取t，条件不满足时取f，否则取u。 例如: If[x>y, a, b, c] 下面，我们来动手练习两题： 练习1（一道竞赛题）： 1998年4月香港数理教育学会主办的初中数学竞赛有这样一道试题：在下面的加法算式中，每个□表示一个数字，任意两个数字都不相同，那么A与B的乘积的最大值是多少？ □＋□□＋□□□＝□□AB
解答：最大值是15。你能给出逻辑推理解法并用计算机加以验证吗？ 练习2（真分数问题）：
□/□□+□/□□＝□/□□ 这是由一位日本数学家在一次数学家大会上提出的一个问题：图3中的加法等式是：两个真分数之和等于第三个真分数（不要求为既约分数），每个□表示一个非零的数字，任意两个数字都不相同。比如：1/32+5/96＝7/84 ，试找出所有可能的解。 本问题利用计算机程序已找到解答，共有 10 个解，但想给出相应的逻辑解法相当困难，可以说计算机模拟解题能力在某些方面确已达到并超过了逻辑推理解题的能力。 转自科研中国

aiweng · 发表于 2012-3-6 10:28

这个问题有点复杂，一起关注，期待高手！

账号		自动登录	找回密码
密码			我要加入

[mathematica] [转帖]mathematica若干使用技巧

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。