计算机仿真——数据处理1

计算机仿真——数据处理

给定数据，已知加速度求解速度问题
【问题提出】
从地面发射一枚火箭，在最初100秒内记录其加速度如下，试求火箭在100秒时的速度。
T(s)=[0 10 20 30 40 50 60 70 80 90 100];
A(m/s2)=[30.00 31.63 33.44 35.47 37.73 40.33 43.29 46.69 50.67 54.01 57.23];
【问题分析】
问题分析：加速度为速度的微分，已知微分求积分，类似于面积问题，可利用梯形法来计算。
【问题解决】
clc;clear;
x=[0 10 20 30 40 50 60 70 80 90 100];
y=[30.00 31.63 33.44 35.47 37.73 40.33 43.29 46.69 50.67 54.01 57.23];
plot(x,y,x,y,‘b.’,‘markersize’,50);
axis([0 100 20 60]);
grid
hold on
s=0:5:100;
z=spline(x,y,s);
plot(s,z,‘rp’);
v=trapz(x,y);
fprintf(‘v=%.2f\n’,v)
【函数分析】
spline函数（三次样条函数插值）
用法：yi=spline(x，y，xi)
式中，x，y为插值点的向量，xi为所求点的横坐标值，yi为所求点的纵坐标值，使用目的是通过三次样条函数插值求函数值。求给定位置横坐标的值（依据三次样条函数）
trapz函数（梯形积分法）
采用最简单的方式求积分。z = trapz(x,y) 其中x表示积分区间的离散化向量，y是与x同维数的向量，表示被积函数，z是返回的积分近似值。
plot函数（绘图平面曲线）
plot(x1,y1,选项1,x2,y2,选项2,…,xn,yn,选项n)
选项内容有一下三类：
1线性选项-实线（默认）、：虚线、-.点画线、–双画线
2颜色选项b蓝（默认）、m、g、y、r、b、c、w
3标记选线.、o、X、+、*、s、d、<、>、p、h
markersize为尺寸标记
【运行结果】

v=4168.75m/s
【反思及拓展】
加入plot(x,v)可以绘制火箭发射时的速度特性曲线。