求一个关于最大公因子的matlab代码

无相天

    求一个matlab程序能满足以下两个条件：
    1、已知有两个整数a和b，要能求出他们的最大公因子c
    我们知道最大公因子c可以写成c=p*a+q*b的形式，p,q也是整数，条件二就是：
    2、要能求出p和q
     求大神帮忙解答，谢谢！                                                   

士心之约

1. MATLAB 有自带的函数gcd

2. 查看函数说明 doc gcd

3. 查看函数代码 edit gcd

4. 调用格式 [c,p,q] = gcd(a,b)

5. 案例： [c,p,q] = gcd(126,66)
    结果： c=6,p=-1,q=2.

附上程序源代码，仅供参考。

1. <FONT color=black size=3>function [g,c,d] = gcd(a,b)
   
2. %GCD    Greatest common divisor.
   
3. %   G = GCD(A,B) is the greatest common divisor of corresponding elements
   
4. %   of A and B.  The arrays A and B must contain integer values and must be
   
5. %   the same size (or either can be scalar). GCD(0,0) is 0 by convention;
   
6. %   all other GCDs are positive integers.
   
7. %
   
8. %   [G,C,D] = GCD(A,B) also returns C and D so that G = A.*C + B.*D.
   
9. %   These are useful for solving Diophantine equations and computing
   
10. %   Hermite transformations.
    
11. %
    
12. %   Class support for inputs A,B:
    
13. %      float: double, single
    
14. %      integer: uint8, int8, uint16, int16, uint32, int32, uint64, int64
    
15. %
    
16. %   See also LCM.
    
17. 
    
18. %   References:
    
19. %   Knuth, Donald, The Art of Computer Programming, Vol. 2, Addison-Wesley:
    
20. %      Reading MA, 1973. Section 4.5.2, Algorithms A and X.
    
21. 
    
22. %
    
23. %   Thanks to John Gilbert for the original version
    
24. %   Copyright 1984-2012 The MathWorks, Inc.
    
25. %   $$$$Revision: 5.14.4.9 $$$$  $$$$Date: 2012/10/15 20:09:56 $$$$
    
26. 
    
27. if ~isequal(size(a),size(b)) && ~isscalar(a) && ~isscalar(b)
    
28.     error(message('MATLAB:gcd:InputSizeMismatch'))
    
29. end
    
30. 
    
31. if ~isscalar(a)
    
32.     siz = size(a);
    
33. else
    
34.     siz = size(b);
    
35. end
    
36. a = a(:); b = b(:);
    
37. 
    
38. if ~isreal(a) || ~isequal(round(a),a) || any(isinf(a)) || ...
    
39.         ~isreal(b) || ~isequal(round(b),b) || any(isinf(b))
    
40.     error(message('MATLAB:gcd:NonIntInputs'))
    
41. end
    
42. 
    
43. if isinteger(a)
    
44.     if ~(strcmp(class(a),class(b)) || (isa(b,'double') && isscalar(b)))
    
45.         error(message('MATLAB:gcd:mixedIntegerTypes'))
    
46.     end
    
47.     classin = class(a);
    
48.     if isa(b,'double') && (b > intmax(classin) || b < intmin(classin))
    
49.         error(message('MATLAB:gcd:outOfRange'));
    
50.     end
    
51.     inttype = true;
    
52. elseif isinteger(b)
    
53.     if ~(isa(a,'double') && isscalar(a))
    
54.         error(message('MATLAB:gcd:mixedIntegerTypes'))
    
55.     end
    
56.     classin = class(b);
    
57.     if a > intmax(classin) || a < intmin(classin)
    
58.         error(message('MATLAB:gcd:outOfRange'));
    
59.     end
    
60.     inttype = true;
    
61. else
    
62.     classin = superiorfloat(a,b);
    
63.     largestFlint = flintmax(classin);
    
64.     if any(abs(a) > largestFlint) || any(abs(b) > largestFlint)
    
65.         warning(message('MATLAB:gcd:largestFlint'));
    
