| 查看: 3482 | 回复: 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» 猜你喜欢
各位大神，目前国内有哪些比较好用的逆合成软件？
已经有13人回复
-
国社科系统bug了，是不是要放榜了？
已经有10人回复
-
现代”学阀”该如何界定
已经有9人回复
-
课题组招2027级博士 上海工程技术大学 激光智能制造方向
已经有4人回复
-
上海工程技术大学 激光智能制造课题组 2027级博士研究生招生
已经有3人回复
-
上海工程技术大学激光智能制造课题组｜2027级博士研究生招生公告
已经有11人回复
-
上海工程技术大学激光智能制造课题组-博士研究生
已经有10人回复
-
申博发邮件
已经有11人回复
-
我的奶奶
已经有3人回复
【答案】应助回帖
|
你这样的程序是根本求解不了的。 fsolve函数是求解非线性方程组数值解的，是不能带符号的。 B=I*r^2/2*(r^2+z^2)^(-3/2)；这里面的I与I1、I2是什么关系？ 把d3Bdz3，d5Bdz5，d7Bdz7写出来，分别带入 I1*d3Bdz3(z=a1)+I2*d3Bdz3(z=a2)=0 I1*d5Bdz5(z=a1)+I2*d5Bdz5(z=a2)=0 I1*d7Bdz7(z=a)+I2*d7Bdz7(z=a2)=0 三个方程，并写入f(x)中，I1、I2、r作为x(4)、x(5)、x(6)，即F内六个表达式。 你先试试吧，不行了把数学模型用图贴出来吧，把已知条件和未知条件描述清楚。 |
17楼2012-05-10 23:25:01
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