Extended Euclidean algorithm(扩展欧几里得算法Matlab实现)

文章目录

• 一、解析
• 二、思路
• • 1、Ini_XY迭代初值
  • 2、Ini_XY迭代矩阵
  • 3、流程
• 三、效果如下
• 四、代码
• • 1、Mian.m
  • 2、GCD.m

一、解析

Main.m
矩阵A为一个n·2的矩阵，每一行存储一对待求解的数据
Ini_XY数组为x,y的初始迭代值
GCD.m
形参：待求解数据矩阵，x,y的初始迭代值
返回值：x,y的迭代值，除数矩阵

二、思路

1、Ini_XY迭代初值

2、Ini_XY迭代矩阵

3、流程

由辗转相除法获取{q1,q2,…,qn}
Qi=[0 1;1 -qi]
若两数互质
则有Ini_XY(n+1)迭代初值=[1 0]
迭代如下：
Ini_XY(n)=Qn*Ini_XY(n+1)
Ini_XY(n-1)=Q(n-1)_Ini_XY(n)
…
Ini_XY(1)=Q1_Ini_XY(2)

三、效果如下

四、代码

1、Mian.m

A=[42 30;1759 550;334 111];
Ini_XY=[1 0];
for i=1:3
    disp("========================================");
    disp("第"+i+"次输入：");
    [X_Y,gcd]=GCD(A(i,:),Ini_XY);
    disp("GCD("+A(i,1)+","+A(i,2)+")="+gcd(2));
    if(gcd(2)==1)
        disp(A(i,1)+","+A(i,2)+"互质");
        disp("于是有");
        disp("GCD("+A(i,1)+","+A(i,2)+")="+A(i,1)+"*x+"+A(i,2)+"*y");
        disp("解得:");
        disp("GCD("+A(i,1)+","+A(i,2)+")="+A(i,1)+"*"+X_Y(1)+"+"+A(i,2)+"*"+X_Y(2));
        disp("x="+X_Y(1)+","+"y="+X_Y(2));
    end
end

2、GCD.m

function [X_Y,gcd] = GCD(A,Initial_XY)
%gcd(a,b)=a*x+b*y
%X_Y=[x y]
%gcd：存储除数的数组
a=A(1);
b=A(2);
gcd=[b];
Q_Array=[];
while(b)
    c=mod(a,b);
    Q_Array=[(a-c)/b Q_Array];
    a=b;
    b=c;
    gcd=[b gcd];
end
%x，y迭代初值
X_Y=Initial_XY';
for i=Q_Array
    Q=[0 1;1 -1*i];
    X_Y=Q*X_Y;
end
X_Y=X_Y';
end

总结

写到这里也结束了，在文章最后放上一个小小的福利，以下为小编自己在学习过程中整理出的一个关于 java开发 的学习思路及方向。从事互联网开发，最主要的是要学好技术，而学习技术是一条慢长而艰苦的道路，不能靠一时激情，也不是熬几天几夜就能学好的，必须养成平时努力学习的习惯，更加需要准确的学习方向达到有效的学习效果。

由于内容较多就只放上一个大概的大纲，需要更及详细的学习思维导图的 点击我的Gitee获取。
还有 高级java全套视频教程 java进阶架构师 视频+资料+代码+面试题！

全方面的java进阶实践技术资料，并且还有技术大牛一起讨论交流解决问题。