Mathematica第二章基本概念下.ppt

科学计算软件 第二章 基本概念 （ 下） 习题解答 2.39 Expand[表达式]：对表达式进行代数展开 分别用=和:=将符号a与b定义成Expand[(x+1)^3]，令x=u+v，再计算a与b(体会两者之间的差别） a=Expand[(x+1)^3] b:=Expand[(x+1)^3] ?b x=u+v a b 习题解答 2.40 Together：把两个分数组合利用通分组成有公分母的一个分数 分别用=和：=把符号y与z定义成Together[a+b]，然后令a=1/x，b=1/(x+1)，再计算y与z。 Clear[a,b,x] y=Together[a+b] z:=Together[a+b] a=1/x; b=1/(x+1); y z 习题解答 2.41 Simplify[表达式]：尝试简化代数表达式 输入a=Simplify[poly]和b:=Simplify[poly]，然后令poly=x2+2x+1，计算a与b的值，并解释结果的差别 Clear[a,b,x] a=Simplify[poly] b:=Simplify[poly] poly=x^2+2x+1 a b 2.6 和与积 Mathematica中和与积的计算相当简单 任何符号都可以做为求和的指标（在下面的描述中使用的是i） 使用步长increment，指标可以负增长 Sum[a[i],{i,imax}]：下标i从1到imax，间隔为1 Sum[a[i],{i,imin,imax}]：下标i从imin到imax Sum[a[i],{i,imin,imax,increment}]：下标i从imin到imax，间隔为increment 2.6 和与积 例46 为了计算前20相邻整数的平方和，输入 Sum[i^2,{i,1,20}] 先生成列表，再对列表中的数据求和 st=Table[i^2,{i,1,20}] Total[st] 2.6 和与积 例47 计算和1/15+1/17+1/19+...+1/51 Sum[1/i,{i,15,51,2}] 注意参与求和的都是有理数形式，结果是精确的！ 例48 给出和1/15+1/17+1/19+......+1/51的近似值 NSum[1/i,{i,15,51,2}] Sum[1./i,{i,15,51,2}] 2.6 和与积 求和限可以是无穷，Mathematica采用复杂的方法计算无穷和。 例49 计算1/1+1/4+1/9+1/16+...... Sum[1/i^2,{i,1,Infinity}] Mathematica的强项在于符号计算！ 2.6 和与积 二重求和可以利用下述语法得到，也可以更方便地两次点击BasicInput模板上的∑符号实现，同样可以扩展到三重和、四重和甚至更多重和。 Sum[a[i,j],{i,imax},{j,jmax}] Sum[a[i,j],{i,imin,imax},{j,jmin,jmax}] Sum[a[i,j],{i,imin,imax,i_increment},{j,jmin,jmax, j_increment}] 2.6 和与积 例50 计算（1/1+1/2+1/3+1/4）+（2/1+2/2+2/3+2/4）+（3/1+3/2+3/3+3/4） Sum[i/j,{i,1,3},{j,1,4}] 在Matlab中，用如下代码实现相同功能 sm=0 for i=1:3 for j=1:4 sm=sm+i/j end end 2.6 和与积 与Sum的语法类似，Mathematica中用Product计算乘积 Product[a[i],{i,imax}] Product[a[i],{i,imin,imax,increment}] 同样可以推广到多重乘积的情况 2.6 和与积 例51 计算从4到9的相邻整数的乘积 Product[i,{i,4,9}] 习题解答 2.45 计算前25个素数的和 Sum[Prime[k],{k,1,25}] 通项表达式为Prime[k] 2.49 给出从1到n相邻整数平方和的一般公式 Sum[k^2,{k,1,n}] 注意n是符号（没有赋予数值） 习题解答 2.52 计算和1+(1+1/2)+(1+1/2+1/3)+...+(1+1/2+1/3+... 1/20) Sum[1/j,{i,1,20},{j,1,i}] 2.53计算(1+1/2)(1+1/2+1/3)…(1+1/2+1/3+…+1/10)的值 NProduct[Sum[1/j,{j,1,i}],{i,2,10}] 习题解答 (1+1/2)(1+1/2+1/3)…(1