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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 物理 » 统计、非线性&软物质 » 【资源】用matlab求解超越方程
51/1返回列表
查看: 2665  |  回复: 13
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【奖励】 本帖被评价6次，作者zxczxc0417增加金币 4.5 个
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zxczxc0417

木虫 (正式写手)

[资源] 【资源】用matlab求解超越方程



function matlabsky

% By LaterComer of MATLAB技术论坛

% See also http://www.matlabsky.com

% Contact me matlabsky@gmail.com

% Modifid 2010-05-14 11:02:15

%

phi=pi/3; % 根据需要修改

[e,fval]=fzero(@(e)objfun(e,phi),rand)



function f=objfun(e,phi)

[K,E] = ellipke(e); % 完整椭圆积分公式

f=(e^2*E-2*(1-e^2)*K)/(2-e^2)/E-cos(phi); % 目标函数
回复此楼

» 收录本帖的淘帖专辑推荐

材料+新能源+催化+外语 condensed matter physics

» 猜你喜欢

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

physicsWL

新虫 (正式写手)
phi=pi/3;
e=linspace(0,1-0.001,1000);
[K,E] = ellipke(e);
F=(e.^2.*E-2*(1-e.^2).*K)./(2-e.^2)./E-cos(phi);
F=abs(F);
[value,index]=min(F);
e(index)
赞一下 回复此楼
7楼2010-05-28 10:29:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答

zxczxc0417

木虫 (正式写手)
e =



    0.9657





fval =



  5.5511e-016
赞一下 回复此楼
2楼2010-05-27 03:05:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ylzhang508

金虫 (小有名气)

★★★ 三星级,支持鼓励

谢谢啊 这个fzero能解复数的方程吗

[ Last edited by ylzhang508 on 2010-5-27 at 16:40 ]
赞一下 回复此楼
4楼2010-05-27 16:32:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tower598

铁杆木虫 (职业作家)

★★★ 三星级,支持鼓励

怎么能求出一系列的解那？
赞一下 回复此楼
5楼2010-05-27 16:41:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 23jxep3nCNZb 2026-09-14 4/200 2026-09-17 13:51 by ObVzqQQh1rrE
[找工作] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-17 13:39 by ObVzqQQh1rrE
[考研] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +6 23jxep3nCNZb 2026-09-14 6/300 2026-09-17 13:29 by ObVzqQQh1rrE
[硕博家园] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-17 13:15 by ObVzqQQh1rrE
[找工作] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 s3fFTmArrBt6 2026-09-14 5/250 2026-09-17 12:39 by ObVzqQQh1rrE
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 mibUvS8DDCwf 2026-09-15 4/200 2026-09-17 06:26 by k06BKNblrGRH
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-17 05:16 by k06BKNblrGRH
[考研] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 Aj1rhIDL5ixY 2026-09-14 4/200 2026-09-17 04:52 by k06BKNblrGRH
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-17 04:28 by k06BKNblrGRH
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-17 03:52 by k06BKNblrGRH
[公派出国] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +6 s3fFTmArrBt6 2026-09-14 6/300 2026-09-17 02:38 by Tql5LQhh5rLK
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-17 02:38 by Tql5LQhh5rLK
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-17 02:02 by Tql5LQhh5rLK
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 LwdutQ8HoqWP 2026-09-13 7/350 2026-09-17 02:02 by Tql5LQhh5rLK
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +6 LwdutQ8HoqWP 2026-09-13 6/300 2026-09-17 01:50 by Tql5LQhh5rLK
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +7 QUjhNVAcOSff 2026-09-13 8/400 2026-09-16 18:37 by rRvxEj3zorCy
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-16 03:14 by DgNGHc3h5tPl
[考博] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +6 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:47 by CgyNCDVNhGVg
[教师之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:26 by C79jjtjAKjEn
[教师之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:46 by C79jjtjAKjEn
信息提示
关闭
请填处理意见
关闭