matlab二分法 matlab二分法代码

matlab求自变量值的问题

3.0000

因为函数是单调的，可以利用二分法来做，给定一个x的初始值，求出F(x)，比较这个F(x)与已知F(x)数%Bisection值的大小， 然后取x/2，再次求出F(x)并比较大小，这样就可以确定出已知F(x)的自变量x取值范围是在0-x/2还是在x/2-x之间，然后再取x/4或者3x/4，求取F(x)值与已知F值比较，知道达到你所要求的精度就可以结束了，这样这个x值实际上是无限的靠近真实解。

matlab二分法 matlab二分法代码

matlab二分法 matlab二分法代码

不知道看明白了没有。

一个二分法的matlab程序，保存后提示Error using mbisec (line 6) Not enough input arguments.

end

你的f函数没有定义啊 建一个函数先

而且你的while abs(feval(f,x))>tol|(b-a)>tol应该改成while abs(feval(f,x))>tol&(b-a)>tol

下面是实验用的

%程序mbisec.m

function [x,k]=mbisec(f,a,b,tol)

%用途：用二分法求非线性方程f(x)=0有根区间[a,b]中的一个根

%格式：[x,k]=mbisec(f,a,b2.0000,tol) f为函数表达式，a,b为区间左右端点，tol为误限

x=(a+b)/2.0;k=0;

while abs(feval(f,x))>tol&(b-a)>tol

b=x;

a=x;

x=(a+b)/2;k=k+1;

f=x^2-4x;

你可以 输入mbisec('test',-2,2,0.01)（你原来的程序在次就出现了0点会无限循环）

证明方程x^3+x-4=0在区间[1,2]内有一个根。如果用二分法求它具有5位有效数字的根,需要二分多少次?用matlab

fb=fc;

%程x^3+x-4=0在区间[1,2]内有且一个根。

%以下是2分法的代码

a(1)=1;b(1)=2;

for i=1:inf

c(i)=(a(i)+b(i))/2;

if c(i)^3+c(i)-4>0

a(i+1)=a(i);

b(i+1)=c(i);

if c(i)^3+c(i)-4<0

a(i+1)=c(i);

b(i+1)=b(i);

if abs(c(i)^3+c(i)-4)<=0.00001 %如果误在10^(-5)以内则跳出循环，我认为控制误合理点

end % if c(i)10^3-floor(c(i)10^3)>0

end %这样做误很大，我试了一下，根为1.4375，此时f(x)=

root=c(i) break %如果只需要根达到5位有效数字的话就把上行代码改为 %0.4080,误很大

f=c(i)^3+c(i)-4

%i为二分的次数 对函数求导，则导函数为X^2+1恒大于0.所以函数在R上单调递增，将1带入函数等于-2，将2带入函数等于6大于0所以在1到2之间，函数必定有根

对函数求导，则导函数为X^2+1恒大于0.所以函数在R上单调递增，将1带入函数等于-2，将2带入函数等于6大于0所以在1到2之间，函数必定有根

帮忙看看二分法的matlab程序，为什么没到精度就输出结果了？？？matlab新手，希望解释清楚点，谢谢

既然是二分法，肯定要循环使用了额

if abs(b-a)>=delta这一句的if改成while

while abs(F(x1))>1e-4ya=yc;

这里面多了一个end，去掉

调用方法：

f=inli0.0020ne('x+1')

a=-2;b=1;delta=0.001;

[c,err,yc]=bisect(f,a,b,delta)

matlab如何解二元二次方程组？

题主给出的方程，用solve函数求解要花费多个小时，为了解决该问题，你可以考虑其他的解方程函数（如vpasolve,二分法，牛顿法等）。虽然solve和vpasolve是属于同一类函数，但vpasolve函数擅长求解方程的数值解，而solve函数擅长求解方程的解析解和数值解，所以vpasolve求解速度要高于solve函数。

用matlab解二元二次方程组，可以直接用solve（）函数来求解。例如：

function f=test(x)

>>[x,y] = solve('x^2y^2 - 2x - 1 = 0','x^2 - y^2 - 1 = 0')

运行结果

对于较复杂的方程，也可以用数值分析的方法，如fsolve（）函数、二分法，牛顿法，求出其数值解。

用MATLAB程序编程：分析方程f(x)=sinx-x/2=0正根的分布情况，并用二分法求正根近似值，使误不超过0.01.

to

clc;clear

a=0;b=1;

fa=1-a-sin(a);

fb=1-b-sin(b);

c=(a+b)/2;

fc=1-c-sin(c);

while abs(fc)>0.510^(-4)

c=(a+b)/2;

fc=1-c-sin(c);

if fbfc>0

b=c;

a=c;

fa=fthec;

format long

fx=fc,x=c

结果：

fx =

-2.414986223420179e-005

x =

0.510986328100

>> x=solve('1-x-sin(x)')

x =

