| 查看: 1608 | 回复: 6 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
monixia木虫 (小有名气)
|
[求助]
求解决matlab解非线性方程组出现的问题 已有1人参与
|
|
|
各位matlab大神，我计算的方程组如下， >> m=0.00040; f=@(x)([3.14159*m^3+1.9635e-7*x(1)-0.002*3.14159*0.0005^2; 6*(m+x(2))*1.9635e-7-3.14159*m^3+1.9635e-7*x(3)+1.9635e-7*x(6)*(m+x(2))/x(4)-0.002*3.14159*0.0005^2; 3.14159*m^3*x(4)/(m+x(2))+1.9635e-7*x(1)*x(5)/(x(1)+x(3))-1/6*1e-6; (6*1.9635e-7*(m+x(2))-3.14159*m^3)*x(4)/(m+x(2))+1.9635e-7*x(3)*x(5)/(x(1)+x(3))+1.9635e-7*x(6)-1/6*1e-6; fix(0.05/(x(1)+x(3)))*(32*0.0010087*x(5)/0.0005^2*x(3)+7.16*0.07275*(3*(1.0087*x(5)/72.75)^(2/3)/0.0005/(1+3.34*(1.0087*x(5)/72.75)^(2/3))))-6455.68*x(6) 0.05/(m+x(2))*(32*0.0010087*x(4)/0.0005^2*x(2)+7.16*0.07275*(3*(1.0087*x(4)/72.75)^(2/3)/0.0005/(1+3.34*(1.0087*x(4)/72.75)^(2/3))))-6455.68*x(6)]); options = optimset('MaxIter',1e30,'MaxFunEvals',1e50,'TolFun',1e-50); x=fsolve(f,[0.005 0.0005 0.001 0.15 0.9 0.3],options) No solution found. fsolve stopped because the relative size of the current step is less than the default value of the step size tolerance squared, but the vector of function values is not near zero as measured by the selected value of the function tolerance. <stopping criteria details> x = 0.0070 -0.0001 0.0030 0.0344 1.0563 0.4352 fsolve stopped because the relative norm of the current step, 2.848645e-15, is less than max(options.TolX^2,eps) = 1.000000e-12. However, the sum of squared function values, r = 3.833839e-15, exceeds sqrt(options.TolFun) = 1.000000e-25. Optimization Metric Options relative norm(step) = 2.85e-15 max(TolX^2,eps) = 1e-12 (default) r = 3.83e-15 sqrt(TolFun) = 1.0e-25 (selected) 如上，初值已经取的很接近实验值了，还是没有满意的结果。将matlab算的带入六个方程， ans = 4.7168e-12 ans = -8.4160e-11 ans = 1.5716e-09 ans = -1.5222e-09 ans = 15.7955 ans = 232.7567 后两项收敛不是很好。我已经尝试把optimset调到很小了，可是还是没有满意的结果。 求助各位大牛，这个该怎么解？还是方程需要修改？ 如果圆满解决，可以加金币~ |
回复此楼» 猜你喜欢
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有80人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab解二元方程组，错误提示，求解决方法
已经有11人回复
- MATLAB解非线性方程组
已经有9人回复
- matlab 解非线性偏微分方程组
已经有6人回复
- MATLAB求解五元非线性方程组，在线等
已经有7人回复
- 非线性二阶微分方程组求解，matlab
已经有12人回复
- Matlab中ode45出现大误差，关于用ODE45解偏微分方程的疑问，有待解决，大家共同讨论！
已经有4人回复
- Matlab用牛顿法求解非线性方程组问题
已经有8人回复
- 用matlab求解非线性方程组说无解，一定是方程组本身无解，还是有可能程序有问题呢？
已经有11人回复
- matlab如何求解一个非线性微分方程组
已经有8人回复
- sos!!1stopt解非线性方程组，一直出现常数定义错误的问题！求大神帮助！
已经有12人回复
- matlab中利用LM算法解非线性方程组的问题
已经有5人回复
- MATLAB求解非线性方程组
已经有5人回复
- 关于matlab求解非线性指数方程组出现问题
已经有3人回复
- 【求助】一些A、B类错误求解
已经有15人回复
- 求Matlab高手解决线性方程组的迭代求解问题
已经有17人回复
- matlab求解非线性方程组，错误提示怎么解决
已经有5人回复
- 求助matlab---fsolve解非线性方程组
已经有6人回复
- 求助！matlab用fsolve函数求解非线性方程组的问题！
已经有19人回复
- matlab报错，求解
已经有5人回复
- matlab求解非线性方程组
已经有16人回复
- matlab问题求助
已经有16人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- 请教Matlab高手解多元非线性方程组中的参数
已经有6人回复
- matlab的fsove 命令求解非线性方程组
已经有6人回复
- 【求助】matlab求解非线性方程组，并画图处理。要求y,z是实数解！
已经有18人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
monixia
木虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 2535.6
- 红花: 3
- 帖子: 80
- 在线: 24.2小时
- 虫号: 1817862
- 注册: 2012-05-15
- 性别: GG
- 专业: 化学反应工程
3楼2014-11-05 19:09:09
monixia
木虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 2535.6
- 红花: 3
- 帖子: 80
- 在线: 24.2小时
- 虫号: 1817862
- 注册: 2012-05-15
- 性别: GG
- 专业: 化学反应工程
6楼2014-11-08 10:42:38