MATLAB及其在电气工程中的应用 教学课件 作者 苏小林 第三章.pptVIP

第3章 数值计算与分析 课堂练习： 1、有一电路的节点电压方程如下，求该电路的各节点电压。 2、下面是某物理量y在几个不同时间点的观察值记录： t y 0.0 0.82 0.3 0.72 0.8 0.63 1.1 0.60 1.6 0.55 2.3 0.50 用衰减指数函数y=c1+c2e-t建立这组数据的模型。 3、对矩阵进行奇异值分解。 A=[1 3;2 4;5 7;6 8]; 课堂练习：1、求矩阵的特征值和特征向量。 A=[2 -4 -1;4 1 -10;-3 15 -10] [t,y]=ode45(@vdp1,[0,20],[2;0]); 数组y的第1列为y1，第2列为y2，每行为它们在对应时间点的解。解法函数的第1个输入参数为定义微分方程的函数句柄。 (4) 用图示观察解的结果 plot(t,y(:,1),-,t,y(:,2), --) title( Solution of van der Pol Equation,\mu=1 ) xlabel( time t ) ylabel( solution y ) legend( y_1 , y_2 ) 图3-2 微分方程的解 【例3-30】将例3-29的范德波耳微分方程中的m取为1000时，微分方程成为刚性的。求该刚性范德波耳微分方程的解。 (1) 编写一阶微分方程组的函数文件vdp1000 function dydt=vdp1000(t,y) mu=1000; dydt=[y(2);mu*(1-y(1)^2)*y(2)-y(1)]; (2) 针对刚性范德波耳微分方程，采用ode15s解法，时间间隔为 [0,3000] tspan=[0 3000]; y0=[2;0]; [t,y]=ode15s(@vdp1,tspan,y0); 课堂练习： 1、求解积分 的数值。 2、求解二重积分 的数值。 3、求解范德波耳(van der Pol)方程 在初始条件 且 情况下的解，并用图示显示解的结果。 4、将上题的范德波耳微分方程中的取为1000时，微分方程成为刚性的。求该刚性范德波耳微分方程的解。 对于线性方程组Ax=b，若A是对称的正定矩阵，其A=RTR，则线性方程组可以转换为RTRx=b，该方程组的解为x=R\(RT\b)，MATLAB命令为x=R\(R’\b)。它与左除x=A\b所得结果完全相同，但对于n×n的矩阵A，命令chol的计算复杂度是o(n3)，而左除命令“\”的计算复杂度o(n2)。 【例3-15】Cholesky分解及其求解线性方程组示例 A=pascal(4); % 线性方程组的系数矩阵； b=[1;3;5;7]; % 线性方程组的常数项矩阵； R=chol(A) % 对系数矩阵进行Cholesky分解 R = % 分解得上三角矩阵 1 1 1 1 0 1 2 3 0 0 1 3 0 0 0 1 x=A\b; % 用左除法求解线性方程组 xf=R\(R’\b); % 用Cholesky法求解线性方程组 可以验证，两种方法所求结果完全相同，但Cholesky分解法的计算复杂度高于左除法。 2. LU分解 LU分解又被称为高斯消去法，它是将一个m×n或n×n的矩阵A分解为一个内在单位下三角矩阵L (形式上不是单位下三角矩阵，但通过行列交换可变为单位下三角矩阵，定义为单位下三角矩阵与变换矩阵的乘积)和上三角矩阵U的乘积，即A=LU。 在MATLAB中，LU分解的常用命令格式如下： · [L,U]=lu(A) 其中A是任意矩阵，L是内在单位下三角矩阵，U是上三角矩阵； · [L,U,P]=lu(A) 其中A是任意矩阵，L是单位下三角矩阵，U是上三角矩阵，P是变换矩阵，满足的条件为PA=LU。 对于