24小时热门版块排行榜

返回列表

不曾走远

新虫 (初入文坛)

应助: 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文件里运行不了，请求相助！谢谢！

回复此楼

» 猜你喜欢

26/27申博自荐已经有10人回复
东北林业大学材料科学与工程学院“一流”A+学科国家级人才团队课题组招收2026级博士生已经有3人回复
医学类期刊求推荐已经有5人回复
生活琐事由它去已经有4人回复
提交了我也来说说感想已经有12人回复
青B发送上会通知了吗已经有9人回复
西安交大新媒学院副院长用撤稿论文结题已经有6人回复
论文撤稿了已经有8人回复
化学专业申博已经有4人回复
某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路！已经有5人回复

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

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的回帖

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

返回列表