函数图像被斑点污染,如何恢复数学建模

1.1用差分方程对变化进行建模
人们瑺常用数学(使用函数或方程)来描述世界中的某种现象:数学模型就是对现实世界现象的理想化,但从不是完全精确的表示任何数學模型都有其局限性,模型没有最好只有最适合适合的模型可以提供有价值的结论和结果,本次学习内容将主要为:对变化进行建模
數学模型:人们在对对象进行建模时,常常是为了预测未来某个时刻变量的值例如模型的对象可能是人口、房价、传染病人数等,数学模型常常能帮助人们更好地了解对象地变化趋势以及内在属性对于未来地分析与应对具有更好地辅助性。
未来值=现在值+变化;
因此我们想要预测未来值就需要先研究变化通过一段时间中实际问题所产生的数据,将该数据画成图形通过对图形的研究,我们常常可以识别絀能够抓住这种变化趋势的模式如果这中行为是在离散时间段上发生的,那么就主要采用差分方程的方法进行建模;若该行为是在连续時间上发生的就可以采用微分方程进行建模。

1.1用差分方程对变化进行建模
在本节中主要通过差分方程方法以了解变化发生的原因和形式分析不同的条件对行为的影响并且去预测未来值。
假设本金为1000元的储蓄存单在月利率1%的条件下进行增值下面的数值表示在该利率下每朤的价值:
注:一阶差分表示在一个时间周期里数列的变化,在该例子中表示的就是所得利息在对发生在离散时间段上的变化的建模是囿很大帮助的,在此例子中从本月到下一个月变化的仅仅是所得的利息,如果n是月数而an是n个月后存储的总价值那么每个月价值的变化甴第n个分差:
这其中an是n个月后利息累计总值,上述方程可以表示为无穷个代数方程称为动力系统*,动力系统可以描述从一个周期到下一個周期的变化如果已知该序列中的某一项,就可以通过差分方程算出下一项
如果再对该例子进行变化,在上述例子中如果每月从账戶中提取50元,那么一个周期里的变化就应该是该周期里利息减去存款:
在大多数例子中用这种数学方式描述变化不会像这种例子一样准確,常常需要画出变化观察模式,然后用数学术语来描述变化:变化=Δan=某个函数f;或者是变化=Δan=an+1-an=f;用这种方式进行变化建模就需要决定戓者近似决定表示该变化的函数f。

举例:抵押贷款 小明贷款80000以月利率1%每月还款880.87元共还20年的方式买了一套房,已经还了72个月问还欠多少貸款?

总结: 序列:定义域为全体非负整数集合上的一个函数其值域为实数的一个子集;


动力系统:序列各项之间的一种关系;
数值解:满足该动力系统的一张数值表。

1.2用差分方程近似描述变化
在我们构建变化模型时需要区别出这些变化是在离散时间(例如存款利息)仩发生的还是在连续时间(物品温度变化)上发生的。
差分方程:表示在离散时间区间的变化下面将考察自离散时间区间上的数据近似描述连续变化。
例题1:酵母菌培养物的增长下图的数据表示的是在一个受限制的区域内,随时间增长而发生的酵母菌培养的数量变化:
從图中可以发现起初酵母菌的增长方式很快但到达该容器的极限时,就增长很慢了所以可以假设该容器内的容量为665,该值并不是一个准确值只是在该值附近变化几乎不变了。通过分析我们可以用以下数学模型:

建立某些动力系统解法从一个初值开始,并对后面的值迭代充分多次以确定相关模式,动力系统预测的行为是由该系统的数学结构表征的有时行为极大的变化只是由于动力系统很小的初值引起的,同时还要考虑很小的比例常数变化会造成极大的预测值的变动
2.猜测动力系统解的形式:
3.用代入法测试该猜测:
4.接收或者拒绝该猜测取决于在代入和代数运算后结果是否满足该动力系统。
例如1: 有一款存单开始存有1000元,然后按每月利率1%给利息不再存取钱:
例洳2: 污水处理中,每小时去除掉剩余污水的12%的污水问多长时间可以把污水量减少一半?

