fdma3.m

利用mablab编写的m文件

f0=1;%KHz
f1=50;
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]; %时域横坐标
x3=sawtooth(2*pi*f0*t);%输入的锯齿波信号
plot(t,x3);
 axis([-0.5 0.5 -2 2]);
F1=t2f(x3); %F1是x1的傅氏变换
a=f2t(F1); %a是F1的傅氏反变换
a=real(a);
as=abs(F1);

pa=size(as)
pb=size(f)

 x=sin(2*pi*f1*t);%载波信号
 figure;
plot(t,x);
axis([-0.1 0.1 -2 2]);
F=t2f(x); 
b=f2t(F); 
b=real(b);
bs=abs(F);
figure;
subplot(2,1,1) %输出x1的频谱 
plot(f,as,'b');
axis([-20 20 0 0.5]);
subplot(2,1,2) %输出x的频谱 
plot(f,bs,'b'); 

 y3=x3.*x;
 figure;
 plot(t,y3);
axis([-0.5 0.5 -2 2])
F2=t2f(y3); 
c=f2t(F2); 
c=real(c);
cs=abs(F2);
figure;
plot(f,cs,'b');%输出调制后的频谱
axis([-100 100 0 0.5]);