DPSK调制解调的Matlab程序

1、Matlab程序Document serial number【UU89WT-UU98YT-UU8CB-UUUT-UUT108】%DPSK调制解调程序清单 %仿真DPSK信号的产生及相干解调过程%将成形滤波器系数量化为10比特后，写入tra_lpf.coe文件中%将接收滤波器系数量化为10比特后，写入rec_lpf.coe文件中 %将DPSK已调数据量化为8比特后，写入Dpsk. txt文件中 ps=l*10飞；码速率为1MHza=0.8;%成形滤波器系数为0.8B=(1+a)*ps; %中频信号处理带宽Fs=8*10飞；采样速率为8MHzfc=2*10飞；载波频率为2MHzN=20000; %

2、仿真数据的长度t=0:1/Fs:(N*Fs/ps-1)/Fs;%产生长度为N,频率为fs的时间序列 s=randint(N,1,2);%产生随机数据作为原始数据,%并将绝对码变换为相对码ds=ones(1,N);for i=2:Nif s(i)=1ds(i)=-ds(i-1);elseds(i)=ds(i-1);end end%对相对码数据以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)%对采样后的数

3、据进行升余弦滤波;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*10八6;%过渡带mag=1 0;%窗函数的理想滤波器幅度dev=0.01 0.01; %纹波 n,wn,beta,ftype=kaiserord(fc,mag,dev,Fs) %获取凯塞窗参数%firpm函数的频段向fpm=0 fc(1)*2/Fs fc(2)*2/Fs %firpm函数的频段向量

4、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

5、_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(成形滤波后信号频谱

6、)；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=l:Len;%产生长度为Len的时间序列x_t二x_t/Fs*106%绘制DPSK解调前后时域波形delay=18;%为便于观察，对解调后的基带波形超前显示delay个点。plot(x_t,s_Ads(1:Len)/2,-,x_t,demod_lpf(