Подсчет количества точных цифр в приближении к истинному значению - PullRequest
Новая
       36

Подсчет количества точных цифр в приближении к истинному значению

0 голосов
/ 19 сентября 2019

Я создал код, который выполняет приближение метода Ньютона.Он печатает в табличном формате аппроксимацию на каждом шаге и соответствующую ошибку.Я хочу добавить столбец, который показывает целочисленное значение, которое представляет количество правильных цифр в приближении против истинного значения.

Я пытаюсь преобразовать каждую ячейку аппроксимации в строку и подсчитывать, сколько цифр является точным,Пример, ок.= 3,14555, правда = 3,1555.Число точных цифр будет 2. Хотя у меня есть эта идея в голове, я все делаю неправильно в своем коде ниже.Знаете ли вы, как создать правильный цикл для достижения этой цели?У меня меньше года опыта в MATLAB;мой умственный инструментарий ограничен. 

% Program Code of Newton's Method to find root
% This program will not produce a result if initial guess is too far from
% true value
clear;clc;format('long','g')
% Can work for various functions
%FUNCTION: 2*x*log(x)-2*log(x)*x^(3)+2*x^(2)*log(x)-x^(2)+1
%INTIAL GUESS: .01
%ERROR: 1.e-8
a=input('Enter the function in the form of variable x:','s');
x(1)=input('Enter Initial Guess:');
error=input('Enter allowed Error:');
% Passing through the function and calculating the derivative
f=inline(a);
dif=diff(str2sym(a));
d=inline(dif);
% Looping through Newton's Method
for i=1:100
x(i+1)=x(i)-((f(x(i))/d(x(i))));
err(i)=abs(x(i+1)-x(i));
% The loop is broken if acceptable error magnitude is reached 
if err(i)<error
break
end
end
root=x(i);
Root = (x(:,1:(end-1)))';
Error = err';

disp('The final approximation is:')
disp(root)
%BELOW IS ALL WRONG, I AM TRYING TO ADD A COLUMN TO 'table' 
%THAT SHOWS HOW MANY DIGITS IN APPROXIMATION IS ACCURATE
iter = 0;
y = zeros(1,length(x)); 
plot(x,y,'+')
zero1 = ('0.327967785331818'); %ACTUAL VALUE
for i = 1:length(Root)
    chr = mat2str(Root(i))
    for j = 1:length(chr(i))
        if chr(i)~=zero1(i)
            iter = 0;
            return
        elseif chr(i)==zero1(i)
            iter = iter + 1;
            acc(i) = iter
        end

    end

end
table(Root, Error)   %ADD ACCURACY COLUMN HERE

1 Ответ

0 голосов
/ 19 сентября 2019

Возможно, что-то вроде умножения обоих чисел на 10 и последующего их определения до тех пор, пока ответы не станут равными: 

approx=3.14555;
truth=3.1555;
approx1=0;
truth1=0;
i=0;
while approx1==truth1
   approx1=floor(approx*10^i);
   truth1=floor(truth*10^i);
   i=i+1;
end
acc=i-1;
...