| 查看: 3721 | 回复: 5 | |||
[求助]
如何用matlab对函数表达式去除虚部 已有2人参与
|
|
请问用matlab在求解函数表达式时，计算得到的表达式运行结果出现虚部，要如何去除？这是我计算用的的代码，想要求通过eq1和eq2求出 x 对 j 的函数表达式，而且保证 x 为正实数。 >> clear >> syms x y j >> eq1=-6.4e7*x+5.6e3*y-2e-11*x*y+j/9.6e-22; %通过两条关系式求解 x 对 j 的表达式 eq2=3.1e7*x-6420*y+3*j/9.6e-22; >> [x,y]=solve(eq1,eq2); >> vpa(x,5) %保留5位小数 ans = - 5.0403e13*j - 0.00010355*(2.3694e35*j^2 + 5.4375e36*j + 3.415e36)^(1/2) - 1.9135e14 0.00010355*(2.3694e35*j^2 + 5.4375e36*j + 3.415e36)^(1/2) - 5.0403e13*j - 1.9135e14 |
回复此楼» 收录本帖的淘帖专辑推荐
 程序 |
» 猜你喜欢
材料考研调剂
已经有25人回复
-
材料工程281还有调剂机会吗
已经有11人回复
-
303求调剂
已经有12人回复
-
290调剂生物0860
已经有14人回复
-
农学0904 312求调剂
已经有3人回复
-
求调剂
已经有7人回复
-
11408。358求调剂
已经有3人回复
-
求助调剂，跨调
已经有8人回复
-
271求调剂
已经有8人回复
-
一志愿厦大0856，306求调剂
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请问：曲线方程如何用Matlab编程~~Matlab菜鸟，比较急！谢谢~~~
已经有3人回复
- 如何用Matlab对函数sin(pi*t)/(pi*t)进行FFT变换
已经有5人回复
- 如何用Matlab拟合大量散点的边界曲线
已经有5人回复
- 请教高手，如何用MATLAB对正弦谐波进行傅里叶分解求出幅值谱和相位谱
已经有2人回复
- 如何用matlab实现分峰拟合？
已经有18人回复
- 求助矩阵高手，如何用Matlab求解矩阵对的不变零点
已经有2人回复
- 如何用matlab对图片中曲线做FFT变换
已经有0人回复
- 如何用matlab对四维数据进行非线性拟合？
已经有5人回复
- 如何用MATLAB实现对曲线扫描的运动仿真？？？
已经有0人回复
- 重金请教-如何用MAtlab对力学模型进行参数识别
已经有12人回复
- 如何用matlab或其他软件批量求解对数方程
已经有6人回复
- 如何用Matlab画出一般多边形内部的Voronoi图
已经有2人回复
- 如何用matlab求一个式子的整数解
已经有4人回复
- 如何用matlab对幂指数函数进行傅里叶变换
已经有1人回复
- 如何用MATLAB软件对细菌多态性进行主成分分析
已经有5人回复
FMStation
至尊木虫 (知名作家)
- 应助: 591 (博士)
- 贵宾: 0.03
- 金币: 18663.7
- 红花: 97
- 帖子: 8891
- 在线: 1491.2小时
- 虫号: 2400059
- 注册: 2013-04-04
- 专业: 计算机应用技术
2楼2016-08-24 21:28:22
送红花一朵 |
还是不行啊，我现在最主要的问题是对下面的函数 a 做出 j （0.01，,100）的图像，因为存在虚部，一直做不出来... >> clear >> a = sym('a','real'); b = sym('b','real'); j = sym('j','real'); >> a=0.0002071*(6.0273e53*j + ((6.0273e53*j - 9.869e53)^2 - 1.0*(1.6667e35*j + 1.0722e36)^3)^(1/2) - 9.869e53)^(1/3) - 1.0081e14*j + 2.0161e-22*((6.0273e53*j + ((6.0273e53*j - 9.869e53)^2 - 1.0*(1.6667e35*j + 1.0722e36)^3)^(1/2) - 9.869e53)^(1/3) + (1.6667e35*j + 1.0722e36)/(6.0273e53*j + ((6.0273e53*j - 9.869e53)^2 - 1.0*(1.6667e35*j + 1.0722e36)^3)^(1/2) - 9.869e53)^(1/3) - 1.3057e18)^2 + (3.4517e31*j + 2.2205e32)/(6.0273e53*j + ((6.0273e53*j - 9.869e53)^2 - 1.0*(1.6667e35*j + 1.0722e36)^3)^(1/2) - 9.869e53)^(1/3) - 2.7041e14; >> ezplot(j,a) 警告: 复数 X 和/或 Y 参数的虚部已忽略 > In ezplot>ezparam (line 392) In ezplot (line 165) In sym/ezplot (line 67) |
3楼2016-08-25 09:59:52
|
即使加上您回答的代码，也不行。依然警告有虚部的存在... >> clear >> a = sym('a','real'); b = sym('b','real'); j = sym('j','real'); >> a=0.0002071*(6.0273e53*j + ((6.0273e53*j - 9.869e53)^2 - 1.0*(1.6667e35*j + 1.0722e36)^3)^(1/2) - 9.869e53)^(1/3) - 1.0081e14*j + 2.0161e-22*((6.0273e53*j + ((6.0273e53*j - 9.869e53)^2 - 1.0*(1.6667e35*j + 1.0722e36)^3)^(1/2) - 9.869e53)^(1/3) + (1.6667e35*j + 1.0722e36)/(6.0273e53*j + ((6.0273e53*j - 9.869e53)^2 - 1.0*(1.6667e35*j + 1.0722e36)^3)^(1/2) - 9.869e53)^(1/3) - 1.3057e18)^2 + (3.4517e31*j + 2.2205e32)/(6.0273e53*j + ((6.0273e53*j - 9.869e53)^2 - 1.0*(1.6667e35*j + 1.0722e36)^3)^(1/2) - 9.869e53)^(1/3) - 2.7041e14; >> b= expand(real(a)); >> ezplot(j,b) 警告: 复数 X 和/或 Y 参数的虚部已忽略 > In ezplot>ezparam (line 392) In ezplot (line 165) In sym/ezplot (line 67) |
4楼2016-08-25 10:03:41
FMStation
至尊木虫 (知名作家)
- 应助: 591 (博士)
- 贵宾: 0.03
- 金币: 18663.7
- 红花: 97
- 帖子: 8891
- 在线: 1491.2小时
- 虫号: 2400059
- 注册: 2013-04-04
- 专业: 计算机应用技术
5楼2016-08-25 12:37:07
Mr__Right
专家顾问 (著名写手)
-
专家经验: +31 - 数学EPI: 1
- 应助: 317 (大学生)
- 金币: 14456.3
- 散金: 500
- 红花: 54
- 帖子: 2716
- 在线: 950.6小时
- 虫号: 1972612
- 注册: 2012-09-04
- 性别: GG
- 专业: 应用数学方法
- 管辖: 程序语言
6楼2016-08-27 15:06:57
