binarymod.m

ook,2psk,2fsk的典型波形和频谱

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功率谱密度');