2实验二及实验三主要命令解析.docVIP

实验二 极限与连续 1．求和命令sum与求积命令prod sum(x) %如果x是向量，则返回x的元素之和，如果x是矩阵，则返回矩阵各列之和。 prod(x) %如果x是向量，则返回x的元素之积，如果x是矩阵，则返回矩阵各列之积。 练习： 求. sum((1:100)) 求的近似值. sum(1./(1:100)) 求的近似值. prod(1:100) 求的近似值. s=0; for n=1:10 s=s+prod(1./(1:n)); end s 2．求极限命令 limit 命令limit用于计算数列或者函数的极限. 其基本形式是 limit(f(x),x,a) 其中 f(x) 是数列或者函数的表达式， a 是自变量x的变化趋势. 如果自变量x趋向于无穷，用 inf代替a. help syms SYMS Short-cut for constructing symbolic objects. SYMS arg1 arg2 ... is short-hand notation for arg1 = sym(arg1); arg2 = sym(arg2); ... Examples: syms x beta real is equivalent to: x = sym(x,real); beta = sym(beta,real); syms k positive is equivalent to: k = sym(k,positive); P14 例1 求极限 . syms n limit((2*n^3+1)/(5*n^3+1),n,inf) 例3 考虑第一个重要极限. ezplot(sin(x)/x,[-pi,pi]) syms x limit(sin(x)/x,x,0) 例4 考虑第二个重要极限. ezplot((1+1/x)^x,[1,100]) syms n limit((1+1/n)^n,n,inf) 试一试下面这句话 eval(1+exp(1)) 例2 求观察当时，函数的变化趋势. 输入 syms x limit(atan(x),x,+inf) limit(atan(x),x,-inf) ezplot(atan(x),[-50,50]) limit(atan(1/x),x,0,right) limit(atan(1/x),x,0,left) 例5 考虑无穷大. ezplot((1+2*x)/(1-x),[-3,4]) syms x limit((1+2*x)/(1-x),x,1,right) limit((1+2*x)/(1-x),x,1,left) 例6 考虑单侧无穷大. ezplot(exp(x^(-1)),[-2,2]) syms x limit(exp(1/x),x,0,left) limit(exp(1/x),x,0,right) 例7 输入 ezplot(x+4*sin(x),[0,20*pi]) 例8 输入 ezplot(x*sin(x),[0,20*pi]) 这个函数无界. 但当时，这个函数不是无穷大. 趋向于无穷大的函数当然无界，而无界函数并不一定是无穷大. 例9 观察可去间断. ezplot(tan(x)/x,[-1,1]) ezplot((sin(x)-x)/x^2,[-pi,pi]) 例10 观察跳跃间断. ezplot(sign(x),[-2,2]) ezplot((exp(1/x)-1)/(exp(1/x)+1),[-2,2]) 例11 观察无穷间断. 分别输入 ezplot((1-x^2)^(-1),[-3,3]) 例12 观察振荡间断. 分别输入 ezplot(sin(1/x),[-pi,pi]) syms x limit(sin(1/x),x,0) MATLAB的输出为ans = -1 .. 1. 这是什么意思？ 例13 有界量乘以无