在高数中的应用.pptVIP

解线性方程组 求极限 求微分 求积分 作图 线性方程组解的理论 若系数矩阵的秩r=n(为未知量的个数),则有唯一解； 若系数矩阵的秩rn,则可能无穷解； 非齐次方程的通解＝ 齐次方程的通解＋非齐次方程的特解 1.1 线性方程唯一解或者特解的解法 [调用格式] X＝A\b; X=inv(A)*b X为方程AX=b的解。 A=[5 6 0 0;1 5 6 0;0 1 5 6; 0 0 1 5 ]; b=[1 0 0 0]; rank(A), X=A\b, inv(A)*b, [a,b]=eig(A)求特征值和特征向量 1.2 齐次线性方程通解的解法 [调用格式] z=null(A,r) Z的列向量为方程AX=0的有理基。 1.3 非齐次线性方程通解的解法 求解步骤： 判断方程AX=b的是否有解，如有则进行下一步计算； 求方程AX=0的通解。 求方程AX=b的特解。 AX=b的通解＝ AX=0通解＋ AX=b的特解 符号计算 Matlab的命令中，如果其操作对象不是数值，而是符号，相应的计算称为符号计算，例如，当a,b是变量时，计算积分 就是符号计算。 基本符号的定义： [调用格式] F=sym(‘a’,‘b’,’c’) Syms a b c 2.1 求极限 [调用格式] g=limit(f,x,a,’left’) G为结果表达式，是函数f在x=a处的左极限，f可以是一个矩阵或者函数，求右极限用’right’,求无穷大处的极限用‘inf’ 3.1 求微分 [调用格式] g=diff(f,x,k) G为结果表达式，是函数f对x处的k阶微分，默认为1阶微分。 4.1 求积分 [调用格式] g1=inf(f,x,) g2=inf(f,x,a,b) g1和g2为结果表达式，g1是函数f对x的不定 积分，g2是函数f对x在[a,b]上的定积分。 6.1 Taylor展开 [调用格式] g=taylor(f,x,a,k) G为结果表达式，将函数f在x=a处展开到第k-1项 MATLAB在高等数学中的应用 MATLAB在高等数学中的应用 §1 Matlab解线性方程组 A=[1 1 -3 -1;3 -1 -3 4;1 5 -9 -8], z=null(A,r),%求有理基 法一： A=[1 1 -3 -1;3 -1 -3 4;1 5 -9 -8], b=[1 4 0]; B=[A,b]; format rat %指定有理数格式输出 if rank(A)==rank(B) %判断是否有解 z=null(A,r),%求有理基 c=A\b, %求特解 end syms k1 k2 %定义符号变量 x=k1*z(:,1)+k2*z(:,2)+c, %通解的表示 法二： A=[1 1 -3 -1;3 -1 -3 4;1 5 -9 -8], b=[1 4 0]; B=[A,b]; C=rref(B), %求增广矩阵的行最简行，可得最简同解方程组 syms k1 k2 %定义符号变量 x=k1*C(:,3)+k2*C(:,4)+ C(:,5) %通解的表示 §2 Matlab求极限 clear; syms x a; R1=a*(1+1/x)^x; R2=sin(x)/x; limit(R1,x,inf), limit(R2,x,inf), clear; syms x; R1=(1+2*x)^(3/sin(x)); limit(R1,x,0), §3 Matlab求微分 clear; syms x ; f=x^3+4*cos(x)-sin(pi/2); diff(f) §4 Matlab求积分 clear; syms x y z; R1=int(x/(1+z^2),z), R2=int(x*log(1+x),0,1), §5 Matlab作图 §6 Matlab在高数中的其他应用 * *