exm052_4.m

北京航空航天大学出版社张志涌、徐彦琴主编的matlab教程附带的mfile例程压缩包

%exm052_4.m

p0=sym('(1+sqrt(5))/2')			%p0为(1+sqrt(5))/2	准确值		  
p1=sym((1+sqrt(5))/2)			%p1是(1+sqrt(5))/2在数值环境下的近似值
pd=double(p0)						%双精度数值
e01=vpa(abs(p0-p1))				%在符号环境32位省缺精度下，观察p0, p1间误差  
e0d=vpa(abs(p0-pd))				%在符号环境32位省缺精度下，观察双精度表达的误差  

p2=vpa(p0)							%p2是p0在32位相对精度下的近似
e02=vpa(abs(p0-p2),64)			%在64位相对精度下，观察p0, p2间的误差  

digits								%验证vpa变精度计算不影响全局的符号数值计算精度