数值试验2.doc

上机实习题（三） SARS的传播及预防问题 非典的爆发和蔓延给我国的经济发展和人民生活带来了很大的影响，下表给出了北京市当年4月份到6月份的疫情数据,通过拟合确诊的累积病人曲线,若延后5天采取严格的预防措施,对疫情的传播所生成的影响做出估计. 日期 已确诊病例累积 现有疑似病例 死亡累积 治愈出院累积 4月20日 297 402 18 33 4月30日 1584 1408 75 90 5月1日 1640 1415 82 100 5月10日 1988 1397 116 175 5月20日 2189 1225 150 395 5月30日 2309 706 176 1006 6月1日 2319 739 181 1124 6月10日 2394 351 184 1747 6月20日 2439 3 191 2189 取拟合曲线的拟合函数为如下非线性函数 试确定拟合函数中的参数：a，b，并推测五天后累积病人数量。 程序如下： t=[1 11 12 21 31 41 42 51 61];t0=66; N=[297 1584 1640 1988 2189 2309 2319 2394 2439]; y=1./N;T=1./t A=[ ones(9,1), T ]; d=A\ y ;a=d(1),b=d(2) N0=1/(a+b/t0); x=0:66;yy=1./(a+b./x); plot(x,yy,t,N,o,66,N0,o) 2．海浪高度数据的拟合问题 海洋水文观测站记录了某海域每天24小时海浪潮高度数据（相对于海堤上的零标尺记号）。工作人员每小时做一次记录，下面是从某年12月1日早上零点开始到12月2日晚上23点结束两天内完整的数据记录。能否根据这些数据预测三天以后（即12月5日）的海浪高度数据。 1~8 2.4000 1.2000 -0.1000 -0.5000 -2.5000 -3.0000 -2.7000 -1.6000 9~16 0.2000 2.1000 3.4000 3.6000 2.9000 1.6000 0.2000 -1.2000 17~24 -2.4000 -3.0000 -3.1000 -2.3000 -0.7000 1.3000 2.9000 3.6000 25~32 3.1000 2.0000 0.6000 0.6000 -2.2000 -3.6000 -3.2000 -2.5000 33~40 -0.9000 -1.1000 2.9000 3.9000 3.6000 2.5000 1.0000 -1.5000 41~48 -2.4000 -3.0000 -3.4000 -3.0000 -1.7000 0.2000 2.2000 3.5000 （1）.问题分析： 要预测三天以后的海浪高度，若能够根据当前数据进行函数拟合，则可以利用拟合曲线来预测三天后的海浪高度情况。 1、由于本题没有给出用何种拟合函数，所以首先对数据特性进行分析，选择一个合适的拟合函数。 绘制海浪高度的散点图和连线图，观察数据分布情况 从图上可以看出数据具有较为明显的周期特性，可以考虑选用正弦函数来作为拟合函数 为了确定数据的周期分量，对数据作频谱分析，观察其频率分布情况。 其频谱分布如下图 从幅值响应可以看出，在归一化的频率里，该信号具有一个很明显的峰值，因此可以考虑选用如下函数对海浪高度进行拟合，这里首先只考虑主峰，并观察精度 （2-1） 其中A为信号的直流分量，B为正弦幅值，为角频率可由信号频响给出可记为已知量，为初始相位由于在sin里边，难于计算，考虑如下公式 （2-2） 故将式2-1进行修改，选择拟合函数如下 （2-3） 式中变量x为从计时起点开始的小时数，参数，其中为归一化数字角频率，由图中取得峰值点距0和1的平均值，利用采样频率为1小时每次估算该海浪高度周期T=12.1363小时，和实际情况基本吻合，同时还可增加FFT点数来更精确寻找频率点。 由上面分析可知，本题化为对如下超定方程组求最小二乘解 根据公式求出最小二乘解 A = 0.0362 B = 0.2087 C = 3.3398 该函数对海浪高度拟合如下图 可以看出数据拟合较为成功 r= 3.7009，34点附近有较大偏差，推测可能纪录有误 （2）、MATLAB源程序 说明： 本程序海浪原始数据放在wave. mat文件48*1的矩阵里， 使用本程序时应放在同一文件夹内。 2、程序在各个绘图点停顿，请按按键进入下一页面。 3、command窗口输出了部分重要参数