幂法 Euler法

news2024/10/6 10:39:49

一、实验原理

根据幂法,Euler法的相关知识和算法编程完成本实验

二、实验内容

     见第7PPT:用规范化幂法计算下面矩阵的按模最大特征值及对应的特征向量

   

见第8PPT:分别用简单欧拉法和改进欧拉法对如下初值问题作验证性计算

    

   

三、实验过程(可以文字说明+运行结果截图)

1)规范化幂法

实验说明:

数据说明

1.矩阵A:

一个 3×3 的系数矩阵,表示为 array[N][N]。

2.初始向量

一个长度为 3 的初始猜测向量,表示为 a0[N]。

3.迭代次数:

count:迭代次数。

收敛条件:

ϵ=0.000001:绝对误差阈值,当特征值的相对变化小于该值时停止迭代。

代码说明

1.GetMax函数:

获取数组中的最大值。

2. matrixMulti函数:

对矩阵 A 和初始向量

进行矩阵相乘,得到迭代后的向量。

3. main函数:

定义矩阵 A 和初始向量

在迭代过程中,对初始向量

进行规范化,然后进行矩阵相乘,得到迭代后的向量。

计算特征值的相对变化,若小于收敛条件 ϵ,则停止迭代。

输出迭代次数和最大特征值。

算法说明

幂法:用于估计矩阵的最大特征值和对应的特征向量。

在每次迭代中,先对当前向量进行规范化,然后进行矩阵与向量的乘积。

通过计算每次迭代后向量中的最大元素,来估计最大特征值的大小。

当最大特征值的相对变化小于给定的绝对误差阈值时,停止迭代。

运行截图:

2)简单欧拉法和改进欧拉法

实验说明:

数据说明

1.初始点

2.步长 h

3.迭代次数 n。

代码说明

1. Euler类:

实现了简单欧拉法。

包括成员函数 fun 用于定义微分方程。

成员函数 Eulerfun 执行简单欧拉法迭代计算。

2.EulerPro类继承自 Euler:

实现了改进的欧拉法。

在 EulerProfun 函数中使用了改进的公式进行迭代计算。

3. main函数:

创建了简单欧拉法和改进欧拉法的对象,并调用相应的函数进行计算。

打印了计算结果。

算法说明

1.简单欧拉法:

使用简单欧拉法逐步逼近微分方程的数值解。

在每一步中,通过计算当前点的斜率,然后使用步长 h进行线性近似,从而得到下一个点的近似值。

2.改进欧拉法:

改进欧拉法在简单欧拉法的基础上,通过使用改进的公式来提高精度。

在每一步中,不仅仅使用当前点的斜率,还使用了下一个点的斜率,然后取两者平均值作为斜率,从而得到下一个点的近似值。

运行截图:

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

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

相关文章

rga_mm: RGA_MMU unsupported Memory larger than 4G!解决

目录 报错完整log如下:解决方案:报错完整log如下: [ 3668.824164] rga_mm: RGA_MMU unsupported Memory larger than 4G! [ 3668.824305] rga_mm: scheduler core[4] unsupported mm_flag[0x0]! [ 3668.824320] rga_mm: rga_mm_map_buffer map dma_buf err

音视频入门基础:H.264专题(4)——NALU Header:forbidden_zero_bit、nal_ref_idc、nal_unit_type简介

音视频入门基础:H.264专题系列文章: 音视频入门基础:H.264专题(1)——H.264官方文档下载 音视频入门基础:H.264专题(2)——使用FFmpeg命令生成H.264裸流文件 音视频入门基础&…

昇思25天学习打卡营第1天|新手上路

这里写自定义目录标题 打卡昇思MindSpore扫盲快速入门 打卡 昇思MindSpore扫盲 第一节基本是一个mindspore的科普扫盲。大概介绍一通mindspore的一些架构,feature,以及其对比于其他同类框架的优势。简单扫读了一遍大概有点印象直接跳过。 快速入门 这…

Spring IOC架构设计解析

