| 查看: 1556 | 回复: 8 | ||
[求助]
matlab求解非线性方程组报错，请各位大神指点！ 已有1人参与
|
|
clear >> syms n beta; >> f1=sym(\'(f01-4*(2-n)*n*f01^(n-1)*sigma1^2*beta*yita1)/(1+2*n*(n-1)*f01^(n-2)*sigma1^2*beta*yita1)=fc1\') f1 = (f01 + 4*f01^(n - 1)*n*sigma1^2*yita1*(n - 2)*beta)/(2*f01^(n - 2)*n*sigma1^2*yita1*(n - 1)*beta + 1) = fc1 >> f2=sym(\'(f02-4*(2-n)*n*f02^(n-1)*sigma2^2*beta*yita2)/(1+2*n*(n-1)*f02^(n-2)*sigma2^2*beta*yita2)=fc2\') f2 = (f02 + 4*f02^(n - 1)*n*sigma2^2*yita2*(n - 2)*beta)/(2*f02^(n - 2)*n*sigma2^2*yita2*(n - 1)*beta + 1) = fc2 >> fo1=4.92; >> fc1=3.76; >> fo2=2.45; >> fc2=2.31; >> sigma1=1.50; >> sigma2=0.635; >> [n,beta]=solve(f1,f2) Warning: Explicit solution could not be found. > In solve at 81 n = [ empty sym ] beta = [] f1和f2是所求的方程组，算法总是报错，不知道该如何求解。。。 |
回复此楼» 猜你喜欢
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
-
材料工程专硕274一志愿211求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- MATLAB解非线性方程组
已经有9人回复
- MATLAB求解五元非线性方程组，在线等
已经有7人回复
- 一个关于有约束非线性方程组的求解问题
已经有7人回复
- matlab 非线性微分方程求解
已经有3人回复
- matlab如何求解带参数的线性方程组
已经有6人回复
- 用matlab求解一个非线性方程组的解
已经有4人回复
- 请教一个matlab求解非线性方程组的问题
已经有9人回复
- 用matlab求解非线性方程组说无解，一定是方程组本身无解，还是有可能程序有问题呢？
已经有11人回复
- 运行matlab出现错误
已经有4人回复
- MATLAB求解非线性方程组
已经有5人回复
- 【求助】非线性方程组的求解问题
已经有6人回复
- 求Matlab高手解决线性方程组的迭代求解问题
已经有17人回复
- 求助matlab 非线性方程组全部解！
已经有6人回复
- matlab求解非线性方程组，错误提示怎么解决
已经有5人回复
- 求助matlab---fsolve解非线性方程组
已经有6人回复
- matlab求解非线性方程组
已经有16人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- 请教Matlab高手解多元非线性方程组中的参数
已经有6人回复
- matlab的fsove 命令求解非线性方程组
已经有6人回复
- 【求助】matlab求解非线性方程组，并画图处理。要求y,z是实数解！
已经有18人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
2楼2014-09-09 09:19:42
NicoleLas
木虫 (正式写手)
- 应助: 66 (初中生)
- 金币: 4288.7
- 红花: 26
- 帖子: 514
- 在线: 177.2小时
- 虫号: 2872321
- 注册: 2013-12-15
- 性别: GG
- 专业: 计算机软件
3楼2014-09-09 10:18:27
4楼2014-09-09 13:05:38
NicoleLas
木虫 (正式写手)
- 应助: 66 (初中生)
- 金币: 4288.7
- 红花: 26
- 帖子: 514
- 在线: 177.2小时
- 虫号: 2872321
- 注册: 2013-12-15
- 性别: GG
- 专业: 计算机软件
【答案】应助回帖
感谢参与，应助指数 +1
|
为图方便，给你用1stOpt编程算了下。 代码： 采用L-M+UGO结果： 可能结果1： n: 1 betta: 0.0879576864290957 yita1: 1.46535105823627 yita2: 0.986840117382106 可能结果2： n: 1 betta: 0.235202634137847 yita1: 0.547990839309009 yita2: 0.369044224009309 上述结果都不一定正确，因为我不知道你的变量的取值范围，而且是否是连续或者整数变量，我假设n为非负整数变量，其余均为实数域连续变量了，你根据你的问题修改下吧，也可以采取其它算法试一下。 |
5楼2014-09-09 18:18:06
|
高手，想跟您请教一下，plot函数画图的时候，为什么x=1:1:10和x=1:0.1:10得到的函数图像走势不同？应该是一样的呀。。。程序代码如下：clear all; >> n=1:1:10; >> y=(4.92-3.72)/(2.45-2.34)-((1.5/0.635)^2)*((4.92/2.45).^(n-2)).*(((n-1)*3.72+(4-2*n)*4.92)./((n-1).*2.34+(4-2*n)*2.45)); >> plot(n,y) 以及 >> clear all; >> n=1:0.1:10; >> y=(4.92-3.72)/(2.45-2.34)-((1.5/0.635)^2)*((4.92/2.45).^(n-2)).*(((n-1)*3.72+(4-2*n)*4.92)./((n-1).*2.34+(4-2*n)*2.45)); >> plot(n,y) 两种情况的函数图像不同，您能给指点一下吗？ |
6楼2014-09-10 20:18:39
NicoleLas
木虫 (正式写手)
- 应助: 66 (初中生)
- 金币: 4288.7
- 红花: 26
- 帖子: 514
- 在线: 177.2小时
- 虫号: 2872321
- 注册: 2013-12-15
- 性别: GG
- 专业: 计算机软件
7楼2014-09-10 21:53:23
8楼2014-09-11 08:09:51
NicoleLas
木虫 (正式写手)
- 应助: 66 (初中生)
- 金币: 4288.7
- 红花: 26
- 帖子: 514
- 在线: 177.2小时
- 虫号: 2872321
- 注册: 2013-12-15
- 性别: GG
- 专业: 计算机软件
9楼2014-09-11 08:46:51
20