源代码:
%2ASK的调制和解调
clear ;
m=[1 1 1000101101];
Lm= length (m);
F =200;
f =800;
A=1;
Q=fF;
if Q> =3
M=500;
t= (0:M-1)M/f;
carry1 =repmat(Asin(2 pi ft),1,Q);
Lcarry1 = length (carry1);
carry2 = kron (ones (size (m)),carry1);
载波信号
ask =kron (m,carry1);
调制后的信号
N= length (ask);
tau= (0:N-1)/(M-1)/f;
时间
Tmin =min(tau);
最小时刻
Tmax = max (tau);
最大时刻
T = ones (size (carry1));
一个数字信号’1’
dsig =kron (m,T);
数字信号波形
figure
subplot(3,1,1);
plot (tau,dsig);
title(‘数字信号)
axis([Tmin Tmax -0.21.21)
各设置坐标范围
subplot(3,1,2);
plot (tau,carry2);
子图分割
title (载波波形’)
画出载波波形
axis ([Tmin Tmax -1.2 A1.2A])
设置坐标范围
subplot(3,1,3);
plot (tau,ask);
各子图分割
tite(经过2ASK调制后的波形”合的
画出调制后的波形
grid on
gxe/iin Tmx -1.2A1.2.A))
L5 sig_ _mul=ask. carry2;
figure
subplot(4.1,1 );plot (tau,sig. jmul);
子图分割
title('信号相乘后的波形’)
画出信号相乘后的波形
axis([Tmin Tmax -0.2 1.21)
(ord,omega_c]=buttord(4pif0.6.4piE0.8,2,30.‘e0;
[num,den]= butter (Ord,omega c,‘s’);
h=tf (num,den);
获得滤波器传递函数滤波
x=Isim(h,sig Jmul,tau);
运用模拟滤波器对信号进行滤波
subplot(4,1,2);
子图分割
plot (tau,x);
画出滤波后的滤形
title (滤波后的波形’)
axis(ITmin Tmax -0.3 0.81);
设置坐标范围
th=0,25;
抽样判决的网值设置
t judge= (0:Lm -1)* Lcarryl +Lcarry12;
抽样判决点的选取
y= (x(化_judge)’;
各抽样判决时刻时的信号值
y_ judge=1* (Yy> =th)+0* (y< =th);
抽样判决信号值的0阶保持
y_value = kron y. _judge ,ones (size (r1));1)抽样判决后的数字信号波形
n_tau=tau+0.5/E;
抽样判决后的信号对应的时间
subplot(4,1,3);
子图分割
plotin tau,y. value);
画出抽样判决后的数字信号波形
title(抽样判决后的数字信号波形’)
axis(min(n tau)max(n. _tau) -0.21.2]1
设置坐标范围
subp1ot(4,1,4);
子图分割
plotitau,dsig);
原信号波形与解调后的信号作对比
title’原始信号波形与2ASK解调后的信号作对比’)号设置坐标范围
axis([TminTmax -0.21.21)
end