一、spring famework 模块解析 二、SpringIOC容器:Core Container: Beans、Core、Context、Expression 三、BeanFactory (该接口定义了容器最基本的功能方法)和 ApplicationContext(对容器提供了更丰富的功能扩展,大部分应用中所使用的基于它实现&…

【机器学习-10】 | Scikit-Learn工具包进阶指南:Scikit-Learn工具包之支持向量机模块研究

🎩 欢迎来到技术探索的奇幻世界👨‍💻 📜 个人主页:一伦明悦-CSDN博客 ✍🏻 作者简介: C软件开发、Python机器学习爱好者 🗣️ 互动与支持:💬评论 &…

高考英语3500词

DAY1 DAY2 DAY3 DAY4 DAY5 DAY6 DAY7 DAY8 DAY9 DAY10 DAY11 DAY12 DAY13 DAY14 DAY15 DAY16 DAY17 DAY18 DAY19 DAY20 DAY21 DAY22 DAY23 DAY24 DAY25 DAY26 DAY27 DAY28 DAY29 DAY30 DAY31 DAY32 DAY33 DAY34 DAY35 DAY36 DAY37 DAY38 DAY39 DAY40

更进一步|瑞数信息位居2023年中国私有云WAF市场份额Top2!

更进一步|瑞数信息位居2023年中国私有云WAF市场份额Top2! 近日,IDC《中国云Web应用防火墙市场份额,2023:WAAP成为共同演进方向》报告正式发布。报告聚焦2023年中国云Web应用防火墙市场的规模、厂商份额以及技术发展变…

Python光束三维二维标量场和算法

🎯要点 🎯矢量计算和绘图:🖊二维,三维 | 🖊转换矢量图 | 🖊矢量和矩阵计算 | 🖊矢量空间变换和计算 | 🖊解矢量线性方程 | 🎯数学和物理矢量计算&#xff1a…

ppdetection-2.7

1, pip install -r requirement.txt 2、 ValueError: paddle.load can not parse the file:C:\Users\HX/.cache/paddle/weights\deepsort_pcb_pyramid_r101.pdparams. 重新下载模型 3、 declarative() got an unexpected keyword argument看一下参数是否写对,比如…

树莓派选型深度指南:架构、性能与应用场景全解析

树莓派作为全球最受欢迎的单板计算机,凭借其低廉的价格、强大的性能和丰富的扩展性,吸引了无数电子爱好者、开发者和教育工作者的青睐。然而,面对琳琅满目的树莓派型号,新手往往不知从何下手。 本文将深入对比几款市面上最热门的…

卸载 ubuntu-wsl2-systemd-script,使用 WSLg 图形用户界面

目录 全新安装 - 以前没有安装 WSL现有 WSL 安装卸载 ubuntu-wsl2-systemd-script使用 Linux GUI参考链接在 Windows 上使用 Linux 开发环境,最好的做法是使用 WSL2。在 WSL 和早期的 WSL2 版本中,并不支持图形用户界面。因此如果想要使用 GUI 程序,需要自行解决。具体方法可…

Notepad++爱国版重磅更新:集成18大高人气插件

Notepad最新版的言论与导向存在价值观问题,经重新编译修正为“统一祖国,振兴中华”版本,传导正确的价值观和爱国情怀。 Notepad功能可以通过各种插件来增强,下面是18大高人气插件 AutoCodepage 功能:在加载、重命名或更…

秋招突击——第七弹——Redis快速入门

文章目录 引言Redis是什么 正文对象String字符串面试重点 List面试考点 压缩列表ZipList面试题 Set面试题讲解 Hash面试重点 HASHTABLE底层面试考点 跳表面试重点 ZSET有序链表面试重点 总结 引言 在项目和redis之间,我犹豫了一下,觉得还是了解学习一下…

ARM服务器虚拟化手机,云手机推流应用案例

大家都知道,ARM 服务器虚拟化手机和云手机推流技术可算是热门话题,不止是企业,个人卖家也会通过云手机推流来获得更多的客源,实现经济自由,但是针对云手机的推流,很多人还是不知道有哪些应用场景~我们可以展…

永洪bi知识点

1、下拉过滤组件和下拉参数组件的区别 下拉过滤组件只能对跟他绑有相同数据集的组件进行过滤 而下拉参数组件是当你设置了筛选条件以后,那么所有的组件,不管你绑定了什么样子的数据集,都能起作用,前提是你这个组件是需要去绑定参…

LED热管理

LED照明系统的热管理 本文提供了用于LED灯具的热管理系统。 包含LED轨道灯具包括照明组件、安装到照明组件上并具有多个孔的夹具壳体,以及将夹具壳体固定到轨道上的安装结构。 照明组件包括具有多个翅片的散热器、安装在所述散热器上的反射器、支撑在所述散热器上…

RabbitMQ的WorkQueues模型

WorkQueues模型 Work queues,任务模型。简单来说就是让多个消费者绑定到一个队列,共同消费队列中的消息。 当消息处理比较耗时的时候,可能生产消息的速度会远远大于消息的消费速度。长此以往,消息就会堆积越来越多,…

80、443端口不能开放也能为IP地址申请SSL证书!

IP地址证书作为一种特定的证书,不同于传统的域名验证证书,IP地址证书是通过验证IP地址来确保安全连接。在证书申请过程中,往往要求短暂开放80或者443端口,如果不能开放,IP地址证书则不能签发。 JoySSL提供的IP地址证书…

如何采集拼多多的商品或店铺数据

怎么使用简数采集器批量采集拼多多的商品或店铺相关信息呢? 简数采集器暂时不支持采集拼多多的商品或店铺相关数据,只能采集页面公开显示的信息,谢谢。 简数采集器采集网站文章资讯等数据特别简单高效:只需输入网站网址&#xf…

文心一言使用笔记

目录 让文心一言提炼已有的内容,模仿给出的案例写一段宣传稿方法例子 发现写出的内容有瑕疵,如何微调?比如文心一言介绍的领导不全如何让文心一言检查语法和表达问题? 如何让文心一言将每个片段用一两句话总结?为了防止…