一段matlab转换为对应的python代码问题

DarkMatter 2009-12-24 07:20:07
有如下一段matlab代码,需要改为python代码实现对应的功能:
(核心部分就两行,一个是包含符号变量的函数声明,一个是函数绘图,这两个功能用sympy和pylab分别可以实现,但无法统一在一起)
syms alfa;
R = 72.5; lamda = R/262;
piston_disp = R*(1-cos(alfa)+lamda/2*sin(alfa)^2);
ezplot(piston_disp);
adc10_disp = subs(piston_disp, alfa, 10/360*2*pi);
本来计划使用python的sympy模块,函数好办,但ezplot对应的Plot函数总出现如下错误,无法运行(windows平台,sympy v0.6.6):
Window initialization failed: Unable to share contexts
请高手指点。
谢谢!
...全文
536 4 打赏 收藏 转发到动态 举报
写回复
用AI写文章
4 条回复
切换为时间正序
请发表友善的回复…
发表回复
DarkMatter 2010-05-26
  • 打赏
  • 举报
回复
[Quote=引用 2 楼 thy38 的回复:]

为什么不用matplotlib?
http://matplotlib.sourceforge.net/
[/Quote]
matplotlib无法进行符号运算,这个问题的难点是要将符号运算和作图放在一起以保证没有对数据流的手工操作,看来实现起来确实比较困难。
DarkMatter 2010-05-26
  • 打赏
  • 举报
回复
[Quote=引用 1 楼 yanqlv 的回复:]

Matlab和.Net,Java都有接口,能用plot,可试试先做成.Net或Java组件
[/Quote]
这个代价太大了吧
thy38 2009-12-28
  • 打赏
  • 举报
回复
为什么不用matplotlib?
http://matplotlib.sourceforge.net/
yanqlv 2009-12-27
  • 打赏
  • 举报
回复
Matlab和.Net,Java都有接口,能用plot,可试试先做成.Net或Java组件
内容概要:本文档聚焦于基于Peng-Robinson状态方程的化工热力学计算研究,系统介绍了利用Matlab代码实现纯组分及多组分系统的压缩因子(z因子)和逸度系数的计算方法,并进一步拓展至泡点压力与露点压力的数值求解。资源涵盖了从状态方程推导、非线性方程求根、相平衡条件判断到迭代算法实现的完整技术路线,提供了可运行的Matlab程序代码,有助于深入理解实际工程中物性参数的计算原理与实现过程; 适合人群:具备化工热力学基础知识和一定Matlab编程能力的高校学生、科研人员及工程技术人员,尤其适合从事过程系统模拟、油气田开发、精馏设计及相关领域研究的专业人士; 使用场景及目标:① 掌握Peng-Robinson方程在真实气体物性计算中的应用;② 实现多组分混合体系的相平衡计算与泡/露点预测;③ 为Aspen Plus等商用流程模拟软件提供底层算法验证与教学参考;④ 支持科研论文复现、课程设计、毕业设计及工业过程优化中的热力学数据支撑需求; 阅读建议:建议读者结合经典化工热力学教材同步学习,重点理解立方型状态方程的根判别、迭代收敛策略与相平衡准则的数学表达,动手调试并运行所提供的Matlab代码以掌握算法细节,还可将其扩展应用于其他状态方程(如SRK、vdW)的对比分析研究。

37,739

社区成员

发帖
与我相关
我的任务
社区描述
JavaScript,VBScript,AngleScript,ActionScript,Shell,Perl,Ruby,Lua,Tcl,Scala,MaxScript 等脚本语言交流。
社区管理员
  • 脚本语言(Perl/Python)社区
  • WuKongSecurity@BOB
加入社区
  • 近7日
  • 近30日
  • 至今

试试用AI创作助手写篇文章吧