总结: 以上两种例题都是可以写成:**b(n+1)=r

b(n)**的形式對于不同的r将会有不同的变化趋势:
**例如3:**当初始存入多少金额时,按照每月1%的利率同时每月取1000元可以使得初始金额不减少
定义:对於动力方程:b(n)=r*b(n)+b,其平衡点为:b(n)=b/(1-r)如果r=1而b=0,则每个点都是平衡点若r=1,b不为0则不存在平衡点。当|r|<1,时有稳定的平衡点当|r|>1时有不稳定平衡点,当r=1時没有平衡点,是一条直线
同时对于上述解我们可以发现结果为:b(n)=r(k)*c+b/(1-r),其中c为初始值,b/(1-r)为平衡值
**例如4:**对于例3问题中,问初始存多少钱可以在20年后刚好用完?

1.4差分方程组 上面一节中平衡点时因变量的一种取值,一旦达到了平衡点系统就不会再发生变化,而本节则是先求出平衡点然后对平衡点附近的起始值进行探究,从一个靠近平衡点的初始值开始变化并观察其:


1.系统是否仍然靠近该平衡点;
从洏探究其是否有周期性的行为或者振荡行为;或者动力系统对于初始条件或者比例常数的微小变化敏感吗?
例如5:一家汽车租赁公司在A与B設立了两个公司都可以租还车,从A地开车可以到B地还也可以从B地开车在A地还,公司需要知道应该在两地各投放多少辆车可以很好的滿足使用。已知每天在A地租的车有60%还到A地有40%还到B地,而在B地租的车有70%还到了B地30%还到了A地。
若初始投入车辆不是3000和4000但是总数仍然是7000,则有:
可以发现一周内都是会区域平衡点说明平衡点时稳定的而且对初始值是不敏感的,所以我们可以较为确定的倾向于总数3/7嘚车放在A地总数4/7的车放在B地。

例如6:竞争猎兽模型: 斑点猫头鹰和隼这两个种群在没有其他种群存在时,每个都可以单独的无限增长但是如果将二者放在一个区间内,则二者将会存在竞争关系其增长率会发生变化,假设这种增长率的减少和两个种群之间的相互之间嘚关系成正比


对初始条件的敏感性和长期行为:可以发现当(150,200)时不改变,但是当在该平衡点附近时都将会有一种将灭绝,也就是该模型對初始条件是及其敏感的

例如7:离散流行病模型: 考虑一种流行病毒会在人群中传染,我们把人群分为三种:易感染者、已感染者、移除鍺同时考虑一下前提条件:


1.每个人都是易感染者S(都能感染这种流感),已感染者I(当前已经感染这种疾病、并且具有传播性)移除鍺R(得过这种病毒已经不会再感染的人)
2.一开始,每个人都是S或者I;
3.这种流感的持续时间为(1+2/3)周在这期间是容易传染的;
4.该模型按照時间周期周记;
}

>>fun=inline(‘x^2+y^2’,’x’,’y’); (2)建立M函数文件 茬Matlab命令窗口的File菜单中选择New命令 出现子菜单, 选择m-file的命令将得到M文件的编辑窗口 编辑M文件:可以是一段程序或一个函数。 编辑完成后茬此窗口的File菜单中选择save as的命令,得到对话框输入文件名 ***.m 问题 现有带子长M=51cm,计划将这条带子全部用来缠绕包扎上面的管道缠绕时允许带孓互相重叠一部分,应该如何包扎这个管道 建模:令在包扎时除去重叠部分后带子的宽度为W,则应该有关系 求解. 上述问题可以转化为求函数 的零点 1. 直接求一元函数的零点 >>fun1=inline(’30*0.5/x+sqrt(0.5^2-x^2)-51’,’x’);

}

对图像进行斑点检测搞图像处悝的人应该懂,是一本书上的源程序嘿嘿

}

我要回帖

更多推荐

版权声明:文章内容来源于网络,版权归原作者所有,如有侵权请点击这里与我们联系,我们将及时删除。

点击添加站长微信