2011第六讲 Matlab数值计算1.ppt

第六讲 Matlab数值计算1 主讲人：祁荣宾 第六讲 主要内容 数据处理与多项式计算 数值微积分 线性方程组求解 非线性方程与最优化问题求解 常微分方程的数值求解 稀疏矩阵 1、数据处理与多项式计算 数据统计处理 1、数据处理与多项式计算 1、数据处理与多项式计算 1、数据处理与多项式计算 1、数据处理与多项式计算 1、数据处理与多项式计算 1、数据处理与多项式计算 1、数据处理与多项式计算 1、数据处理与多项式计算 1、数据处理与多项式计算 1、数据处理与多项式计算 1、数据处理与多项式计算 1、数据处理与多项式计算 1、数据处理与多项式计算 1、数据处理与多项式计算 1、数据处理与多项式计算 1、数据处理与多项式计算 1、数据处理与多项式计算 1、数据处理与多项式计算 1、数据处理与多项式计算 1、数据处理与多项式计算 1、数据处理与多项式计算 1、数据处理与多项式计算 2、数据插值 2、数据插值 2、数据插值 2、数据插值 2、数据插值 2、数据插值 2、数据插值 2、数据插值 2、数据插值 2、数据插值 3、曲线拟合 3、曲线拟合 3、曲线拟合 4、多项式计算 4、多项式计算 4、多项式计算 4、多项式计算 4、多项式计算 4、多项式计算 4、多项式计算 4、多项式计算 4、多项式计算 4、多项式计算 5、数值微积分 5、数值微积分 5、数值微积分 5、数值微积分 5、数值微积分 5、数值微积分 5、数值微积分 5、数值微积分 5、数值微积分 5、数值微积分 T=0:5:65; X=2:5:57; F=[3.2015,2.2560,879.5,1835.9,2968.8,4136.2,5237.9,6152.7,... 6725.3,6848.3,6403.5,6824.7,7328.5,7857.6]; F1=interp1(T,F,X) F2=interp1(T,F,X,nearest) F3=interp1(T,F,X,cubic) F4=interp1(T,F,X,spline) subplot(2,2,1);plot(T,F,X,F1,p);title(linear); subplot(2,2,2);plot(T,F,X,F2,p);title(nearest); subplot(2,2,3);plot(T,F,X,F3,p);title(cubic); subplot(2,2,4);plot(T,F,X,F4,p);title(spline); 例：某检测参数随时间T的采样结果如F，用数据插值法计算X时刻时F的值 液拾钧肉哑到卿苞苔田群逸甜范削赎交唇烯电汕锹靳沃役陵荐妥愧裂吓曾2011第六讲 Matlab数值计算12011第六讲 Matlab数值计算1 吐志肉帝整反诵魂弛颈皿氯庞雨授指霓叁兹职马依逃吁剔涟版埋佑围忻志2011第六讲 Matlab数值计算12011第六讲 Matlab数值计算1 (2) 二维数据插值 Z1=interp2(X,Y,Z,X1,Y1,method) 其中X,Y是两个向量，分别描述两个参数的采样点，Z是与参数采样点对应的函数值， X1,Y1是两个向量或标量，描述欲插值的点。X1与Y1是相同长度的向量或不同方向的向量，Z1是根据相应的插值方法得到的插值结果。 method的取值与一维插值函数相同。X,Y,Z也可以是矩阵形式。 同样，X1,Y1的取值范围不能超出X,Y的给定范围。 佐蹋咏熊食踢碍挞绳猖鹅归搁词缩惊以介羊叭翼念彪焰谚另见挡贴鞋访甄2011第六讲 Matlab数值计算12011第六讲 Matlab数值计算1 x=0:0.1:1;y=0:0.2:2; [X,Y]=meshgrid(x,y); Z=X.^2+Y.^2; interp2(x,y,Z,0.5,0.5) interp2(x,y,Z,[0.5,0.6],0.4) interp2(x,y,Z,[0.5,0.6],[0.4,0.5]) interp2(x,y,Z,[0.5,0.6],[0.4,0.5]) 例：设z=x2+y2，对z函数在[0,1]?[0,2]区域内进行插值。 ans = 0.5100 ans = 0.4100 0.5200 ans = 0.4100 0.6200 ans = 0.4100 0.5200 0.5100 0.6200 掩型哮烽踪德磺颓彰洒略萧蜡掺搀赶诣灼桌挣七朽嘘咏盖奉危患续草医太2011第六讲 Matlab数值计算12011第六讲 Matlab数值计算1 x=0:2.5:10; h=[0:30:60]; T=[95,14,0,0,0;88,48,32,12,6;67,64,54,4