| 查看: 3132 | 回复: 19 | ||
[求助]
求助！matlab用fsolve函数求解非线性方程组的问题！
|
|
我要用fsolve函数求解一个非线性方程组，程序代码如下： function p=fsolveDemon1 x0=unifrnd(0.1,0.1); options=optimset('MaxFunEvals',20000,'MaxIter',2000); p=fsolve(@f,x0,options); function F=f(x) syms I1 I2 F=[I1*x(1)^2*(3*x(3)*(x(1)^2+x(3)^2)^(-7/2)-7*x(3)^3*(x(1)^2+x(3)^2)^(-9/2))+I2*x(2)^2*(3*x(3)*(x(2)^2+x(3)^2)^(-7/2)-7*x(3)^3*(x(2)^2+x(3)^2)^(-9/2)); I1*x(1)^2*(-5*x(3)*(x(1)^2+x(3)^2)^(-9/2)+30*x(3)^3*(x(1)^2+x(3)^2)^(-11/2)-33*x(3)^5*(x(1)^2+x(3)^2)^(-13/2))+I2*x(2)^2*(-5*x(3)*(x(2)^2+x(3)^2)^(-9/2)+30*x(3)^3*(x(2)^2+x(3)^2)^(-11/2)-33*x(3)^5*(x(2)^2+x(3)^2)^(-13/2)); I1*x(1)^2*(35*x(3)*(x(1)^2+x(3)^2)^(-11/2)-385*x(3)^3*(x(1)^2+x(3)^2)^(-13/2)+1001*x(3)^5*(x(1)^2+x(3)^2)^(-15/2)-715*x(3)^7*(x(1)^2+x(3)^2)^(-17/2))+I2*x(2)^2*(35*x(3)*(x(2)^2+x(3)^2)^(-11/2)-385*x(3)^3*(x(2)^2+x(3)^2)^(-13/2)+1001*x(3)^5*(x(2)^2+x(3)^2)^(-15/2)-715*x(3)^7*(x(2)^2+x(3)^2)^(-17/2))]; end end 可是运行总是出现错误，请求各位高手帮帮我，我刚开始学习matlab,发现这个论坛特别好，刚刚注册，还没有太多的金币送给大家！不胜感激！！！ |
回复此楼» 本帖已获得的红花（最新10朵）
521438» 猜你喜欢
求调剂
已经有23人回复
-
301求调剂
已经有15人回复
-
304求调剂（085602，过四级，一志愿985）
已经有17人回复
-
302分求调剂 一志愿安徽大学085601
已经有12人回复
-
288环境专硕,求调材料方向
已经有23人回复
-
环境专硕调剂
已经有6人回复
-
22408 调剂材料
已经有6人回复
-
285求调剂
已经有12人回复
-
求调剂
已经有6人回复
-
085600材料与化工301分求调剂院校
已经有19人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 如何使用matlab求解非线性方程组的所有整数解？
已经有9人回复
- 求助解6变量的指数方程组的方法
已经有4人回复
- 求助有关Matlab有约束非线性最优化问题
已经有9人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- matlab的fsove 命令求解非线性方程组
已经有6人回复
- MATLAB非线性方程组求解~~~反应釜模拟
已经有9人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】matlab求解非线性方程组，并画图处理。要求y,z是实数解！
已经有18人回复
- 【求助】MATLAB 如何解矩阵方程（达人救命）
已经有6人回复
- 【求助】求助解决这样一个非线性方程组
已经有8人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
- 【求助】非线性方程组求解
已经有3人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
- 【求助】用fsolve求解非线性一元方程
已经有5人回复
2楼2012-05-07 20:49:56
jinpeng_6118
木虫 (正式写手)
- 应助: 28 (小学生)
- 金币: 5027
- 散金: 25
- 红花: 2
- 帖子: 856
- 在线: 147.9小时
- 虫号: 1019189
- 注册: 2010-05-15
- 性别: GG
- 专业: 微生物遗传学
3楼2012-05-08 16:01:50
4楼2012-05-08 17:15:11
jinpeng_6118
木虫 (正式写手)
- 应助: 28 (小学生)
- 金币: 5027
- 散金: 25
- 红花: 2
- 帖子: 856
- 在线: 147.9小时
- 虫号: 1019189
- 注册: 2010-05-15
- 性别: GG
- 专业: 微生物遗传学
5楼2012-05-08 17:40:14
6楼2012-05-08 18:54:23
★
xiegangmai: 金币+1, 谢谢参与！ 2012-05-08 23:09:32
xiegangmai: 金币+1, 谢谢参与！ 2012-05-08 23:09:32
|
本帖内容被屏蔽 |
7楼2012-05-08 21:05:46
jinpeng_6118
木虫 (正式写手)
- 应助: 28 (小学生)
- 金币: 5027
- 散金: 25
- 红花: 2
- 帖子: 856
- 在线: 147.9小时
- 虫号: 1019189
- 注册: 2010-05-15
- 性别: GG
- 专业: 微生物遗传学
8楼2012-05-08 21:46:20
|
其实我是把自己的方程简化了下，我这里有6个未知数(I1，I2,r1,r2,a1,a2.)3个方程，求不出来具体的数值的。具体方程如下： F=[I1*r1^2*(3*a1*(r1^2+a1^2)^(-7/2)-7*a1^3*(r1^2+a1^2)^(-9/2))+I2*r2^2*(3*a2*(r2^2+a2^2)^(-7/2)-7*a2^3*(r2^2+a2^2)^(-9/2)); I1*r1^2*(-5*a1*r1^2+a1^2)^(-9/2)+30*a1^3*(r1^2+a1^2)^(-11/2)-33*a1^5*(r1^2+a1^2)^(-13/2))+I2*r2^2*(-5*a2*(r2^2+a2^2)^(-9/2)+30*a2^3*(r2^2+a2^2)^(-11/2)-33*a2^5*(r2^2+a2^2)^(-13/2)); I1*r1^2*(35*a1*(r1^2+a1^2)^(-11/2)-385*a1^3*(r1^2+a1^2)^(-13/2)+1001*a1^5*(r1^2+a1^2)^(-15/2)-715*a1^7*(r1^2+a1^2)^(-17/2))+I2*r2^2*(35*a2*(r2^2+a2^2)^(-11/2)-385*a2^3*(r2^2+a2^2)^(-13/2)+1001*a2^5*(r2^2+a2^2)^(-15/2)-715*a2^7*(r2^2+a2^2)^(-17/2))]; 我加了约束条件a1=a2，这样还有5个未知数于是我就先不管I1,I2。只要能求出r1,与r2,以及r1与a1之间的关系就行！即r1/r2=具体的数值，r1/a1=具体的数值就行。看看你会么？ |
9楼2012-05-09 10:49:36
10楼2012-05-09 12:31:19
