二叉树Morris Traversal

常见的二叉树遍历有递归和栈循环两种方式,其实还有另一种更为巧妙的遍历方式Morris Traversal。

Morris Traversal的空间复杂度为O(1),时间复杂度为O(n)

我们知道,在深度搜索遍历的过程中,之所以要用递归或者是用非递归的栈方式,都是因为其他的方式没法记录当前节点的parent,而如果在每个节点的结构里面加个parent 分量显然是不现实的,那么Morris是怎么解决这一问题的呢?好吧,他用得很巧妙,实际上是用叶子节点的空指针来记录当前节点的位置,然后一旦遍历到了叶子节点,发现叶子节点的右指针指向的是当前节点,那么就认为以当前节点的左子树已经遍历完成。

以inorder为例,初始化当前节点为root,它的遍历规则如下:

  • 如果当前节点为空,程序退出。
  • 如果当前节点非空,
    • 如果当前节点的左儿子为空,那么输出当前节点,当前节点重置为当前节点的右儿子。
    • 如果当前节点的左儿子非空,找到当前节点左子树的最右叶子节点(此时最右节点的右儿子有两种情况,一种是指向当前节点,一种是为空,你也许感到奇怪,右节点的右儿子怎么可能非空,注意,这里的最右叶子节点只带的是原树中的最右叶子节点。),若其最右叶子节点为空,令其指向当前节点,将当前节点重置为其左儿子,若其最右节点指向当前节点,输出当前节点,将当前节点重置为当前节点的右儿子,并恢复树结构,即将最右节点的右节点再次设置为NULL

morris 中序遍历

[cpp]  view plain  copy
  1. void bst_morris_inorder(struct bst_node *root)  
  2. {  
  3.     struct bst_node *p = root, *tmp;  
  4.   
  5.     while (p) {  
  6.         if (p->left == NULL) {  
  7.             printf("%d ", p->key);  
  8.             p = p->right;  
  9.         }  
  10.         else {  
  11.             tmp = p->left;  
  12.             while (tmp->right != NULL && tmp->right != p)  
  13.                 tmp = tmp->right;  
  14.             if (tmp->right == NULL) {  
  15.                 tmp->right = p;  
  16.                 p = p->left;  
  17.             }  
  18.             else {  
  19.                 printf("%d ", p->key);  
  20.                 tmp->right = NULL;  
  21.                 p = p->right;  
  22.             }  
  23.         }  
  24.     }  
  25. }  
代码相当简单,局部变量也只需要两个。以下简单讲述其实现方式。


morris traversal 原理很简单,利用所有叶子结点的right 指针,指向其后继结点,组成一个环,在第二次遍历到这个结点时,由于其左子树已经遍历完了,则访问该结点

如下图为morris 的一个示例



morris 前序遍历

[cpp]  view plain  copy
  1. void bst_morris_preorder(struct bst_node *root)  
  2. {  
  3.     struct bst_node *p = root, *tmp;  
  4.   
  5.     while (p) {  
  6.         if (p->left == NULL) {  
  7.             printf("%d ", p->key);  
  8.             p = p->right;  
  9.         }  
  10.         else {  
  11.             tmp = p->left;  
  12.             while (tmp->right != NULL && tmp->right != p)  
  13.                 tmp = tmp->right;  
  14.             if (tmp->right == NULL) {  
  15.                 printf("%d ", p->key);                  
  16.                 tmp->right = p;  
  17.                 p = p->left;  
  18.             }  
  19.             else {  
  20.                 tmp->right = NULL;  
  21.                 p = p->right;  
  22.             }  
  23.         }  
  24.     }  
  25. }  
比较简单,只要注意访问结点的顺序即可

morris 后序遍历

