| 查看: 3412 | 回复: 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» 猜你喜欢
2026国自然函评费到账
已经有17人回复
-
让我中一个面上吧！
已经有15人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有5人回复
-
跳槽后在研项目怎么办？
已经有10人回复
-
什么时候开奖？
已经有10人回复
-
今天放榜吗？
已经有16人回复
-
93BebMhtakh前后11位开头都是大写
已经有8人回复
-
范进中举一文的中心思想
已经有4人回复
-
只有每年这种时候来逛逛小木虫
已经有26人回复
-
filecode，4个jtjc了
已经有16人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 如何使用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人回复
|
其实我是把自己的方程简化了下，我这里有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
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
40