fdma2.m

利用mablab编写的m文件

f0=1;%KHz,输入信号频率
f1=30;%载波频率
N=2^10;
dt=0.001; %ms
df=1/(N*dt); %KHz
T=N*dt; %截短时间
Bs=N*df/2; %系统带宽
f=[-Bs+df/2:df:Bs] %频域横坐标
t=[-T/2+dt/2:dt:T/2]; %时域横坐标
x2=square(2*pi*f0*t);%输入的方波信号
plot(t,x2);
axis([-0.6 0.6 -2 2]);
F1=t2f(x2); %F1是x1的傅氏变换
a=f2t(F1); %a是F1的傅氏反变换
a=real(a);
as=abs(F1);

 x=sin(2*pi*f1*t);%载波信号
 figure;
plot(t,x);
axis([-0.2 0.2 -2 2]);
F=t2f(x); 
b=f2t(F); 
b=real(b);
bs=abs(F);
figure;
subplot(2,1,1) %输出x1的频谱 
plot(f,as,'b');
subplot(2,1,2) %输出x的频谱 
plot(f,bs,'b');
 
 y1=x2.*x;
 figure;
 plot(t,y1);
 axis([-0.5 0.5 -2 2])
F2=t2f(y1); 
c=f2t(F2); 
c=real(c);
cs=abs(F2);
figure;
plot(f,cs,'b');%输出调制后的频谱