文章目录
一、要求二、设计三、代码四、结果五、结论
一、要求
画出MSK信号产生的信息、I路、Q路信号对比图以及接收端解调的信息恢复的对比图。
二、设计
1.MSK调制 2.MSK解调
三、代码
clear all
;clc;
N=20
;
T=1
;
fc=10
;
Fs=100
;
bitstream=randi
([0
,1
],1
,N
);
bitstream1=2
*bitstream
-1
;
b0=1
;
for i=1:N
encode_output
(i
)=b0
*bitstream1
(i
);
b0=encode_output
(i
);
end
I=
[];Q=
[];
for i=1:N
if mod
(i
,2
)~=0
I=
[I,encode_output(i)];
else
Q=
[Q,encode_output(i)];
end
end
bit_data=
[];
for i=1:N
bit_data=
[bit_data,encode_output(i)*ones(1,T*Fs)];
end
I_data=
[];Q_data=
[];base_wave=
-T:1
/Fs:T
-1
/Fs
;
for i=1:N
/2
I_data=
[I_data,I(i)*cos(pi*base_wave/(2*T))];
Q_data=
[Q_data,Q(i)*cos(pi*base_wave/(2*T))];
end
% Q delay
number_delay=length
(base_wave
)/2
;
Q_data1=
[zeros(1,number_delay),Q_data(1:length(Q_data)-number_delay)];
% plot
figure
();
t=0:1
/Fs:N
*T
-1
/Fs
;
subplot
(3
,1
,1
)
plot
(t
,bit_data
);legend
('Bitstream')
subplot
(3
,1
,2
)
plot
(t
,I_data
);legend
('I Bitstream')
subplot
(3
,1
,3
)
plot
(t
,Q_data1
);legend
('Q Bitstream')
% carrier signal
bit_t=0:1
/Fs:N
*T
-1
/Fs
;
I_carrier=cos
(2
*pi
*fc*bit_t
);
Q_carrier=cos
(2
*pi
*fc*bit_t
+pi
/2
);
% transmit signal
MSK_signal=I_data
.*I_carrier
+Q_data1
.*Q_carrier
;
snr=1
;
MSK_receive=awgn
(MSK_signal
,snr
);
% demodulate
I_output=MSK_receive
.*I_carrier
;
Hd=myfilter
;
I_filter_ouput=
filter(Hd
,I_output
);
Q_output=MSK_receive
.*Q_carrier
;
Q_filter_ouput=
filter(Hd
,Q_output
);
for i=1:N
/2
if I_filter_ouput
((2
*i
-1
)*number_delay
)>0
I_recover
(i
)=1
;
else
I_recover
(i
)=
-1
;
end
if Q_filter_ouput
(2
*i
*number_delay
)>0
Q_recover
(i
)=1
;
else
Q_recover
(i
)=
-1
;
end
end
bit_recover=
[];
for i=1:N
if mod
(i
,2
)~=0
bit_recover=
[bit_recover,I_recover((i-1)/2+1)];
else
bit_recover=
[bit_recover,Q_recover(i/2)];
end
end
% decode
for i=1:N
if i==1
bit_recover1
(i
)=bit_recover
(i
);
else
bit_recover1
(i
)=bit_recover
(i
)*bit_recover
(i
-1
);
end
end
recover_data=
[];
for i=1:N
recover_data=
[recover_data,bit_recover1(i)*ones(1,T*Fs)];
end
% plot
figure
();
t=0:1
/Fs:N
*T
-1
/Fs
;
bit_stream=
[];
for i=1:N
bit_stream=
[bit_stream,bitstream1(i)*ones(1,T*Fs)];
end
subplot
(2
,1
,2
)
plot
(t
,bit_stream
);legend
('original Bitstream')
subplot
(2
,1
,1
)
plot
(t
,recover_data
);legend
('recover Bitstream')
四、结果
myfilter是通过matlab命令行调用filter designer设计生成的函数
五、结论
在MSK的调制与解调中，产生和接收到的信号能够进行还原，即调制产生的I、Q以及MSK信号在经过解调之后能够还原出原来的信号。尽管过滤器设计不同，但是仍然不影响最后的结果。
