版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
样条是绘图员用于描绘光滑曲线的一种机械器件,它是一些易弯曲材料制成的窄条或棒条.在绘制需要通过某点的光滑曲线时,对它在这些点的位置上“压铁”,它就被强制通过或接近图表上确定的描绘点.“样条函数”这个术语意在点出这种函数的图象与机械样条画出的曲线很象.第五节样条插值样条是绘图员用于描绘光滑1一、三次样条插值问题的提法一、三次样条插值问题的提法2数值分析(15)样条插值课件3数值分析(15)样条插值课件4数值分析(15)样条插值课件5二、三转角方程求解法二、三转角方程求解法6数值分析(15)样条插值课件7数值分析(15)样条插值课件8数值分析(15)样条插值课件9数值分析(15)样条插值课件10数值分析(15)样条插值课件11数值分析(15)样条插值课件12数值分析(15)样条插值课件13数值分析(15)样条插值课件14(2)构造三弯矩方程(2)构造三弯矩方程15(2)构造三弯矩方程(2)构造三弯矩方程16数值分析(15)样条插值课件17数值分析(15)样条插值课件18数值分析(15)样条插值课件19即:即:20数值分析(15)样条插值课件21四、样条插值的收敛性四、样条插值的收敛性22数值分析(15)样条插值课件23Chzh1.mChzh1.m24三次样条插值常用命令(1)y0=interp1(x,y,x0,'spline')y0=interp2(x,y,x0,'spline')(2)y0=spline(x,y,x0)(一维)(3)pp=spline(x,y);(结构数组)
y0=fnval(pp,x0);poly(x0,y0),pp.coefs(系数)drp=fnder(pp,n),drv=fnval(drp,a)求导数itp=fnint(pp),itv=fnval(itp,[a,b])*[-1,1]求积分fnplt(pp),fnplt(drp,n),fnplt(itp)绘图
(4)pp=csape(x,y,'边界条件选项');
(5)pp=csapi(x,y);
演示程序sptest.m三次样条插值常用命令25【例】根据连续时间函数的采样数据,利用spline重构该连续函数,并检查重构误差。t=-5:0.5:5;w=exp(-abs(t));N0=length(t);tt=linspace(t(1),t(end),10*N0);ww=spline(t,w,tt);error=max(abs(ww-exp(-abs(tt))))plot(tt,ww,'b');holdonstem(t,w,'filled','r');holdofferror=0.0840【例】根据连续时间函数的采样数据,利用spline重构该连续26【例】对于函数,很容易求得。本例将借此演示样条函数求数值不定积分、导函数的能力。(1)不定积分样条函数、导数样条函数的求取和精度分析。x=(0:0.1:1)*2*pi;y=sin(x);pp=spline(x,y);int_pp=fnint(pp);der_pp=fnder(pp);%xx=(0:0.01:1)*2*pi;err_yy=max(abs(ppval(pp,xx)-sin(xx)))err_int=max(abs(ppval(int_pp,xx)-(1-cos(xx))))err_der=max(abs(ppval(der_pp,xx)-cos(xx)))
err_yy=0.0026err_int=0.0010err_der=0.0253【例】对于函数,很容易求得27(2)不定积分样条函数、导数样条函数的使用%计算y(x)在区间[1,2]上的定积分DefiniteIntegral.bySpline=ppval(int_pp,[1,2])*[-1;1]; DefiniteIntegral.byTheory=(1-cos(2))-(1-cos(1));%计算dy(3)/dxDerivative.bySpline=fnval(der_pp,3);Derivative.byTheory=cos(3);Derivative.byDiference=(sin(3.01)-sin(3))/0.01;DefiniteIntegral,Derivative
DefiniteIntegral=bySpline:0.9563byTheory:0.9564Derivative=bySpline:-0.9895
byTheory:-0.9900byDiference:-0.9907(2)不定积分样条函数、导数样条函数的使用28(3)绘制三个样条函数的图形fnplt(pp,'b-');holdonfnplt(int_pp,'m:'),fnplt(der_pp,'r--');holdofflegend('y(x)','S(x)','dy/dx')(3)绘制三个样条函数的图形29习题五
P208-----14数值实验题五
P210-----4ginput(10)习题五P208-----1430
样条是绘图员用于描绘光滑曲线的一种机械器件,它是一些易弯曲材料制成的窄条或棒条.在绘制需要通过某点的光滑曲线时,对它在这些点的位置上“压铁”,它就被强制通过或接近图表上确定的描绘点.“样条函数”这个术语意在点出这种函数的图象与机械样条画出的曲线很象.第五节样条插值样条是绘图员用于描绘光滑31一、三次样条插值问题的提法一、三次样条插值问题的提法32数值分析(15)样条插值课件33数值分析(15)样条插值课件34数值分析(15)样条插值课件35二、三转角方程求解法二、三转角方程求解法36数值分析(15)样条插值课件37数值分析(15)样条插值课件38数值分析(15)样条插值课件39数值分析(15)样条插值课件40数值分析(15)样条插值课件41数值分析(15)样条插值课件42数值分析(15)样条插值课件43数值分析(15)样条插值课件44(2)构造三弯矩方程(2)构造三弯矩方程45(2)构造三弯矩方程(2)构造三弯矩方程46数值分析(15)样条插值课件47数值分析(15)样条插值课件48数值分析(15)样条插值课件49即:即:50数值分析(15)样条插值课件51四、样条插值的收敛性四、样条插值的收敛性52数值分析(15)样条插值课件53Chzh1.mChzh1.m54三次样条插值常用命令(1)y0=interp1(x,y,x0,'spline')y0=interp2(x,y,x0,'spline')(2)y0=spline(x,y,x0)(一维)(3)pp=spline(x,y);(结构数组)
y0=fnval(pp,x0);poly(x0,y0),pp.coefs(系数)drp=fnder(pp,n),drv=fnval(drp,a)求导数itp=fnint(pp),itv=fnval(itp,[a,b])*[-1,1]求积分fnplt(pp),fnplt(drp,n),fnplt(itp)绘图
(4)pp=csape(x,y,'边界条件选项');
(5)pp=csapi(x,y);
演示程序sptest.m三次样条插值常用命令55【例】根据连续时间函数的采样数据,利用spline重构该连续函数,并检查重构误差。t=-5:0.5:5;w=exp(-abs(t));N0=length(t);tt=linspace(t(1),t(end),10*N0);ww=spline(t,w,tt);error=max(abs(ww-exp(-abs(tt))))plot(tt,ww,'b');holdonstem(t,w,'filled','r');holdofferror=0.0840【例】根据连续时间函数的采样数据,利用spline重构该连续56【例】对于函数,很容易求得。本例将借此演示样条函数求数值不定积分、导函数的能力。(1)不定积分样条函数、导数样条函数的求取和精度分析。x=(0:0.1:1)*2*pi;y=sin(x);pp=spline(x,y);int_pp=fnint(pp);der_pp=fnder(pp);%xx=(0:0.01:1)*2*pi;err_yy=max(abs(ppval(pp,xx)-sin(xx)))err_int=max(abs(ppval(int_pp,xx)-(1-cos(xx))))err_der=max(abs(ppval(der_pp,xx)-cos(xx)))
err_yy=0.0026err_int=0.0010err_der=0.0253【例】对于函数,很容易求得57(2)不定积分样条函数、导数样条函数的使用%计算y(x)在区间[1,2]上的定积分DefiniteIntegral.bySpline=ppval(int_pp,[1,2])*[-1;1]; DefiniteIntegral.byTheory=(1-cos(2))-(1-cos(1));%计算dy(3)/dxDerivative.bySpline=fnval(der_pp,3);Derivative.byTheory=cos(3);Derivative.byDiference=(sin(3.01)-sin(3))/0.01;DefiniteIntegral,Derivative
DefiniteIntegral=bySpline:0.9563byTheory:0.9564Derivative=bySpline:-0.9895
byTheory:-0.9900byDife
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 青岛版小学信息技术三年级上册教案
- 2024年玉米淀粉行业发展现状分析:我国玉米淀粉产量为3644万吨
- 2025届高考地理一轮复习选择题题型滚动练 地球上的大气
- hro服务协议模板
- 2024-2025学年高中英语高一下册牛津译林版教学设计合集
- 2024-2025学年初中信息技术(信息科技)八年级 第12册滇人版(旧版)教学设计合集
- 2024-2025学年小学数学五年级上册冀教版(2024)教学设计合集
- 2024年幼儿园资助工作计划
- 2024年6月5日环境日总结
- 2024-2025学年初中美术九年级上册赣美版(2024)教学设计合集
- 二年级上册《劳动》全册教学设计教案
- 财经素养知识考试题及答案
- 培训学校岗位薪酬绩效管理制度与方案(共三篇)
- 2024年秋季学期新外研版英语三年级上册课件 Unit 6 第2课时(Speed up)
- 国开2024年秋《机电控制工程基础》形考任务1答案
- 第二单元《我们自己》(共七课教案)科学一年级上册教科版
- 国企2024年上半年纪检监察重点工作汇报
- 2024年欧洲生物气溶胶实时在线监测系统市场主要企业市场占有率及排名
- 2024年典型事故案例警示教育手册15例
- 《金刚石、石墨和C60》第一课时名师课件
- 中国天眼完整版本
评论
0/150
提交评论