matlab解多元二次方程组,我算出来的结果误差好大?怎么办?急,跪求! 200
functionF=myfun(x)F=[x(1)-x(7)-0.25*x(4)+0.003*x(3);x(2)-x(6)+0.23*x(3)+0.003*x(4);1....
function F = myfun(x)
F = [x(1)-x(7)-0.25*x(4)+0.003*x(3);
x(2)-x(6)+0.23*x(3)+0.003*x(4);
1.5-x(5)-1.51*x(3);
x(6)-x(5)+0.07*x(3);
x(7)-1.51*x(4);
x(6)*x(4)+x(7)*x(3)+0.02*x(3)*x(4)-2;
x(1)-0.02*x(3)+0.4*x(4)-x(8);
x(2)-0.02*x(4)-0.4*x(3)-x(9);
x(1)*x(3)+x(2)*x(4)-0.1;
x(2)*x(3)-x(1)*x(4)-0.15;
x(10)-x(14)*x(13);
x(14)-10000-0.2744*x(13);
10000*x(13)+x(11);
0.00000001*(x(3)*x(8)+x(4)*x(9))-x(10);
0.00000001*(x(5)*x(7)+x(6)*x(8))-x(11);
0.02*x(17)-0.4*x(18)+x(15)-1;
0.02*x(18)+0.4*x(17)+x(16);
x(18)-0.1];
x0 = [1;3;3;1;5;1;8;1;1;1;1;1;1;10000;1;1;1;1];
options=optimset('TolFun',1e-30);
[x,fval]=fsolve(@myfun,x0,options)
初值是自己随便给的,有没有方法算的精确一点 这个算出来的答案比较无语
fsolve stopped because it exceeded the function evaluation limit,
options.MaxFunEvals = 1800 (the default value).
是不是因为这个提示,我应该怎么修改? 展开
F = [x(1)-x(7)-0.25*x(4)+0.003*x(3);
x(2)-x(6)+0.23*x(3)+0.003*x(4);
1.5-x(5)-1.51*x(3);
x(6)-x(5)+0.07*x(3);
x(7)-1.51*x(4);
x(6)*x(4)+x(7)*x(3)+0.02*x(3)*x(4)-2;
x(1)-0.02*x(3)+0.4*x(4)-x(8);
x(2)-0.02*x(4)-0.4*x(3)-x(9);
x(1)*x(3)+x(2)*x(4)-0.1;
x(2)*x(3)-x(1)*x(4)-0.15;
x(10)-x(14)*x(13);
x(14)-10000-0.2744*x(13);
10000*x(13)+x(11);
0.00000001*(x(3)*x(8)+x(4)*x(9))-x(10);
0.00000001*(x(5)*x(7)+x(6)*x(8))-x(11);
0.02*x(17)-0.4*x(18)+x(15)-1;
0.02*x(18)+0.4*x(17)+x(16);
x(18)-0.1];
x0 = [1;3;3;1;5;1;8;1;1;1;1;1;1;10000;1;1;1;1];
options=optimset('TolFun',1e-30);
[x,fval]=fsolve(@myfun,x0,options)
初值是自己随便给的,有没有方法算的精确一点 这个算出来的答案比较无语
fsolve stopped because it exceeded the function evaluation limit,
options.MaxFunEvals = 1800 (the default value).
是不是因为这个提示,我应该怎么修改? 展开
2个回答
展开全部
可以去matlab贴吧看一看,大神很多的
已赞过
已踩过<
评论
收起
你对这个回答的评价是?
富港检测东莞有限公司
2024-12-25 广告
2024-12-25 广告
正弦振动多用于找出产品设计或包装设计的脆弱点。看在哪一个具体频率点响应最大(共振点);做正弦振动找富港,富港工业检测技术有限公司是一家专业的第三方检测机构,拥有完善的质量管理体系,先进的检测设备,优秀的技术人才;已取得CNAS、CMA、IS...
点击进入详情页
本回答由富港检测东莞有限公司提供
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询