24小时热门版块排行榜    

查看: 1062  |  回复: 2

jackie309

银虫 (正式写手)

[求助] 请教函数作为一个参数的用法,谢谢

请教一个问题:
如何调用这样的一个m函数,谢谢

bisection(f,a,b)
function p = bisection(f,a,b)

% provide the equation you want to solve with R.H.S = 0 form.
% Write the L.H.S by using inline function
% Give initial guesses.
% Solves it by method of bisection.
% A very simple code. But may come handy

if f(a)*f(b)>0
    disp('Wrong choice bro')
else
    p = (a + b)/2;
    err = abs(f(p));
    while err > 1e-7
   if f(a)*f(p)<0
       b = p;
   else
       a = p;         
   end
    p = (a + b)/2;
   err = abs(f(p));
    end
end

比如我现在需要的f函数是 y = 3*x^3 + x*x -2


我的做法是在一m文件,或者command windows:

ff=@(x)  x.*x+3.*x+2;
p = bisection(ff,-3,4);

运行,提示的出错信息:
Wrong choice bro
Error in bisection (line 9)
if f(a)*f(b)>0

Output argument "p" (and maybe others) not assigned during call to "D:\Program
Files\MATLAB\R2013a\bisection.m>bisection".

Error in f (line 4)
p = bisection(ff,-3,4);
回复此楼

» 猜你喜欢

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

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

月只蓝

主管区长 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
用inline函数:

ff=inline(' 3*x.^3 + x.*x -2','x');
bisection(ff,-3,4)

ans =

    0.7754
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
2楼2013-08-21 23:39:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军

【答案】应助回帖

感谢参与,应助指数 +1
提示说明是返回参数p没有指定,
在disp('Wrong choice bro')下面加一行p = 0;
matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com
3楼2013-08-21 23:49:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jackie309 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 时间戳又变了8-15 +7 archvillain 2026-08-15 11/550 2026-08-15 15:06 by Yeuchan
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 k0dTPqJtl0jt 2026-08-14 5/250 2026-08-15 14:58 by QXtWNz7PI7wZ
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +26 医学老男孩 2026-08-13 57/2850 2026-08-15 14:25 by JChen2024
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +7 majunge000 2026-08-11 9/450 2026-08-15 13:34 by kiz2011
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 7lpolszZVXgi 2026-08-14 8/400 2026-08-15 11:53 by KxMI1BYBxWX1
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 HFw0lei2R37i 2026-08-14 5/250 2026-08-15 11:41 by KxMI1BYBxWX1
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 HFw0lei2R37i 2026-08-14 9/450 2026-08-15 11:33 by KxMI1BYBxWX1
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 7lpolszZVXgi 2026-08-14 6/300 2026-08-15 09:28 by sunzitan
[基金申请] filecode +7 cratir 2026-08-14 11/550 2026-08-15 06:19 by 学员8dgXkO
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 k0dTPqJtl0jt 2026-08-14 4/200 2026-08-15 02:21 by 4wMiSEwB6436
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23 Lanmanbaby 2026-08-09 37/1850 2026-08-14 13:40 by Equinoxhua
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
信息提示
请填处理意见