A Mecânica Quântica Simplética Relativística (RSQM) considera a álgebra não comutativa de funções no espaço de fase G =(q,p)  e num espaço de Hilbert associado H_{G} para construir uma representação unitária da álgebra de Lie do Grupo de Poincaré. Em consequência uma equação de Klein-Gordon e uma Equação de Dirac são obtidas em termos de variáveis  no espaço de fase e a função de Wigner pode ser derivada diretamente a partir de soluções destas equações. Nesta formula\c{c}\~ao os operadores Q e P são dados por Q=q*=q + (ih/2)(d/dp) e P=p*=p - (ih/2)(d/dq) ; consequentemente as soluções para as equações de Klein-Gordon ou Dirac com potenciais requerem um tratamento especial, para isto estendemos neste trabalho o método de Supersimetria na Mecânica Quântica (SUSYQM) para RSQM: a Hierarquia de Hamiltonianos é gerada e soluções exatas  são obtidas para cada potencial. Tomando uma expressão geral para potencial vetor e potencial escalar em termos de parâmetros apropriados nós analisamos o Oscilador Harmônico de Dirac, um potencial tipo (-A/Q) externo e um potencial que altera a massa. O espectro de auto-valores, as correspondentes auto-funções (amplitudes de quase-

probabilidade no espaço de fase) e funções de Wigner são derivadas e discutidas.