数学建模竞赛用盗版matlab

2、已知速度曲线v(t) 上的四个数据点下表所示

t=[0.15,0.16,0.17,0.18];

v=[3.5,1.5,2.5,2.8];

x=0.15:0.001:0.18

y=i n t e r p1(t,v,x,'s p l i n e')

S=t r a p z(x,y)

p=p o l y f i t(x,y,5);

d p=p o l y d

e r(p);

d p x=p o l y v a l(d p,0.18)

运行结果

S=

0.0687

Dpx=-

3、计算图片文件tu.bmp 给出的两个圆A，B 的圆心，和两个圆的两条外公切线和两条内公切线的切点的坐标。

（1）计算A 圆的圆心坐标

I=imread('tu.bmp');

[m,n]=size(I)

BW=im2bw(I)

BW(:,200:512)=1;

figure, imshow(BW)

ed=edge(BW);

[y,x]=find(ed);

x0=mean(x), y0=mean(y)

r1=max(x)-min(x),r2=max(y)-min(y)

r=(r1+r2)/4

x0 =109.7516

y0 =86.7495

r1 =162

r2 =158

r =80

（2）B圆的圆心坐标和半径

I=imread('tu

1 蒙特卡洛 2 数据处理 2.1 数据拟合

多项式曲线拟合：p = polyfit(x, y, m);

m为拟合多项式的次数。从高次到低次将系数返回到p中。

求多项式在x0处的值y0：y0 = polyval(p, x0);

非线性曲线最小二乘法拟合：[x, resnorm] = lsqcurvefit(fun, x0, xdata, ydata); fun为给定的函数，x0为初值。返回fun中的系数向量x和残差的平方和resnorm。

非线性曲线最小二乘法拟合：[x, resnorm] = lsqnonlin(fun, x0, LB, UB, option, para1, para2, …);

fun为给定的函数，x0为初值，LB为系数下限，UB为系数上限，para为函数fun所需要的参数（依序）。返回fun中的系数向量x和残差的平方和resnorm。

设置选项option：optimset(‘MaxIter’, 300, ‘TolX’, 1e-10, ‘TolFun’, 1e-10);

MaxIter为最大允许的迭代次数，TolX为x的终止公差，TolFun为函数值的终止公差。

非线形回归：[beta, r, j] = nli

MATLAB数学建模习题1

一、单项选择题（将选择答案写在答题纸上，每小题2分共20分）

1．在MATLAB命令窗口中键入命令，Vname=prod(7:9)/prod(1:3)，可计算组合数

如果省略了变量名Vname，MATLAB表现计算结果将用下面的哪一变量名做缺省变量名

A）ans； B）pi； C）NaN； D）eps

2．宝石切割问题中，石料左右长度、前后长度、上下高度分别为a1、a2、a3，即a1×a2×a3(cm3)，而精品尺寸为b1×b2×b3(cm3)。操作时，同向切割连续两次再旋转刀具。某一切割方案的切割面积依次为：2a1a2? 2a1b3 ? 2b2b3，则这一切割方案为

A）左右?前后?上下； B）上下?前后?左右； C）前后?上下?左右； D）前后? 左右?上下 3．机场指挥塔位置：北纬30度35.343分，东经104度2.441分，在MATLAB中用变量B=[30 35.343]表达纬度，L=[104 2.441]表达经度。将数据转化为以度为单位的实数，下面正确的语句是

A） P=B(1)+B(2)/60，Q=L(1)+L(2)； B) P =

数统 应数 20121323003 王妍 数统 应数 20121323022 胡可旺

1 生成5阶矩阵，使其元素满足均值为1，方差为4的正态分布； 代码：y=1+sqrt(4)*randn(5) 结果:

2，生成一个20行5列矩阵A，其元素分别以概率0.7、0.1、0.2取自0、1、4这三个整数， 然后将“删除A的全零行”之后的结果赋给矩阵B； for i = 1:20 for j = 1:5

p = rand(); if p<=0.7

A(i,j) = 0;

elseif p>0.7&&p<=0.8 A(i,j) = 1; else

A(i,j) = 4; end end end B = []; for i = 1:20 flag = true; for j = 1:5

if A(i,j)~=0 flag=false; break; end end if ~flag

B = [B;A(i,:)]; end end B

3, （

2016年周口师范学院第四届数学建模竞赛题目

（请先阅读“数学建模竞赛论文格式规范”）

A题 养老金发放问题

伴随着我国老龄化程度不断提高，老年人在我国人口中所占的比例越来越高，我国已经步入了老龄化社会。老龄化人口的增多对我国现行的社会保障制度带来了严峻的挑战，为应对这种挑战，备受关注的延迟退休年龄政策的步伐渐行渐明朗。党的十八届三中全会已敲定：实行渐进式延长退休年龄，直至达到65岁退休。在新的制度下大家对退休后的养老金能否满足生活需求存在疑问。

问题1：收集关于养老金发放的相关政策建立不同人群养老金发放数量的数学模型；部分企事业单位为提高岗位竞争力，为职工建立职业年金（职业年金是退休金的补充），建立含有企业年金的养老金发放数学模型；

问题2：收集周口市养老金发放相关资料，使用问题1建立的数学模型计算不同人群的养老金发放数量。

参考资料：

【1】周口市统计局网：http://www.zktj.gov.cn/； 【2】周口统计年鉴等。

2016年周口师范学院第四届数学建模竞赛题目

（请先阅读“数学建模竞赛论文格式规范”）

B题 周口市区住宅价格问题

随着周口市城市化进程的推进，城市住宅投资规模不断扩大，住宅销售面积逐年增加，住宅价格下

