24小时热门版块排行榜

>论坛更新日志 (3248)
>虫友互识 (354)
>文献求助 (203)
>休闲灌水 (180)
>导师招生 (154)
>考博 (78)
>博后之家 (50)
>硕博家园 (48)
>考研 (48)
>论文投稿 (40)
>基金申请 (34)
>论文道贺祈福 (31)
>公派出国 (30)
>找工作 (26)
>教师之家 (25)
>招聘信息布告栏 (23)

返回列表

不曾走远

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1
帖子: 8
在线: 3.2小时
虫号: 2062635
注册: 2012-10-15
专业: 制造系统与自动化

[求助] matlab 问题求助，真心求大侠相助！

function [r,n]=half(fun,a,b)
n=0;
while abs(fun(a)-fun(b))>1e-5
c=0.5*(a+b);
if fun(a)*fun(c)<0
b=c;
elseif fun(b)*fun(c)<0
a=c;
end
n=n+1;
end
r=c;
clear;
clc;
fun=inline('x^2-x-1','x');
[r,n]=half(fun,0,10)

在M文件里运行不了，请求相助！谢谢！

回复此楼

» 猜你喜欢

为什么中国大学教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？已经有3人回复
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急已经有4人回复
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急已经有3人回复
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急已经有4人回复
“人文社科而论，许多学术研究还没有达到民国时期的水平” 已经有5人回复
过年走亲戚时感受到了所开私家车的鄙视链已经有11人回复
什么是人一生最重要的？已经有4人回复
版面费该交吗已经有3人回复
今年春晚有几个节目很不错，点赞！已经有12人回复
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低已经有12人回复

» 本主题相关价值贴推荐，对您同样有帮助:

1楼 2012-10-21 12:54:30

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

csl_2012

捐助贵宾 (初入文坛)

应助: 6 (幼儿园)
金币: 838.3
散金: 300
红花: 1
帖子: 50
在线: 113小时
虫号: 1734989
注册: 2012-04-03
性别: GG
专业: 动力学与控制

【答案】应助回帖

★
感谢参与，应助指数 +1
jjdg: 金币+1, 辛苦了 2012-10-22 00:54:24

我试了一下，程序没有错误。计算结果：r = 1.618034839630127,n = 22
是不是你对MATLAB还不太熟悉，m文件名要和函数名相同，而且上面代码中从clear开始的代码要出现在命令窗口中，不能在m文件里。
另外求方程根，二分法效率较低，可以采用如牛顿迭代法等，详见数值分析。
希望对你有用。

赞一下

回复此楼

2楼2012-10-21 23:06:51

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主不曾走远的主题更新

返回列表