matlab 求解黄金分割率

求解黄金分割率

clear;
clc;
%% 本程序求解黄金分割率的比值
% r^2 - r - 1 = 0的解就是比值。
p = [1 -1 -1];
% 此数组代表了上式的二次项系数、一次项系数和常数项。
r = roots(p);
print_str = sprintf('r^2 - r - 1 = 0的结果是:%f和%f\n', r);
disp(print_str);


%% 第二种计算的方法。
clear;
clc;
syms r; r2 = solve(r^2-r-1,r);
print_str = sprintf('r*r - r - 1 = 0的结果是:%f和%f\n', r2);
disp(print_str);

%% 提高精确度的方法。
clear;
clc;
syms r; r2 = solve(r^2-r-1,r);
print_str = sprintf('r*r - r - 1 = 0的结果是:%f和%f\n', r2);
disp(print_str);
precise_res=vpa(r2,50000)
disp(precise_res);