| 查看: 742 | 回复: 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) 系统提示'"对无效变量名或定义数字的字符串的支持将在将来的版本中删除。创建符号表达式，首先创建符号变量，然后对其使用操作。"" 我不知道该怎么办了，求助，非常感谢！ |
回复此楼» 猜你喜欢
求取一些关于纳米材料和纳米技术相关的英文PPT。
已经有0人回复
-
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有216人回复
-
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有22人回复
-
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
26申博推荐：南京航空航天大学国际前沿科学研究院光学方向招收博士生！
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
新西兰Robinson研究所招收全奖PhD
已经有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
- 金币: 68120.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