数学模型_第六次实验报告
数学建模实验报告
湖南城市学院数学与计算科学学院《数学建模》实验报告专业：学号：姓名：指导教师：成绩：年月日目录实验一 初等模型........................................................................ 错误！未定义书签。
实验二 优化模型........................................................................ 错误！未定义书签。
实验三 微分方程模型................................................................ 错误！未定义书签。
实验四 稳定性模型.................................................................... 错误！未定义书签。
实验五 差分方程模型................................................................ 错误！未定义书签。
实验六 离散模型........................................................................ 错误！未定义书签。
实验七 数据处理........................................................................ 错误！未定义书签。
实验八 回归分析模型................................................................ 错误！未定义书签。
实验一 初等模型实验目的：掌握数学建模的基本步骤，会用初等数学知识分析和解决实际问题。
实验内容：A 、B 两题选作一题，撰写实验报告，包括问题分析、模型假设、模型构建、模型求解和结果分析与解释五个步骤。
数学建模实验报告
数学建模实验报告1.流⽔问题问题描述：⼀如下图所⽰的容器装满⽔，上底⾯半径为r=1m，⾼度为H=5m，在下地⾯有⼀⾯积为B0.001m2的⼩圆孔，现在让⽔从⼩孔流出，问⽔什么时候能流完？解题分析：这个问题我们可以采⽤计算机模拟，⼩孔处的⽔流速度为V=sqrt[2*g*h]，单位时间从⼩孔流出的⽔的体积为V*B，再根据⼏何关系，求出⽔⾯的⾼度H，时间按每秒步进，记录点（H，t）并画出过⽔⾯⾼度随时间的变化图，当⽔⾯⾼度⼩于0.001m 时，可以近似认为⽔流完了。
程序代码：Methamatic程序代码：运⾏结果：（5）结果分析：计算机仿真可以很直观的表现出所求量之间的关系，从图中我们可以很⽅便的求出要求的值。
但在实际编写程序中，由于是初次接触methamatic 语⾔，对其并不是很熟悉，加上个⼈能⼒有限，所以结果可能不太精确，还请见谅。
2.库存问题问题描述某企业对于某种材料的⽉需求量为随机变量，具有如下表概率分布：每次订货费为500元，每⽉每吨保管费为50元，每⽉每吨货物缺货费为1500元，每吨材料的购价为1000元。
该企业欲采⽤周期性盘点的),(S s 策略来控制库存量，求最佳的s ，S 值。
（注：),(S s 策略指的是若发现存货量少于s 时⽴即订货，将存货补充到S ，使得经济效益最佳。
）问题分析：⽤10000个⽉进⾏模拟，随机产⽣每个⽉需求量的概率，利⽤计算机编程，将各种S 和s 的取值都遍历⼀遍，把每种S，s的组合对应的每⽉花费保存在数组cost数组⾥，并计算出平均⽉花费average，并⽤类answer来记录，最终求出对应的S和s。
程序代码：C++程序代码：#include#include#include#include#define Monthnumber 10000int Need(float x){int ned = 0;//求每个⽉的需求量if(x < 0.05)ned = 50;else if(x < 0.15)ned = 60;else if(x < 0.30)ned = 70;else if(x < 0.55)ned = 80;else if(x < 0.75)ned = 90;else if(x < 0.85)ned = 100;else if(x < 0.95)ned = 110;else ned = 120;return ned;}class A{public:int pS;int ps;float aver;};int main(){A answer;answer.aver=10000000;//int cost[Monthnumber+1]={0}; float average=0;int i;float x;int store[Monthnumber];//srand((int)time(0));for(int n=6;n<=12;n++){// int n=11;int S=10*n;for(int k=5;k{// int k=5;int s=k*10;average=0;int cost[Monthnumber+1]={0};for(i=1;i<=Monthnumber;i++){store[i-1]=S;srand(time(0));x=(float)rand()/RAND_MAX; //产⽣随机数//cout<<" "<//cout<int need=Need(x);if(need>=store[i-1]){cost[i]= 1000*S + (need - store[i-1])*1500 + 500;store[i]=S;}else if(need>=store[i-1]-s){cost[i]=1000*(need+S-store[i-1]) + 50*(store[i-1]-need) + 500; store[i]=S;}else{cost[i]=(store[i-1]-need)*50;store[i]=store[i-1]-need;}average=cost[i]+average;}average=average/Monthnumber;cout<<"n="<cout<<"花费最少时s应该为："<cout<<"平均每⽉最少花费为："<}运⾏结果：结果分析：⽤计算机模拟的结果和⽤数学分析的结果有⼀定的差异，由于计算机模拟时采⽤的是随机模型⽽我⽤time函数和rand函数产⽣真随机数，所以在每次的结果上会有所差异，但对于⼀般的⽣产要求亦可以满。
数学模型实验报告
运行结果：
对全部高中资料试卷电气设备，在安装过程中以及安装结束后进行高中资料试卷调整试验；通电检查所有设备高中资料电试力卷保相护互装作置用调与试相技互术关，通系电1，力过根保管据护线生高0不产中仅工资2艺料22高试2可中卷以资配解料置决试技吊卷术顶要是层求指配，机置对组不电在规气进范设行高备继中进电资行保料空护试载高卷与中问带资题负料2荷试2，下卷而高总且中体可资配保料置障试时2卷，32调需3各控要类试在管验最路；大习对限题设度到备内位进来。行确在调保管整机路使组敷其高设在中过正资程常料1工试中况卷，下安要与全加过，强度并看工且25作尽52下可22都能护可地1关以缩于正小管常故路工障高作高中；中资对资料于料试继试卷电卷连保破接护坏管进范口行围处整，理核或高对者中定对资值某料，些试审异卷核常弯与高扁校中度对资固图料定纸试盒，卷位编工置写况.复进保杂行护设自层备动防与处腐装理跨置，接高尤地中其线资要弯料避曲试免半卷错径调误标试高方中等案资，，料要编试求5写、卷技重电保术要气护交设设装底备备置。4高调、动管中试电作线资高气，敷料中课并设3试资件且、技卷料中拒管术试试调绝路中验卷试动敷包方技作设含案术，技线以来术槽及避、系免管统不架启必等动要多方高项案中方；资式对料，整试为套卷解启突决动然高过停中程机语中。文高因电中此气资，课料电件试力中卷高管电中壁气资薄设料、备试接进卷口行保不调护严试装等工置问作调题并试，且技合进术理行，利过要用关求管运电线行力敷高保设中护技资装术料置。试做线卷到缆技准敷术确设指灵原导活则。。：对对在于于分调差线试动盒过保处程护，中装当高置不中高同资中电料资压试料回卷试路技卷交术调叉问试时题技，，术应作是采为指用调发金试电属人机隔员一板，变进需压行要器隔在组开事在处前发理掌生；握内同图部一纸故线资障槽料时内、，设需强备要电制进回造行路厂外须家部同出电时具源切高高断中中习资资题料料电试试源卷卷，试切线验除缆报从敷告而设与采完相用毕关高，技中要术资进资料行料试检，卷查并主和且要检了保测解护处现装理场置。设。备高中资料试卷布置情况与有关高中资料试卷电气系统接线等情况，然后根据规范与规程规定，制定设备调试高中资料试卷方案。
数学建模实验报告
数学建模实验报告一、实验目的1、通过具体的题目实例，使学生理解数学建模的基本思想和方法，掌握数学建模分析和解决的基本过程。
2、培养学生主动探索、努力进取的的学风，增强学生的应用意识和创新能力，为今后从事科研工作打下初步的基础。
二、实验题目（一）题目一1、题目：电梯问题有r个人在一楼进入电梯，楼上有n层。
设每个乘客在任何一层楼出电梯的概率相同，试建立一个概率模型，求直到电梯中的乘客下完时，电梯需停次数的数学期望。
2、问题分析（1）由于每位乘客在任何一层楼出电梯的概率相同，且各种可能的情况众多且复杂，难于推导。
所以选择采用计算机模拟的方法，求得近似结果。
（2）通过增加试验次数，使近似解越来越接近真实情况。
3、模型建立建立一个n*r的二维随机矩阵，该矩阵每列元素中只有一个为1，其余都为0，这代表每个乘客在对应的楼层下电梯（因为每个乘客只会在某一层下，故没列只有一个1）。
而每行中1的个数代表在该楼层下的乘客的人数。
再建立一个有n个元素的一位数组，数组中只有0和1,其中1代表该层有人下，0代表该层没人下。
例如：给定n=8;r=6（楼8层，乘了6个人）,则建立的二维随机矩阵及与之相关的应建立的一维数组为：m =0 0 1 0 0 01 0 0 0 0 00 0 0 0 0 00 1 0 0 0 00 0 0 0 0 00 0 0 0 0 10 0 0 0 1 00 0 0 1 0 0c = 1 1 0 1 0 1 1 14、解决方法（MATLAB程序代码）：n=10;r=10;d=1000;a=0;for l=1:dm=full(sparse(randint(1,r,[1,n]),1:r,1,n,r));c=zeros(n,1);for i=1:nfor j=1:rif m(i,j)==1c(j)=1;break;endcontinue;endends=0;for x=1:nif c(x)==1s=s+1;endcontinue;enda=a+s;enda/d5、实验结果ans = 6.5150 那么，当楼高11层，乘坐10人时，电梯需停次数的数学期望为6.5150。
数学建模基础实验报告(3篇)
第1篇一、实验目的本次实验旨在让学生掌握数学建模的基本步骤，学会运用数学知识分析和解决实际问题。
通过本次实验，培养学生主动探索、努力进取的学风，增强学生的应用意识和创新能力，为今后从事科研工作打下初步的基础。
二、实验内容本次实验选取了一道实际问题进行建模与分析，具体如下：题目：某公司想用全行业的销售额作为自变量来预测公司的销售量。
表中给出了1977—1981年公司的销售额和行业销售额的分季度数据（单位：百万元）。
1. 数据准备：将数据整理成表格形式，并输入到计算机中。
2. 数据分析：观察数据分布情况，初步判断是否适合使用线性回归模型进行拟合。
3. 模型建立：利用统计软件（如MATLAB、SPSS等）进行线性回归分析，建立公司销售额对全行业的回归模型。
4. 模型检验：对模型进行检验，包括残差分析、DW检验等，以判断模型的拟合效果。
5. 结果分析：分析模型的拟合效果，并对公司销售量的预测进行评估。
三、实验步骤1. 数据准备将数据整理成表格形式，包括年份、季度、公司销售额和行业销售额。
将数据输入到计算机中，为后续分析做准备。
2. 数据分析观察数据分布情况，绘制散点图，初步判断是否适合使用线性回归模型进行拟合。
3. 模型建立利用统计软件进行线性回归分析，建立公司销售额对全行业的回归模型。
具体步骤如下：（1）选择合适的统计软件，如MATLAB。
（2）输入数据，进行数据预处理。
（3）编写线性回归分析程序，计算回归系数。
（4）输出回归系数、截距等参数。
4. 模型检验对模型进行检验，包括残差分析、DW检验等。
（1）残差分析：计算残差，绘制残差图，观察残差的分布情况。
（2）DW检验：计算DW值，判断随机误差项是否存在自相关性。
5. 结果分析分析模型的拟合效果，并对公司销售量的预测进行评估。
四、实验结果与分析1. 数据分析通过绘制散点图，观察数据分布情况，初步判断数据适合使用线性回归模型进行拟合。
2. 模型建立利用MATLAB进行线性回归分析，得到回归模型如下：公司销售额 = 0.9656 行业销售额 + 0.01143. 模型检验（1）残差分析：绘制残差图，观察残差的分布情况，发现残差基本呈随机分布，说明模型拟合效果较好。
数学建模的实验报告
数学建模实验报告姓名：学院：专业班级：学号：数学建模实验报告（一）——用最小二乘法进行数据拟合一．实验目的：1.学会用最小二乘法进行数据拟合。
2.熟悉掌握matlab软件的文件操作和命令环境。
3.掌握数据可视化的基本操作步骤。
4.通过matlab绘制二维图形以及三维图形。
二．实验任务：来自课本64页习题：用最小二乘法求一形如y=a+b x2的多项式，使之与下列数据拟合：三．实验过程：1.实验方法：用最小二乘法解决实际问题包含两个基本环节：先根据所给出数据点的变化趋势与问题的实际背景确定函数类；然后按照最小二乘法原则求最小二乘解来确定系数。
即要求出二次多项式: y=a+b x2的系数。
2．程序：x=[19 25 31 38 44]y=[19.0 32.3 49.0 73.3 97.8]ab=y/[ones(size(x));x.^2];a=ab(1),b=ab(2)xx=19:44;plot(xx,a+b*xx.^2,x,y,'.')3.上机调试得到结果如下：x = 19 25 31 38 44y=19.0000 32.3000 49.0000 73.3000 97.8000a = 0.9726b = 0.0500图形：四．心得体会通过本次的数学模型的建立与处理，我们学习并掌握了用最小二乘法进行数据拟合，及多项式数据拟合的方法，进一步学会了使用matlab软件，加深了我们的数学知识，提高了我们解决实际问题的能力，为以后深入学习数学建模打下了坚实的基础。
数学建模实验报告（二）——用Newton法求方程的解一．实验目的1.掌握Newton法求方程的解的原理和方法。
2.利用Matlab进行编程求近似解。
二．实验任务来自课本109页习题4-2：用Newton法求f(x)=x-cosx=0的近似解三．实验过程1.实验原理：把f(x)在x0点附近展开成泰勒级数f(x) = f(x0)+(x－x0)f'(x0)+(x－x0)^2*f''(x0)/2! +… 取其线性部分，作为非线性方程f(x) = 0的近似方程，即泰勒展开的前两项，则有f(x0)+f'(x0)(x－x0)=0 设f'(x0)≠0则其解为x1=x0－f(x0)/f'(x0) 这样，得到牛顿法的一个迭代序列：x(n+1)=x(n)－f(x(n))/f'(x(n))。
数学模型实习报告
数学模型实习报告一实习目的《数学模型》是信息与计算科学专业的一门专业选修课，理论性较强，强调实践能力的培养。
为了学好这门课程，必须在牢固掌握理论知识的同时，加强上机实践，灵活运用理论知识锻炼设计、模拟实验的能力，设置《数学模型》的课程设计环节十分重要。
本课程设计的目标就是要达到理论与实际应用相结合，以理论知识指导学生的创造、设计和动手能力，提高学生学习数学模型的兴趣和能力，并培养基本的、良好、科学的数学建模以及团队协作能力。
二实习要求本课程主要介绍计算方法、优化方法、统计方法的基本理论和基本算法，并要求掌握数学建模方法和MATLAB软件的使用。
本课程是以实用为最终目的。
要求学生能综合运用数学基础知识，进行数据的分析和处理、并利用MATLAB软件进行计算机求解。
课程的实践性比较强，强调培养学生的动手动脑能力、开创与创新意识以及解决实际问题的能力。
设计中要求综合运用所学知识，上机解决一些与实际应用结合紧密的、规模较大的问题，通过数据分析、处理等各环节的训练，使学生深刻理解、牢固掌握数据建模的方法，掌握分析、解决实际问题的能力。
三实习内容教师薪金的确定（一），问题的提出某地人事部门为研究中学教师的薪金与他们的资历、性别。
教育程度及培训情况等因素之间的关系，要建立一个数学模型，分析人事策略的合理性，特别是考察女教师是否受到不公正的待遇，以及她们的婚姻状况是否会影响收入。
为此，从当地教师中随机选了3414位进行观察，然后从中保留了90个观察对象，得到了下表给出的相关数据。
尽管这些数据具有一定的代表性，但是仍有统计分析的必要。
现将表中数据的符号介绍如下：Z ~月薪（元）；1X~工作时间（月）；2X=1~男性，2X=0~女性，3X=1~男性或单身女性，3X=0~已婚女性；4X~学历（取值0~6，值越大表示学历越高）；5X=1~受雇于重点中学，5X =0~其它；6X =1~受过培训的毕业生，6X =0~为受过培训的毕业生或受过培训的肄业生；7X =1~已两年一上未从事教学工作，7X =0~其它。
数学建模优秀实验报告
一、实验背景与目的随着科学技术的不断发展，数学建模作为一种解决复杂问题的有力工具，在各个领域都得到了广泛应用。
本实验旨在通过数学建模的方法，解决实际问题，提高学生的数学思维能力和解决实际问题的能力。
二、实验内容与步骤1. 实验内容本实验选取了一道具有代表性的实际问题——某城市交通拥堵问题。
通过对该问题的分析，建立数学模型，并利用MATLAB软件进行求解，为政府部门提供决策依据。
2. 实验步骤（1）问题分析首先，对某城市交通拥堵问题进行分析，了解问题的背景、目标及影响因素。
通过查阅相关资料，得知该城市交通拥堵的主要原因是道路容量不足、交通信号灯配时不当、公共交通发展滞后等因素。
（2）模型假设为简化问题，对实际交通系统进行以下假设：1）道路容量恒定，不考虑道路拓宽、扩建等因素；2）交通信号灯配时固定，不考虑实时调整；3）公共交通系统运行正常，不考虑公交车运行时间波动；4）车辆行驶速度恒定，不考虑车辆速度波动。
（3）模型构建根据以上假设，构建以下数学模型：1）道路容量模型：C = f(t)，其中C为道路容量，t为时间；2）交通流量模型：Q = f(t)，其中Q为交通流量；3）拥堵指数模型：I = f(Q, C)，其中I为拥堵指数。
（4）模型求解利用MATLAB软件，对所构建的数学模型进行求解。
通过编程实现以下功能：1）计算道路容量C与时间t的关系；2）计算交通流量Q与时间t的关系；3）计算拥堵指数I与交通流量Q、道路容量C的关系。
（5）结果分析与解释根据求解结果，分析拥堵指数与时间、交通流量、道路容量之间的关系。
针对不同时间段、不同交通流量和不同道路容量，提出相应的解决方案，为政府部门提供决策依据。
三、实验结果与分析1. 结果展示通过MATLAB软件求解，得到以下结果：（1）道路容量C与时间t的关系曲线；（2）交通流量Q与时间t的关系曲线；（3）拥堵指数I与交通流量Q、道路容量C的关系曲线。
2. 结果分析根据求解结果，可以得出以下结论：（1）在高峰时段，道路容量C与时间t的关系曲线呈现下降趋势，说明道路容量在高峰时段不足；（2）在高峰时段，交通流量Q与时间t的关系曲线呈现上升趋势，说明交通流量在高峰时段较大；（3）在高峰时段，拥堵指数I与交通流量Q、道路容量C的关系曲线呈现上升趋势，说明拥堵指数在高峰时段较大。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验报告（二）院（系）理学院课程名称：数学模型日期2012年3月14 班级A0921 学号01 实验室文理楼209 专业信息与计算科学姓名计算机号12实验名称拟合成绩评定所用软件Matlab 教师签名实验目的或要求1.直观了解拟合基本内容2.掌握用数学软件求解拟合问题实验步骤、心得体会实验题目：1.通过测量10个成年男性的身高与体重，得到如下数据身高x/m1.65 1.70 1.72 1.75 1.78 1.80 1.82 1.79 1.68 1.76体重y/kg61 70 72 74 76 79 80 77 66 73试求身高与体重的关系，并作出关系图实验步骤：（1）以身高为横轴，体重为纵轴，将上述数据对(,)(1,2,,10)i ix y i 在直角坐标系中作出散点图，程序如下：x=[1.65 1.70 1.72 1.75 1.78 1.80 1.82 1.79 1.68 1.76]; %x代表身高y=[61 70 72 74 76 79 80 77 66 73]; %y代表体重plot(x,y,'r*'); %画出数据的散点图运行程序后出来下图：（2）所给的数据点(,)(1,2,,)i i x y i n = 的分布大致呈直线状态，可用直线来拟合这些数据点的曲线。
x=[1.65 1.70 1.72 1.75 1.78 1.80 1.82 1.79 1.68 1.76]; %x 代表身高 y=[61 70 72 74 76 79 80 77 66 73]; %y 代表体重p=polyfit(x,y,1) %对(x,y)做一阶线性拟合运行结果为： p =103.3868 -107.6100即用最小二乘法求得的拟合直线方程为103.3868107.6100y x =-（3）绘制出拟合出的曲线和散点图的对比x=[1.65 1.70 1.72 1.75 1.78 1.80 1.82 1.79 1.68 1.76]; %x 代表身高 y=[6170 72 74 76 79 80 77 66 73]; %y 代表体重p=polyfit(x,y,1) %对(x,y)做一阶线性拟合z=polyval(p,x); %计算拟合后的多项式p 中x 点对应的函数值 plot(x,y,'r*',x,z) %绘制散点图2.已知热敏电阻数据：温度t (ºC) 20.5 32.7 51.0 73.0 95.7 电阻R (Ω) 765 8268739421032求60ºC 时的电阻R ． 实验步骤：(1) 以温度为横轴，电阻为纵轴，将上述数据对(,)(1,2,,5)i i t R i = 在直角坐标系中作出散点图，程序如下：t=[20.5 32.7 51.0 73.0 95.7]; %t 代表温度 R=[765 826 873 942 1032]; %R 代表电阻 plot(t,R,'r*'); %画出数据的散点图 运行程序后出来下图：（2）所给的数据点(,)(1,2,,)i i t R i n = 的分布大致呈直线状态，可用直线来拟合这些数据点的曲线。
t=[20.5 32.7 51.0 73.0 95.7]; %t 代表温度 R=[765 826 873 942 1032]; %R 代表电阻 p=polyfit(t,R,1) %做一阶线性拟合 运行结果为：p =3.3987 702.0968即用最小二乘法求得的拟合直线方程为=+R 3.3987t702.0968（3）绘制出拟合出的曲线和散点图的对比t=[20.5 32.7 51.0 73.0 95.7]; %t代表温度R=[765 826 873 942 1032]; %R代表电阻p=polyfit(t,R,1) %对(t,R)做一阶线性拟合z=polyval(p,t); %计算拟合后的多项式p中t点对应的函数值plot(t,R,'r*',t,z) %绘制散点图R0= polyval(p,60); %求60ºC时的电阻R运行结果：R0 =906.02123.设有数据如表所列，请拟合2~6、10次数时的多项式，并作出数据对比图x 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1y -0.447 1.978 3.280 6.160 7.080 7.340 7.660 9.560 9.480 9.300 11.2 实验步骤：1）将上述数据对(,)(1,2,,11)i i x y i = 在直角坐标系中作出散点图，程序如下：x=[0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1]; %x 的数据y=[-0.447 1.978 3.280 6.160 7.080 7.340 7.660 9.560 9.4809.300 11.200]; %y 的数据 plot(x,y ,'r*'); %画出数据的散点图运行程序后出来下图：2）所给的数据点(,)(1,2,,)i i x y i n = ，用2次的多项式来拟合这些数据的曲线。
程序如下： clearx=[0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1]; %x 的数据y=[-0.447 1.978 3.280 6.160 7.080 7.340 7.660 9.560 9.480 9.300 11.200]; %y 的数据p=polyfit(x,y ,2) %求2次拟合多项式的系数xi=linspace(0,1,100); %建立拟合多项式曲线的横坐标数列yi=polyval(p,xi); %计算拟合后的多项式对应xi 点的函数值yi plot(x,y ,'r*',xi,yi) 运行结果： p =-9.8108 20.1293 -0.0317 即2次拟合多项式为：229.810820.12930.0317p x x =-+-运行程序后出来下图：所给的数据点(,)(1,2,,)i i x y i n ，用3~6次的多项式来拟合这些数据的曲线。
程序如下： clearx=[0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1]; %x 的数据y=[-0.447 1.978 3.280 6.160 7.080 7.340 7.660 9.560 9.480 9.300 11.200]; %y 的数据%-------------------3次拟合曲线--------------------------------------- p3=polyfit(x,y ,3) %求3次拟合多项式的系数xi=linspace(0,1,100); %建立拟合多项式曲线的横坐标数列yi=polyval(p3,xi); %计算拟合后的多项式对应xi 点的函数值yi subplot(2,2,1),plot(x,y ,'r*',xi,yi) title('3次拟合曲线')%-------------------4次拟合曲线--------------------------------------- p4=polyfit(x,y ,4) %求4次拟合多项式的系数xi=linspace(0,1,100); %建立拟合多项式曲线的横坐标数列yi=polyval(p4,xi); %计算拟合后的多项式对应xi 点的函数值yi subplot(2,2,2),plot(x,y ,'r*',xi,yi) title('4次拟合曲线')%-------------------5次拟合曲线--------------------------------------- p5=polyfit(x,y ,5) %求3次拟合多项式的系数 xi=linspace(0,1,100); %建立拟合多项式曲线的横坐标数列yi=polyval(p5,xi); %计算拟合后的多项式对应xi 点的函数值yi subplot(2,2,3),plot(x,y ,'r*',xi,yi) title('5次拟合曲线')%-------------------6次拟合曲线--------------------------------------- p6=polyfit(x,y ,6) %求3次拟合多项式的系数xi=linspace(0,1,100); %建立拟合多项式曲线的横坐标数列yi=polyval(p6,xi); %计算拟合后的多项式对应xi 点的函数值yi subplot(2,2,4),plot(x,y ,'r*',xi,yi) title('6次拟合曲线') 运行结果： p3 =16.0758 -33.9245 29.3246 -0.6104 p4 =22.0571 -28.0385 -6.3531 23.8103 -0.4516 p5 =11.6827 -7.1496 -2.5312 -15.4072 24.9093 -0.465 p6 =1.0e+003 *0.6980 -2.0822 2.3469 -1.2208 0.2678 0.0018 -0.0004 即3次拟合多项式为：6104.03246.299245.330758.1623-+-=x x x y 4次拟合多项式为：4516.08103.233531.60385.280571.22234-+--=x x x x y 5次拟合多项式为：465.09093.244072.155312.21496.76827.112345-+---=x x x x x y6次拟合多项式为：0.0004-0.0018x 0.2678x1.2208x-2.3469xx 2.0822- x 0.698023456+++=y运行程序后出来下图：clearx=[0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1]; %x 的数据y=[-0.447 1.978 3.280 6.160 7.080 7.340 7.660 9.560 9.480 9.300 11.200]; %y 的数据p=polyfit(x,y ,10) %求10次拟合多项式的系数xi=linspace(0,1,100); %建立拟合多项式曲线的横坐标数列yi=polyval(p,xi); %计算拟合后的多项式对应xi 点的函数值yi plot(x,y ,'r*',xi,yi) 运行结果： p =1.0e+006 *Columns 1 through 9-0.4644 2.2965 -4.8773 5.8233 -4.2948 2.0211 -0.6032 0.1090 -0.0106Columns 10 through 110.04 -0.0000即10次拟合多项式为：10987654100.4644 2.2965 4.8773 5.8233 4.2948 2.02110.60320.p xx x x x x x =-+-+-+-+运行程序后出来下图：4.用下列一组数据拟合0.02=+中的参数a,b,kc t a be-()ktt100 200 300 400 500 600 700 800 900 1000 j310c⨯ 4.54 4.99 5.35 5.65 5.90 6.10 6.26 6.39 6.50 6.59 j1）编写M文件curvefun1.mfunction f=curvefun1(x,tdata)f=x(1)+x(2)*exp(-0.02*x(3)*tdata)%其中x(1)=a; x(2)=b；x(3)=k;2）输入命令tdata=100:100:1000cdata=1e-03*[4.54,4.99,5.35,5.65,5.90,6.10,6.26,6.39,6.50,6.59];x0=[0.2,0.05,0.05]; %拟合参数的初始值x=lsqcurvefit ('curvefun1',x0,tdata,cdata) %拟合相应的参数值，放入到x当中f= curvefun1(x,tdata) %求拟合后数据点对应的函数plot(tdata,cdata,'r*',tdata,f) %画出拟合后函数与原始数据对比图3）运算结果为：f =0.0043 0.0051 0.0056 0.0059 0.00610.0062 0.0062 0.0063 0.0063 0.0063x = 0.0063 -0.0034 0.2542得到图像：4）结论：a=0.0063, b=-0.0034, k=0.2542即拟合的表达式为：0.0051=-()0.00630.0034tc t e-备注：本实验报告用于各学科与计算机应用相关课程的实验，务必按时完成。