- 1、本文档共66页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
数值分析第二章 插值.ppt
§5 Cubic Spline Sketch of the Algorithm: Cubic Spline ① 计算 ?j , ? j , gj ; ② 计算 Mj (追赶法等) ; ③ 找到 x 所在区间 ( 即找到相应的 j ) ; ④ 由该区间上的 S[j](x) 算出 f(x) 的近似值。 插值法小结 ? Lagrange : 给出 y0 … yn,选基函数 li(x),其次数为 节点数 –1。 ? Newton ? Ln(x),只是形式不同;节点等距或渐增节点时方便处理。 ? Hermite: 给出 yi 及 yi ’,选 ?i(x) 及 ? i(x) 。 ? Spline:分段低次, 自身光滑, f 的导数只在边界给出。 ? §6 快速傅立叶变换 /* Fast Fourier Transform */ ? 问题的背景 /* background */ 傅立叶变换 —— 函数展开为三角级数 设 f(x) 周期为2?,在 [0, 2? ] 上展开为三角级数 , 其中Cj 为复系数, ,则实际计算时要取级数的前 N 项,并要求在区间的N 个等分点上与f(x)重合。 即:给定[0, 2? ] 上N 个等分点 上的函数值 ,令 满足插值条件 。 N 个未知数 N 个方程 Discrete Fourier Transform Inverse of DFT 总之要进行形如 的计算,其中 已知, §6 Fast Fourier Transform ? Fast Fourier Transform 快速计算 ( j = 0, 1, …, N?1),其中 直接计算需复数乘法 次 N 2 ? 降到 N · logN 由于W 的周期性W qN+s = W s , W k j 实际上只有 这 个不同的值。若 N 为偶数,则W k j只有 个不同值。 1 0 ... - N W W N N / 2 先合并同类项,再做乘法。 I’m gonna need some magic here … §6 Fast Fourier Transform 例:N = 23 = 8, 计算 , j = 0, 1, 2, 3, 4, 5, 6, 7 技巧:将 k , j 先记为二进制数 /* binary numbers */ = ? + ? + ? = = ? + ? + ? = ) ( 2 2 2 ) ( 2 2 2 0 1 2 0 0 1 1 2 2 0 1 2 0 0 1 1 2 2 j j j j j j j k k k k k k k k j W ) 2 2 2 ( ) 2 2 2 ( 0 0 1 1 2 2 0 0 1 1 2 2 ? + ? + ? ? ? + ? + ? = k k k j j j W ) ( ) 2 2 2 ( ) 2 2 2 ( 0 1 2 0 1 0 2 1 3 2 1 2 0 3 1 4 2 2 k k k j k k k j k k k j W + ? + ? + ? + ? + ? + ? = ) ( ) 0 ( ) 0 0 ( 0 1 2 0 0 1 1 0 2 k k k j k k j k j W + + = 2?3次乘法 全部计算需要2?3 ?8次乘法 一般地:取 N = 2p , 每个Cj 用 2p 次乘法,共用 2N log2N 次乘法。 利用 ,还可以进一步化简到 N(p?1)/2 次乘法。 HW: p.138 #1 §4 Newton’s Interpolation ? 等距节点公式 /* Formulae with Equal Spacing */ 向前差分 /* forward difference */ i i i f f f - = ? + 1 i k i k i k i k f f f f 1 1 1 1 ) ( - + - - ?
您可能关注的文档
最近下载
- 2025年浙江安防职业技术学院单招职业技能测试题库(名校卷).docx
- 2024年广州中考英语听说考试——“信息询问” 情景训练课件.pptx
- 四年级下册单元作文《热爱生命》.ppt
- 骆驼奶营养价值课件.pptx
- 统编版高中语文教材中杜甫四首诗《登高》《登岳阳楼》《蜀相》《客至》联读.docx
- 附件3-2:12千伏一二次融合环网柜(箱)及配电自动化终端(DTU)标准化设计方案(2021版).pdf VIP
- 第6课 戊戌变法 课件(共34张PPT).pptx
- 2023年上海杉达学院工程管理专业《管理学》科目期末试卷A(有答案).docx VIP
- 分包考核管理办法和考核表.pdf VIP
- 送教上门康复教案.docx
文档评论(0)