DPSK调制解调的Matlab程序

%DPSK调制解调程序清单%仿真DPSK信号的产生及相干解调过程%将成形滤波器系数量化为10比特后，写入tra_lpf.coe文件中%将接收滤波器系数量化为10比特后，写入rec_lpf.coe文件中%将DPSK已调数据量化为8比特后，写入Dpsk.txt文件中ps=1*106; %码速率为1MHza=0.8; %成形滤波器系数为0.8B=(1+a)*ps; %中频信号处理带宽Fs=8*106; %采样速率为8MHzfc=2*106; %载波频率为2MHzN=20000; %仿真数据的长度t=0:1/Fs:(N*Fs/ps-1)/Fs;%产生长度为N,频率为fs的时间序列s=randint(N,1,2); %产生随机数据作为原始数据,%并将绝对码变换为相对码ds=ones(1,N);for i=2:N if s(i)=1 ds(i)=-ds(i-1); else ds(i)=ds(i-1); endend%对相对码数据以Fs频率采样Ads=upsample(ds,Fs/ps);%设计平方升余弦滤波器n_T=-2 2;rate=Fs/ps;T=1;Shape_b = rcosfir(a,n_T,rate,T);%figure(4);freqz(Shape_b)%对采样后的数据进行升余弦滤波;rcos_Ads=filter(Shape_b,1,Ads);%产生载频信号f0=sin(2*pi*fc*t); %产生DPSK已调信号dpsk=rcos_Ads.*f0; %与相干载波相乘，实现相干解调demod_mult=dpsk.*f0;%设计接收端低通滤波器fc=ps 3.1*106; %过渡带mag=1 0; %窗函数的理想滤波器幅度dev=0.01 0.01; %纹波n,wn,beta,ftype=kaiserord(fc,mag,dev,Fs) %获取凯塞窗参数fpm=0 fc(1)*2/Fs fc(2)*2/Fs 1; %firpm函数的频段向量magpm=1 1 0 0; %firpm函数的幅值向量rec_lpf=firpm(n,fpm,magpm); %firpm函数返回的最优滤波器系数%对乘法运算后的数据进行低通滤波，输出解调后的基带信号demod_lpf=filter(rec_lpf,1,demod_mult);figure(1)%绘制成形滤波后信号频谱、DPSK信号频谱、DPSK信号时域波形m_rcos_Ads=20*log10(abs(fft(rcos_Ads,1024);m_rcos_Ads=m_rcos_Ads-max(m_rcos_Ads);m_dpsk=20*log10(abs(fft(dpsk,1024);m_dpsk=m_dpsk-max(m_dpsk);%设置幅频响应的横坐标单位为MHzx_f=0:(Fs/length(m_dpsk):Fs/2;x_f=x_f/106;%只显示正频率部分的幅频响应mrcos_Ads=m_rcos_Ads(1:length(x_f);mdpsk=m_dpsk(1:length(x_f);%设置时域波表的横坐标单位为usLen=100;%设置时域波形显示的点数x_t=1:Len;%产生长度为Len的时间序列x_t=x_t/Fs*106;%显示所需的频谱及时域波形subplot(311); plot(x_f,mrcos_Ads);legend(成形滤波后信号频谱);xlabel(频率(MHz);ylabel(幅度(dB);grid on;subplot(312); plot(x_f,mdpsk); legend(DPSK已调信号频谱);xlabel(频率(MHz);ylabel(幅度(dB);grid on;subplot(313); plot(x_t,dpsk(101:Len+100);grid on;legend(DPSK时域信号波形);xlabel(时间(us);ylabel(幅度(V);grid on;figure(2)%对相对码重复Fs/ps倍采样，便于绘图比较s_Ads=rectpulse(ds,Fs/ps);Len=500;%设置时域波形显示的点数x_t=1:Len;%产生长度为Len的时间序列x_t=x_t/Fs*106;%绘制DPSK解调前后时域波形delay=18;%为便于观察，对解调后的基带波形超前显示delay个点。plot(x_t,s_Ads(1:Len)/2,-,x_t,demod_lpf(de