| 查看: 1577 | 回复: 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是所求的方程组，算法总是报错，不知道该如何求解。。。 |
回复此楼» 猜你喜欢
267求调剂
已经有13人回复
-
292求调剂
已经有3人回复
-
322求调剂
已经有8人回复
-
348求调剂
已经有7人回复
-
材料调剂
已经有5人回复
-
349求调剂
已经有4人回复
-
环境工程 085701，267求调剂
已经有16人回复
-
一志愿北交材料工程总分358
已经有5人回复
-
材料0856 英一数二 323 求调剂
已经有8人回复
-
318求调剂
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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
