应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 求助！matlab用fsolve函数求解非线性方程组的问题！
201/2返回列表12下一页
查看: 3011  |  回复: 19
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

521438

铜虫 (初入文坛)

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-05-07 17:19:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

521438

铜虫 (初入文坛)
送鲜花一朵
有人看到我的帖子么？大家知道么？帮我看一下。十分感谢！
赞一下 回复此楼
2楼2012-05-07 20:49:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jinpeng_6118

木虫 (正式写手)
LZ你好，请问下 MatLab  5.2 这个软件哪里能有下载的？谢谢
赞一下 回复此楼
人生百态原为海，看破红尘方为岸！
3楼2012-05-08 16:01:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

521438

铜虫 (初入文坛)
引用回帖:
3楼: Originally posted by jinpeng_6118 at 2012-05-08 16:01:50:
LZ你好，请问下 MatLab  5.2 这个软件哪里能有下载的？谢谢

为什么要5.2的呢？我用的是7.0的！需要的话，可以发给你！我也不知道哪里有5.2的，不好意思啊！
赞一下 回复此楼
4楼2012-05-08 17:15:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jinpeng_6118

木虫 (正式写手)
引用回帖:
4楼: Originally posted by 521438 at 2012-05-08 17:15:11:
为什么要5.2的呢？我用的是7.0的！需要的话，可以发给你！我也不知道哪里有5.2的，不好意思啊！

好的 ，谢谢你了，麻烦你把软件发到我邮箱吧：jinpeng_6118@126.com
赞一下 回复此楼
人生百态原为海，看破红尘方为岸！
5楼2012-05-08 17:40:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

521438

铜虫 (初入文坛)
引用回帖:
5楼: Originally posted by jinpeng_6118 at 2012-05-08 17:40:14:
好的 ，谢谢你了，麻烦你把软件发到我邮箱吧：jinpeng_6118@126.com

已经发给你了，你查收一下吧!
赞一下 回复此楼
6楼2012-05-08 18:54:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiakai003

禁虫 (小有名气)
★
xiegangmai: 金币+1, 谢谢参与！ 2012-05-08 23:09:32
本帖内容被屏蔽
7楼2012-05-08 21:05:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jinpeng_6118

木虫 (正式写手)
引用回帖:
6楼: Originally posted by 521438 at 2012-05-08 18:54:23:
已经发给你了，你查收一下吧!

已经收到了，非常感谢你的帮助
赞一下 回复此楼
人生百态原为海，看破红尘方为岸！
8楼2012-05-08 21:46:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

521438

铜虫 (初入文坛)
引用回帖:
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=具体的数值就行。看看你会么？
赞一下 回复此楼
9楼2012-05-09 10:49:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiegangmai

版主 (职业作家)

我没头衔

优秀版主优秀版主优秀版主

【答案】应助回帖

★
感谢参与，应助指数 +1
臭水沟: 金币+1, 谢谢交流~~ 2012-05-09 15:27:07
你的函数里包含符号I1,I2,肯定解不出数值解的。
不用简化，把数学模型贴出来看看。

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
明德厚学、求是创新
10楼2012-05-09 12:31:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 521438 的主题更新
201/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