bpsk15.m

语音端点检测

a=1; 
fc=2; 
ns=8; 
n=500; 
ts=1; 
T=300 
dt=ts/fc/ns; 
t=0:dt:n*ts-dt; 
lt=length(t); 
d=sign(randn(1,n)); 
dd=sigexpand((d+1)/2,fc*ns); 
gt=ones(1,fc*ns);a=1; 
fc=2; 
ns=8; 
n=500; 
ts=1; 
T=300 
dt=ts/fc/ns; 
t=0:dt:n*ts-dt; 
lt=length(t); 
d=sign(randn(1,n)); 
dd=sigexpand((d+1)/2,fc*ns); 
gt=ones(1,fc*ns); 
figure(1) 
subplot(2,2,1); 
dnrz=conv(dd,gt); 
plot(t,dnrz(1:length(t))); 
axis([0 10 0 1.2]); 
ylabel('输入信号'); 
subplot(2,2,2); 
[f,dnrzf]=t2f(t,dnrz); 
plot(f,10*log10(abs(dnrzf).^2/T)); 
axis([-2 2 -50 10]); 
ylabel('输入信号功率谱密度'); 
ht=a*cos(2*pi*fc*t); 
s_2ask=dnrz(1:lt).*ht; 
subplot(2,2,3) 
plot(t,s_2ask); 
axis([0 10 -1.2 1.2]); 
ylabel('ask'); 
[f,s_2askf]=t2f(t,s_2ask); 
subplot(2,2,4) 
plot(f,10*log10(abs(s_2askf).^2/T)); 
axis([-fc-4 fc+4 -50 10]); 
ylabel('ask功率谱密度'); 
figure(2) 
subplot(2,2,1); 
dnrz=conv(dd,gt); 
plot(t,dnrz(1:length(t))); 
axis([0 10 0 1.2]); 
ylabel('输入信号'); 
subplot(2,2,2); 
[f,dnrzf]=t2f(t,dnrz); 
plot(f,10*log10(abs(dnrzf).^2/T)); 
axis([-2 2 -50 10]); 
ylabel('输入信号功率谱密度'); 
d_2psk=2*dnrz-1; 
s_2psk=d_2psk(1:lt).*ht; 
subplot(2,2,3) 
plot(t,s_2psk); 
axis([0 10 -1.2 1.2]); 
ylabel('2psk'); 
subplot(2,2,4) 
[f,s_2pskf]=t2f(t,s_2psk); 
plot(f,10*log10(abs(s_2pskf).^2/T)); 
axis([-fc-4 fc+4 -50 10]); 
ylabel('2psk功率谱密度'); 
sd_2fsk=2*dnrz-1; 
s_2fsk=a*cos(2*pi*fc*t+2*pi*sd_2fsk(1:length(t)).*t); 
figure(3) 
subplot(2,2,1); 
dnrz=conv(dd,gt); 
plot(t,dnrz(1:length(t))); 
axis([0 10 0 1.2]); 
ylabel('输入信号'); 
subplot(2,2,2); 
[f,dnrzf]=t2f(t,dnrz); 
plot(f,10*log10(abs(dnrzf).^2/T)); 
axis([-2 2 -50 10]); 
ylabel('输入信号功率谱密度'); 
subplot(2,2,3) 
plot(t,s_2fsk); 
axis([0 10 -1.2 1.2]); 
ylabel('2fsk'); 
subplot(2,2,4) 
[f,s_2fskf]=t2f(t,s_2fsk); 
plot(f,10*log10(abs(s_2fskf).^2/T)); 
axis([-fc-4 fc+4 -50 10]); 
ylabel('2fsk功率谱密度');