66.     end
    
67.     inttype = false;
    
68. end
    
69. 
    
70. if nargout <= 1
    
71.     % intmin in signed integers requires special handling
    
72.     iminIndex = [];
    
73.     if inttype
    
74.         imin = intmin(classin);
    
75.         if imin < 0
    
76.             iminIndex = xor(a == imin, b == imin);
    
77.         end
    
78.     end
    
79.     u = max(abs(a),abs(b));
    
80.     v = min(abs(a),abs(b));
    
81.     u(iminIndex) = u(iminIndex)/2;
    
82.     vnz = v>0;
    
83.     while any(vnz)
    
84.         t = rem(u,v);
    
85.         u(vnz) = v(vnz);
    
86.         v(vnz) = t(vnz);
    
87.         vnz = v>0;
    
88.     end
    
89.     g = reshape(u,siz);
    
90. else
    
91.     if inttype
    
92.         if intmin(classin) == 0    % unsigned integers not supported
    
93.             error(message('MATLAB:gcd:unsupportedType'));
    
94.         end
    
95.     end
    
96.     len = prod(siz);
    
97.     if issparse(a) || issparse(b)
    
98.         u = spalloc(len,3,nnz(a)+len);
    
99.     else
    
100.         u = zeros(len,3,classin);
     
101.     end
     
102.     u(:,1) = 1;
     
103.     u(:,3) = a;
     
104.     if issparse(b)
     
105.         v = spalloc(len,3,nnz(b)+len);
     
106.     else
     
107.         v = zeros(len,3,classin);
     
108.     end
     
109.     v(:,2) = 1;
     
110.     v(:,3) = b;
     
111.     vnz = v(:,3)~=0;
     
112.     while any(vnz)
     
113.         if inttype
     
114.             q = idivide(u(:,3),v(:,3),'fix');
     
115.         else
     
116.             q = fix( u(:,3)./v(:,3));
     
117.         end
     
118.         t = u - bsxfun(@times, v, q);
     
119.         u(vnz,:) = v(vnz,:);
     
120.         v(vnz,:) = t(vnz,:);
     
121.         vnz = v(:,3)~=0;
     
122.     end
     
123.    
     
124.    
     
125.     g = reshape(u(:,3),siz);
     
126.     c = reshape(u(:,1),siz).*sign(g);
     
127.     d = reshape(u(:,2),siz).*sign(g);
     
128.     g = abs(g);
     
129.     % correct overflow conditions in signed integers
     
130.     if inttype
     
131.         overflow1 = reshape(a == intmin(classin) & b == -1, siz);
     
132.         overflow2 = reshape(a == -1 & b == intmin(classin), siz);
     
133.         g(overflow1 | overflow2) = 1;
     
134.         c(overflow1) = 0;
     
135.         d(overflow1) = -1;
     
136.         c(overflow2) = -1;
     
137.         d(overflow2) = 0;
     
138.     end
     
139. end
     
140. </FONT>

复制代码

森之张卫东

        同学，要学会使用Matlab中的帮助文档！！！

无相天

森之张卫东 发表于 2015-10-23 18:50
同学，要学会使用Matlab中的帮助文档！！！

我懂得打开这个，但看不懂这个代码，满足我的条件不需要这么复杂吧

无相天

士心之约 发表于 2015-10-23 18:10
1. MATLAB 有自带的函数gcd

2. 查看函数说明 doc gcd

我懂得打开这个，但看不懂这个代码，满足我的条件不需要这么复杂吧

林威123

我是来混禁言的

无相天

没有我想要的答案，帮助文档谁不懂打开，我看不懂matlab的帮助文档，看不出其中编程的思路，才在此求助。就贴个matlab文档在这，不值300金币。

hzlhm

function GCD=Euclid(m,n)
%欧几里德算法
r=mod(m,n);
while r~=0
    r=mod(m,n);
    m=n;
    n=r;
end
GCD=m;

>> m=319;n=377;Euclid(m,n)
ans =
    29