%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(A*sin(2 *pi* f*t),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.2*A]) 设置坐标范围 subplot(3,1,3); plot (tau,ask); 各子图分割 tite(经过2ASK调制后的波形”合的 画出调制后的波形 grid on gxe/iin Tmx -1.2*A1.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(4*pi*f*0.6.4*pi*E*0.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
源代码:
