MATLAB教程4方程组ppt课件.pptVIP

1.4.4 解方程与方程组 解方程可分为两种 （1）对符号变量方程求解析解； （2）对非符号变量方程求数值解。 A=[3 5 -7;2 -12 3;-1 9 8]; b=[34;-56;27]; x1=A\b x2=linsolve(A,b) ②超定方程（方程个数多于自变量个数）此时不存在唯一解。求解超定方程一般采用最小二乘法。 * * 欢 迎 学 习 微机应用软件 一、 求符号变量方程的解析解 调用格式：solve(eqn1, eqn2 ,… var1 , var2 …) 其中eqni为方程，vari为求解变量，若不指定，则为默认值。 1、求方程的解 例1： syms a b c x f=a*x^2+b*x+c solve(f) ans = [ 1/2/a*(-b+(b^2-4*a*c)^(1/2))] [ 1/2/a*(-b-(b^2-4*a*c)^(1/2))] 也可对指定变量求解 b=solve(f,b) b = -(a*x^2+c)/x 例2：可直接将方程代入 s=solve(cos(2*x)+sin(x)=1) s = [ 0] [ pi] [ 1/6*pi] [ 5/6*pi] 例3：若符号表达式不含等号，则函数solve会自动将表达式转成等号右端为0的符号表达式 solve(p*tan(x)-r) ans = atan(r/p) 例4：符号运算将数值结果用精确值表示，不化为小数。 s=solve(x^3-2*x^2=x-1) %得到精确数值解 vpa(s,10) %化为10位数的小数 又： x=solve(2*sin(3*x-pi/4)=1) x = 5/36*pi 2、使用指令solve解方程组 [x,y]=solve(1/x^3+1/y^3=28,1/x+1/y=4,x,y) [x y]=solve(x+y-98,x^(1/3)+y^(1/3)-2) [u v]=solve(u^3+v^3-98,u+v-2) 或： syms u v [e f]=solve(u^3+v^3-98,u+v-2) [x,y]=solve(1/x^3+1/y^3=28,1/x+1/y=4) syms x y z r1 r2 r3 [x y z]=solve(1/x+1/(y+z)-1/r1, 1/y+1/(z+x)-1/r2,1/z+1/(x+y)-1/r3) syms x y alpha [x,y]=solve(x^2*y^2-1,x-(y/2)-alpha) 二、 求线性方程组的解析解与数值解 对形式为 Ax＝b 的线性方程 组（这里 A 是n×m的矩阵，b 是已知矢量，x 是未知矢量） ， 可分为三种： n=m且A非奇异 恰定方程 可尝试计算精确解 nm 超定方程 可尝试计算最小二乘解 nm 欠定方程 可尝试计算含有最少m的 基解 2、举例： ①恰定方程 1、调用格式： ① linsolve(A,b) ％求解析解 ② x =A \ b ％求数值解 注意：A,b可以是符号矩阵或数值矩阵 = A x = b syms a11 a12 a13 a21 a22 a23 a31 a32 a33 b1 b2 b3; A= [a11 a12 a13; a21 a22 a23 ;a31 a32 a33 ]; b= [b1 ; b2 ; b3]; X=linsolve(A,b) XX=A\b A为符号矩阵 例: x1+2x2=1 2x1+3x2=2 3x1+4x2=3 a=[1 2;2 3;3 4];b=[1;2;3]; x=a\b x = 1.00 0 = AA x = b 例：有一组实验数据，在时间 t 测得数值y，如表所示: 0.50 0.55 0.60 0.63 0.72 0.82 y 2.2