输入自定义f(x)，用matlab编程f(x)=0的根的指令，[x,f）＝root b(fn,a,b),算法为二分法，要具体的编程！！

1.1.31000

function [x,f]=rootb(fn,a,b)

tol=1e-6; %精度

cntmax=1000;%迭代次数

cnt=0; %迭代次数

fa=fn(a);

fb=fn(b);

if (fafb>0)

error('a,b的函数值同号')

while cnt

if (b-a

break;

x=(a+b)/2;

f=fn(x);

%[x,f] %调试

if (faf<=0)

b=x;

fb=f;

a=x;

fa=f;

cnt=cnt+1;

%%%%以上是函数文件，可以这样试一下：

f=@(x)x^3-x-2;

[x,ff]=rootb(f,1,2)

你需要二分法的函数，你应该把rootb()函数传上来

matlab怎么解方程

for i=1:12000001

matlab怎么解方程方法如下：

一般用matlab求解方程有两种方法：用软件自带的solve（）函数、fsolve（）函数、roots（）函数等；利用数值分析的方法（如二分法，牛顿法等）自行编程其数值函数。

拓展资料：

方程（equation）是指含有未知数的等式。是表示两个数学式（如两个数、函数、量、运算）之间相等关系的一种等式，使等式成立的未知数的值称为“解”或“根”。求方程的解的过程称为“解方程”。

通过方程求解可以免去逆向思考的不易，直接正向列出含有欲求解的量的等式即可。方程具有多种形式，如一元一次方程、二元一次方程、一元二次方程等等，还可组成方程组求解多个未知数。

解分数方程的方法如下：

1、看等号两边是否可以直接计算。

2、如果两边不可以直接计算，就运用和积商的公式对方程进行变形。

3、对可以相加减的项进行通分。

4、两边同时除以一个不为零的数。

注意：

（1）都含有未知数的项才能相加减，或者都不含有未知数的项才能相加减。

（2）除以一个数等于乘以这个数的倒数。

乘法分配律的应用：

1、加法交换律：a+b=b+a加法结合律：(a+b)+c=a+(b+c)。

2、乘法交换律：a×b=b×a乘法结合律：(a×b)×c=a×(b×c)。

3、乘法分配律：（a+bif F((x0+x1)/2) >= 0）×c=a×c+b×c。

4、减法的性质：a－b－c=a－(b+c)。

5、除法的性质：a÷b÷c=a÷(b×c)。

（注意：去括号时，括号前面是减号的，去掉括号，括号里的每一项要变号，也就是括号里的加号要变减号，减号要变成加号。这是运用了减法的性质）

matlab怎么计算方程组

of

matlab中解方程组还是很方便的，例如，对于代数方程组Ax=b(A为系数矩阵，非奇异)的求解，MATLAB中有两种方法：

else

(1)x=inv(A)b — 采用求逆运算解方程组；

(2)x=AB — 采用左除运算解方程组

PS：使用左除的运算效率要比求逆矩阵的效率高很多~

例:

x1+2x2=8

2x1+3x2=13

>>A=[1,2;2,3];b=[8;13];

>>x=inv(A)b

x =

2.00

3.00

>>x=AB

x =

2.00

3.00；

即二元一次方程组的解x1和x2分别是2和3。

对于同学问到的用matlab解多次的方程组，有符号解法，方法是：先解出符号解,然后用vpa(F,n)求出n位有效数字的数值解.具体步骤如下：

syms x1 x2 x3 x4 x5 f1 f2 f3 f4 f5

f1=(x1+x2+x3+x4+x5+22.55x1x4+7.63x3x4+148.222x1^3x4+5.8610^13x1^12x4^7+67.03x1x4^2+212.17x1x4^6+338.08x1x5+4948932.908x1^3x5^2+3522485477x1^4x5^3+3.46x4x5+4.109x3x5+6.03x3x5^2+4.847x3^2x5+135498.74x1^3x2x4^3+7.5510^14x1^11x2x4^7)-1;

f2=0.5653(1/3x2+135498.74x1^3x2x4^3+7.5510^14x1^11x2x4^7)-0.2946(0.5x1+22.55x1x4+444.666x1^3x4+7.0310^14x1^12x4^7+67.03x1x4^2+212.17x1x4^6+338.08x1x5+14846798.72x1^3x5^2+1.4110^10x1^4x5^3+406496.24x1^3x2x4^3+8.30510^15x1^11x2x4^7);

f4=0.05(x4+22.55x1x4+7.63x3x4+148.222x1^3x4+4.10210^14x1^12x4^7+134.06x1x4^2+1273.02x1x4^6+3.46x4x5+406496.24x1^3x2x4^3+5.28510^15x1^11x2x4^7)-0.3456(0.5x3+7.63x3x4+4.109x3x5+6.03x3x5^2+9.694x3^2x5);

f5=0.3456(x5+338.08x1x5+9897865.816x1^3x5^2+1.05610^10x1^4x5^3+3.46x4x5+4.109x3x5+12.06x3x5^2+4.847x3^2x5)-0.05(x4+22.55x1x4+7.63x3x4+148.222x1^3x4+4.10210^14x1^12x4^7+134.06x1x4^2+1273.02x1x4^6+3.46x4x5+406496.24x1^3x2x4^3+5.28510^15x1^11x2x4^7);

[solx1 solx2 solx3 solx4 solx5]=solve(f1==0,f2==0,f3==0,f4==0,f5==0,x1,x2,x3,x4,x5);

matlab 的solve函数求解等式中的未知数，未知数形式为N^-小数，如下，求解多个小时，怎么办？

对于本例，只要vpasolve替代solve函数，就可以顺利得到其数值解。

。6.0000。。。。。

symif fafb>0,break,ends N

eqn=((1286-Sm)/107000)(2N)^-0.0875 + 0.15(2N)^-0.6594 ==e

Ni=vpasolve(eqn,N)