《数学软件》实验指导书.docVIP

《数学建模与数学实验》实验指导书（3+1） 实验一：matlab编程 学时：2学时 实验目的：熟悉matlab编程，掌握用matlab进行函数拟合的方法。 实验内容： 写一个函数rs=f(s)，对传进去的字符串变量s，删除其中的小写字母，然后将原来的大写字母变为小写字母，得到rs返回。例如s=”aBcdE,Fg?”，则rs=”be,f?”。提示：可利用find函数和空矩阵。 f(x)的定义如下： 写一个函数文件f(x)实现该函数，要求参数x可以是向量。 求[100,999]之间能被23整除的数的个数。提示：可利用find和length函数。 一个自然数是素数，且它的各位数字位置经过任意对换之后仍为素数，则成为绝对素数。例如113是绝对素数。试求所有三位的绝对素数。 根据美国人口从1790年到1990年间的人口数据（如下表），确定人口指数增长模型（Logistic模型）中的待定参数，估计出美国2010年的人口，同时画出拟合效果的图形。 美国人口统计数据 年 份 1790 1800 1810 1820 1830 1840 1850 人口(×106) 3.9 5.3 7.2 9.6 12.9 17.1 23.2 年 份 1860 1870 1880 1890 1900 1910 1920 人口(×106) 31.4 38.6 50.2 62.9 76.0 92.0 106.5 年 份 1930 1940 1950 1960 1970 1980 人口(×106) 123.2 131.7 150.7 179.3 204.0 226.5 提示： 选择的一般形式是： if expression {commands} end 如果在表达式(expression)里的所有元素为真，就执行if和end语句之的命令串{commands}. 允许一组命令以固定的和预定的次数重复的形式是： for x=array {commands} end 在for和end语句之间的命令串{commands}按数组（array）中的每一列执行一次. 在每一次迭代中，x被指定为数组的下一列，即在第n次循环中，x=array(：，n) 一些常用函数： sum(a) ：对数组a求和；若a为一矩阵，对a的每一列求和（得到一行向量）；sum(a,2)对每一行求和。 max(a),min(a)的用法和sum一样。 find(a)：找到a中不为0的元素的下标；find(a==2)：找到a中等于2的元素的下标 length(a)：数组a的长度 字符串连接：s=[s1, s2,num2str(1234)] 函数如下定义： function 返回值=函数名（自变量名） 文件名.m必须和函数名一样，如果不一样，函数以文件名为主。 人口模型： 指数增长模型： Logistic模型： 可参考拟合函数：a=lsqcurvefit(example_curvefit_fun,a0,x,y); 实验二：Lingo求解线性规划问题 学时：4学时 实验目的：掌握用Lingo求解线性规划问题的方法，能够阅读Lingo结果报告。 实验内容：（选做两题以上） 1、求解书本上P130的习题1： 某银行经理计划用一笔资金进行有价证券的投资，可供购进的证券以及其信用等级、到期年限、收益如下表1所示，按照规定，市政证券的收益可以免税，其他证券的收益需按50%的税率纳税，此外还有以下限制： 1）政府及代办机构的证券总共至少要购进400万元； 2）所购证券的平均信用等级不超过1.4（信用等级数字越小，信用程序越高）； 3）所购证券的平均到期年限不超过5年。 表 1 证券名称 证券种类 信用等级 到期年限 到期税前收益（%） A 市政 2 9 4.3 B 代办机构 2 15 5.4 C 政府 1 4 5.0 D 政府 1 3 4.4 E 市政 5 2 4.5 （1）若该经理有1000万元资金，应如何投资？ （2）如果能够以2.75%的利率借到不超过100万元资金，该经理应如何操作？ （3）在1000万元资金情况下，若证券A的税前收益增加为4.5%，投资应否改变？若证券C的税前收益减少为4.8%，投资应否改变？ 列出线性规划模型，然后用Lindo求解，根据结果报告得出解决方案。 2、建立模型并求解P130页第3题。（建立线性规划模型的技巧：问什么假设什么，如何雇用即雇用多少全时服务员以12:00-1:00为午餐, 雇用多少全时服务员以1:00-2:00为午餐，雇佣多少从9:0