matlab线性多部法求常微分方程数值解

用Adamas内差二步方法,内差三步方法,外差二步方法,外差三步方法这四种方法计算。
在这里插入图片描述
在这里插入图片描述
中k为1和2.

在这里插入图片描述
k为2和3

代码

function chap1_adams_methodu0 = 1;
T = 2;
h = 0.1;
N= T/h;
t = 0:h:T;
solu = exact1(t);f = @f1;
u_inter_2s = adams_inter_2steps(f,u0,t,h,N);
u_extra_2s = adams_extra_2steps(f,u0,t,h,N);
figure(1)
plot(t,u_inter_2s,'*',t,u_extra_2s, 'o',t, solu,'r')
legend('Adams-inter-2s', 'Adams-extra-2s','Exact-soln')u_inter_3s = adams_inter_3steps(f,u0,t,h,N);
u_extra_3s = adams_extra_3steps(f,u0,t,h,N);
figure(2)
plot(t, u_inter_3s,'*', t, u_extra_3s, 'o', t, solu, 'r')
legend('Adams-inter-3s', 'Adams-extra-3s', 'Exact-soln')endfunction u = adams_inter_2steps(f, u0,t,h,N)
u = zeros(N+1,1);
u(1) = u0;
% u(2) = u(1) + h*f(t(1), u(1)) ;
u(2) = exact1(1*h);
eps_in = 1e-6;
K_in = 6;
for n = 2:Nf_nm1 = f(t(n-1),u(n-1));f_n = f(t(n), u(n));s1 = u(n);du = 1;k = 1;
while abs(du)>eps_in & k<K_ins2 = u(n) + h*( 5*f(t(n+1), s1)+ 8*f_n - f_nm1 )/12;du = s2- s1;s1 = s2;k = k + 1;
end
u(n+1) = s2;
end
endfunction u = adams_inter_3steps(f, u0,t,h,N)
u = zeros(N+1,1);
u(1) = u0;
% u(2) = u(1) + h*f(t(1), u(1)) ;
% u(3) = u(2) + h*f(t(2), u(2)) ;
u(2) = exact1(1*h);
u(3) = exact1(2*h);
eps_in = 1e-6;
K_in = 6;
for n = 3:Nf_nm2 = f(t(n-2), u(n-2));f_nm1 = f(t(n-1), u(n-1));f_n = f(t(n), u(n)); s1 = u(n);du = 1;k = 1;
while abs(du)>eps_in & k<K_ins2 = u(n)+ h*(9*f(t(n+1),s1) + 19*f_n - 5*f_nm1 + f_nm2 ) /24;du = s2 - s1;s1 = s2;
k = k + 1;end
u(n+1) = s2;
end
endfunction u = adams_extra_2steps(f,u0,t,h,N)
u = zeros(N+1,1);
u(1) = u0;
% u(2) = u(1) + h*f(t(1), u(1)) ;
u(2) = exact1(h);
for n = 2:N
f_nm1 = f(t(n-1), u(n-1));
f_n = f(t(n), u(n));
u(n+1) = u(n)+ h*( 3*f_n - f_nm1 )/2;
end
endfunction u = adams_extra_3steps(f, u0,t, h,N)
u = zeros(N+1,1);
u(1) = u0;
% u(2) = u(1) + h*f(t(1), u(1)) ;
% u(3) = u(2) + h*f(t(2), u(2));
u(2) = exact1(1*h);
u(3) = exact1(2*h);
eps_in = 1e-6;
for n = 3:N
f_nm2 = f(t(n-2), u(n-2));
f_nm1 = f(t(n-1), u(n-1));
f_n = f(t(n), u(n));
u(n+1) = u(n)+ h*( 23*f_n - 16*f_nm1 + 5*f_nm2 )/12;
end
endfunction f = f1(t, u)
f = -5*u ;
endfunction f = exact1(t)
f = exp(-5*t);
end

结果:
在这里插入图片描述
在这里插入图片描述

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/855794.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

B类IP介绍

