《数学实验》报告

实验名称      Matlab基础知识      

20##年4月

一、        【实验目的】

1.       了解Matlab软件，熟悉软件界面和基本操作。

2.       了解Matlab基本运算，掌握Matlab数据类型以及矩阵和数组基本运算。

3.       掌握Matlab程序设计，掌握M文件程序控制语句等。

二、        【实验任务】

P16 第1题：计算表达式e^12+23^3log2(5)/tan21的值

P27 第2题：矩阵A=[123;456;789],B=[468;556;322],计算A*B，A.*B，并比较两者的区别。

P27 第3题：已知A=[52;91]，B=[12;92]，做简单的关系运算A>B,A==B，A<B，

并做逻辑运算（A==B）&（A<B），（A==B）&（A>B）

P34 第1题：用pi/4=1-1/3+1/5-1/7+…公式求pi近似值，直到某一项的绝对值小于10^-6为止。

P34 第3题：编写函数，计算1！+2！+…+50！。

三、        【实验程序】

P16 第1题：

exp(12)+23^3*log2(5)/tan(21)

P27 第2题：

A=[1 2 3;4 5 6;7 8 9];

B=[4 6 8;5 5 6;3 2 2];

C=A*B

D=A.*B

P27 第3题：

A=[5 2;9 1];

B=[1 2;9 2];

C=(A>B)

D=(A==B)

E=(A<B)

F=(A==B)&(A<B)

G=(A==B)&(A>B)

P34 第1题：

m=1;

sum=0;

n=1;

while 1/m>1.0e-6

 sum=sum+n*(1/m);

 m=m+2;

 n=-n;

end

pi=4*sum

P34 第3题：

x=sym('1');

for i=1:50

x=x*sym(i);

end

x

四、        【实验结果】

P16 第1题：

P27 第2题：

P27 第3题：

P34 第1题：

P34 第3题：

五、        【实验总结】

 

第二篇：MATLAB基础性实验报告6

生 验

报

学 实 告

了解插值与拟合的基本原理和方法；掌握用MATLAB计算插值与作最小二乘多项式拟合和曲线拟合的方法；通过范例展现求解实际问题的初步建模过程；

通过动手作实验学习如何用插值与拟合方法解决实际问题，提高探索和解决问题的能力。这对于学生深入理解数学概念，掌握数学的思维方法，熟悉处理大量的工程计算问题的方法具有十分重要的意义。

二、实验仪器、设备或软件： 电脑,MATLAB软件

三、实验内容

1．编写插值方法的函数M文件；

2．用MATLAB中的函数作函数的拟合图形；

3．针对实际问题，试建立数学模型，并求解。

四、实验步骤

1．开启软件平台——MATLAB，开启MATLAB编辑窗口；

2．根据各种数值解法步骤编写M文件；

3．保存文件并运行；

4．观察运行结果(数值或图形)；

5．写出实验报告，并浅谈学习心得体会。

五、实验要求与任务

根据实验内容和步骤，完成以下具体实验，要求写出实验报告（实验目的→问题→数学模型→算法与编程→计算结果→分析、检验和结论→心得体会）。

1．天文学家在19xx年8月的7次观测中，测得地球与金星之间距离（单位：米），并取得常用对数值，与日期的一组历史数据如下表：

由此推断何时金星与地球的距离（米）的对数值为9.93518？

解：输入命令

days=[18 20 22 24 26 28 30];

distancelogs=[9.96177 9.95436 9.94681 9.93910 9.93122 9.92319 9.91499]; t1=interp1(distancelogs,days,9.93518) %线性插值

t2=interp1(distancelogs,days,9.93518,'nearest') %最近邻点插值

t3=interp1(distancelogs,days,9.93518,'spline') %三次样条插值

t4=interp1(distancelogs,days,9.93518,'cubic')

计算结果：

t1 =

24.9949

t2 =

24

t3 =

25.0000

t4 =

%三次插值

25.0000

综上所得，可推断25日金星与地球的距离（米）的对数值为9.93518。

