多符号表达式的共同子表达式提取教程

生成的符号表达式,可能会存在过于冗长的问题,且多个符号表达式中,有可能存在相同的计算部分,如果不进行处理,计算过程中会导致某些算式计算多次,从而影响计算效率。

那么多个符号表达式生成函数时,如何将多个符号表达式中相同算子提取出来,以便实现计算效率的提升?

一种方案是:共同子表达式提取(Common Subexpression Elimination,CSE) ,即找到多个表达式中相同的子表达式,并将其提取出来。这可以通过建立一个表达式树或图的数据结构,然后检测重复的子树来实现。一旦发现相同的子表达式,你可以将其计算结果保存起来,而不是多次计算。这通常需要对表达式进行遍历和分析。

1. Python实现(推荐)

在Python中,可以利用sympy中的cse实现,案例如下:

from sympy import symbols, cse, expand# 定义符号变量
x, y, z = symbols('x y z')# 创建一些符号表达式
expr1 = x**2 + y**3 + x**z
expr2 = x**2 - y**2 + 5*x**z
expr1 = expand(expr1)
expr2 = expand(expr2)combined_expr = [expr1, expr2]# 进行共同子表达式提取
common_subexpr, simplified_expr = cse(combined_expr)# 打印提取出的共同子表达式
print("共同子表达式:", common_subexpr)# 打印简化后的表达式
print("简化后的表达式:", simplified_expr)"""
Output:
>> 共同子表达式: [(x0, x**2), (x1, x**z)]
>> 简化后的表达式: [x0 + x1 + y**3, x0 + 5*x1 - y**2]
"""

2. MATLAB实现

  • 法1:利用subexpr,优点:可读性强;缺点:不够简洁,原因是该函数一次只能提取1个变量,为此,可以考虑进行多次迭代提取,可参考MATLAB加速计算:符号表达式——子表达式法,代码传送门
    clear
    A=str2sym('[a b; c d]')
    [V,D]=eig(A)subexpr([V;D])
    

在这里插入图片描述

  • 法2:将符号表达式保存成.m文件,利用MATLAB Coder生成C,生成的C中会自动进行CSE处理。但可读性差

在 MATLAB 中,共同子表达式提取(Common Subexpression Elimination,CSE)通常由 MATLAB 编译器自动处理。MATLAB 的 JIT(Just-In-Time)编译器会尝试优化你的代码,包括检测和消除共同的子表达式以提高执行效率。


MATLAB 中 JIT 编译器的工作原理如下:

  • 解释执行和分析: 初始时,MATLAB 会解释执行你的代码,并进行一些分析以了解哪些表达式是重复计算的。

  • JIT 编译: 一旦 MATLAB 确定了可能的优化机会,它会对代码进行 JIT 编译。这时,MATLAB 将生成优化后的机器代码。

  • 执行优化后的代码: 在 JIT 编译之后,MATLAB 将执行优化后的代码,其中可能包括共同子表达式提取等优化。


    虽然 MATLAB 的 JIT 编译器会尝试进行一些优化,但并不保证在所有情况下都能实现最优的共同子表达式提取。在某些特殊情况下,你可能需要手动考虑一些优化策略,如使用局部变量来存储重复计算的结果。


    总的来说,MATLAB 通常会在后台自动处理共同子表达式提取,而无需显式的用户干预。如果你对具体的代码片段有疑虑,可以使用 MATLAB 的 Profiler 工具来分析代码性能,并查看是否存在潜在的优化机会。

参考链接:

  • https://docs.sympy.org/latest/modules/simplify/simplify.html#sympy.simplify.cse_main.cse
  • https://zhuanlan.zhihu.com/p/673052435

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

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

相关文章

[机器学习]KNN——K邻近算法实现

一.K邻近算法概念 二.代码实现 # 0. 引入依赖 import numpy as np import pandas as pd# 这里直接引入sklearn里的数据集,iris鸢尾花 from sklearn.datasets import load_iris from sklearn.model_selection import train_test_split # 切分数据集为训练集和测试…

2024年数学建模美赛 分析与编程

