24小时热门版块排行榜

>论坛更新日志 (3358)
>虫友互识 (271)
>文献求助 (251)
>导师招生 (165)
>硕博家园 (105)
>招聘信息布告栏 (76)
>休闲灌水 (74)
>考博 (72)
>论文道贺祈福 (61)
>博后之家 (51)
>考研 (46)
>教师之家 (44)
>公派出国 (43)
>论文投稿 (42)
>基金申请 (40)
>找工作 (28)

返回列表

wtisunday

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 1875.3
散金: 30
红花: 1
帖子: 936
在线: 568.7小时
虫号: 444484
注册: 2007-10-27
性别: GG
专业: 高分子材料的加工与成型

[求助] m文件调用问题

大侠们，我依据教科书编了一个二分法的程序（标准的啊），然后调用的时候总是出现如下情况，还请赐教！本人乃初学菜鸟级的
f=inline('x^3-x-1');
>> bisect(f,1,2)
??? Error using ==> feval
Undefined function or method 'f' for input arguments of type 'double'.

Error in ==> bisect at 10
ya=feval('f',a)

回复此楼

» 猜你喜欢

第一性原理计算方向2026级博士申请 PRB*1，四级484 已经有1人回复
求助VISSIM破解版软件已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有232人回复
求2026年在台湾举行的物理和材料领域国际学术会议信息已经有0人回复
求国际会议网站已经有1人回复
求取一些关于纳米材料和纳米技术相关的英文PPT。已经有0人回复
【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有17人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复

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

如何用MATLAB 实现化学反应方程式（写程序代码）? 已经有20人回复
matlab中M文件的调用已经有12人回复
论文引用的参考文献问题已经有9人回复
归一化问题已经有11人回复
请教一个调用数据文件的程序问题，谢谢已经有4人回复
关于M文件的一个问题已经有4人回复
【讨论】关于一个m文件中函数调用及画图问题已经有7人回复
【求助】matlab 中，几个m文件调用已经有14人回复
【求助】matlab生成旋转体？已经有9人回复
【求助】CuBTC MOF + CO2 MuSiC 问题已经有101人回复
【问题已经解决】怎样使用endnote按照文献引用的先后顺序自动编号已经有7人回复
【求助】EndNote引用出现“文件未找到”的问题已经有8人回复

1楼2011-12-01 09:51:19

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

dbb627

荣誉版主 (著名写手)

专家经验: +4
计算强帖: 12
应助: 289 (大学生)
贵宾: 0.589
金币: 24640.4
散金: 551
红花: 61
沙发: 1
帖子: 1246
在线: 1794.8小时
虫号: 149791
注册: 2005-12-29
性别: GG
专业: 污染控制化学
管辖: 计算模拟

【答案】应助回帖

wtisunday(金币+10): 版主?涓Π� 2011-12-01 12:04:32

CODE:

function [approx_root,n]=bisect(f,a,b)

 %a,b为根的区间

fa=feval(f,a)

fb=feval(f,b)

 n=0;

 while(abs(b-a)> 0.0001)

     c=(a+b)/2;

     approx_root=c;

         fc=feval(f,c)

     [a,c,b;fa,fc,fb];

     n=n+1;

if ( sign(fb) * sign(fc)<=0)

            a=c;

            fa=fc;

          else

            b=c;

            fb=fc;

          end

        end

>> f=inline('x^3-x-1');
>> bisect(f,1,2)

ans =

1.3248

>> f(1.3248)

ans =

3.4991e-004

赞一下

回复此楼

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

2楼2011-12-01 10:54:52

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

wtisunday

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 1875.3
散金: 30
红花: 1
帖子: 936
在线: 568.7小时
虫号: 444484
注册: 2007-10-27
性别: GG
专业: 高分子材料的加工与成型

function [c,err,yc]=bisect(f,a,b,delta)
% f是所要求解的函数
% a和b分别为有根区间的左右限
% delta是允许的误差界
% c为所求近似值
% yc为函数f在c上的函数值
% err为c的误差估计
if nargin<4 delta=1e-5;
end
ya=feval('f',a)
yb=feval('f',b)
if yb==0,c=b,return,
end
if ya*yb>0,
disp('(a,b)不是有根区间');
return
end
max1=1+round((log(b-a)-log(delta))/log(2));
for k=1:max1
c=(a+b)/2;
yc=feval('f',c);
if yc==0 a=c;b=c;break,
elseif yb*yc>0
b=c;yb=yc;
else
a=c;ya=yc;
end
if (b-a) end
k,c=(a+b)/2,err=abs(b-a),yc=feval('f',c)

我这个哪出错了吗

赞一下

回复此楼

3楼2011-12-01 11:23:23

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