test.m

完整的幅相调制信号可k阶非线性功率谱分析

clear all;
clc
R=100000;
fc=200000;
fs=800000;
rf=0.5;
N=100000;
shape=2;%升余弦脉冲成型
mode1=1;%测试BPSK信号生成
mode2=2;%测试QPSK信号生成
mode3=3;%测试OQQPSK信号生成
mode4=4;%测试pi4QPSK信号生成
mode5=5;%测试18PSK信号生成
mode6=6;%测试16QAM信号生成
mode7=7;%测试32QAM信号生成
mode8=8;%测试64QAM信号生成
mode9=9;

[y1,x1]=modla(N,fs,fc,R,shape,mode1,rf);%BPSK矩形脉冲成型信号生成
[y2,x2]=modla(N,fs,fc,R,shape,mode2,rf);%BPSK矩形脉冲成型信号生成
[y3,x3]=modla(N,fs,fc,R,shape,mode3,rf);%BPSK矩形脉冲成型信号生成
[y4,x4]=modla(N,fs,fc,R,shape,mode4,rf);%BPSK矩形脉冲成型信号生成
[y5,x5]=modla(N,fs,fc,R,shape,mode5,rf);%BPSK矩形脉冲成型信号生成
[y6,x6]=modla(N,fs,fc,R,shape,mode6,rf);%BPSK矩形脉冲成型信号生成
[y7,x7]=modla(N,fs,fc,R,shape,mode7,rf);%BPSK矩形脉冲成型信号生成
[y8,x8]=modla(N,fs,fc,R,shape,mode8,rf);%BPSK矩形脉冲成型信号生成
[y9,x9]=modla(N,fs,fc,R,shape,mode9,rf);%BPSK矩形脉冲成型信号生成

[P1,f1]=kthspec(y1,2,fs);
[P2,f2]=kthspec(y2,4,fs);
[P3,f3]=kthspec(y3,4,fs);
[P4,f4]=kthspec(y4,8,fs);
[P5,f5]=kthspec(y5,8,fs);
[P6,f6]=kthspec(y6,4,fs);
[P7,f7]=kthspec(y7,4,fs);
[P8,f8]=kthspec(y8,2,fs);
[P9,f9]=kthspec(y9,4,fs);
subplot(3,3,1)
plot(3,t(f1,10*log10(P1))
title('BPSK平方谱');

subplot(33,2)
plot(f2,10*log10(P2))
title('QPSK平方谱');

subplot(3,3,3)
plot(f3,10*log10(P3))
title('OQPSK四次方谱');

subplot(3,3,4)
plot(f4,10*log10(P4))
title('pi4QPSK八次方谱');

subplot(3,3,5)
plot(f5,10*log10(P5))
title('8PSK平方谱');

subplot(3,3,6)
plot(f6,10*log10(P6))
title('16QAM平方谱');

subplot(3,3,7)
plot(f7,10*log10(P7))
title('32QAM平方谱');

subplot(3,3,8)
plot(f8,10*log10(P8))
title('64QAM平方谱');

subplot(3,3,9)
plot(f8,10*log10(P8))
title('128QAM平方谱');