2024年数学建模美赛 分析与编程 1、本专栏将在2024年美赛题目公布后,进行深入分析,建议收藏; 2、本专栏对2023年赛题,其它题目分析详见专题讨论; 2023年数学建模美赛A题(A drought stricken plant communi…

python脚本实现浏览器驱动chromedriver的版本自动升级

chromedriver的版本号与chrome浏览器版本不匹配时在运行程序时就会报错 用下面的脚本可以自动安装chromedriver的最新版本到指定路径 from webdriver_manager.utils import get_browser_version_from_os from webdriver_manager.chrome import ChromeDriverManager import re…

Spring JPA与Hibernate学习使用

EntityManager 在使用持久化工具的时候,一般都有一个对象来操作数据库,在原生的Hibernate中叫做Session,在 JPA 中叫做EntityManager,在MyBatis中叫做SqlSession,通过这个对象来操作数据库。 EntityManager是 JPA 中…

比VS Code快得多

Zed 是一款支持多人协作的代码编辑器,底层采用 Rust,且默认支持 Rust,还自带了 rust-analyzer,主打“高性能”。1 月 24 日,备受关注的 Zed 项目宣布正式开源。 Zed 代码库将采用 Copyleft 许可证,其中编辑…

JavaSE——运算符、运算符优先级、API、Scanner

目录 基本的算术运算符 自增自减运算符 赋值运算符 关系运算符 逻辑运算符 三目运算符 运算符优先级 API Scanner 基本的算术运算符 符号作用加-减*乘/除%取余 基本与C语言的基本算术运算符一致 注意:两个整数相除结果还是整数 public static void main…

C++PythonC# 三语言OpenCV从零开发(7):图像的阈值

文章目录 相关链接前言阈值阈值使用代码PythonCCsharpcsharp代码问题 总结 相关链接 C&Python&Csharp in OpenCV 专栏 【2022B站最好的OpenCV课程推荐】OpenCV从入门到实战 全套课程(附带课程课件资料课件笔记) OpenCV一个窗口同时显示多张图片 …

C 变量

目录 1. C变量 2. C变量定义 2.1 变量初始化 2.2 C中的变量声明 3. C中的左值(Lvalues)和右值(Rvalues) 1. C变量 在C语言中,变量可以根据其类型分为以下几种基本类型: 整型变量:用…

自然语言nlp学习 三

4-8 Prompt-Learning--应用_哔哩哔哩_bilibili Prompt Learning(提示学习)是近年来在自然语言处理领域中,特别是在预训练-微调范式下的一个热门研究方向。它主要与大规模预训练模型如GPT系列、BERT等的应用密切相关。 在传统的微调过程中&a…

Mac 上 Ln 指令使用说明

ln指令 ln 是 Mac 中又一个非常重要命令,它的功能是为某一个文件在另外一个位置建立一个同不的链接,这个命令最常用的参数是 -s,具体用法是:ln –s 源文件 目标文件。 实例: /Users/apple/Desktop git:(master) ✗&g…

将vite项目(vue/react)使用vite-plugin-pwa配置为pwa应用,只需要3分钟即可

将项目配置为pwa模式,就可以在浏览器里面看到安装应用的选项,并且可以将web网页像app一样添加到手机桌面或者pad桌面上,或者是电脑桌面上,这样带来的体验就像真的在一个app上运行一样。为了实现这个目的,我们可以为vue…

Ubuntu18搭建Kubernate集群

1.基本环境 版本介绍 ubuntu 18.04 docker 20.10 kube-apiserver 1.21.3 kube-controller-manage 1.21.3 kube-scheduler 1.21.3 kube-proxy 1.21.3 kubernetes 21.10 2.ubuntu环境必备 Docker环境 参考:Ubuntu18配置Docker 关闭防火墙 ufw disable关闭…

算法设计与分析实验:滑动窗口与二分查找

目录 一、寻找两个正序数组的中位数 1.1 具体思路 1.2 流程展示 1.3 代码实现 1.4 代码复杂度分析 1.5 运行结果 二、X的平方根 2.1 具体思路 2.2 流程展示 2.3 代码实现 2.4 代码复杂度分析 2.5 运行结果 三、两数之和 II-输入有序数组 3.1 采用二分查找的思想 …

【InternLM 大模型实战】作业与笔记汇总

笔记1:书生浦语大模型全链路开源体系 笔记2:轻松玩转书生浦语大模型趣味Demo 作业2:https://github.com/xiaomile/InternLM-homework/tree/main/%E4%BD%9C%E4%B8%9A1 笔记3:基于 InternLM 和 LangChain 搭建你的知识库 作业3&…

LeetCode —— 43. 字符串相乘

😶‍🌫️😶‍🌫️😶‍🌫️😶‍🌫️Take your time ! 😶‍🌫️😶‍🌫️😶‍🌫️😶‍🌫️…

Bio-Rad(Abd serotec)独特性抗体

当一种抗体与另一种抗体的独特型结合时,它被称为抗独特型抗体。抗体的可变部分包括独特的抗原结合位点,称为独特型。独特型(即独特型)内表位的组合对于每种抗体都是独特的。 如今开发的大多数治疗性单克隆抗体是人的或人源化的,用于诱导抗药…

华为策略路由+NQA配置

---NQA--- [RouterA] nqa test-instance admin NQA [RouterA-nqa-admin-vlan10] test-type icmp [RouterA-nqa-admin-vlan10] destination-address ipv4 对方地址 [RouterA-nqa-admin-vlan10] frequency 10 [RouterA-nqa-admin-vlan10] probe-count 2 [RouterA-nqa-admin-vlan…

【国产MCU】-认识CH32V307及开发环境搭建

认识CH32V307及开发环境搭建 文章目录 认识CH32V307及开发环境搭建1、CH32V307介绍2、开发环境搭建3、程序固件下载1、CH32V307介绍 CH32V307是沁恒推出的一款基于32位RISC-V设计的互联型微控制器,配备了硬件堆栈区、快速中断入口,在标准RISC-V基础上大大提高了中断响应速度…

【HTML 基础】标签的属性

文章目录 1. src - 图像和脚本的来源2. alt - 替代文本3. href - 超链接目标4. class - 类5. id - 标识符总结 在 HTML 中,标签不仅定义了文档的结构和内容,还可以通过属性为标签提供附加信息,以调整标签的行为或外观。以下是一些常见的 HTML…

java 社区资源管理系统Myeclipse开发mysql数据库web结构java编程计算机网页项目

一、源码特点 java Web社区资源管系统是一套完善的java web信息管理系统 ,对理解JSP java编程开发语言有帮助,系统具有完整的源代码和数据库,系统主要采用B/S模式开发。开发环境为TOMCAT7.0,Myeclipse8.5开发,数据库为Mysql5.…