fdma1.m

利用mablab编写的m文件

%================================BEGIN====================================%
f0=1;%KHz
f1=10;
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]; %时域横坐标
% t=[-1:0.001:1];
% f=[-512*0.1:0.1:0.1*511];
x1=sin(2*pi*f0*t);%输入的正弦信号
plot(t,x1);
axis([-0.5 0.5 -2 2]);
F1=t2f(x1); %F1是x1的傅氏变换
a=f2t(F1); %a是F1的傅氏反变换
a=real(a);
as=abs(F1);

pa=size(f);
pb=size(as);
 
 x=sin(2*pi*f1*t);%载波信号
 figure;
plot(t,x);
axis([-1 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([-4 4 0 1]);
subplot(2,1,2) %输出x的频谱 
plot(f,bs,'r');
axis([-100 100 0 1]);
 
 y1=x1.*x;
 figure;
 plot(t,y1);
 axis([-1 1 -2 2])
F2=t2f(y1); 
c=f2t(F2); 
c=real(c);
cs=abs(F2);
figure;
plot(f,cs,'r');%输出调制后的频谱
axis([-30 30 0 1]);