第一章_MATLAB快速入门.pptVIP

图2-1 实例1的图形曲线 2.函数作图 实例2作y=sin(1/x)在[-2，2]上的图形。 fplot(sin(1/x),[-2,2]) 实例3作函数z=x^2+y在([0,8],[0,8])上的图形。 D=[0:0.5:8]; [x,y]=meshgrid(D); surf(x,y,x.^2+y) 实例4作函数z=exp(-(x^2+y^2)/2)在([0,3],[0,3])上的图形。 D=[0:0.1:3]; [x,y]=meshgrid(D); surf(x,y,exp(-(x.^2+y.^2)/2)) 图2-2 实例2的图形曲线 图2-3 实例3的图形曲线 图2-4 实例4的图形曲线 实例5 绘制y=exp(-x^2)在[-3，3]上的条形图 x=-3:.3:3; bar(x,exp(-x.^2)) 图2-5 实例5的图形曲线 2.4 符号微积分初步 表2-9 MATLAB的函数分析函数 limit(F,x,a) 求极限值 dif(F,’x’) 求F对x的一阶导数 dif(F) 求微分 int(F) 求不定积分 int(F,a,b) 求定积分 symsum(S) 对无穷级数表达式S中的符号变量k从0到k-1求和 taylor(f) 将函数表达式展开成泰勒展开式 solve(S) 求解代数方程的解 solve(‘S1’,’S2’) 求解由代数方程S1=0和S2=0组成的方程组的解 dsolve(‘de’,’con’) 求解微分方程在初始条件下的解 实例1求极限 syms x a limit((x^2-1)/(x-1),x,1) ans = 2 实例2 求极限 limit(sin(x)/x,x,a) ans = sin(a)/a 实例3 求函数y=cos(ax^2-1)的一阶导数。 diff(cos(a*x^2),x) ans = -2*sin(a*x^2)*a*x 实例4 求积分 int(x^2/sqrt(1+x^2)) %求解不定积分 ans = 1/2*x*(1+x^2)^(1/2)-1/2*asinh(x) int(x^3*atan(x),0,1) %求定积分 ans = 1/6 实例8 将函数exp(1/2x)展开成幂级数. taylor(exp(1/2*x)) ans = 1+1/2*x+1/8*x^2+1/48*x^3+1/384*x^4+1/3840*x^5 实例9 求解高次方程x^4-2x^3+1=0的解 solve(x^4-2*x^3+1) ans = 1 1/3*(19+3*33^(1/2))^(1/3)+4/3/(19+3*33^(1/2))^(1/3)+1/3 -1/6*(19+3*33^(1/2))^(1/3)-2/3/(19+3*33^(1/2))^(1/3)+1/3+1/2*i*3^(1/2)*(1/3*(19+3*33^(1/2))^(1/3)-4/3/(19+3*33^(1/2))^(1/3)) -1/6*(19+3*33^(1/2))^(1/3)-2/3/(19+3*33^(1/2))^(1/3)+1/3-1/2*i*3^(1/2)*(1/3*(19+3*33^(1/2))^(1/3)-4/3/(19+3*33^(1/2))^(1/3)) 实例11 求微分方程y’’-3y’-4y=3t-1的通解。 syms x y t dsolve(D2y-3*Dy-4*y=3*t-1) ans = exp(4*t)*C2+exp(-t)*C1+13/16-3/4*t 2.5 线性系统与矩阵的特征值 表2-10 MATLAB简单矩阵函数 Det(A) 求行列式 A’ 求矩阵A的转置 rank (A) 求秩 inv (A) 求逆矩阵 pinv (A) 求矩阵A的伪逆 trace (A) 求矩阵A的迹，对角线元素之和 eig (A) 求特征值 1.矩阵的基本运算 举例： A=[1 3 5;4 5 6;5 6 7]; B=[2 4 6;9 6 3;4 3 7]; A*B ans = 49 37 50 77 64 81 92 77 97 B*A ans = 48 62 76 48 75 102 51