tdma.m

利用mablab编写的m文件

%==============================Begin=================================%
t0=0.38;                             %信号持续时间
ts=0.0001;                           %采样时间间隔
f0=25;                               %信号频率
% Fs=1/ts;
% df=0.3;

t1=[0:ts:(t0/4+t0/40)];              %第一个时间段
x1=sin(2*pi*f0*t1);                  %第一个时间段上传输正弦信号
% [X1,x1,df1]=fft_seq(x1,ts,df);
% X1=X1/Fs;
% f=[0:df1:df1*(length(x1)-1)]-Fs/2;

t2=[(t0/4+t0/40):ts:(t0/2-t0/16)];   %第二个时间段
x2=square(2*pi*f0*t2);               %第二个时间段上传输方波信号

t3=[(t0/2-t0/16):ts:t0];             %第三个时间段
x3=sawtooth(2*pi*f0*t3);             %第三个时间段上传输锯齿波信号

figure;                              %原始传输信号图1
subplot(3,1,1);plot(t1,x1);axis([0 0.15  -2 2]);
subplot(3,1,2);plot((t2-t2(1,1)),x2);axis([0 0.1 -2 2]);
subplot(3,1,3);plot((t3-t3(1,1)),x3);axis([0 0.3 -2 2]);

t=[t1,t2,t3];
x=[x1,x2,x3];
figure;                            %时分复用后的信号图
plot(t,x);
axis([0 0.4 -2 2]);

% figure;
% plot(f,abs(fftshift(X1)))

%==================以下是从服用信号中提取出正弦、方波和锯齿波================%
t11=t(1,1:length(t1));
x11=x(1,1:length(x1));
figure;                          %分离出的信号图，并与前面的图1进行比较
subplot(3,1,1);
plot(t11,x11,'bla');
axis([0 .15 -2 2]);

t12=t(1,length(t1):length(t1)+length(t2));
x12=x(1,length(x1):length(x1)+length(x2));
t12=t12-t12(1,1);
subplot(3,1,2);
plot(t12,x12,'bla');
axis([0 .1 -2 2]);

t13=t(1,length(t1)+length(t2):length(t1)+length(t2)+length(t3));
x13=x(1,length(x1)+length(x2):length(x1)+length(x2)+length(x3));
t13=t13-t13(1,1);
subplot(3,1,3);
plot(t13,x13,'bla');
axis([0 .3 -2 2]);
%==============================End=========================================%