24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

521438

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 82.9
帖子: 19
在线: 13.2小时
虫号: 1801287
注册: 2012-05-07
性别: MM
专业: 电磁场与电路

[求助] 求助！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

» 猜你喜欢

回收溶剂求助已经有7人回复
职称评审没过，求安慰已经有40人回复
硝基苯如何除去已经有3人回复
A期刊撤稿已经有4人回复
垃圾破二本职称评审标准已经有17人回复
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿已经有22人回复
申请26博士已经有5人回复
EST投稿状态问题已经有7人回复
毕业后当辅导员了，天天各种学生超烦已经有4人回复
求助文献已经有3人回复

» 本主题相关价值贴推荐，对您同样有帮助:

如何使用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人回复

1楼 2012-05-07 17:19:44

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

521438

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 82.9
帖子: 19
在线: 13.2小时
虫号: 1801287
注册: 2012-05-07
性别: MM
专业: 电磁场与电路

引用回帖:

7楼: Originally posted by jiakai003 at 2012-05-08 21:05:46:
弱弱问下，I1 I2没有定义？

其实我是把自己的方程简化了下，我这里有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=具体的数值就行。看看你会么？