非线性方程求解.ppt

MATLAB求解非线性方程函数非线性方程非线性方程组非线性方程多项式函数rootsfzerofsolve第31页，共48页，星期日，2025年，2月5日多项式求根函数roots多项式的表达式约定如下：对于多项式，用以下行向量表示：这样就把多项式问题转化为向量问题Matlab提供了多种多项式计算函数，如多项式求根函数roots，求多项式的值，polyval；多项式乘法，conv；多项式除法，deconv；多项式微分，polyder；多项式拟合，polyfit第32页，共48页，星期日，2025年，2月5日函数rootsr=roots(c)，用于求解多项式的根其中，行向量c的元素是多项式的系数，按多项式次数降序排列如果c中含有n+1个元素，则多项式为n次roots可以获得多项式的所有根第33页，共48页，星期日，2025年，2月5日例题：求方程的根c=[1-10-1];r=roots(c)r=1.4656-0.2328+0.7926i-0.2328-0.7926ipolyval(c,r(1))ans=-2.5535e-015第34页，共48页，星期日，2025年，2月5日关于非线性方程求解第1页，共48页，星期日，2025年，2月5日课堂练习某蒸馏釜的操作压力位106.7kpa，其中溶液含苯摩尔分数为0.2，甲苯为0.8，求此溶液的泡点和平衡的气相组成。苯甲苯溶液可以看做理想溶液，组分蒸汽压为：a为苯，b为甲苯第2页，共48页，星期日，2025年，2月5日数值计算的误差误差的分类模型误差观测误差截断误差舍入误差计算机对所储存的数据位数有限制计算时只截取有限项第3页，共48页，星期日，2025年，2月5日误差的分类绝对误差相对误差定义：设X*是数X的近似值，如果X*的绝对误差限时他的某一位数位的半个单位，并且从X*左起第一个非零数字到该数位共有N为，则称这个N个数字为X*的有效数字，也称X*近似X时有N为有效数字有效数字第4页，共48页，星期日，2025年，2月5日算法设计的原则避免两个数量级相差很大的数运算算法中应尽量减少运算次数避免两个相近的数相减递推运算中，防止误差的积累避免用绝对值对过小的数作除数第5页，共48页，星期日，2025年，2月5日非线性方程求解非线性方程数值求解基本原理Matlab求解非线性方程函数－roots，fzero，fsolve第6页，共48页，星期日，2025年，2月5日引言在945.36kPa（9.33atm）、300.2K时，容器中充以2mol氮气，试求容器体积。已知此状态下氮气的P-V-T关系符合范德华方程，其范德华常数为a＝4.17atm?L/mol2，b＝0.0371L/mol数学模型：范德华方程变形可得关于V的非线性方程第7页，共48页，星期日，2025年，2月5日非线性方程非线性方程包括：高次代数方程、超越方程及其它们的组合与线性方程相比，非线性方程求解问题无论从理论上还是从计算公式上都要复杂得多对于高次代数方程，当次数4时，则没有通解公式可用，对于超越方程既不知有几个根，也没有同样的求解方式。实际上，对于n≥3代数方程以及超越方程都采用数值方法求近似根。第8页，共48页，星期日，2025年，2月5日非线性方程数值求解原理第9页，共48页，星期日，2025年，2月5日逐步扫描法逐步扫描法效率较低，常用于求根的初始近似值开始结束输入初值A，步长hX1=A，Y1=f(X1)X2=X1+hY2=f(X2)打印A，h，X1Y1*Y20X1=X2Y1=Y2是否第10页，共48页，星期日，2025年，2月5日寻找根区间functionmasearch(fun,a,b,h)n=(b-a)/h;x1=zeros(1,n);x2=zeros(1,n);a1=a;b1=a1+h;k=1;while(b1b)iffeval(fun,a1)*feval(fun,b1)0x1(k)=a1;x2(k)=b1;elsea1=b1;b1=a1+h;continue;enda1=b1;b1=a1+h;k=k+1;endfori=1:kif