| 查看: 296 | 回复: 0 | |||
[求助]
多元函数f=f(t,u,P1,P2,P3)的求解和数值提取和作图的问题
|
|
syms t u P1 P2 P3; a12=21.98*10^7; a111=2.57*10^9; a112=6.95*10^9; a123=13.13*10^9; S11=11*10^-12; S12=-4.1*10^-12; S44=27.5*10^-12; Q11=0.084; Q12=-0.025; Q44=0.035; >> A13=a12+Q12*(Q11+Q12)/(S11+S12); >> A12=a12-[(Q11^2+Q12^2)*S12-2*Q11*Q12*S11]/(S11^2-S12^2)+Q44^2/2/S44; >> eq1=sym('A1=(3.255*t-1960)*10^5-u*(Q11+Q12)/(S11+S12)'); eq2=sym('A3=(3.255*t-1960)*10^5-2*u*Q11/(S11+S12)'); eq3=sym('A11=(-0.3775*t-53.92)*10^5-0.5*[(Q11^2+Q12^2)*S11-2*Q11*Q12*S12]/(S11^2-S12^2)'); eq4=sym('A33=(-0.3775*t-53.92)*10^5+Q12^2/(S11+S12)'); >> eq5=sym('0=2*A1*(P1+P2)+4*A11*(P1^3+P2^3)+2*A13*(P1*P3^2++P2*P3^2+P3*P1^2+P3*P2^2)+2*A3*P3+4*A33*P3^3+2*A12*(P1*P2^2+P1^2*P2)+6*a111*(P1^5+P2^5+P3^5)+2*a112*(2*P1^3*P2^2+2*P1^3*P3^2+P1^4*P2+P3^4*P2+P1^4*P3+2*P3^3*(P1^2+P2^2)+P2^4*P3)+2*a123*(P1*P2^2*P3^2+P1^2*P2*P3^2+P1^2*P2^2*P3);'); >> [t u P1 P2 P3]=solve(eq1,eq2,eq3,eq4,eq5,t,u,P1,P2,P3); Warning: Explicit solution could not be found. > In solve at 160 主要想做出来f=f(t,u,P1,P2,P3)=0时，的取值 谢谢啊 |
回复此楼» 猜你喜欢
假如你的研究生提出不合理要求
已经有12人回复
-
实验室接单子
已经有7人回复
-
全日制（定向）博士
已经有5人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
对氯苯硼酸纯化
已经有3人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有3人回复
找到一些相关的精华帖子，希望有用哦~
f(x,y)和g(x,y)连续,有关f(x,y)/g(x,y)连续的问题？
已经有7人回复
- 求助，有关于matlab作图的问题
已经有7人回复
- 请问一下matlab这个是什么问题
已经有10人回复
- 在origin中如何实现y=ax^b曲线拟合
已经有4人回复
- 如何从函数图像中提取数值，在拟合原函数？
已经有3人回复
- 请问如何画出下面函数f的图形
已经有4人回复
- 如何用mathematica 9做一个关于极坐标的函数（如f(r,theta)）的contour图？
已经有6人回复
- 求助：Matlab中如何绘制三维隐函数图像
已经有20人回复
- 各位大侠，vasp下计算完能带如何作图
已经有12人回复
- 求助数值计算结果和仿真结果有什么不一样的
已经有6人回复
- 请教大虾，一个分段函数的 matlab 程序实现
已经有4人回复
- Ansys载荷步计算结果的保存问题
已经有8人回复
- 求大虾指导 一个关于atoi函数的问题
已经有3人回复
- 关于SIESTA中如何计算功函数的一种做法探讨
已经有8人回复
- 【求助】请高手帮我解答个比较疑惑的问题（UNIFAC估算）
已经有4人回复
- 【求助】matlab的fit的结果读出问题，cfit的类型
已经有7人回复
- 【求助】二元函数拟合
已经有12人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人