matlab secant method

最新推荐文章于 2023-09-08 01:49:22 发布

weixin_30919235

最新推荐文章于 2023-09-08 01:49:22 发布

阅读量184

点赞数

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/yxwkf/p/4553880.html

本文介绍了一个使用Matlab编写的割线法脚本，该脚本用于求解非线性方程，具体应用为寻找指定数值的平方根。用户需输入目标数值及两个初始猜测值，通过迭代计算直至达到所需精度。

% Matlab script to illustrate the secant method
% to solve a nonlinear equation

% this particular script finds the square root of a number M
% (input by the user)

% note that the function we are trying to zero is f(x) = x^2 - M.
% this function is hard-coded in the script.

g=9.8065;
k=0.00341;
% f(x)=log(cosh(t*srt(g*k)))/k;

format long

% get user input
M = input('Please enter the number whose square root you want: ')
t0 = input('Please enter the first of two starting guesses: ')
t1 = input('Please enter the second of two starting guesses: ')

% iteration counter
k = 1
% compute first secant iterate to enter loop
s = (((log(cosh(t1*sqrt(g*k)))/k)-M)-((log(cosh(t0*sqrt(g*k)))/k)-M) )/(t1-t0);
% s = ( (x1^2-M) - (x0^2-M) ) / (x1 - x0);
t = t1 - (((log(cosh(t1*sqrt(g*k)))/k)-M))/s
% x = x1 - (x1^2-M)/s
disp('Hit return to continue')
pause

while abs(t-t1) > eps*abs(t),
% reset guesses
t0 = t1;
t1 = t;
% increment iteration counter
k = k + 1
% compute and display secant iterate
s = (((log(cosh(t1*sqrt(g*k)))/k)-M)-((log(cosh(t0*sqrt(g*k)))/k)-M) )/(t1-t0);
% s = ( (x1^2-M) - (x0^2-M) ) / (x1 - x0);
% x = x1 - (x1^2-M)/s
t = t1 - (((log(cosh(t1*sqrt(g*k)))/k)-M))/s
disp('Hit return to continue')
pause
end

转载于:https://www.cnblogs.com/yxwkf/p/4553880.html