arma=ARMAProcess[{(Subscript[a, 1]),(Subscript[a, 2])},{(Subscript[b, 1]),(Subscript[b, 2])},σ];
ws=WeakStationarity[arma]

RegionPlot[ws,{(Subscript[a, 1]),-2,2},{(Subscript[a, 2]),-2,2},GridLines->Automatic]

aspt=ProcessParameterAssumptions[arma];
rest=(Subscript[a, 1])≠0&&(Subscript[a, 2])≠0&&(Subscript[b, 1])≠0&&(Subscript[b, 2])≠0&&σ≠0;
wsARMA=arma/.FindInstance[ws&&rest&&aspt,{(Subscript[a, 1]),(Subscript[a, 2]),(Subscript[b, 1]),(Subscript[b, 2]),σ}][[1]]

WeakStationarity[wsARMA]

in=TimeSeriesInvertibility[ARMAProcess[{(Subscript[a, 1]),(Subscript[a, 2])},{(Subscript[b, 1]),(Subscript[b, 2])},(σ^2)]]

RegionPlot[in,{(Subscript[b, 1]),-2,2},{(Subscript[b, 2]),-2,2},GridLines->Automatic]

aspt=ProcessParameterAssumptions[arma];
rest=(Subscript[a, 1])≠0&&(Subscript[a, 2])≠0&&(Subscript[b, 1])≠0&&(Subscript[b, 2])≠0&&σ≠0;
invARMA=arma/.FindInstance[in&&rest&&aspt,{(Subscript[a, 1]),(Subscript[a, 2]),(Subscript[b, 1]),(Subscript[b, 2]),σ}][[1]]

TimeSeriesInvertibility[invARMA]