[cpp]  view plain  copy
  1. static void bst_morris_reverse(struct bst_node *node, struct bst_node *last)  
  2. {  
  3.     struct bst_node *p = node, *x, *y;  
  4.     if (p == last) {  
  5.         printf("%d ", last->key);  
  6.         return;  
  7.     }  
  8.   
  9.     /* change right to parent pointer */  
  10.     x = p->right;  
  11.     for (;;) {  
  12.         if (x == last) {  
  13.             x->right = p;  
  14.             break;  
  15.         }  
  16.         y = x->right;  
  17.         x->right = p;  
  18.         p = x;  
  19.         x = y;  
  20.     }  
  21.   
  22.     /* visit each */  
  23.     x = last;  
  24.     for (;;) {  
  25.         printf("%d ", x->key);  
  26.         if (x == node)  
  27.             break;  
  28.         x = x->right;  
  29.     }   
  30.   
  31.     /* revert right pointer */  
  32.     p = last;  
  33.     x = last->right;  
  34.     for (;;) {  
  35.         if (x == node) {  
  36.             x->right = p;  
  37.             break;  
  38.         }  
  39.         y = x->right;  
  40.         x->right = p;  
  41.         p = x;  
  42.         x = y;  
  43.     }  
  44. }  
  45.   
  46.   
  47. void bst_morris_postorder(struct bst_node *root)  
  48. {  
  49.     struct bst_node dummy;  
  50.     struct bst_node *p, *tmp;  
  51.   
  52.     dummy.left = root;  
  53.     dummy.right = NULL;  
  54.     p = &dummy;  
  55.   
  56.     while (p) {  
  57.         if (p->left == NULL) {  
  58.             p = p->right;  
  59.         }  
  60.         else {  
  61.             tmp = p->left;  
  62.             while (tmp->right != NULL && tmp->right != p)  
  63.                 tmp = tmp->right;  
  64.             if (tmp->right == NULL) {  
  65.                 tmp->right = p;  
  66.                 p = p->left;  
  67.             }  
  68.             else {  
  69.                 bst_morris_reverse(p->left, tmp);  
  70.                 tmp->right = NULL;  
  71.                 p = p->right;  
  72.             }  
  73.         }  
  74.     }  
  75. }  
后序遍历比较复杂,它的思路是利用中序遍历,所以首先产生了一个假的根结点,其左子树为原来的二叉树,从假的根结点开始中序遍历

它和中序遍历有所不同,在发现当前结点左子树为空时,不访问此结点(后序遍历需要保证访问完右子树后才能访问根结点),直接访问右子树。

第二次遍历到某个结点时,将该结点左子树的最右路径反序输出即可,对应函数为bst_morris_reverse



  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
在信号处理领域,DOA(Direction of Arrival)估计是一项关键技术,主要用于确定多个信号源到达接收阵列的方向。本文将详细探讨三种ESPRIT(Estimation of Signal Parameters via Rotational Invariance Techniques)算法在DOA估计中的实现,以及它们在MATLAB环境中的具体应用。 ESPRIT算法是由Paul Kailath等人于1986年提出的,其核心思想是利用阵列数据的旋转不变性来估计信号源的角度。这种算法相比传统的 MUSIC(Multiple Signal Classification)算法具有较低的计算复杂度,且无需进行特征值分解,因此在实际应用中颇具优势。 1. 普通ESPRIT算法 普通ESPRIT算法分为两个主要步骤:构造等效旋转不变系统和估计角度。通过空间平移(如延时)构建两个子阵列,使得它们之间的关系具有旋转不变性。然后,通过对子阵列数据进行最小二乘拟合,可以得到信号源的角频率估计,进一步转换为DOA估计。 2. 常规ESPRIT算法实现 在描述中提到的`common_esprit_method1.m`和`common_esprit_method2.m`是两种不同的普通ESPRIT算法实现。它们可能在实现细节上略有差异,比如选择子阵列的方式、参数估计的策略等。MATLAB代码通常会包含预处理步骤(如数据归一化)、子阵列构造、旋转不变性矩阵的建立、最小二乘估计等部分。通过运行这两个文件,可以比较它们在估计精度和计算效率上的异同。 3. TLS_ESPRIT算法 TLS(Total Least Squares)ESPRIT是对普通ESPRIT的优化,它考虑了数据噪声的影响,提高了估计的稳健性。在TLS_ESPRIT算法中,不假设数据噪声是高斯白噪声,而是采用总最小二乘准则来拟合数据。这使得算法在噪声环境下表现更优。`TLS_esprit.m`文件应该包含了TLS_ESPRIT算法的完整实现,包括TLS估计的步骤和旋转不变性矩阵的改进处理。 在实际应用中,选择合适的ESPRIT变体取决于系统条件,例如噪声水平、信号质量以及计算资源。通过MATLAB实现,研究者和工程师可以方便地比较不同算法的效果,并根据需要进行调整和优化。同时,这些代码也为教学和学习DOA估计提供了一个直观的平台,有助于深入理解ESPRIT算法的工作原理。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值