使用octave符号运算求解不定积分、微分方程等(兼容matlab)
1.求解1/(1+cos(x))^2的不定积分。
在和学生讨论一道物理竞赛题的时候,出现了这个函数的积分求解需求。查积分表也可写出答案。但是可以使用octave的符号运算工具箱来做。
syms x; y = 1/(1+cos(x))^2; int(y)
既可以得到结果:
ans = (sym)3/x\ /x\tan |-| tan|-|\2/ \2/------- + ------6 2
octave中的符号工具箱实际上是调用了sympy的核心库。所以看自来结果有符号艺术的感觉。
2.求解微分方程
在一些场合,需要求解微分方程的解析解。
举个例子,求解谐振子的微分方程。hamonic.m
syms x(t) DE = diff(x,2) + w^2*x == 0 dsolve(DE)
得到的结果是:
>> hamonicDE = (sym)22 dw *x(t) + ---(x(t)) = 02dtans = (sym)-I*t*w I*t*wx(t) = C1*e + C2*e
这里面octave的微分方程表达方法和matlab不大一样,二阶导函数不是写作D2(x),而是使用 diff(x,2)来表达。
转载于:https://www.cnblogs.com/xiaohezyd/p/10799030.html
使用octave符号运算求解不定积分、微分方程等(兼容matlab)相关推荐
- matlab中符号运算求解结果出现的是1i不是li
说实话,也就我这种粗心大头才研究了半天为啥是li,早看到1i就能早点明白这不过是代表复数 因为 i 容易和变量搞混了,所以matlab在使用符号函数求解时,其结果会推荐你使用 1i .而且当你可以将 ...
- matlab符号函数绘图法_《MATLAB符号运算及其应用》(黄忠霖著).pdf
<MATLAB符号运算及其应用>(黄忠霖著).pdf [General Information] 书名=MATLAB符号运算及其应用 作者=黄忠霖,黄京编著 页数=425 出版社=北京市: ...
- matlab常见符号运算( 因式分解、展开、合并、简化及通分等,计算导数,积分,符号求和、代数方程和微分方程求解)
常见符号运算 因式分解.展开.合并.简化及通分等 因式分解factor(f) factor 也可用于正整数的分解 函数展开 expand(f) 多项式展开 三角函数展开 合并同类项 collect(f ...
- 【MATLAB】符号数学计算(七):符号微积分、符号微分方程求解、符号代数方程求解
一.符号表达式的极限 limit(F,x,a):求当时,符号表达式F的极限. limit(F,a):符号表达式F采用默认自变量(可由函数findsym求得),该函数求F的自变量趋于a时的极限值. li ...
- matlab整理符号表达式,[2018年最新整理]MATLAB符号运算与符号方程求解.ppt
[2018年最新整理]MATLAB符号运算与符号方程求解 MATLAB符号计算 1 符号对象 2 符号微积分 3 级 数 4 符号方程求解 9.1 符号对象 9.1.1 建立符号对象 1.建立符号变量 ...
- 「SymPy」符号运算(3) (非)线性方程(组)求解、数列求和、连乘、求极限
目录 导言 解方程(组) solve函数 solveset函数 求和 ∑ \sum ∑ 连乘 ∏ \prod ∏ 求函数极限 求数列极限 导言 在前两篇文章中,我们学习了SymPy的输入输出.基本符号 ...
- 数学建模学习笔记(19)Matlab符号运算
% Matlab的符号运算% 1.符号变量的创建 % syms 符号变量名 % 注释:可以用一个syms语句同时创建多个变量 syms y; syms a b;% 2.符号方程的创建 % 符号方程可以 ...
- matlab数值与符号运算
matla数值与符号运算 1. 多项式计算 主函数 clc clear close %多项式 x^4-12*x^3+25*x+116 p=[1 -12 0 25 116]%多项式求值函数 polyva ...
- 「SymPy」符号运算(4) 微积分与有限差分
目录 导言 积分 不定积分 定积分 多重积分 求导 一阶导数 高阶导数 偏导数 有限差分 常微分差分 差分系数 高阶差分 偏微分差分 导言 在前几篇中,我们学习了SymPy的基本语法.方程求解等基础知 ...
- MATLAB08:符号运算
pdf版本笔记的下载地址: MATLAB08_符号运算(访问密码:3834) MATLAB08:符号运算 创建符号变量 创建符号数字 创建符号变量 符号运算 符号表达式的化简与代入 符号表达式的化简 ...
最新文章
- 上交大博士生:博三经历5轮审稿被拒后,如今发10篇论文成为审稿人
- Spring中的@ Component,@ Repository和@Service批注有什么区别?
- 华为鸿蒙王者荣耀视频,王者荣耀鸿蒙2.0版
- c语言中栈区运用原理形象图,C语言实现使用动态数组来构造栈结构
- 爱奇艺深度语义表示学习的探索与实践
- colab加载google drive并且调试时跳过不运行
- maven建立webapp项目时显示Cannot change version of project facet Dynamic web module to 2.5
- 【CodeForces - 608C】Chain Reaction (二分 或 dp ,思维)
- git checkout 和 git reset
- PAP CHAP认证
- 央行宣布降准0.5个百分点 释放长期资金约8000多亿元
- 测试思路系列:《谷歌的软件测试之道》读书笔记
- 【Go语言】【16】GO语言的并发
- 一个工具类实现自定义Tablayout的下划线宽度
- dos中进入其他盘中的方法
- Python 入门演示 1
- UNIX网络编程 第1卷 数据分享
- 武汉工程大学matlab,Lorenz系统动力学行为的MATLAB仿真与分析[1]
- 吉林大学珠海学院计算机录取分数线,大学介绍 | 吉林大学珠海学院(附录取分数线,重点专业)...
- win10 设备管理器显示便捷设备
热门文章
- 拓端tecdat|使用Python中Keras的LSTM递归神经网络进行时间序列预测
- 拓端tecdat|R语言在BRFSS数据中可视化分析探索糖尿病的影响因素
- (10)数据分析-变量分析
- 基于MATLAB 关于圆的研究
- libSM.so.6: cannot open shared object file: No such file or directoryapt-file search libSM.so.6
- Linux学习笔记(8)
- 机器学习---分类、回归、聚类、降维的区别
- python LDA实践入门学习
- 周志华机器学习西瓜书速记第一章绪论
- Hyperledger Fabric教程(10)-- peer命令-链码chaincode