| 查看: 1074 | 回复: 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) 系统提示'"对无效变量名或定义数字的字符串的支持将在将来的版本中删除。创建符号表达式，首先创建符号变量，然后对其使用操作。"" 我不知道该怎么办了，求助，非常感谢！ |
回复此楼» 猜你喜欢
青基已中
已经有1人回复
-
散金币祈福
已经有186人回复
-
物理学I论文润色/翻译怎么收费?
已经有223人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有109人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有64人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有4人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有26人回复
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
- 金币: 68175.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.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
10