matlab.m

mpsk生成代码,各种方式生成的代码

%-----------------mpsk生成程序--------------------------------%
clear all;close all;
clc;
M=2;%bpsk
%-------------------------------------------------------------%
%M=4;%Qpsk
%M=8;%8psk
%M=16;%16psk
%--------------------------------------------------------------%
fs=1;%归一化采样频率
fc=0.2;%载波频率即中心频率
Tc=1/fs;%载波周期
fo=0.02;%码元速率
To=1/fo;%码元宽度
N=100;%码元个数
dt=1/fs;%采样间隔
t=0:dt:(N*To-dt);%时间长度
mm=length(t);%总采样点数
d=ceil(rand(1,N)*M);%码元个数个M进制数
for k=1:N
    for m=1:50
       p(50*(k-1)+m)=d(k);%采样总点数个M进制数
   end;
end;
A=10;%幅度
Q=ones(1,length(p));
B=A*sin((2*pi*fc*t+2*pi*(p-Q)/M));
%-----------------------------------------------------------------------%
%B=A*cos((2*pi*fc*t+2*pi*(p-Q)/M));%初相位不同的调制
%B=A*exp(j*(2*pi*fc*t+2*pi*(p-Q)/M));%不同生成方法
%-----------------------------------------------------------------------%
signal=real(B);%MPSK调制信号
%-----------------------------------------------------------------------%
%-----------调用函数法生成-----------------------------------------------%
%p=randint(1,mm,M);
%signal=pskmod(p,M);
%-----------------------------------------------------------------------%
figure(1);
plot(t(1:5000),signal(1:5000));
%plot(t(1:5000)/350,signal(1:5000));%公频为70hz的mpsk仿真波形
title('BPSK signal');
xlabel('time');ylabel('magnitude');