not sure if all correct,

* 即户主取得收入1时，配偶取得取得收入1-4的比例各自为多少
mat SP_INC=J(4,4,.)
forv i=1/4 {
  egen x`i'=sum(mainr*income`i'), by(householdid)
  gen sp_income`i' = x`i'*(mainr==0)
  forv j=1/4 {
          su income`j' if x`i' == 1 & mainr == 0
          mat SP_INC[`i',`j'] = r(mean)
  }
}
drop x*
svmat SP_INC
gen HD_INC = _n in 1/4
gr bar SP_INC*, over(HD_INC)

*将户主的变量信息直接赋予配偶
forv i=1/4 {
  egen hdsp_income`i'=sum(mainr*income`i'), by(householdid)
}

*直接保存户主或者配偶中至少有1个人取得过income1的收入的那部分数据drop _all
forv i=1/4 {
  egen x`i'=sum(income`i'), by(householdid)
}
keep if x1 >0
drop x*

* 如果要保存户主或配偶中个人收入较高的那个作为唯一的观测值
egen double x = max(pincome), by(householdid)
keep if x == pincome
drop x

it does work,thank you!