CEBusPL.m

f=[0:359]/359*1.75e5+0.75e5;
%f=4e5;
t=[0:359]/359*1e-4;
chirp1=sin(2*pi*f.*t);
chirp=3.5*[chirp1(131:360) chirp1(1:130)];
figure(1)
plot(t,chirp)
xlabel('Time (Sec.)')
ylabel('Magnitude')
F=fft(chirp);
figure(2)
py=10*log10((abs(F(1:100))).^2/360/10/0.001/3.6e6*2);
plot([1:100]/360*3.6e6, py)
grid
xlabel('Frequency (Hz)')
ylabel('PSD (dBm/Hz)')