山野莽夫

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

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

2016年6月14日 3864点热度 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来减少垃圾评论。了解我们如何处理您的评论数据。

标签聚合
wordpress onedrive c语言 ppt 宝塔面板 虚拟机 模板 地震学程序
最新 热点 随机
最新 热点 随机
Azure Student 微软云 学生订阅 免费12个月用量避坑注意点集合 MP3音频文件格式详细解析 python按固定采样点个数分割wav格式音频 愉快使用谷歌免费人工智能平台colab,训练你的神经网络模型,为你的学术生活添砖加瓦 华为云版轻量应用服务器-云耀云服务器简单体验评测 Cloudflare 免费CDN自定义节点ip之自选cloudflare 高速节点ip工具分享
直立长方体重力异常及其一阶偏导数 wordpress更换域名 windows系统下最新版QQ、QQ轻聊版7.9(14314)TIM 3.0防撤回补丁以及使用方法 简单网站源码(应付建站公选课作业) 谷歌云平台gcp300美元试用失败经历。 本站提供免费office365教育版A1子号(5Tonedrive网盘),可以自助申请

COPYRIGHT © 2021 shanyemangfu.com. ALL RIGHTS RESERVED.

Theme Kratos Made By Seaton Jiang

蜀ICP备15031791号-2