东华大学Matlab-ppt-第4章-函数和方程.ppt

MATLAB数学实验;第四章函数和方程;4.1预备知识：零点;4.1预备知识：极值;4.1预备知识：极值;4.1预备知识：最小二乘拟合;4.2函数零点MATLAB指令;MATLAB中一个多项式用系数降幂排列向量来表示。;Fun=inline(‘funstr’,’var’)定义一个inline

函数,其中funstr是函数的表达式,

var是变量名;x=fzero(Fun,x0)返回一元函数Fun的一

个零点,其中Fun为函数句柄、

inline函数或匿名函数。

x0为标量时,

返回函数在x0附近的零点；

x0为向量[a,b]时,

返回在[a,b]中的零点;[x,f,h]=fsolve(Fun,x0)

x:返回多元函数Fun在x0附近的一个零点，其中x0为迭代初值向量;

f:返回Fun在x的函数值,应该接近0;

h:返回值如果大于0，说明计算结果可靠，否那么计算结果不可靠。;例3求函数y=xsin(x2-x-1)在(-2,-0.1)内的零点

fun=inline(x*sin(x^2-x-1),x)

fzero(fun,[-2-0.1]);例4求方程组在原点附近的解;min(y)返回向量y的最小值

max(y)返回向量y的最大值;例5.求二元函数f(x,y)=5-x4-y4+4xy在原点附近的极大值。

解：maxf?min(-f)

x?x(1),y?x(2);4.2最小二乘拟合MATLAB指令;c=lsqnonlin(Fun,c0)使用迭代法搜索最

优参数c.其中Fun是以参数c(可以是向量)

为自变量的函数,表示误差向量

y-f(c,x)(x,y为数据)，

c0为参数c的近似初值(与c同维向量);用lsqnonlin再解例4.2;用lsqcurvefit再解例4.2;迭代法是从解的初始近似值x0(简称初值)开始，利用某种迭代格式xk+1=g(xk),

求得一近似值序列x1,x2,…,xk,xk+1,…

逐步逼近于所求的解?(称为不动点)。

最常用的迭代法是牛顿迭代法，其迭代格式为;牛顿法程序newton1.m;牛顿法程序newton.m;例6求方程x2-3x+ex=2的正根(要求精度?=10-6);2．线性化拟合;不难算出，你向银行总共借了25.2万，

30年内共要还51.696万，

这个案例中贷款年利率是多少呢？;解设xk为第k个月的欠款数，

a为月还款数,r为月利率。;常识上，r应比当时活期存款月利率略高一些。我们用活期存款月利率0.0198/12作为迭代初值，用fzero求解

clear;

fun=inline(25.2*(1+r)^360-((1+r)^360-…

1)/r*0.1436,r)

r=fzero(fun,0.0198/12);

R=12*r

得年利率为5.53%.;月还款计算公式;每次订货需要收取一定量的生产准备费。

没用完的配件，要在仓库里储存一段时间，为此要付出储存费。

假设订货量很小，那么需频繁定货，造成生产准备费的增加；

反之，假设订货量很大，定货周期延长而使生产准备费减少但会造成储存费的增加。

如何确定适宜的订货量？;解先作一些必要的假设将问题简化

1〕汽车工厂对配件的日需求量是恒定的，

每日为r件；

2〕所订配件按时一次性交货，

生产准备费每次k1元；

3〕储存费按当日实际储存量计算，

储存费每日每件k2元;

4〕你的工厂不允许缺货。

设一次订货x件，那么订货周期为T=x/r,

第t天的储存量为q(t)=x-rt,0tT;第t天的储存费为k2q(t)

一个周期的总储存费为

一个周期总费用C(x)=k1+k2x2/(2r)

优化目标是使单位产品费用

f(x)=C(x)/x=k1/x+k2x/(2r)到达最小

由f’(x)=0,即-k1/x2+k2/(2r)=0

得;线性迭代xk+1=axk+b要么收敛于它的不动点x=b/(1-a)，要么趋于无穷大。

不收敛的非线性迭代可能会趋于无穷大，也可能趋于一个周期解，

但也有可能在一个有限区域内杂乱无章地游荡，这类由确定性运动导致的貌似随机的现象称为混沌现象;*昆虫数量的Logistic模型

xk+1=