matlab系统仿真课设.docx

《系统仿真与matlab》综合试题题 目： 香烟过滤嘴问题编 号： 11 姓 名：班 级：学 号：联系方式：成 绩： 目录1.1、问题描述………………………………………………………………………………………11.2、模型假设及要求………………………………………………………………………………12.1、模型构建………………………………………………………………………………………22.2、模型实现………………………………………………………………………………………33.1、运行指南………………………………………………………………………………………43.2、实例分析………………………………………………………………………………………54.程序代码………………………………………………………………………………………….71.1问题描述 香烟制造商既要满足瘾君子的需要，又要顺应减少吸烟危害的潮流．还要获取丰厚的利润．于是普遍地在香烟上安装了过滤嘴．过滤嘴的作用倒底有多大，与使用的材料和过滤嘴的长度有什么关系，要从定量的角度回答这些问题就要建立一个描述吸烟过程的数学模型．吸烟时毒物吸入人体的过程大致是这样的：毒物基本上均匀地分布在烟草中．吸烟时点燃处的烟草大部分化为烟雾，毒物由烟雾携带着一部分直接进入空中，另一部分沿香烟穿行．在穿行过程中又部分地被未点燃的烟草和过滤嘴吸收而沉积下来，剩下的进入人体．被烟草吸收而沉积下来的那一部分毒物，当香烟燃烧到那里的时候又通过烟雾部分进入空气，部分沿香烟穿行，这个过程一直继续到香烟燃烧至过滤嘴处为止．于是我们看到，原来分布在烟草中的毒物除了进入空气和被过滤嘴吸收的一部分外，剩下的全都被人体吸入。实际的吸烟过程非常复杂并且因人而异．为了能建立一个初步的模型，可以认为毒物随烟雾进入空气和沿香烟穿行的数量比例、烟雾穿行的速度、过滤嘴和烟草对毒物的吸收率等在吸烟过程中都是常数．1.2模型假设及要求模型假设： 1．烟草和过滤嘴的长度分别是l1和l2,香烟总长l=l1+l2．毒物M(毫克)均匀分布在烟草中．密度为w0=M/l1.2．点燃处毒物随烟雾进入空气和沿香烟穿行的数量比例是a’:a，a’+a=1. 3．未点燃的烟草和过滤嘴对随烟雾穿行的毒物的吸附率（单位时间内毒物被吸收的比例)分别是常数b和β．4．烟雾沿香烟穿行的速度是常数v，香烟燃烧速度是常数u，且vu.5．将一支烟吸完后毒物进入人体的总量(不考虑从空气的烟雾中吸入的)记作Q。仿真要求： 系统输入为烟草和过滤嘴的长度l1和l2,毒物总质量M，点燃处毒物随烟雾进入空气和沿香烟穿行的数量比例a’:a，点燃的烟草和过滤嘴对随烟雾穿行的毒物的吸附率b和β，烟雾沿香烟穿行的速度v，香烟燃烧速度u．系统输出为毒物进入人体的量Q。要求有输入、输出界面及仿真过程。1.2.1模型构建 设在t=0的时刻点燃香烟，在过滤嘴末端截面单位时间的流量为q(t)，则由题意易得：Q=∫q(t)dt（0tl1/u），故而以下对q（t）进行推算。 设距离烟草前端距离为x的界面为q（x,t）,设各点烟雾毒物的密度为p（x,t）,在点燃的瞬间有如下： b∫p(x,0)dx 0xl1 积分均为由x至x+△x q(x,0)-q(x+△x，0)= 此式由物质守恒定律得 β∫p(x,0)dx l1xl取△x趋近于0并将q(x,0)=v*p(x,0)带入可得： bp（x）/v 0xl1 dp/dx= βp(x)/v l1xl设各点毒物密度为w（x,t）由题意知点燃瞬间存在q(0,0)=auw(0,0),故而有p（0,0)=auw(0,0)/v。由上式微分以及p(0,0)可解得在0xl1时p(x,0)的表达式： p(x,0)=(auw(x,0)/v)*exp(-bx/v) 0xl1由函数连续性可求得p(l1,0)=(auw(x,0)/v)*exp(-bx/v),继而由p（l1,0）与上式微分求l1xl时的表达式： p(x,0)= (auw(x,0)/v)*exp((((β-b)*l1-β*x))/v) l1xl故而有过滤嘴端口处:q(l,0)=v* p(l,0)= auw(l,0)*exp((((β-b)*l1-β*l))/v)继而求q(l,t): 先设h（t）为t时刻点燃时单位时间内烟草释放的有毒物质的量，h(t)=uw(x,t)=uw(ut,t),2.由p(x,0)的表达式变换可得： (ah(t)/v)*exp(（-bx+but）/v) 0xl1 p(x,t)= (ah(t)/v)*exp(((β-b)*l1+but-βx)/v) l1xl故而q(l,t)= ah(t) *exp(((β-b)*l1+but-βl)/v),因此，求出h(t)即可得