第八章matlab解方程.pptVIP

线性方程组求解 非线性方程的根 非线性方程的根 例： [x,f,h]=fsolve(f,x0)返回一元或者多元函数f在x0附近的一个零点，其中x0为迭代初值，f返回f在x0的函数值，应该接近0；h返回值如果大于0，说明计算结果可靠，否则计算结果不可靠。 例 求函数 例 求下列方程组在原点附近的解 Matlab 符号方程求解器 Matlab 符号方程求解器 求解方程函数小结 * 线性方程组求解 linsolve(A,b)：解线性方程组 例：解方程组 A=[1 2 –1; 1 0 1; 1 3 0]; b=[2;3;8]; x=linsolve(A,b) b是列向量！ X=A\b Matlab 非线性方程的数值求解 fzero(f,x0)：返回一元函数 f在 x0 附近的根。 方程可能有多个根，但 fzero 只给出距离 x0 最近的一个 fzero 先找出一个包含 x0 的区间，使得 f 在这个区间两个端点上的函数值异号，然后再在这个区间内寻找方程 f=0 的根；如果找不到这样的区间，则返回 NaN。 x0 是一个标量，不能缺省 由于 fzero 是根据函数是否穿越横轴来决定零点，因此它无法确定函数曲线仅触及横轴但不穿越的零点，如 |sin(x)| 的所有零点。 fzero 的另外一种调用方式 fzero(f,[a,b]) 方程在 [a,b] 内可能有多个根，但 fzero 只给出一个 求方程 f=0 在 [a,b] 区间内的根。 参数 f 可通过以下方式给出： fzero(x^3-3*x+1,2); f=inline(x^3-3*x+1); fzero(f,2) fzero(@(x)x^3-3*x+1,2); f 不是方程！也不能使用符号表达式！ fzero(sin(x),10) fzero(@sin,10) fzero(x^3-3*x+1,1) fzero(x^3-3*x+1,[1,2]) fzero(x^3-3*x+1=0,1) X fzero(x^3-3*x+1,[-2,0]) f=inline(x^3-3*x+1); fzero(f,[-2,0]) fun=inline(‘x*sin(x^2-x-1)’,’x’) fplot(fun,[-2,-0.1]);grid on x1=fzero(fun,[-1,-1.2]),x2=fzero(fun,[-1.2,-0.1]) 或x1=fzero(fun,-1.6),x2=fzero(fun,-0.6) 或[x1,f1,h1]=fsolve(fun,-1.6), [x2,f2,h2]=fsolve(fun,-0.6) fun=inline(‘[-x(1)+4*x(2)+x(1)^2/8,4*x(1)-x(2)+exp(x(1))/10-1]’) [x,f,h]=fsolve(fun,[0,0]) fzero ,fsolve只能求实根 s=solve(f,v)：求方程关于指定自变量的解； s=solve(f)：求方程关于默认自变量的解。 f 可以是用字符串表示的方程，或符号表达式； 若 f 中不含等号，则表示解方程 f=0。 solve 例：解方程 x^3-3*x+1=0 syms x; f=x^3-3*x+1; s=solve(f,x) s=solve(x^3-3*x+1,x) s=solve(x^3-3*x+1=0,x) solve 也可以用来解方程组 solve( f1 , f2 , ... , fN , v1 , v2 , ... , vN) 求解由 f1 , f2 , ... , fN 确定的方程组关于 v1 , v2 , ... , vN 的解 例：解方程组 [x,y,z]=solve(x+2*y-z=27,x+z=3, ... x^2+3*y^2=28,x,y,z) 输出变量的顺序要书写正确！ solve 在得不到解析解时，会给出数值解。 例：解方程组 [x,y,z]=solve(x+2*y-z=27,x+z=3, ... x^5+3*y^2=28,x,y,z) roots(p)：多项式的所有零点，p 是多项式系数向量。 fzero(f,x0)：求 f=0 在 x0 附近的根，f 可以使用 inline、字符串、或 @，但不能是方程或符号表达式！ solve(f,v)：求方程关于指定自变量的解，f 可以是用字符串表示的方程、符号表达式或符号方程； solve 也可解方程组(包含非线性)； 得不到解析解时，给出数值解。 linsolve(A,b)：解线性方程组。 [x,f,h]=fsolve