山野莽夫

  • 归档
    • 随笔
    • 建站资源
    • 分享
    • 代码
  • 地球物理学
    • 专业课
    • 概念解释
  • 计算机
  • 互联网
  • 教程
  • 规划
  • 实验室
    • 珍藏的软件
    • 贴吧云签到
    • A1账号自助申请
山野莽夫
小学生的挣扎的点点滴滴
  1. 首页
  2. 地球物理学
  3. 正文

计算方法---复化科斯特公式

2016年6月14日 4094点热度 1人点赞 0条评论
  1. 实验四 复化科斯特公式求积
    1. 实验要求(含有算法说明和程序说明)

    复化公式就是把一个区间分成若干个小区间,然后在每一个小区间分别应用梯形公式,科斯特公式等计算结果,然后把得到的结果进行累加,从而提高公示的精度。

    复化科斯特公式比复化辛普生公式有更高的精度,把一个大区间分成n份,然后把每个小区间分成四份计算。这是插值型求积公式,为四次插值函数积分的近似值,对于科斯特公式的精度,是可以用于常规运算的。一开始公式没加入精度控制,会无法判断需要计算几次,需要加入精度控制,也就是误差估计,但是误差估计公式中的五阶倒数比较难实现,diff函数试过多次没有成功,所以我就采取了后一项减去前一项的方式配合默认显示精度进行循环次数控制,此法比较简单有效,

    1. 源代码

    function kst(a,b)

    n=2;a=0;b=1;

    for i=1:500

    x= a:(b-a)/(4*n):b;

    m=4*n+1;

    h=(b-a)/n;

    c=(h/90)*(sum(7*f(x(1:4:m-1)))+sum(32*f(x(2:4:m)))+sum(12*f(x(3:4:m)))+sum( 32*f(x(4:4:m)))+ sum(7*f(x(5:4:m))));

    fprintf('s(%d)=%f\n',n,c);

    t(i+1)=c;


    if(i>=1)

    t(i)=c;


    end


    if(i>=2)


    if(abs(t(i)-t(i-1))<=1e-7)

    fprintf('总共计算了%d次',i)


    break;


    end


    end

    n=n*2;

    end

  1. 实验结果


>> kst2

s(2)=0.909605

s(4)=0.909604

s(8)=0.909604

总共计算了3次>>

标签: 复化 科斯特公式 计算方法
最后更新:2018年4月6日

小菜菜

菜鸟

打赏 点赞
< 上一篇
下一篇 >

文章评论

razz evil exclaim smile redface biggrin eek confused idea lol mad twisted rolleyes wink cool arrow neutral cry mrgreen drooling persevering
取消回复

此站点使用Akismet来减少垃圾评论。了解我们如何处理您的评论数据。

标签聚合
宝塔面板 虚拟机 c语言 wordpress onedrive 模板 ppt 地震学程序
最新 热点 随机
最新 热点 随机
Azure Student 微软云 学生订阅 免费12个月用量避坑注意点集合 MP3音频文件格式详细解析 python按固定采样点个数分割wav格式音频 愉快使用谷歌免费人工智能平台colab,训练你的神经网络模型,为你的学术生活添砖加瓦 华为云版轻量应用服务器-云耀云服务器简单体验评测 Cloudflare 免费CDN自定义节点ip之自选cloudflare 高速节点ip工具分享
地震矩张量计算 地震小科普之地震发生观点(成因) 直立长方体重力异常及其一阶偏导数 关于全球付globalcash万事达虚拟信用卡的安全和使用问题的提示和下车方法 腾讯云cos搭建个人ftp服务器 回到家乡-陈宁

COPYRIGHT © 2021 shanyemangfu.com. ALL RIGHTS RESERVED.

Theme Kratos Made By Seaton Jiang

蜀ICP备15031791号-2