一、B类ip给谁用&#xff1a; 城域网用的大多就是B类地址。 二、基本介绍 ip地址范围 - 理论范围 128.0.0.0 ~ 191.255.255.255&#xff1a;10000000 00000000 00000000 00000000 ~ 10111111 11111111 11111111 11111111 - 实际&#xff08;正常情况&#xff09;可用地址范围…

2024-6-19(沉默springboot)

1.spring开启事务支持 事务在逻辑上是一组操作&#xff0c;要么执行&#xff0c;要不都不执行。主要是针对数据库而言的&#xff0c;比如说 MySQL。 业务场景eg: public void savePosts(PostsParam postsParam) {// 保存文章save(posts);// 处理标签insertOrUpdateTag(posts…

【C++】生成随机数

生成 l , r l,~r l, r 范围内的随机数 C o d e Code Code #include <iostream> #include <random> using namespace std;int get_rand(int l, int r) {random_device rd; // 用于获取随机数种子的设备mt19937 gen(rd()); // 以随机数种子初始化随机数生成器unif…

Leetcode 3187. Peaks in Array

Leetcode 3187. Peaks in Array 1. 解题思路2. 代码实现 题目链接&#xff1a;3187. Peaks in Array 1. 解题思路 这一题算是一个套路题&#xff0c;基本就是典型的segment tree的题目。 对于任意一个query&#xff0c;segment tree可以直接获得对应范围内的peak的数目&…

python使用多进程加快处理数据

最近需要用python循环遍历xlsx文件&#xff0c;但总感觉特慢&#xff0c;于是记录一下使用多进程让cpu性能拉满&#xff0c;加快速度。 python使用多进程加快处理数据 1.原始代码2.使用多进程改写代码 1.原始代码 下面展示 原始代码。 这个举例代码的是用来循环遍历我文件夹下…

2023年全国大学生数学建模竞赛A题定日镜场的优化设计(含word论文和源代码资源)

文章目录 一、部分题目二、部分论文基于数学算法研究定日镜场的优化设计一、 问题重述1.1 问题背景1.2 需要解决的问题 二、 模型假设三、 问题分析3.1 问题一的分析3.2 问题二的分析 四、 符号说明五、 模型的建立与求解5.1 问题一模型的建立与求解 三、完整word版论文和源代码…

MYSQL基础查询

示例&#xff1a;user_profile iddevice_idgenderageuniversityprovince12138male21北京大学Beijing23214male复旦大学Shanghai36543female20北京大学Beijing42315female23浙江大学Zhejiang55432male25山东大学Shandong 查询所有列 select * from user_profile&#xff1b;查询…

各种文件类型

一、pkl结尾的文件 .pkl结尾的文件通常是指使用Python的pickle模块序列化后的文件。Pickle是Python中的一个模块&#xff0c;用于将Python对象转换为字节流&#xff0c;以便保存到文件中或者通过网络传输。这个过程被称为“序列化”&#xff0c;而将字节流重新转换回Python对象…

React最新面试攻略

一.简单聊聊React中的组件生命周期及其作用解&#xff1a; React组件的生命周期可以分为三个阶段&#xff1a;挂载阶段、更新阶段和卸载阶段。 挂载阶段包括constructor、componentWillMount、render和componentDidMount等生命周期方法&#xff0c;用于初始化组件和挂载后的操…

项目训练营第二天

项目训练营第二天 用户登录逻辑 1、账户名不少于4位 2、密码不少于8位 3、数据库表中能够查询到账户、密码 4、密码查询时用同样加密脱敏处理手段处理后再和数据库中取出字段进行对比&#xff0c;如果账户名未查询到&#xff0c;直接返回null 5、后端设置相应的脱敏后用户的s…

就因为没在大屏项目加全屏按钮,早上在地铁挨了领导一顿骂

“嗯嗯”&#xff0c;“嗯嗯”&#xff0c;“那产品也没说加呀”&#xff0c;“按F11不行吗&#xff1f;”&#xff0c;“嗯嗯”&#xff0c;“好的”。 早上在4号线上&#xff0c;我正坐在地铁里&#xff0c;边上站着的妹子&#xff0c;我看他背着双肩包&#xff0c;打着电话…

