MATLAB实验课第1讲基本操作-1.1-1.8分解.ppt

* * 求导数 diff(S,v) 求表达式S对变量v的一阶导数。 diff(S,v,n) 求表达式S对变量v的n阶导数。 如求sin(x)+ex的三阶导数，键入命令 diff(sin(x)+x*exp(x),3) 得 ans = -cos(x)+3*exp(x)+x*exp(x) * * sym a t x; f=[a,t^3;t*cos(x);log(x)]; df=diff(f); dfdt2=diff(f,t,2); dfdxdt=diff(diff(f,x),t); * * 求和 symsum(S) 对通项S求和，其中k为变量。且从0变到k-1。 symsum(S,v) 对通项S求和，指定其中v为变量。且v从0变到v-1。 symsum(S,a,b) 对通项S求和，其中k为变量。且从a变到b。 symsum(S,v,a,b) 对通项S求和，指定其中v为变量。且v从a变到b。 * * 例：键入k=sym(‘k’); symsum(k) 得 ans = 1/2*k^2-1/2*k 又例如：键入 symsum(k^2,0,10)得 ans = 385 又例如：键入 symsum(x^k/sym(k!),k,0,inf)得 ans = exp(x) 这最后的一个例子是无穷项求和。 * * 求极限 limit(P) 表达式P中自变量趋于零时的极限。 limit(P,a) 表达式P中自变量趋于a时的极限。 limit(P,x,a,left) 表达式P中自变量x趋于a时的左极限。 limit(P,x,a,right) 表达式P中自变量x趋于a时的右极限。 * * 例如：键入 P=sym(sin(x)/x); limit(P) 得 ans = 1 键入 P=sym(1/x); limit(P,x,0,right) 得 ans = inf 键入 P=sym((sin(x+h)-sin(x))/h);h=sym(h); limit(P,h,0) 得 ans = cos(x) 键入 v=sym([(1+a/x)^x,exp(-x)]); limit(v,x,inf,left) 得 ans = [ exp(a), 0] * * 线性方程组的求解 线性方程组的形式为A*X=B；其中A至少行满秩。 X=linsolve(A,B) 输出方程的特解X。 例如：键入 A=sym([cos(t),sin(t);sin(t),cos(t)]); B=sym([1;1]); c=linsolve(A,B) c = [ 1/(sin(t)+cos(t))] [ 1/(sin(t)+cos(t))] * * 代数方程的求解 solve(P,v) 对方程P中的指定变量v求解。v可省略。 solve(P1,P2,…，Pn,v1,v2,…，vn) 对方程P1，P2，…Pn中的指定变量v1, v2…vn求解。 * * 例：可输入 solve(p+sin(x)=r) 得： ans =-asin(p-r) 例输入：P1=x^2+x*y+y=3;P2=x^2-4*x+3=0; [x,y]=solve(P1,P2) 得： x = [ 1; 3] y = [1; -3/2] 可输入：P1=a+u^2+v^2=0;P2=u-v=1; [u,v]=solve(P1,P2,u,v) 得： u = [ 1/2+1/2*(-1-2*a)^(1/2); 1/2-1/2*(-1-2*a)^(1/2)] v =[ -1/2+1/2*(-1-2*a)^(1/2) ;-1/2-1/2*(-1-2*a)^(1/2)] * * 对于有些无法求出解析解的非线性方程组，MATLAB只给出一个数值解。这一点可以 从表示解的数字不被方括号括住而确定。例如：键入： [x,y]=solve(sin(x+y)-exp(x)*y=0,x^2-y=2) 得： x = -6.0173272500593065641097297117905 y = 34.208227234306296508646214438330 由于这两个数字没有被[ ]括住，所以它们是数值解。 * * 另外，可利用solve来解线性方程组的通解。例如：键入 P1=2*x1+7*x2+3*x3+x4=6; P2=3*x1+5*x