tiaoxiangxinhao.m

来自「关于模拟调制解调的源代码」· M 代码 · 共 44 行

clear
echo on
t0=0.25;
tz=0.0005;
fc=200;
kf=50;
fz=1/tz;
t=[0:tz:t0];
df=0.25;
%定义信号序列
m=zeros(1,501);
for i=1:1:125
    m(i)=i;
end
for i=126:1:375
    m(i)=m(125)-i+125;
end
for i=376:1:501
    m(i)=m(375)+i-375;
end
m=m/50;
[M,m,df1]=fftseq(m,tz,df);
M=M/fz;
f=[0:df1:df1*(length(m)-1)]-fz/2;
for i=1:length(t)
    mn(i)=m(i);
end
u=cos(2*pi*fc*t+mn);
[U,u,df1]=fftseq(u,tz,df);
U=U/fz;
pause;       
subplot(2,1,1)
plot(t,m(1:length(t)))
axis([0 0.25-3 3]);           xlabel('时间')
title('信号波形');             subplot(2,1,2)
plot(t,u(1:length(t)))
axis([0 0.15-2.1 2.1]);           xlabel('时间')
%title('调相信号的时域波形')；     pause
subplot(2,1,1)
plot(f,abs(fftshift(M)))
%xlabel('频率')；                 title('信号频谱')
subplot(2,1,2)
plot(f,abs(fftshift(U)))
%title('调相信号的频谱')；     xlabel('频率')