第讲方程的图形法迭代法直接法(第次课)教材.ppt

第二讲 方程及方程组解法(下) 内容：本讲主要讲解方程(组)求解的MATAB直接 解法，补充三种迭代数列图示方法 目的：掌握MATLAB方程（组）求解的相关函数 要求：能够处理带应用背景的方程问题 MATLAB软件直接求解函数： solve fsolve fzero roots 迭代数列图示法： 线性联结图 蛛网图 费根鲍姆图 迭代的复杂行为：分岔 倍周期 混沌 方程(组)直接求解函数：solve 利用MATLAB内置函数，我们可以直接对一些方程或方程组进行求解，免去书写代码的时间，优点是快速高效，缺点是缺乏自编函数的灵活性 solve 对单变量方程f(x)=0求解(解析解) ： 例1 求解方程 ax2+bx+c=0 x=solve(a*x^2+b*x+c) 或者 x=solve(a*x^2+b*x+c=0) pretty(x) 方程(组)直接求解函数：solve 方程(组)直接求解函数：solve 方程(组)直接求解函数：fsolve 方程(组)直接求解函数：fsolve 方程(组)直接求解函数：fsolve 方程(组)直接求解函数：fsolve 方程直接求解函数：fzero 方程直接求解函数： roots 第3章(补充)收敛与混沌 ( chaos ) 混沌现象的最佳注解-蝴蝶效应 混沌现象的最佳注解-蝴蝶效应 分形(Fractal)部分与整体的相似 分形几何学又称大自然几何学 迭代数列图示法-线性联结图 迭代数列图示法-线性联结图 迭代数列图示法-蛛网图 迭代数列图示法-蛛网图 迭代数列图示法-费根鲍姆图 迭代数列图示法-费根鲍姆图 迭代产生的分岔与混沌观察实验 分形的世界——分形图案欣赏 实验一: 编程实现课本引例实验二: 编写二分法和牛顿法通用程序That’s all~3Q! * * 第二讲 方程的图形法 迭代法直接法 solve 对单变量方程f(x)=0求解(数值解) ： 例2 求解方程x3-2x2=x-1 x=solve(x^3-2*x^2=x-1) double(x) %实值化处理 fplot(x^3-2*x^2-x+1,[-5,5]); set(findobj(gca,type,line,color,b),linewidth,2); %加粗 hold on; grid on; axis([-1 3 -10 10]); line([-5 5],[0 0],color,r,linewidth,2); %作x轴 solve 对单变量方程f(x)=0求解(无穷解?) ： 例3 求解方程tan(x)=sin(x) x=solve(tan(x)=sin(x)) fplot( tan(x)-sin(x) ,[-10*pi,10*pi]); set(findobj(gca,type,line,color,b),linewidth,2); %加粗 hold on; grid on; line([-40 40],[0 0],color,r,linewidth,2); %作x轴 solve 对多变量方程fn(xn)=0求解(解析解)： 例4 求解方程组 x2y2=0 x-1/2y=b s=solve(x^2*y^2,x-y/2=b); [s.x, s.y] %输出结果整理 solve 对多变量方程f(xn)=0求解(数值解)： 例5 求解方程组 x2y2-2x-1=0 x2-y2-1=0 s=solve(x^2*y^2-2*x-1,x^2-y^2-1=0); xy=[s.x, s.y], double(xy)%输出结果整理+数值化 ezplot(x^2*y^2-2*x-1); hold on; set(findobj(gca,type,line,color,b), color,r); %改色 ezplot(x^2-y^2-1); grid on; set(findobj(type,line), linewidth,2); %加粗 注意: 须小心使用solve，实际中最好结合多种方法 fsolve 对非线性方程组的求解（以下为标准语法） 例6 求解非线性方程组 解法1: 直接用solve函数求解 s=solve(sin(x)+y^2+log(z)-7,3*x+2^y-z^3+1,x+y+z-5); [s.x,s.y,s.z] 解法2写法1: 编写被调函数，供fsolve调用求解 Step1：编写独立被调函数 nxxf.m function eq=nxxf(x) eq(1)= sin(x(1))+x(2)^2+log(