yhkn.net
当前位置:首页 >> mAtlAB画图:用二分法求方程x^3Cosx+6sin^5x%3=0 >>

mAtlAB画图:用二分法求方程x^3Cosx+6sin^5x%3=0

首先,自定义函数 dichotomy_fun=inline('x^3-3*x^2-x+3','x') 然后,执行下列命令 dichotomy(dichotomy_fun,0,1,1e-6) ...

用二分法,你得先找出一个根所在的区间 function [a,b]=findbracket(f,x0) %f is the function evaluated %x0 is the starting point %a is the left boundary %b is the right boundary %[a,b] is the interval around x0 wanted a=x0; b=x0; d...

这个定义一个函数。 f =@(x)x^3-3*x^2-x+3; 的意思是说定义一个函数f(x),它只有一个自变量。 使用时,直接可以用f(1)它就是x=1时的函数值。 另外多参数可以是: f=@(x,y)sqrt(x.^2+y.^2),求点到原点的距离。

先建立二分法的fun.m文件,代码如下: function fun(a,b,e) %f是自定义的函数 %a为隔根区间左端点,b为隔根区间右端点,e为绝对误差限 if nargin==2 e=1.0e-6; elseif nargin=b input('隔根区间输入错误!'); return; end a1=a; b1=b; c1=(a1+b1...

%问题出在第7行,x3=(x1+x1)/2,结果肯定是进入死循环了,改为x3=(x1+x2)/2应该就可以了。

Matlab里面采用的算法我不知道是不是二分法,但是结果都很精确的。这里你可以尝试fzero命令解,当然你也可用solve命令。 syms x f=@(x)x^4-3*x+1; x=fzero(f,0.3) 运行结果: x = 0.337666765642802 如果用solve的话:x=solve('x^4-3*x+1','x') ...

%我想证明很简单,令f(x)=x^3+x-4,则f'(x)=3*x^2+1>0;f(x)单调递增,f(1)=-20,所以方 %程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)=...

源代码: 在matlab中保存为:demimethod.m function [x,k]=demimethod(a,b,f,emg) % a,b :求解区间的两个端点 % f :所求方程的函数名 % emg :精度指标 % x:所求近似解 % k: 循环次数 fa=feval(f,a); fab=feval(f,(a+b)/2); k=0; while abs(b-a)>...

举个简单例子: 问题1:求f(x)=1-x-sinx=0在【0,1】的根 误差不超过0.5*10^(-4) 解答: 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); if fa*fb>0,break,end while abs(fc)>0.5*10^(-4) c=(a+b)/2; fc=1-c-sin(...

在Matlab下输入:edit,然后将下面两行百分号之间的内容,复制进去,保存 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function y=zhidao_fx(x) y=x.^3-sin(x)-12*x+1; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 在Matlab下输...

网站首页 | 网站地图
All rights reserved Powered by www.yhkn.net
copyright ©right 2010-2021。
内容来自网络,如有侵犯请联系客服。zhit325@qq.com