2．在某海域测得一些点(x,y)处的水深z由下表给出，船的吃水深度为5英尺，在矩形区域（75，200）×（-50，150）里的哪些地方船要避免进入。

（1） 输入插值基点数据；

（2） 在矩形区域（75，200）×（-50，150）作二维插值；

（3） 作海底曲面图；

（4） 作出水深小于5的海域范围，即z=5的等高线。

解：

程序：

%输入插值基点数据

x=[129 140 103.5 88 185.5 195 105 157.5 107.5 77 81 162 162 117.5];

y=[7.5 141.5 23 147 22.5 137.5 85.5 -6.5 -81 3 56.5 -66.5 84 -33.5];

z=[4 8 6 8 6 8 8 9 9 8 8 9 4 9];

z=-z;

%在矩形区域（75，200）×（-50，150）作二维插值

cx=75:0.5:200;

cy=-50:0.5:150;

cz=griddata(x,y,z,cx,cy','cubic');

%作海底曲面图

subplot(1,2,1),meshz(cx,cy,cz)

xlabel('x'),ylabel('y'),zlabel('z')

%作出水深小于5的海域范围，即z=5的等高线

subplot(1,2,2),[c,h]=contour(cx,cy,cz);

clabel(c,h,-5)

插值后作出的海底曲面图及等高线图如下：

若船的吃水深度为5英尺，在矩形区域（75，200）×（-50，150）里如上图等高线-5m内的地方船要避免进入。

3．用电压V=10伏的电池给电容器充电，电容器上t时刻的电压为

t

v(t)?V?(V?V0)e?，其中V0是电容器的初始电压，?是充电常数。试由下面一组（t，v）数据确定V0和 ?。

?

解一：

（1）用命令lsqcurvefit。

1)编写M文件curvefun1.m

function f=curvefun1(x,tdata)

f=10-(10-x(1))*exp(-tdata/x(2));

2）主程序xitithree1.m如下

tdata=[0.5 1 2 3 4 5 7 9];

cdata=[6.36 6.48 7.26 8.22 8.66 8.99 9.43 9.63];

x0=[0.4316,1];

x=lsqcurvefit('curvefun1',x0,tdata,cdata)

f=curvefun1(x,tdata)

3）运行主程序，得结果为

x =

5.5577 3.5002

f =

6.1490 6.6616 7.4913 8.1147 8.5832 8.9353 9.3987

9.6604

即拟合得V0=5.5577， ?=3.5002。

（2）用命令lsqnonlin。

1)编写M文件curvefun2.m

function f=curvefun2(x)

tdata=[0.5 1 2 3 4 5 7 9];

cdata=[6.36 6.48 7.26 8.22 8.66 8.99 9.43 9.63]; f=cdata-10+(10-x(1))*exp(-tdata/x(2));

2）主程序xitithree2.m如下

x0=[0.4316,1];

x=lsqnonlin('curvefun2',x0)

f=curvefun2(x)

3）运行主程序，得结果为

x =

5.5577 3.5002

f =

0.2110 -0.1816 -0.2313 0.1053 0.0768 0.0547 0.0313 -0.0304

结果同上，即拟合得V0=5.5577， ?=3.5002。 解二：

（1）对将要拟合的非线性模型v(t)?V?(V?V0)e?，建立M文件volum.m如下：

function yhat=volum(beta,t)

yhat=10-(10-beta(1))*exp(-t./beta(2)); ?t

（2）输入数据：

t=[0.5 1 2 3 4 5 7 9];

y=[6.36 6.48 7.26 8.22 8.66 8.99 9.43 9.63];

beta0=[5 3]';

（3）求回归系数：

[beta,r,J]=nlinfit(t',y','volum',beta0);

beta

得结果：beta =

5.5577

3.5002 即得回归模型为：v(t)?10?(10?5.5577)e

（5） 预测及作图：

[YY,delta]=nlpredci('volum',t',beta,r,J);

plot(t,y,'k+',t,YY,'r') ?t3.5002

七、指导教师评语及成绩：