matlab实现二分法解方程

function [x] = halfSolve(a, b, f ,tol)  % tol 是 tolerance 的缩写，表示绝对误差c = (a + b) / 2; k = 1;m = 1 + round((log(b - a) - log(2 * tol)) / log(2)); % <1>while k <= m + 10if (f(c) == 0)x = c;break;elseif f(a) * f(c) < 0b = (a + b) / 2;elsea = (a + b) / 2;endc = (a + b) / 2; k = k + 1;endx = c; % 这里加分号是为了不再命令行中输出k % 不加分号就会在控制台输出c
endclc;clear all;
v_y = 10;
R = 565.98;      %光源轨迹半径syms lamda
fun = @(lamda)(lamda + (R/v_y) * sin(lamda) - pi);
a = 0;
b = pi;
tol = 1.0e-15;         %控制误差精度
result = halfSolve(a, b,fun,tol);%% 求交点
syms lamda
eqn1 =  lamda + (R/v_y) * sin(lamda) - pi  == 0;
solx1 = vpasolve(eqn1 , lamda,[0,pi])

参考链接：
https://www.cnblogs.com/fanlumaster/p/14551258.html

matlab实现二分法解方程相关推荐

二分法求方程根matlab,matlab用二分法求方程的正根,要求误差小于0.0005
matlab编程题:用二分法求方程x^3-3*x-1=0的根先建立二分法的fun.m文件,代码如下:functionfun(a,b,e)%f是自定义的函数%a为隔根区间左端点,b为隔根区间右端点,e ...
matlab中solve解方程,怎么用Matlab解方程？
Matlab是一款功能强大的数学软件,我们常常会使用它来解决一些数学难题,但是有时候我们刚装上这个软件的时候,就有一些问题不懂得怎么处理,比如说怎么用Matlab解方程,为了帮助大家解决这个问题,小编 ...
【Matlab 控制】解方程solve() 或代入函数求值 eval() subs()
解方程用 solve() syms x eqn=sin(x)==1; solve(eqn,x) 得出方程结果 ans =pi/2 From: MATLAB的solve函数代入求值用 eval() 配 ...
解题思考F. 解方程（二分法解方程）
题目描述 Problem Description 给定方程 8x^4 + 7x^3 + 2x^2 + 3x + 6 == Y,请计算x在[0,100]范围内的解. Input 输入数据首先是一个正整数 ...
matlab中solve解方程,matlab解方程组.我的每一个方程都很长,如果用solve,solve里的参数就会很长.有什么办法呢?...
问题描述: matlab解方程组.我的每一个方程都很长,如果用solve,solve里的参数就会很长.有什么办法呢? s= (((22*r)/5 - 479/50)*(7*r^2 + (149*r)/ ...
matlab用solve解方程错误提示,MATLAB中使用solve解决方程组的问题
希望使用MATLAB的solve函数解出一个带有虚数的方程组,但是一直提示计算错误,要么就是算不出来结果,希望大佬们能帮帮忙程序如下: syms a1 a2; a=[a1 a2]; C11=3.06 ...
matlab匿名函数解方程,matlab 数学建模匿名函数@
匿名函数的基本用法: handle = @(arglist)anonymous_function 其中handle为调用匿名函数时使用的名字.arglist为匿名函数的输入参数,可以是一个,也可以是多 ...
二分法解具有单调性的方程
解方程的手段有很多,如二分法.牛顿迭代法等等,本次介绍的是采用二分法,在使用二分法解方程时,有一点应该特别注意,就是那种"具有单调性"的函数才可以,否则是会有问题的 ** HDU ...
学习Matlab强大的符号计算（解方程）
学习Matlab强大的符号计算(解方程) 分类: 学习其它语言 2011-09-19 10:34 3135人阅读评论(0) 收藏举报 matlab [plain] view plaincopyp ...

matlab实现二分法解方程

matlab实现二分法解方程相关推荐

最新文章

热门文章