2014-09-27 23 views
1

我必須導入一個大的整數(1000位)到matlab中,對其運行一些計算。但是,當我輸入它時,由於matlab使用科學記數法這一事實,我似乎失去了準確性。 有什麼辦法可以得到實際的整數?在matlab中導入大整數

這裏的實際數據我不得不進口:

73167176531330624919225119674426574742355349194934 
96983520312774506326239578318016984801869478851843 
85861560789112949495459501737958331952853208805511 
12540698747158523863050715693290963295227443043557 
66896648950445244523161731856403098711121722383113 
62229893423380308135336276614282806444486645238749 
30358907296290491560440772390713810515859307960866 
70172427121883998797908792274921901699720888093776 
65727333001053367881220235421809751254540594752243 
52584907711670556013604839586446706324415722155397 
53697817977846174064955149290862569321978468622482 
83972241375657056057490261407972968652414535100474 
82166370484403199890008895243450658541227588666881 
16427171479924442928230863465674813919123162824586 
17866458359124566529476545682848912883142607690042 
24219022671055626321111109370544217506941658960408 
07198403850962455444362981230987879927244284909188 
84580156166097919133875499200524063689912560717606 
05886116467109405077541002256983155200055935729725 
71636269561882670428252483600823257530420752963450 
+4

你有沒有簽出[VPI](http://www.mathworks.com/matlabcentral/fileexchange/22725-variable-precision-integer-arithmetic)? – 2014-09-27 23:24:56

+1

John D'Errico爲Matlab做出了一些很大的貢獻,其中很多都非常穩固,所以我一定會檢查上面Bruce的鏈接。否則,我懷疑你將不得不編寫一些'split/process/merge'算法,或者編寫一個處理128或256位整數的自定義類,因爲這肯定比任何Matlab數據類型可以處理的位數多。 – Hoki 2014-09-27 23:33:59

回答

1

這樣一個大的整數不能在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