| 查看: 1064 | 回复: 6 | |||
[求助]
Matlab 求解非线性方程组 已有1人参与
|
|
使用matlab求解带有椭圆积分的非线性方程组，总是有错误，请各位大牛帮忙看看如何解决。 其中 B，f，l, h都是常数。 m文件和运行结果见附件。谢谢  1.png  2.jpg |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : 2.jpg
2015-05-07 15:01:40, 59.48 K
» 猜你喜欢
国自然面上和省基金B类撒花
已经有17人回复
-
有没有学校收留
已经有3人回复
-
312求调剂
已经有3人回复
-
华师大读博
已经有5人回复
-
又一批高校组建人工智能学院 师资行吗 不是骗人吗
已经有5人回复
-
急需审稿人！！！
已经有3人回复
-
申博/考博
已经有8人回复
-
295分求调剂
已经有6人回复
-
085600材料与化工调剂
已经有6人回复
-
有没有接收比较快的sci期刊呀，最好在一个月之内的，研三孩子求毕业
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Matlab如何解同余方程组
已经有3人回复
- Matlab用fsolve求解非线性方程，卡了三四天了跪求帮助！
已经有35人回复
- MATLAB解非线性方程组
已经有9人回复
- 请教高手，如何解非线性方程组！！！
已经有13人回复
- MATLAB求解五元非线性方程组，在线等
已经有7人回复
- 非线性二阶微分方程组求解，matlab
已经有12人回复
- matlab解方程组求助攻
已经有4人回复
- matlab 求解线性方程组Ax=b
已经有10人回复
- matlab 非线性微分方程求解
已经有3人回复
- 求助：非线性方程组的求解（郁闷中）
已经有26人回复
- 用matlab解非线性方程组
已经有4人回复
- 非常简单！！MATLAB求解非线性一元方程
已经有5人回复
- Matlab求解多元多次方程组问题
已经有7人回复
- 用matlab求解非线性方程组说无解，一定是方程组本身无解，还是有可能程序有问题呢？
已经有11人回复
- matlab如何求解一个非线性微分方程组
已经有8人回复
- MATLAB求解非线性方程组
已经有5人回复
- 求Matlab解非线性方程代码
已经有12人回复
- matlab求解非线性方程组，错误提示怎么解决
已经有5人回复
- 求助 MATLAB解方程组-fslove
已经有7人回复
- 求助matlab---fsolve解非线性方程组
已经有6人回复
- matlab求解非线性方程组
已经有16人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- 请教Matlab高手解多元非线性方程组中的参数
已经有6人回复
- 【求助】非线性方程组求解
已经有3人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
wurongjun
专家顾问 (职业作家)
-
专家经验: +831 - 数学EPI: 9
- 应助: 791 (博后)
- 贵宾: 0.308
- 金币: 24609
- 散金: 310
- 红花: 75
- 帖子: 3004
- 在线: 881.4小时
- 虫号: 1368482
- 注册: 2011-08-14
- 性别: GG
- 专业: 计算数学与科学工程计算
- 管辖: 数学
2楼2015-05-07 18:32:34
3楼2015-05-07 21:42:56
wurongjun
专家顾问 (职业作家)
-
专家经验: +831 - 数学EPI: 9
- 应助: 791 (博后)
- 贵宾: 0.308
- 金币: 24609
- 散金: 310
- 红花: 75
- 帖子: 3004
- 在线: 881.4小时
- 虫号: 1368482
- 注册: 2011-08-14
- 性别: GG
- 专业: 计算数学与科学工程计算
- 管辖: 数学
4楼2015-05-07 22:10:43
|
fun.m文件： function f=fun(x) B=2; ff=0.5; l=20; hh=1; N=(x(2)/2-ff*sin(x(1)))/cos(x(1)); t=N*sin(x(1))-ff*cos(x(1)); hh=sqrt((x(2)/2)^2+t^2); k=sqrt((ff+hh)/(2*hh)); e=sqrt((hh-t)/(hh+t)); v=asin(sqrt(((t+ff)/(hh+t))/k^2)); f(1)=sqrt(B/hh^5)*(-t*hh*mfun('EllipticF',v,k)+2*t*hh*mfun('EllipiticE',v,k)-sqrt(2*h*(ff+t))*(x(2)/2+t*e))-l; f(2)=sqrt(B/hh^5)*(x(2)*hh*mfun('EllipticF',v,k)-2*x(2)*hh*mfun('EllipticE',v,k)+sqrt(2*hh*(ff+t))*(2*t+x(2)*e))-2*h; 命令窗口： [x,f,h]=fsolve(@fun,[0.1,1]) 非常感谢！ |
5楼2015-05-08 10:16:15
wurongjun
专家顾问 (职业作家)
-
专家经验: +831 - 数学EPI: 9
- 应助: 791 (博后)
- 贵宾: 0.308
- 金币: 24609
- 散金: 310
- 红花: 75
- 帖子: 3004
- 在线: 881.4小时
- 虫号: 1368482
- 注册: 2011-08-14
- 性别: GG
- 专业: 计算数学与科学工程计算
- 管辖: 数学
6楼2015-05-08 10:56:04
7楼2015-05-08 12:15:27
5