Need Help with Matlab (Signal Processing)

Hello, I have a Matlab project for signal analysis (Modified Periodgram )
And since I have no experience with matlab i thought I can ask for a lil bit of help .

Q:

X = A1 sin (ω1 nT+Φ1) + A2 sin (ω2 nT+Φ2) + V

Where:
A1=0.1 , A2=1 ; amplitudes in V
ω1= 2Π0.1, ω2= 2Π0.15 , angular freq. in rad/s
T=1 sampling time in s
Φ1,Φ2, statistically independent phase processes distributed between [0,2Π]
V is a zero-mean white gaussian noise process of σ²= 0.01

how to generate 50 different realizations of the process X whith N=128 samples each , estimate power spectrum using the modified periodgram , and plot the magnitude/frequency using a dB /rad scales

I already have a draft of this matlab code but I believe I have tons of mistakes , I'm looking forward to getting help from you guys , thanks a lot .