数学建模matlab_课件_11、12、方程组(更新).pptVIP

求解方程组 1. 矩阵求逆解法 利用求系数矩阵A的逆阵A-1，我们可以得到矩阵求逆解法。对于线性代数方程组Ax=b，等号两侧各左乘A-1，有： A-1Ax=A-1b 由于A-1A=I，故得： x=A-1b 【例1】试用矩阵求逆解法求解矩阵A为系数矩阵的线性代数方程组Ax=b的解。 A=[1 -1 1;5 -4 3;2 1 1]; b=[2;-3;1]; x=inv(A)*b x = -3.8000 1.4000 7.2000 2. 直接解法 对于线性代数方程组Ax=b，我们可以运用左除运算符“\”象解一元一次方程那样简单地求解： x=A\b 当系数矩阵A为N*N的方阵时，MATLAB会自行用高斯消去法求解线性代数方程组。若右端项b为N*1的列向量，则x=A\b可获得方程组的数值解x（N*1的列向量）；若右端项b为N*M的矩阵，则x=A\b可同时获得同一系数矩阵A、M个方程组数值解x（为N*M的矩阵），即x(:,j)=A\b(:,j)，j=1,2,…M。 解法1：分别解方程组 (1)Ax=b1；(2)Ay=b2 A=[1 -1 1;5 -4 3;2 1 1]; b1=[2;-3;1]; b2=[3;4;-5]; x=A\b1 x = -3.8000 1.4000 7.2000 解法2：将两个方程组连在一起求解：Az=b b=[2 3;-3 4;1 -5] z=A\b z = -3.8000 -3.6000 1.4000 -2.2000 7.2000 4.4000 很明显，这里的解z的两个列向量便是前面分别求得的两组解x和y 解： A=[2 3 1;1 -2 4;3 8 -2;4 -1 9]; b=[4;-5;13;-6]; B=[A,b]; rref(B) ans = 1 0 2 -1 0 1 -1 2 （B的行最简形） 0 0 0 0 0 0 0 0 既得 得通解： 习题 1、解方程组Ax＝b，分别用求逆解法与直接解法求其解。 习题 2、求方程组 的通解 到此为止 谢谢? 二、非线性方程组 1、求非线性方程及方程组的解：solve 命令调用格式为： solve(’eqn1’,’eqn2’,...,’eqnN’)：对 N 个方程的默认变量求解； solve(’eqn1’,’eqn2’,...,’eqnN’,’var1,var2,...,varN’)：对 N 个方程的 var1,var2,...,varN 变量求解。要注意变量的英文字母顺序，且在变量前不可有空格； S=solve(’eqn1’,’eqn2’,...,’eqnN’,’var1’,’var2’,...,’varN’)：对 N 个方程的’var2’,...,’varN‘ 变量求解；S 是一个结构数组; [x1,x2,…,xn]=solve(‘eqn1’,‘eqn2’,...,‘eqnN’,‘var1’,‘var2’,...,‘varN’)：对变量 var1,var2,..., varN 求解，求解的结果分别赋给 x1,x2,…,xn(按照变量var1,var2,...varN 在英文字母中的顺序给 x1,x2,…,xn赋值)。 例题 1、求解方程 2、求解方程组 3、求解方程组 2、求解常微分方程 函数 dsolve 用来求常微分方程的符号解。 在符号方程中，用符号表达式中包含的字母“ D”来代替微分运算，符号 D2、D3、…DN 分别对应于第二、第三、…第 N 阶导数。因变量是位于 D后面的变量，缺省的自变量为 t。 dsolve 的调用格式为： S=dsolve(’eqn1’,’eqn2’, ...) 输入参数包含三部分内容：微分方程、初始条件和指定独立变量。输出S 是结构数组。 初始条件或边界条件写成‘y(a)=b’或‘Dy(a)=b’等。a、b 可以是变量使用符以外的其他字符。 当初始条件少于微分方程数时，在所得解中将出现任意常数符 C1，C2，….，解中任意常数符的数目等于所缺少的初始条件数。 当无输出参数时，MATLAB 工作内存中在 y1，y2，…定义的输出参数中保存计算结果。 例题、 1、求解微分方程 2、求解微分方程 ，其中初始条件为y(0)=1,y’(0)=0 3、求解微分方程组