第三讲matlab在线性代数和高等数学中的应用幻灯片.pptVIP

例. 求一个正交变换将二次型 化为标准型 A=[1 1 0 -1;1 1 -1 0;0 -1 1 1;-1 0 1 1] [P,D]= eig(A) syms x1 x2 x3 x4; X=[x1;x2;x3;x4] Y=P*X 2. 矩阵的相似对角化 例. 化上述矩 阵A为对角阵。 D=V-1AV 2.4 线性方程组的直接求解法在MATLAB中的实现 齐次线性方程组的直接求解法 例. 用基础解系表示齐次方程组 的通解。 A=[1 1 1 1 1;3 2 1 1 -3;0 1 2 2 6;5 4 3 3 -1]; B=null(A,r); syms k1 k2 k3; x=k1*B(:,1)+k2*B(:,2)+k3*B(:,3) * 2. 非齐次线性方程组的直接求解法 （1） 求逆法 对AX=B，若|A|不等于0，则解由inv(A)*B获得 （2） 左除法与右除法 例. 求非齐次线性方程组 的通解。 A=[2 3 1 4;1 -2 4 -5;3 8 -2 13; 4 -1 9 -6]; rref(A) …… 3. 符号方程组的求解方法 （1） 线性方程组AX=B的符号解 X=linesolve(A,B) 只给出特解 例. 求 的解。 syms a b c; A=[a 0 0;0 b 0]; B=[1;c] X=linesolve(A,B) （2） 非线性方程组的符号解 [x1,x2,x3,…]=solve(e1,e2,e3,…) 例 解 e1=sym(a+b+x=y); e2=sym(2*a*x-b*y=-1); e3=sym(2*(a+b)=x+y); e4=sym(a*y+b*x=4); [a,b,x,y]=solve(e1,e2,e3,e4) 3 常微分方程的符号解 3 常微分方程的符号解 r=dsolve(eqn1,eqn2,...,var) r=dsolve(eqn1,eqn2,... ,var) 功能：求微分方程或微分方程组的通解 r=dsolve(eqn1,eqn2,..., ‘cond1, cond2,...,var) r=dsolve(eqn1,eqn2,..., ‘cond1,cond2,..., var) 功能：求微分方程(组)满足初始条件的特解 在Matlab中，约定D1表示一次微分，D2表示二次微分，依此类推， Dn表示n次微分，符号Dy相当于Dy/Dt，默认自变量为t。 例1 求解一阶微分方程 及x=0时的特解。 dsolve(Dy=1+y^2, x) dsolve(Dy=1+y^2, y(0)=1 ,x) 例2 求解二阶微分方程 dsolve(x^2*D2y+x*Dy+(x^2-1/2)*y=0,y(pi/2)=2,Dy(pi/2)=-2/pi,x) 4 概率论基础知识及其在MATLAB中的实现 §3.1 随机事件及其概率 rand(N) 返回一个N×N的随机矩阵 rand(N,M) 返回一个N×M的随机矩阵 rand(P1,P2,…,Pn) 返回一个P1×P2×…×Pn的随机矩阵 round(X) 对向量或矩阵X的每个分量四舍五入取整 fix(x) 对x朝零方向取整 floor(x) 对x朝负无穷大方向取整 ceil(x) 对x朝正无穷大方向取整 二项分布 binopdf(X,N,P) 计算二项分布的密度函数 binocdf(X,N,P) 计算二项分