第6讲(2)数据读入、数据插值.pptVIP

例: 向量t和p表示从1900～1990年的每隔10年的美国人口普查数据： t=1900:10:1990; p=[75.995 91.972 105.711 123.203 131.669 150.697 179.323 203.212 226.505 249.633]; 根据人口普查数据估计1975年的人口，并利用插值估计从1900～2000年每年的人口数。 interp1(t, p, 1975) %估计1975年的人口 再利用插值估计1900～2000年每年的人口数，并利用画图命令得出曲线分布: x=1900:1:2000; y=interp1(t, p, x, spline); %样条插值 plot(t, p, :o, x, y, -r) 例：在一 天24小时内，从零点开始每间隔2小时测得的环境温度数据分别为 12，9，9，1，0，18 ，24，28，27，25，20，18，15，13， 推测中午12点（即13点）时的温度． x=0:2:24; y=[12 9 9 10 18 24 28 27 25 20 18 15 13]; a=13; y1=interp1(x,y,a,spline) 若要得到一天24小时的温度曲线： xi=0:1/3600:24; yi=interp1(x,y,xi, spline); plot(x,y,o ,xi,yi) 例: 利用插值运算对峰值函数peaks插入更多的栅格。 [X, Y]=meshgrid(-3:.25:3); Z=peaks(X, Y); [XI, YI]=meshgrid(-3:.125:3); ZI=interp2(X, Y, Z, XI, YI); mesh(X, Y, Z), hold, mesh(XI, YI, ZI+15) hold off axis([-3 3 -3 3 -5 20]) 例: 给定雇员数据如下： years=1950:10:1990; %工作年份 service=10:10:30; %服役时间，即在职时间 wage=[150.697 199.592 187.625 179.323 195.072 250.287 203.212 179.092 322.767 226.505 153.706 426.730 249.633 120.281 598.243]; %工资 利用二维数据插值计算一个雇员在工作15年后在1975年所获得的工资。 数据读入、数据插值 1. 利用函数读取数据 MATLAB 提供了在主窗口（命令窗口、编辑窗口等）直接输入和利用xlsread、textread等函数直接读取数据等方式。由于大部分的金融数据来自各种机构提供的数据库，这些数据库数据的输出大多会支持Excel和文本文件（txt）的输出格式。因此本节介绍如何利用MATLAB 提供的xlsread 函数来读取*.xls 数据文件。 一 数据的读入 Excel数据文件的读取 格式 num=xlsread(filename) [num,txt,raw]=xlsread(filename) [num,txt,raw]=xlsread(filename,sheet,range) filename：Excel文件名（*.xls）； sheet：工作表名，可用字符串名称或直接用数字表示； range：单元格区域； num：返回xls中的数值型数据； txt：返回xls中的文本（字符串）， raw：返回未处理的形式（字符串）。 【例】 建立2012年6月4日至6月15日上证指数文件名为上证指数历史数据(2011.11).xls的Excel文件，如图所示，并将此文件存放在MATLAB自动搜索文件夹下，比如：D:\用户目录\Documents\MATLAB\上证指数历史数据(2011.11).xls。 上证指数Excel数据表（一） （1）在MATLAB 命令窗口输入： num=xlsread(上证指数历史数据(2011.11)) （2）若在Excel文件“上证指数历史数据(2011.11)”中打开sheet2表 由于sheet2表排在第2张表的位置，所以输入命令： [num,txt,raw]=xlsread(上证指数历史数据(2011.11),2) （3）若只选取数值，可输入命令： num=xlsread(shzhindex,2,B2:C11) （4）如果Excel表中某张表的名字给定，如将“Sheet1”改为名称“全部数据”，则可直接运行命令： [num,txt]=xlsread(上证指数历史数据(2011.11)