24小时热门版块排行榜

>论坛更新日志 (3346)
>导师招生 (1255)
>虫友互识 (624)
>休闲灌水 (518)
>考博 (235)
>招聘信息布告栏 (88)
>论文投稿 (48)
>硕博家园 (46)
>教师之家 (25)
>考研 (17)
>学术会议 (13)
>公派出国 (13)
>文献求助 (10)
>博后之家 (3)
>科研资料 (3)
>找工作 (3)

返回列表

wtisunday

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 1698.3
散金: 30
红花: 1
帖子: 939
在线: 570.5小时
虫号: 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)

回复此楼

» 猜你喜欢

深圳大学应用技术学院招聘凝聚态物理博士后已经有124人回复
求助中医院临床中医中标的青年标书学习，谢谢已经有5人回复
物理学I论文润色/翻译怎么收费? 已经有79人回复
浙大物院表面物理课题组长期招收普博及直博学生已经有259人回复
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry 已经有47人回复
25届硕士申请，材料计算方向已经有2人回复
课题组招1-2名硕士已经有1人回复
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上已经有0人回复
斯坦福大学与Elsevier世界前2%学者新排名查询已经有0人回复
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！已经有19人回复

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

如何用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 (幼儿园)
金币: 1698.3
散金: 30
红花: 1
帖子: 939
在线: 570.5小时
虫号: 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的回帖