图灵机

用程序语言来模拟图灵机的运作过程实现对一个数乘以2

1.题目分析
用程序代码来模拟图灵机的工作步骤，将输入的一个十进制数转化为二进制，利用书上的拓展原理将二进制码转化为图灵机能运作的拓展序列（磁带），再将拓展序列进行XNx2的变化，得到“x2”后的拓展序列，转化为二进制，最后转化为十进制。
二进制转化为拓展码的过程实际上是在一串零里,根据二进制的形式加1，所有的0,1都是储存在数组里，所有的变化是在数组的结构基础上实现。
为了可以更直观的感受图灵机的工作步骤，把输入的数字的拓展码进行XNx2的每一步变化后的序列都输出。
2.算法构造
（1）实现十进制到二进制
十进制整数转换为二进制整数采用"除2取余，逆序排列"法。具体做法是：用2整除十进制整数，可以得到一个商和余数；再用2去除商，又会得到一个商和余数，如此进行，直到商为小于1时为止，然后把先得到的余数作为二进制数的低位有效位，后得到的余数作为二进制数的高位有效位，依次排列起来。将这些余数存储在数组里。
（2）实现二进制到拓展序列
拓展序列就是磁带上的数列。将磁带tape[]数组的后五位数初始化为11000，人工设置结束标志。这样只需将二进制数转化为它的拓展码。从存储二进制的数组a[]中倒序读取0/1，读到1就从tape[]的倒数第六位开始往后两个零，假设读到a[i],在a[i]存入1，然后i–,反之如果在a[]中读到零，就让i=i-2。这样一直读完a[]中的二进制码。
（3）实现XNx2
模拟图灵机的这六条指令：00->00R; 01->10R; 10->01R;11->100R; 100>111R; 110->01STOP;
将这六个指令看成六个条件，图灵机执行这六条指令与内态和在磁条上读取的数字有关。例如00->00R，当内态为0且此时读到磁条上的数也为0时，内态和这时磁条上的数字都不变，再读磁条的下一个数，当读到tape[]的倒数第六位时结束读取，这时tape[]中的拓展序列已被更新到XNx2后的结果。
（4）实现拓展序列到二进制
创建一个新数组b[j]来存储得到的二进制码，从tape[]的倒数第六位开始读起，从b[]的最末开始存储。当读到tape[i],这一位的前后两位都为1时，将1存到b[j]中，i=i-2;j–;不满足此条件时，存入0到b[j]里，i–;j–;直到读完tape[]。
（5）实现二进制到十进制
把二进制数首先写成加权系数展开式，然后按十进制加法规则求和。从b[j]的最后一位读起，计算每一位对应的十进制数字（m=(b[i])*pow(2,j-i-2)；）然后求和，就得到此二进制码的对应的十进制数字。
3.算法实现

#include <stdio.h>
#include<stdlib.h>
#include<math.h>
#include<iostream.h>
#define MAX 22
void translate(int y)
   {  
	int x,a[34],i=0,count=0,k;
    x=y;
	while(y!=0)//实现十进制数进行XNx2之后的十进制结果
	{
		a[i++]=y%2; //将结果存放到数组中
		y=y/2;
	}
	printf("十进制整数%d转换为二进制数是: ",x);
	for(i--;i>=0;i--)
	{	printf("%d",a[i]);
     	count++;
	}//通过数组输出

	cout<<endl;
	int tape[MAX]={0};
	for(i=MAX-5;i<MAX;i++)
	{
	if(i>MAX-4)
	{tape[i]=0;}
	else tape[i]=1;
	}
	i=0;
	int j=MAX-6;
	int in=0;
    while(i<=(count-1))
	{
       if(a[i]==1)
	   {tape[j-1]=1;j=j-2;i=i+1; }
	   if(a[i]==0)
	   {j=j-1;i=i+1;}
	}
    printf("转化为拓展码为：  ");
		for(j=0;j<20;j++)
    {
		printf("%d",tape[j]);
	}
    cout<<endl;
	printf("以下为模拟图灵机工作的步骤变化\n");
	for(i=0;i<MAX;)//模拟图灵机的工作步骤
	{
        if(in==0&&tape[i]==0)
		{ 	in=0;tape[i]=0;i++;	
		for(k=0;k<MAX;k++)
		{
			printf("%d",tape[k]);
		}
		printf("\n");
		}
		if(in==0&&tape[i]==1)
		{in=1;tape[i]=0;i++;    }
		if(in==1&&tape[i]==0)
		{in=0;tape[i]=1;i++;
        for(k=0;k<MAX;k++)
		{
			printf("%d",tape[k]);
		}
		printf("\n");
		}
		if(in==1&&tape[i]==1)
		{	in=10;tape[i]=0;i++;
       for(k=0;k<MAX;k++)
		{
			printf("%d",tape[k]);
		}
		printf("\n");
		}
		if(in==10&&tape[i]==0)
		{in=11;tape[i]=1;i++;
        for(k=0;k<MAX;k++)
		{
			printf("%d",tape[k]);
		}
		printf("\n");
		}
		if(in==11&&tape[i]==0)
		{in=0;tape[i]=1;i=MAX;
        for(k=0;k<MAX;k++)
		{
			printf("%d",tape[k]);
		}
		printf("\n");
		}
	}
    printf("XNx2后的拓展码为：");
	for(i=2;i<22;i++)
	{
		printf("%d",tape[i]);
	}	
    cout<<endl;
    int b[MAX];
	j=0;
    for(i=0;i<MAX-3;) //将拓展码收缩为二进制
	{
		if((tape[i]==1&&tape[i+1]==0)||(tape[i]==0&&tape[i+1]==1))
		{b[j]=1;j++;i=i+2;}
		else 
		{b[j]=0;j++;i=i+1;}
	}
    printf("拓展码转化为二进制：  ");
    for(i=0;i<j-1;i++)
	{
		while(b[i]!=0)
		{  for(;i<j-1;i++)
			{printf("%d",b[i]);}
		break;
		}
	}
    cout<<endl;
	int sum=0,m=0;
	for(i=j-2;i>=0;i--)
	{
        m=(b[i])*pow(2,j-i-2);//计算当前二进制所在位的权值
		sum=sum+m;//累加权值得到十进制数
	}
	cout<<"转化为十进制数为："<<sum<<endl;    
}
void main()
{	
	int x;
	printf("请输入一个十进制整数：");
	cin>>x;
	translate(x);
}

4.调试、测试及运行结果
4.1调试

4.2测试
（1）十进制到二进制

（2）二进制到拓展序列

（3）XNx2

（4）拓展序列到二进制

（5）二进制到十进制

4.3完整运行截图

5.经验归纳
经过此次上机实验，我对图灵机的工作步骤更加掌握。将图灵机的工作步骤用代码实现是需要具体分成几部分来实现的，我觉得最难的是拓展码变为二进制码，这需要我熟练的运用数组，变化数组，对数组中的元素进行操纵。
不能只想到怎么实现算法，要能将自己及所想的东西用程序语言实现才表示自己真的掌握了。这次程序极大部分用到了数组，对数组中的元素进行来回的操纵，很容易把自己搞混，这就要求自己的脑中很明确算法步骤，这一步做完，下一步应该怎么做，有些变量的值也在不断变化，自己都要将这些考虑到。