| 查看: 3581 | 回复: 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 |
回复此楼» 收录本帖的淘帖专辑推荐
 程序 |
» 猜你喜欢
请问有评职称，把科研教学业绩算分排序的高校吗
已经有4人回复
-
孩子确诊有中度注意力缺陷
已经有12人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
-
2026博士申请-功能高分子，水凝胶方向
已经有6人回复
-
论文投稿，期刊推荐
已经有4人回复
-
硕士和导师闹得不愉快
已经有13人回复
-
请问2026国家基金面上项目会启动申2停1吗
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
请问：曲线方程如何用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
- 金币: 18061.7
- 红花: 97
- 帖子: 8891
- 在线: 1485.3小时
- 虫号: 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
- 金币: 18061.7
- 红花: 97
- 帖子: 8891
- 在线: 1485.3小时
- 虫号: 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