2016年周口师范学院第四届数学建模竞赛题目

（请先阅读“数学建模竞赛论文格式规范”）

A题 养老金发放问题

伴随着我国老龄化程度不断提高，老年人在我国人口中所占的比例越来越高，我国已经步入了老龄化社会。老龄化人口的增多对我国现行的社会保障制度带来了严峻的挑战，为应对这种挑战，备受关注的延迟退休年龄政策的步伐渐行渐明朗。党的十八届三中全会已敲定：实行渐进式延长退休年龄，直至达到65岁退休。在新的制度下大家对退休后的养老金能否满足生活需求存在疑问。

问题1：收集关于养老金发放的相关政策建立不同人群养老金发放数量的数学模型；部分企事业单位为提高岗位竞争力，为职工建立职业年金（职业年金是退休金的补充），建立含有企业年金的养老金发放数学模型；

问题2：收集周口市养老金发放相关资料，使用问题1建立的数学模型计算不同人群的养老金发放数量。

参考资料：

【1】周口市统计局网：http://www.zktj.gov.cn/； 【2】周口统计年鉴等。

2016年周口师范学院第四届数学建模竞赛题目

（请先阅读“数学建模竞赛论文格式规范”）

B题 周口市区住宅价格问题

随着周口市城市化进程的推进，城市住宅投资规模不断扩大，住宅销售面积逐年增加，住宅价格下

matlab 数学建模

1、设A、B两方案的净现金流量（单位：万元）如下表所示： （1）设折现率为10%，计算两个方案的净现值；

（2）计算两个方案的内部收益率。

2、某厂生产的一种电器的销售量y与竞争对手的价格x1和本厂的价格x2有关。下表是该厂在二十个城市的销售记录。

matlab 数学建模

（1） 根据这些数据建立本厂的需求函数模型；

设方程为 y=a + b*x1 + c*x2 Matlab 程序设计如下：

x1=[ 100 110 105 115 118 145 112 165 170 190 156 268 150 200 230 160 270 250 280 286];

x2=[120 140 138 130 136 148 110 150 165 160 130 145 166 145 150 140 180 145 200 220]; y=[102 100 110 115 105 98 95 93 90 89 77 69 92 60 85 82 65 69];

X=[ones(20,1) x1’ x2’];

[b,bint,r,rint,stats]

matlab入门知识

数学建模与数学实验MATLAB入门 入门

后勤工程学院数学教研室

matlab入门知识

●

MATLAB作为线性系统的一种分析和仿真

工具，是理工科大学生应该掌握的技术工具，它 作为一种编程语言和可视化工具，可解决工程 工程、 工程 科学计算和数学学科 科学计算 数学学科中许多问题。 数学学科 向量、数组 矩阵的基 ● MATLAB建立在向量 数组 矩阵 向量 数组和矩阵 础上，使用方便，人机界面直观，输出结果可 视化。 ● ● 矩阵是MATLAB的核心 矩阵 MATLAB的进入与运行方式(两种)

matlab入门知识

MATLAB入门 入门一、变 量 与 函 数 二、数 三、 矩 组 阵

四、 MATLAB编程 编程 五、 实 验 作 业

matlab入门知识

一、变 量 与 函 数1、变量 、 MATLAB中变量的命名规则 命名规则是： 命名规则 (1)变量名必须是不含空格的单个词； (2)变量名区分大小写； (3)变量名最多不超过19个字符； (4)变量名必须以字母打头，之后可以是 任意字母、数字或下划线，变量名中 不允许使用标点符号.

matlab入门知识

特殊变量表特殊变量 ans pi eps flops inf NaN i,j na

WORD整理版

§12.5 灰色预测

我们通常所说的系统是指：由客观世界中相同或相似的事物和因素按一定的秩序相互关联、相互制约而构成的一个整体．例如：工程技术系统、社会系统、经济系统等．如果一个系统中具有充足的信息量，其发展变化的规律明显、定量描述方便、结构与参数具体，则这种系统通常称为白色系统．如果一个系统的内部特征全部是未知的，则称此系统为黑色系统．如果系统内部信息和特征是部分已知的，另一部分是未知的，这种系统称为灰色系统．例如：社会系统、农业系统、经济系统、气象系统、生物系统等．对于这类系统，内部因素难以辨识，相互之间的关系较为隐蔽，人们难以准确了解这类系统的行为特征．因此，对于这类问题进行定量描述，即建立模型难度较大．区别白色系统与灰色系统的重要标志是系统内各因素之间是否具有确定的关系．

灰色系统分析方法主要是根据具体灰色系统的行为特征数据，充分利用数量不多的数据和信息寻求相关因素自身与各因素之间的数学关系，建立相应的数学模型．目前，灰色系统理论在实际中已得到了广泛的应用，例如：在工程技术、经济管理、气象预报以及政治、社会、工业、农业等领域都取得了一定的应用成果

《Matlab与数学建模》综合练习

1．按顺序进行如下的操作：

（1）产生一个5阶魔术方阵A；并计算A'与A-1(即inv(A))； >> A=magic(5) A =

17 24 1 23 5 7 4 6 13 10 12 19 11 18 25

>> A'

ans =

17 23 4 24 5 6 1 7 13 8 14 20 15 16 22 >> inv(A)

ans =

-0.0049 0.0512 0.0431 -0.0373 -0.0303 0.0031 0.0047 -0.0065 0.0028 0.0050

（2）求A的特征值；>> P=poly(A) P =

1.0e+006 *

0.0000 -0.0001

>> roots(P)

8 15 14 16 20 22 21 3 2