软件无线电实验matlab程序Word文档格式.doc
《软件无线电实验matlab程序Word文档格式.doc》由会员分享，可在线阅读，更多相关《软件无线电实验matlab程序Word文档格式.doc（4页珍藏版）》请在冰点文库上搜索。
%频域信号
Sn=abs（Sn_fft）;
%figures
figure
（1）;
plot（t,sn）;
title（'
s（n）的时域波形'
）;
gridon;
xlabel（'
t'
xlim（[00.021]）;
ylabel（'
幅度'
figure
（2）;
xax=[-N/2:
N/2-1]/N*fs/1000;
%将f轴单位变为kHz
plot（xax,20*log10（fftshift（Sn）））;
s（n）的频域波形'
f（kHz）'
幅度（dB）'
ylim（[-200100]）;
（2）、SDR带通采样理论
f=1.3e3;
%基频
fo=100e6;
%载波频率
fs=4e3;
%采样频率
t=0:
0.001:
2;
st=cos（2*pi*（f+fo）*n/fs）;
%带通抽样信号
St1=fft（st）;
St=abs（St1）;
st0=cos（2*pi*（f+fo）*t）;
St0=fft（st0）;
St2=abs（St0）;
plot（n,st）;
带通抽样信号时域图像'
n'
xlim（[01026]）;
plot（t,st0）;
带通原始信号时域图像'
figure（3）;
xax=[-N/2:
plot（xax,20*log10（fftshift（St）））;
带通信号抽样后频域图像'
（3）、频率调制信号
%频率调制
fs=1e3;
%抽样频率
ts=1/fs;
%采样率
n=0:
1:
N-1;
t=n/fs;
kf=100.1;
A=100;
Aw=10;
fc=10000;
fm=10;
%调制信号频率
mf=kf*Aw/（2*pi*fm）;
s1=A*cos（2*pi*fc*t）.*cos（mf*sin（2*pi*fm*t））;
s2=A*sin（2*pi*fc*t）.*sin（mf*sin（2*pi*fm*t））;
s=s1-s2;
%频率调制信号
S0=fft（s）;
S1=fftshift（S0）;
S_0=abs（S1）;
plot（t,s）;
频率调制信号时域图像'
xlim（[0,0.5]）;
t（s）'
plot（xax,20*log10（S_0））;
频率调制信号频域图像'
（4）、幅度调制信号、DSB信号、USB信号
clearall;
clc;
closeall;
N=10000;
N;
ma=0.5;
fc=200;
x=cos（2*pi*fc*t）;
%生成载波
fm=1;
m=cos（2*pi*fm*t）;
%基带信号
s=（1+ma*m）.*x;
%幅度调制信号
s1=m.*x;
%双边带信号DSB
s2=cos（2*pi*（fc+fm）*t）;
%SSB信号中的USB
S1=fft（s）;
S_1=abs（S1）;
S2=fft（s1）;
S_2=abs（S2）;
S3=fft（s2）;
S_3=abs（S3）;
幅度调制信号时域图像'
xlim（[01.5]）;
plot（xax,20*log10（fftshift（S_1）））;
幅度调制信号频域图像'
xlim（[0.190.21]）;
plot（t,s1）;
双边带调制信号时域图像'
figure（4）;
xax1=[-N/2:
plot（xax1,20*log10（fftshift（S_2）））;
双边带调制信号频域图像'
figure（5）;
plot（t,s2）;
USB调制信号时域图像'
xlim（[01.0]）;
figure（6）;
plot（xax1,20*log10（fftshift（S_3）））;
USB调制信号频域图像'
（5）、2ASK信号
%2ASK信号的时域和频域分析
%清除一切。
m=[111000101101];
%数字信号的序列（自己随便产生的）
Lm=length（m）;
%序列的长度
F=200;
%数字信号的带宽（'
1'
或'
0'
的带宽是200Hz）
f=800;
%正弦载波信号的频率
A=1;
%载波的幅度
Q=f/F;
%频率比，即一个码元宽度中的正弦周期个数
M=500;
%一个正弦周期内的采样点数
t=（0:
M-1）/M/f;
%一个正弦信号周期内的时间
carry1=repmat（A*sin（2*pi*f*t）,1,Q）;
%一个码元宽度内的正弦载波信号
Lcarry1=length（carry1）;
%一个码元宽度内的信号长度
carry2=kron（ones（size（m））,carry1）;
%载波信号
ask=kron（m,carry1）;
%调制后的ASK信号
N=length（ask）;
%ASK信号的长度
tau=（0:
N-1）/（M-1）/f;
%时间
Tmin=min（tau）;
%最小时间
Tmax=max（tau）;
%最大时间
T=ones（size（carry1））;
%一个数字信号‘1’
dsig=kron（m,T）;
Ask=fft（ask）;
Ask0=abs（Ask）;
%images
plot（tau,dsig）;
axis（[TminTmax-0.21.2]）;
数字信号序列的图像'
plot（tau,carry2）;
xlim（[00.06]）;
载波时域图像'
plot（tau,ask）;
axis（[TminTmax-1.2*A1.2*A]）;
2ASK信号的时域图像'
N/2-1]/N*Lm*F;
plot（xax,20*log10（fftshift（Ask0）））;
xlim（[-300300]）;
2ASK信号的频域图像'
（6）、16路信道化发射机
a=1;
%幅度
fs=50;
%内插前采样频率
k=1:
16;
%16路信号
I=16;
%内插系数
fk=1:
sk=zeros（16,N）;
%创建16行N列的零矩阵用来存放16路时域信号
Sk=zeros（16,N）;
%创建16行N列的零矩阵用来存放16路频域信号
fork=1:
16
sk（k,:
）=a*（1+0.5*exp（j*2*pi*k*fk/fs））;
%输入16路时域信号
end
Sk（k,:
）=fft（sk（k,:
））;
%16路信号频域
N/2-1]/N*fs;
plot（xax,20*log10（abs（fftshift（Sk（:
1:
N）））））;
16原始复信号时域叠加图'
f（Hz）'
xlim（[-525]）;
h0=fir1（265,0.5）;
HW=fft（h0,N）;
plot（xax1,10*log10（abs（fftshift（HW））））;
滤波器响应'
%I倍内插，即16倍内插，此时fs变大
N_1=N*16;
st=zeros（16,N_1）;
St=zeros（16,N_1）;
forg=0:
N-1
st（k,g*I+k）=sk（k,g+1）;
end
sum=zeros（1,N_1）;
sum=sum+st（k,1:
N_1）;
St（k,:
）=fft（sum（1,1:
N_1））;
St（k,1:
N）=St（k,1:
N）.*HW;
xax3=[-N_1/2:
N_1/2-1]/N_1*fs*16;
plot（xax3,sum（1,:
内插后信号时域图像'
xlim（[-6060]）;
xax2=[-N/2:
N/2-1]/N*fs*16;
plot（xax2,db（abs（fftshift（St（:
内插后信号频域图像'
xlim（[-10200]）;
升级会员 