二分法matlab程序

functiontestclearclcf0=@(x)x-cos(x);f2fen(f0,0,pi/2)%下面是二分法的函数部分functionx2=f2fen(f,a,b)x1=[a;b];p1=f

请楼主参考下面的代码.>>f=inline('exp(-x)-sin(pi*x/2)');%随便给个函数.自己输入吧.a=0;b=1;%输入区间.也是按照题目要求输入.err=(b-a)/2;k=0;

%f是要求根的函数,f=inline('','');ya=f(a)可能输入的f不能用f(a)求函数值,换成feval(f,a);试试

f=inline('x^2-x-2');>>[c,err,yc]=bisect(f,0,3,0.01)c=2.0010err=0.0059yc=0.0029-----------%使用二分法求解上面超

functiontesta=1;b=2;tol=1e-4;N=10000;k=0;fa=f(a);fork=1:Np=(a+b)/2;fp=f(p);if(fp==0||(b-a)/2再问：粘贴了您的

在Matlab下输入：edit,然后将下面两行百分号之间的内容,复制进去,保存%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%functiony=zh

兄弟,西理工的?

这是源代码:在matlab中保存为:bisection.mfunction rtn=bisection(fx,xa,xb,n,delta)% 二分法解方程% fx是由方程

可以,只要在二分法外面,加一个循环就可以了.比如,以i为变量从1到1000,只要f(x(i-1))*f(xi)

这个有例子区间二分法：与对分查找法相同1区间二分法求出的仅仅是方程的一个单根,如果方程有重根或者多个根时,在做区间二分法时就会出现分叉,这样方程有几个根,就会产生几个实数序列,每一个实数序列的极限便是

二分法数学方面：一般地,对于函数f(x),如果存在实数c,当x=c时f(c)=0,那么把x=c叫做函数f(x)的零点.解方程即要求f(x)的所有零点.先找到a、b,使f(a),f(b)异号,说明在区间

主体用c写的#includeusingnamespacestd;intg（inta）{intx；x=a*a*a-a-1if（x>0）{return(1);}elsereturn(0);}voidmai

#include"stdio.h"#defineMAX100typedefstructnode{intkey;}NODE;intbinsearch(NODEr[MAX],intk,intn){in

建议楼主遇到关于matlab的问题就到mathworks网站的fileexchange里找.下面是二分法的函数文件,你直接设置输入参数就可以了function[c,err,yc]=bisect(f,a

程序如下：clear,clc;a=0;%a=input('inputa:');b=1;%b=input('inputb:');err=10^-5;y1=a*exp(a)

那调用两次不就ok了.区间你是用a、b定义的嘛.

建议你在第三行插入一个end,以完成函数的定义而且函数的定义只有在m文件里面才能做到.（这很重要!）还有,好像没有与while相应的end对了,你的函数在[-pi/2,pi/2]里面只有一个解0,所以

functiontestclearclc%实验方程：3*x.^2+x+2*exp(x)=0%原函数f=@(x)3*x.^2+x-2*exp(x);%导函数df=@(x)6*x+1-2*exp(x);%

A=[123-4-2;-3-4-12135;2100-310;4149-137];%A为系数矩阵fori=2:4forj=i:5fork=1:(i-1)A(i,j)=A(i,j)-A(i,k)*A(k

%二分法symsx;fun=input('(输入函数形式)fun=');a=input('（输入二分法下限）a=');b=input('（输入二分法上限）b=');%二分法求根%f=inline(x^