一、MSK的附加相位图程序.docVIP

一、MSK的附加相位图程序 MATLAB程序如下： data_len =7; %码元个数 sample_number = 10; %采样个数 Rb = 2400; %码元速率 fc = 9600; %载波频率 data=[-1,1,1,1,-1,1,-1]; for i = 1:sample_number data_sample(i:sample_number:data_len*sample_number)= data; end phase = zeros(1,data_len*sample_number); phase(1) = data_sample(1) * pi/2/sample_number; for i = 2:data_len*sample_number phase(i)=phase(i-1)+data_sample(i-1)*pi/2/sample_number; end plot(phase) title(MSK附加相位图) grid 结果如图 二、MSK的功率谱 其matlab程序代码如下： fb=2400; Tb=1/fb; fc=9600; f1=4800;f2=6400; fs=fc/fb; f=-3*fb:0.001*fb:3*fb; Mf=(1/16)*Tb*((cos(2*pi*(f-fs)*Tb).^2))./(pi.^2*(1-(4*(f-fs)*Tb).^2)); Ff=(1/16)*Tb*(abs((sin(pi*(f+f1/fb)*Tb))./(pi*(f+f1/fb)))).^2+... (1/16)*Tb*(abs((sin(pi*(f-f1/fb)*Tb))./(pi*(f-f1/fb)))).^2+... (1/16)*Tb*(abs((sin(pi*(f+f2/fb)*Tb))./(pi*(f+f2/fb)))).^2+... (1/16)*Tb*(abs((sin(pi*(f-f2/fb)*Tb))./(pi*(f-f2/fb)))).^2+... (1/16)*(imp(f,-f1)+imp(f,f1)+imp(f,-f2)+imp(f,f2)) plot(f,Mf,b-);title(MSK的功率谱);ylabel(幅度); hold on plot(f,Ff*10.^5,r.) hold off 结果如图： 三、MSK调制解调 （1）时延函数 function out = delay(data,n,sample_number) out = zeros(1,length(data)); out(n*sample_number+1:length(data)) = data(1:length(data)-n*sample_number); （2）差分函数 function [data_diff] = difference(data) %差分编码 %************************************************************************** %data 输入信号 %data_diff 差分编码后信号 %************************************************************************** %-------------------------------------------------------------------------- data_diff = zeros(1,length(data)); data_diff(1) = 1 * data(1); %1为差分编码的初始参考值 for i = 2:length(data) data_diff(i) = data_diff(i-1) * data(i); end (3)基带调制函数 function [signal_out,I_out,Q_out] = mod_msk(data,data_len,sample_number,Rb) %MSK基带调制 %************************************************************************** % data 调制信号 % data_len 码元个数 % sample_number 每个码元采样点数 % Rb