MATLAB求解一个比较复杂的方程 - 仿真模拟 - MATLAB - 第2页小木虫论坛-学术科研互动平台

明德厚学、求是创新

11楼2011-07-01 17:20:48

cyzhixue

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 9.5
红花: 1
帖子: 25
在线: 7.7小时
虫号: 1260375
注册: 2011-04-09
专业: 高分子材料结构与性能

引用回帖:

Originally posted by 信彼南山 at 2011-07-01 16:38:26:
运行结果：
a =

-4.9837156546668308559867719464884-.34823945996364925341340692246585*i
-4.9837156546668308559867719464884+.34823945996364925341340692246585*i
4.983715654666830855986771946 ...

哇塞，，你们都好强啊。。我这个菜鸟真是太佩服了。。向你们学习！！谢谢你们啦。。

12楼2011-07-03 21:59:24

cyzhixue

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 9.5
红花: 1
帖子: 25
在线: 7.7小时
虫号: 1260375
注册: 2011-04-09
专业: 高分子材料结构与性能

引用回帖:

Originally posted by xiegangmai at 2011-06-30 11:15:59:
呵呵，应该可以了，换用'levenberg-marquardt'算法，增加迭代次数。

用fsolve求解，初值选取有很大关系，当初值选为[ 100; 100; 0 ]时，
[ X, Fval, exitflag ] = fsolve( @cyzhixue, [ 100; 100; 0 ], optio ...

你好，谢谢你回答我的问题，我还想问一下就是这个初值的选择有什么依据吗？？我又换了个方程，选了好多都不行啊。。数学没学好呵呵。。

13楼2011-07-03 22:07:29

Luckner

新虫 (初入文坛)

应助: 1 (幼儿园)
金币: 31
帖子: 24
在线: 8.7小时
虫号: 1516658
注册: 2011-11-30
性别: GG
专业: 燃烧学

引用回帖:

3楼: Originally posted by xiegangmai at 2011-06-30 11:03:09:
尝试了一下用符号求解，找不到解析解。

[code]syms a b c;
solve('a/(1.0+exp(-b*(0.003-c)))-a/(1.0+exp(b*c))=1.460','a/(1.0+exp(-b*(0.012-c)))-a/(1.0+exp(b*c))=4.388','a/(1.0+exp(-b*(0.030- ...

请问这样的问题怎么解决呀？谢谢
函数是：
function f=myhs1(x)
f=[(2*x(1)+3*x(2)+x(3))^2*(x(1)+x(2)-x(3))/((0.4125-x(1)-x(2))*(0.21*0.825-0.5*x(1))^0.5*(1.65+1.5*x(1)+2*x(2))^1.5)-3.389*exp(11);(2*x(1)+3*x(2)+x(3))^3*(x(1)+x(2)-x(3))/((0.4125-x(1)-x(2))*(0.4125-x(2)-x(3))*(1.65+1.5*x(1)+2*x(2))^2)-1.155;(2*x(1)+3*x(2)+x(3))*x(3)/((x(1)+x(2)-x(3))*(0.4125-x(2)-x(3)))-8.469];

运行如下：
>> x0=[0.1 0.1 0.1];
>> [r,val]=fsolve(@myhs1,x0,opt)

结果为：
Solver stopped prematurely.
r =

0.3465 0.0658 0.2101

val =

1.0e+003 *

0.4843
1.8293
-0.0001

我试了好多组初值都不行，要么出现上面的情况，要么是“No solution found.”
请问这是怎么回事呀？还是哪儿出问题了？
不胜感激！

物换星移，沧海桑田，人心却没有多少改变。

14楼2011-12-19 22:01:00

lylpiao

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 6761
红花: 1
帖子: 391
在线: 233.6小时
虫号: 1758002
注册: 2012-04-16
性别: MM
专业: 微/纳机械系统

送红花一朵

想问一下楼主，会不会解这个方程：
cos(v-t)*(-k*sin(h)+c*i*cos(h))-(i*sin(h)*(n*t+n*pi-b+k*v)+i*a*cos(h))*sin(v-t)=0;其中a,b,c,h,i都是已知的，只有v和t未知，求v=f(t)的形式，要怎么求？

把一切交给时间！

15楼2015-12-14 11:08:26