零停机部署——特征开关(Feature Toggles)的应用

引言 在现代软件开发和部署中&#xff0c;零停机部署技术是实现高可用性和无缝用户体验的关键。本文将讨论功能开发开关(Feature Toggles)的类型并分析它们的优缺点&#xff0c;同时提供相关的例子和演示。 PS: https://github.com/WeiXiao-Hyy/blog整理了后端开发的知识网络…

即时到账支付系统源码第四方支付平台源码(支付宝/QQ钱包/微信二维码收款+附配套软件)

即时到账支付系统源码第四方支付平台源码价值10万&#xff0c;支付宝/QQ钱包/微信二维码收款&#xff0c;附配套软件 开发语言&#xff1a;phpmysql 这个是一个可以跟码支付一样用自己的二维码收款的网站 还可以作为即时到账 代收款 或者易支付使用后台配置好就行&#xff…

python 微信群发工具

效果如下 python 微信消息群发工具 1、环境安装依据python手机app自动化测试环境搭建-CSDN博客 2、元素获取 启动uiautomatorviewer.bat工具,如下图 点击获取手机界面,如下图 3、根据获取的界面元素id或class名称编写代码 python简易代码如下 from appium import webdrive…

2024.06.19 刷题日记

41. 缺失的第一个正数 这个题目的通过率很低&#xff0c;是一道难题&#xff0c;类似于脑筋急转弯&#xff0c;确实不好想。实际上&#xff0c;对于一个长度为 N 的数组&#xff0c;其中没有出现的最小正整数只能在 [1,N1] 中。 这个结论并不好想&#xff0c;举个例子&#x…

移动硬盘在苹果电脑上无法识别的诊断与恢复策略

一、问题描述 在数字时代&#xff0c;移动硬盘已成为我们存储和传输数据的重要工具。然而&#xff0c;当我们将移动硬盘插入苹果电脑时&#xff0c;有时会遇到无法识别的情况&#xff0c;这让我们感到十分困扰。本文将详细探讨移动硬盘插苹果电脑后读不出来的现象&#xff0c;…

浸没式液冷服务器的换热效率及节能潜力分析

服务器浸没式液冷的换热效率及节能潜力 摘要&#xff1a;我们针对服务器浸没式液冷实验台进行了深入测试&#xff0c;探究了不同室外温度和服务器发热功率对系统制冷PUE的影响。实验数据显示&#xff0c;该系统的制冷PUE值介于1.05至1.28之间&#xff0c;高效节能特点显著。 在…

坚持刷题|反转链表

文章目录 题目思考实现1. 迭代方式实现链表翻转2. 递归方式实现链表翻转 Hello&#xff0c;大家好&#xff0c;我是阿月。坚持刷题&#xff0c;老年痴呆追不上我&#xff0c;今天继续链表&#xff1a;反转链表 题目 LCR 024. 反转链表 思考 翻转链表是一个常见的算法问题&a…

为什么带上符号位计算,仍然可以算出正负数对应的补码

文章目录 &#x1f680;前言&#x1f680; 为什么要引入原反补码✈️ 利用加法表示减法✈️ 关于数字在计算机中的存储&#x1f680; 数学解释 &#x1f680; 为什么带上符号位计算&#xff0c;仍然可以算出正负数对应的补码 &#x1f680;前言 为应付期末&#xff0c;速成数电…

Linux系统之配置Nginx负载均衡

Linux系统之配置Nginx负载均衡 一、Nginx介绍1.1 Nginx简介1.2 Nginx反向代理1.3 相关概念二、本次实践介绍2.1 本次实践简介2.2 本次实践环境规划三、部署两台web服务器3.1 运行两个Docker容器3.2 编辑测试文件四、配置负载均衡4.1 安装nginx软件4.2 编辑nginx配置文件4.3 启动…