Matlab中所有 m 序列的特征多项式怎么求？

题目：寻找与特征多项式8进制表示的1021的 m 序列构成优选对的所有 m 序列的特征多项式？
由题目转化1021的八进制转二进制得[1,0,0,0,0,0,1,0,0,0,1]，r=9，
我只知道pr=primpoly(9,'all')是输出所有9级的本原序列，我认为思路是求出1021的相关函数绝对值的最大值，然后用输出所有9级的本原序列的相关函数绝对值的最大值去和1021的作比较然后输出，可能要用for循环，可是鄙人不会写，课上也不好意思提问，所以在这里请教各位学长，望能帮助鄙人。