這樣一個大的整數不能在IEEE表示浮點標準。檢出this答案爲可以表示不失精度的最大雙倍數(其1.7977e+308
)。這可以通過在MATLAB中輸入realmax
來獲得。
您可以使用vpi
(可在此處獲得,如評論中所述),或者您可以使用MATLAB內置的vpa
。
這是你如何使用vpa
R=vpa('7316717653133062491922511967442657474235534919493496983520312774506326239578318016984801869478851843858615607891129494954595017379583319528532088055111254069874715852386305071569329096329522744304355766896648950445244523161731856403098711121722383113622298934233803081353362766142828064444866452387493035890729629049156044077239071381051585930796086670172427121883998797908792274921901699720888093776657273330010533678812202354218097512545405947522435258490771167055601360483958644670632441572215539753697817977846174064955149290862569321978468622482839722413756570560574902614079729686524145351004748216637048440319989000889524345065854122758866688116427171479924442928230863465674813919123162824586178664583591245665294765456828489128831426076900422421902267105562632111110937054421750694165896040807198403850962455444362981230987879927244284909188845801561660979191338754992005240636899125607176060588611646710940507754100225698315520005593572972571636269561882670428252483600823257530420752963450');
您可以檢查以下內容:
vpa('R+1000-R')
以上的答案是1000預期。不要忘記把你的表情放在引號中。否則,您將無邊界傳遞給vpa
而不是1000位數字。
如果你想使用vpi
,它是一個漂亮的工具箱,請繼續下載。走進它的根目錄,並運行以下命令:
a=vpi('7316717653133062491922511967442657474235534919493496983520312774506326239578318016984801869478851843858615607891129494954595017379583319528532088055111254069874715852386305071569329096329522744304355766896648950445244523161731856403098711121722383113622298934233803081353362766142828064444866452387493035890729629049156044077239071381051585930796086670172427121883998797908792274921901699720888093776657273330010533678812202354218097512545405947522435258490771167055601360483958644670632441572215539753697817977846174064955149290862569321978468622482839722413756570560574902614079729686524145351004748216637048440319989000889524345065854122758866688116427171479924442928230863465674813919123162824586178664583591245665294765456828489128831426076900422421902267105562632111110937054421750694165896040807198403850962455444362981230987879927244284909188845801561660979191338754992005240636899125607176060588611646710940507754100225698315520005593572972571636269561882670428252483600823257530420752963450')
那麼,與vpi
的優點如下:
的vpi
輸出:
a=vpi(<<Your 1000 digit number in quotes>>); %output prints 1000 digits on screen.
的vpa
輸出:
R=vpa(<<Your 1000 digit number in quotes>>);
此打印:
R =
7.3167176531330624919225119674427e999
而且,與vpi
,你可以做這樣的事情:
a=vpi('7316717653133062491922511967442657474235534919493496983520312774506326239578318016984801869478851843858615607891129494954595017379583319528532088055111254069874715852386305071569329096329522744304355766896648950445244523161731856403098711121722383113622298934233803081353362766142828064444866452387493035890729629049156044077239071381051585930796086670172427121883998797908792274921901699720888093776657273330010533678812202354218097512545405947522435258490771167055601360483958644670632441572215539753697817977846174064955149290862569321978468622482839722413756570560574902614079729686524145351004748216637048440319989000889524345065854122758866688116427171479924442928230863465674813919123162824586178664583591245665294765456828489128831426076900422421902267105562632111110937054421750694165896040807198403850962455444362981230987879927244284909188845801561660979191338754992005240636899125607176060588611646710940507754100225698315520005593572972571636269561882670428252483600823257530420752963450')
b=a+1
b-a %output of this yields 1.
我有點不能做的b-a
的操作vpa
並獲得答案1
。
你有沒有簽出[VPI](http://www.mathworks.com/matlabcentral/fileexchange/22725-variable-precision-integer-arithmetic)? – 2014-09-27 23:24:56
John D'Errico爲Matlab做出了一些很大的貢獻,其中很多都非常穩固,所以我一定會檢查上面Bruce的鏈接。否則,我懷疑你將不得不編寫一些'split/process/merge'算法,或者編寫一個處理128或256位整數的自定義類,因爲這肯定比任何Matlab數據類型可以處理的位數多。 – Hoki 2014-09-27 23:33:59