matlab实验资料报告材料.doc

word
word
1 / 32
word
MATLAB
数学实验报告
指导教师：
班 级：
小组成员：
时间：201_/_/_
Matlab第二次实验报告
小组成员：
称性可以简便运算。
3题目：实验八,河流流量估计与数据插值
目的：由一些测量数据经过计算处理，解决一些生活实际问题。
问题：实验八上机练习题第三题：瑞士地图如下列图，为了算出他的国土面积，做以下测量，由西向东为x轴，由南向北为y轴，从西边界点到东边界点划分为假如干区域，测出每个分点的南北边界点y1和y2，得到以下数据(mm)。比例尺1:2222，计算瑞士方公里。
word
word
8 / 32
word
测量数据如下：
x=[ 34 48 56 61 91 96 101 104 106 118 142 146 150 157 158]
;
y1=[44 45 47 50 50 38 30 30 34 36 34 41 45 46 43 37 33 28 32 65 55 54 52 50 66 66 68];
y2=[44 59 70 72 93 100 110 110 110 117 118 116 118 118 121 124 121 121 121 122 116 83 81 82 86 85 68];
问题分析：先由题目给定的数据作出瑞士地图的草图，再根据梯形法，使用trapz语句，来估算瑞士国土的面积。
编程：x=[ 34 48 56 61 91 96 101 104 106 118 142 146 150 157 158];
word
word
9 / 32
word
y1=[44 45 47 50 50 38 30 30 34 36 34 41 45 46 43 37 33 28 32 65 55 54 52 50 66 66 68];
y2=[44 59 70 72 93 100 110 110 110 117 118 116 118 118 121 124 121 121 121 122 116 83 81 82 86 85 68];
plot(x,y1,'r.','markersize',15);
plot(x,y2,'r.','markersize',15);
axis([0 160 0 135])
grid;hold on
t=7:158;
u1=spline(x,y1,t);
u2=spline(x,y2,t);
plot(t,u1)
plot(t,u2)
s1=trapz(t,u1);
s2=trapz(t,u2);
s=(s2-s1)*2222*22222/10000000;
fprintf('S=%.0f',s)
word
word
10 / 32
word
结果分析：
上图为由所给数据绘制出的瑞士地图。
上图为运算结果，计算出瑞士的方公里，与准确值41288较为接近。
word
word
11 / 32
word
纪录：使用梯形分割的方法，trapz语句可以方便计算不规如此图形面积，但存在一定误差。
4题目：实验七：圆周率的计算与数值积分
目的：将数值积分最根本的原理应用于matlab之中，解决一些与积分有关的问题。
问题：实验七上机练习题第一题：〔排洪量〕，为了计算最大排洪量，需要计算其断面积，试根据所给数据(m)用梯形法计算其断面积。
问题分析：河床断面可近似分割成假如干曲边梯形，近似处理把它们当做梯形来计算面积可使问题得到简化。
编程：
clc;clear;
x=[0 4 10 12 15 22 28 34 40];
y=[0 1 3 6 8 9 5 3 0];
word
word
12 / 32
word
y1=10-y;
plot(x,y1,'k.','markersize',15);
axis([0 40 0 10]);
grid;
hold on
t=0:40;
u=spline(x,y1,t);
plot(t,u);
s=40*10-trapz(t,u);
fprin