利用matlab软件确定模型中的参数问题.pdfVIP

利用Matlab软件确定 模型中的参数问题 数学与信息科学学院 孔祥庆 数学建模实验项目1 （1） 一、实验名称：利用Matlab软件确定模型中的参数问题 二、实验目的： 掌握利用数学软件确定数学模型中的参数方法 三、实验内容 1、熟悉Matlab软件中数据拟合命令 2、进一步熟悉用初等数学建模的一些方法 3、掌握数学模型中参数的确定方法 四、实验步骤 1、学习Matlab软件中的命令polyfit() 若数据(x ,y ), (x ,y )要拟合是m次多项式函数， 1 1 n n m m -1 y a x a x ... a x a 1 2 m m 1 数学建模实验项目1 （2 ） 则可用命令格式：a=polyfit(x,y,m) 其中X，Y是数据向量，m是多项式的次数。 返回的a是系数构成的向量： a [a 1 a 2 ...a m a m 1 ] 例1 有一只对温度敏感的电阻,已经测得了一组温 度t 和电阻R的数据，试建立R和t关系。 0 t( C) 20.5 32.7 51.0 73.0 95.7 R(W) 765 826 873 942 1032 解：先看一下温度与电阻大 存在什么关系： 输入如下命令： 数学建模实验项目1 （3 ） t=[20.5 32.7 51 73 95.7]; R= [765 826 873 942 1032]; plot(t,R,’r*’) xlabel(‘t’),ylabel(‘R’) aa=polyfit(t,R,1); 运行后得到如下图形： 说明R和t的关系是线性关系，则可用polyfit命令来拟合它的系 数 aa=polyfit(t,R,1); a=aa (1) 运行得：a＝3.3987 b=aa (2) b＝702.0968 数学建模实验项目1 （4 ） 这样得到R与t的关系为： R＝3.3987t+ 702.0968 它们的图形为： 数学建模实验项目1 （5 ） 例2 一册书的成本费y与印刷的册数x有关，统计数据 如下，试用 y=a+b/x去拟合以上数据。 x (千册) 1 2 3 4 5 6 7 8 9 10 i yi （元 10.15 5.52 4.08 2.85 2.11 1.62 1.41 1.30 1.21 1.15 ） Matlab程序： x=1:10; x1=1./x; y=[10.15 5.52 4.08 2.85 2.11 1.62 1.41 1.30 1.21 1.15]; A=polyfit (x1,y,1); 数学建模实验项目1 （6 ） 得结果： 一册书的成本费y 与印刷的册数x 的关系式： 10.2487