MATLAB应用实例分析1 - matlab.docxVIP

.20-0.2y-0.4-0.6-0.8-0.8-0.6-0.4-0.200.20.40.60.8x用plot()函数），得到直观的表示。例：在化学反应中，为研究某化合物的浓度随时间的变化规律，测得一3000);yi=polyval(p,xi);plot(xi,yi)x0=[0.1,0.1];3ff方程组：

.20-0.2y-0.4-0.6-0.8-0.8-0.6-0.4-0.200.20.40.60.8x

用plot()函数），得到直观的表示。例：在化学反应中，为研究某化合物的浓度随时间的变化规律，测得一

3000);yi=polyval(p,xi);plot(xi,yi)x0=[0.1,0.1];3ff

方程组：3x14x1x25x29x26x37x32x37.5412分析：对于线性方程组求解，常用线性

仅举一些运用MATLAB的例子，这些问题在数学建模中时常遇到，希望能帮助同学们在短时间内方便、快捷的使用MATLAB解决数学建模中的问题，并善用这一工具。

常用控制命令：

clc：%清屏；clear：%清变量；save：%保存变量；load：%导入变量

一、利用公式直接进行赋值计算

本金P以每年n次，每次i%的增值率（n与i的乘积为每年增值额的百分比）增加，当增加到r×P

时所花费的时间T为：（利用复利计息公式可得到下式）rPP(10.01i)nTTlnrnln(10.01i)（r2,i0.5,n12）

MATLAB的表达形式及结果如下：

r=2;i=0.5;n=12;%变量赋值

T=log(r)/(n*log(1+0.01*i))

计算结果显示为：

T=11.5813

即所花费的时间为T=11.5813年。

分析：上面的问题是一个利用公式直接进行赋值计算问题，实际中若变量在某个范围变化取很多值时，使用MATLAB，将倍感方便，轻松得到结果，其绘图功能还能将结果轻松的显示出来，变量之间的变化规律将一目了然。

若r在[1,9]变化，i在[0.5,3.5]变化；我们将MATLAB的表达式作如下改动，结果如图1。

r=1:0.5:9;

i=0.5:0.5:3.5;

n=12;

p=1./(n*log(1+0.01*i));

p=1./(n*log(1+0.01*i));T=log(r)*p;

plot(r,T)

xlabel(r)%给x轴加标题

ylabel(T)%给y轴加标题

q=ones(1,length(i));

text(7*q-0.2,[T(14,1:5)+0.5,T(14,6)-0.1,T(14,7)-0.9],num2str(i))

T

40

35

30

25

20

15

10

5

0.5

0.5

1

1.5

2

2.5

3

3.5

0

2

3

4

5

r

图1

6

7

8

9

1

绘制其曲线。MATLAB的表达形式如下：ezplot(1/y-log(y)+log(-1+y)29310.8972-0.6236%

绘制其曲线。MATLAB的表达形式如下：ezplot(1/y-log(y)+log(-1+y)

29310.8972-0.6236%建立系数矩阵%建立常数项矩阵%求方程组的解在数学建模竞赛中，我们

.11,0.16,0.35,0.48,0.61,0.71,0.85];plot(x,y,o)[p

vefit()函数得到的关系为：ε=-3.5810+0.5344σresidresid2分别代表运用

1

从图1中既可以看到T随r的变化规律，而且还能看到i的不同取值对T—r曲线的影响（图中的六条曲线分

别代表i的不同取值）。

二、方程组的求解

8x

求解下面的方程组：3x

1

4x

1

x

2

5x

2

9x

2

6x

3

7x

3

2x

3

7.5

4

12

分析：对于线性方程组求解，常用线性代数的方法，把方程组转化为矩阵进行计算。axbxa1bxa\b

MATLAB的表达形式及结果如下：

a=[816;357;492];b=[7.5;4;12];x=a\b

计算结果显示为：

x=

1.2931

0.8972

-0.6236

%建立系数矩阵

%建立常数项矩阵

%求方程组的解

三、数据拟合与二维绘图

三、数据拟合与二维绘图

在数学建模竞赛中，我们常会遇到这种数据表格问题，如果我们仅凭眼睛观察，很难看到其中的规律，也就更难写出有效的数学表达式从而建立数学模型。因此可以利用MATLA拟合函数，即polyfit()函数，并