specture.m

本程序实现随机数的产生和在相同调制方式下功率谱比较

%－－分别求两信号序列的功率谱密度函数
A=fft(sym0_fc,200);
Paa=A.*conj(A)/200;               %－－第一组信号序列的功率谱密度函数
A1=fft(a1,200);
Paa1=A1.*conj(A1)/200;            %－－第二组信号序列的功率谱密度函数
f=1000*(0:100)/200;               %－－频域采样，分辨率5Hz
figure(3)
subplot(211);
plot(f,Paa(1:101),'r');
title('Frequency content for signal1')
xlabel('frequency (Hz)')
grid on
subplot(212)
plot(f,Paa1(1:101),'r');
title('Frequency content for signal2');
xlabel('frequency (Hz)');
grid on;