| 查看: 884 | 回复: 6 | ||||
[求助]
MATLAB求解复杂方程组，求助！！！
|
|
MATLAB新手。想用MATLAB求解复杂二元多次方程组。方程式如下所述： d/x=-1.131*[log(a/x)]^3+13.635*[log(a/x)^2]-30.594[[log(a/x)]]+29.267； b/ac=(-1.40557y^3+0.77526y^2+0.15830y-0.06831)[ log(a/x)]^3+(17.93006y^3-9.22091y^2-2.37733n+0.86295[[log(a/x)] ^2+(-79.99715y^3+40.55620y^2+9.00157y-2.54543[[log(a/x)]+(122.65069y^3-63.88418y^2-9.58936y+6.20045) 其中，a,b,c,d为参量，他们的数值在Excel表格中，x,y为未知数，我想利用MATLAB，求解出上述方程的解。在a,b,c,d参量变化的情况下，求解出x,y的值。 我看到网上说可以使用m文件自定义函数，也可以使用solve语句，但是我用solve语句写出来如下， syms x y a b c d a= xlsread('20190718.xlsx','sheet2','B2:B6'); b= xlsread('20190718.xlsx','sheet2','C2:C6'); c= xlsread('20190718.xlsx','sheet2','A2:A6'); d= xlsread('20190718.xlsx','sheet2','D2  6');[x,y]=solve('-1.131*[log(a/x)]^3+13.635*[log(a/x)^2]-30.594[[log(a/x)]]+29.267==d/x,b/ac==(-1.40557y^3+0.77526y^2+0.15830y-0.06831)[ log(a/x)]^3+(17.93006y^3-9.22091y^2-2.37733n+0.86295[[log(a/x)] ^2+(-79.99715y^3+40.55620y^2+9.00157y-2.54543[[log(a/x)]+(122.65069y^3-63.88418y^2-9.58936y+6.20045)',x,y) x=simple(x),y=simple(y) 系统提示'"对无效变量名或定义数字的字符串的支持将在将来的版本中删除。创建符号表达式，首先创建符号变量，然后对其使用操作。"" 我不知道该怎么办了，求助，非常感谢！ |
回复此楼» 猜你喜欢
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有140人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有6人回复
-
【修正版】单环用磁约束低速电子实验方案（简化版）
已经有0人回复
-
桂林理工大学物理学专业招收调剂，还有三个名额！！！
已经有22人回复
-
考博自荐
已经有2人回复
-
山东大学第二批博士研究生招生
已经有0人回复
-
中国科学院东莞材料科学与技术研究所-2026年博士招生-吴昊研究员-磁学与自旋电子学
已经有0人回复
-
《电磁学》教材推荐
已经有1人回复
-
【急招】合肥工大核聚变材料计算方向2026级工程博士生
已经有0人回复
2楼2019-07-22 13:25:48
3楼2019-07-22 14:00:00
4楼2019-07-22 14:00:19
5楼2019-07-22 14:02:01
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68137.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
|
你的a b c d每个都是包含5个数的数组 你是想对于每组a b c d都解一次方程组 一共是解5次 对吧？ 提供几点思路： 1 不要用solve函数 这个函数不适合求复杂方程的数值解 改用fsolve函数 2 由于a b c d有5组 如果想一次求出所有的解 需要写个循环 比如用for i=1:5 上述两点 我以前写过类似的代码回在别人的帖子里 具体翻一下我回复记录 3 MATLAB的fsolve函数需要提供初值 如果求解困难 可以试试1stOpt软件来解 发自小木虫Android客户端 |
» 本帖已获得的红花（最新10朵）
Q莹儿
6楼2019-07-22 17:12:18
7楼2019-07-22 23:57:41