完整word版软件无线电习题Word下载.docx
- 文档编号:17584378
- 上传时间:2022-12-07
- 格式:DOCX
- 页数:12
- 大小:46.32KB
完整word版软件无线电习题Word下载.docx
《完整word版软件无线电习题Word下载.docx》由会员分享,可在线阅读,更多相关《完整word版软件无线电习题Word下载.docx(12页珍藏版)》请在冰豆网上搜索。
268.5
300
三、对一圆形阵列的波束形成进行仿真,假设发射信号载频为0.8GHz,圆形阵列半径为0.8米,在圆周上均匀布置28个阵元。
1.画出指向10度的方向图;
2.如果目标在0度,有一不相干的干扰信号在-20度,干扰噪声功率比为35dB。
请用自适应波束形成方法画出方向图;
3.采用旁瓣对消的方法(选取两个阵元作为辅助天线),计算对消比
第一题源程序:
%=====================================================================%
%软件无线电作业题1
%对调制信号进行正交分解、滤波、抽取、解调
%====================================================================%
clc;
closeall;
clearall;
fs=40e6;
%信号时钟
Ts=1/fs;
f0=10.7e6;
%载频
f1=1000;
f2=2000;
D1=10;
%CIC滤波抽取因子
D2=8;
%3级半带滤波抽取因子
D3=5;
%低通滤波抽取因子
%===================================================================%
%CIC滤波器设计
S1_cic=ones(1,D1);
%一级CIC
[H1,F1]=freqz(S1_cic,1,1024,fs);
S2_cic=conv(S1_cic,S1_cic);
%两级CIC
[H2,F2]=freqz(S2_cic,1,1024,fs);
S3_cic=conv(S2_cic,S1_cic);
%三级CIC
[H3,F3]=freqz(S3_cic,1,1024,fs);
figure;
plot(F1/(fs/2),20*log10(abs(H1))-max(20*log10(abs(H1))),'
b'
),grid;
holdon;
plot(F2/(fs/2),20*log10(abs(H2))-max(20*log10(abs(H2))),'
g'
);
plot(F3/(fs/2),20*log10(abs(H3))-max(20*log10(abs(H3))),'
m'
xlabel('
\fontsize{12}\bf归一化频率(\times\pirad/sample)'
ylabel('
\fontsize{12}\bf幅值(dB)'
),title('
\fontsize{12}\bfCICF幅频响应'
boxon;
legend('
\bf单级CIC'
'
\bf两级CIC'
\bf三级CIC'
Location'
SouthWest'
set(gca,'
FontWeight'
bold'
FontSize'
12);
holdoff;
%=======================================================================%
%HB滤波器设计
B2=firhalfband(8,blackman(9));
%9阶HBF
B2=conv(conv(B2,B2),B2);
%三级级联
[H4,F4]=freqz(B2,1,1024,fs);
figure
(2),plot(F4/(fs/2),20*log10(abs(H4))),grid;
title('
\fontsize{12}\bfHBF幅频响应'
%FIRF设计
fs1=fs/(8*10);
[N,Fo,Ao,W]=firpmord([20e335e3],[10],[10^(-100/20)10^(-100/20)],fs1);
%%低通滤波器
b1=firpm(127,Fo,Ao,W);
[H5,F5]=freqz(b1,1,1024,fs1);
figure(3),plot(F5/(fs1/2),20*log10(abs(H5))),grid;
\fontsize{12}\bfFIFR幅频响应'
%鉴频FIR设计
fs2=fs1/5;
[N,Fo,Ao,W]=firpmord([20e330e3],[10],[10^(-80/20)10^(-80/20)],fs2);
b2=firpm(63,Fo,Ao,W);
[H6,F6]=freqz(b2,1,1024,fs2);
figure(4),plot(F6/(fs2/2),20*log10(abs(H6))),grid;
\fontsize{12}\bf鉴频FIFR幅频响应'
%========================================================================%
%滤波器总的幅频响应
HH=2.*H3.*H4.*H5;
figure(5),plot((1:
1024)*fs/1024,20*log10(abs(HH))),grid;
\fontsize{12}\bf频率/Hz'
\fontsize{12}\bf总的滤波器幅频响应'
%=========================================================================%
%数据处理,给出每级数据经滤波抽取后的结果
N=400000;
n=0:
400000-1;
phi=0.7*pi*sin(2*pi*f1*n*Ts+pi/6)+0.7*pi*sin(2*pi*f2*n*Ts+pi/3);
Sn=cos(2*pi*f0*n*Ts+phi);
%+1/sqrt(20)*randn(1,N);
%%%SNR=20dB%%%
I1=Sn.*cos(2*pi*f0*n*Ts);
Q1=Sn.*sin(2*pi*f0*n*Ts);
%正交变换
II2=2*filter(S3_cic,1,I1);
QQ2=2*filter(S3_cic,1,Q1);
%CIC滤波
I2=II2(1:
D1:
N);
Q2=QQ2(1:
%10倍抽取
II3=filter(B2,1,I2);
QQ3=filter(B2,1,Q2);
%HB滤波
I31=II3(1:
2:
N/D1);
Q31=QQ3(1:
%2倍抽取
I32=filter(B2,1,I31);
Q32=filter(B2,1,Q31);
I33=I32(1:
N/(2*D1));
Q33=Q32(1:
I34=filter(B2,1,I33);
Q34=filter(B2,1,Q33);
I3=I34(1:
N/(2*2*(D1)));
Q3=Q34(1:
I41=filter(b1,1,I3);
Q41=filter(b1,1,Q3);
%低通滤波
I4=I41(1:
D3:
N/(D1*D2));
Q4=Q41(1:
%5倍抽取
%%%%%%%画图显示%%%%%%%%%
figure(6);
plot((0:
Ts:
1000*Ts),Sn(1:
1001)),gridon;
\fontsize{12}\bf时间(s)'
),ylabel('
\fontsize{12}\bf幅值(v)'
\fontsize{12}\bf原始信号'
figure(7);
plot((1:
40000),I1(1:
40000),'
r-'
(1:
40000),Q1(1:
b-'
),legend('
\bfI1'
\bfQ1'
\fontsize{12}\bfI1与Q1路时域信号'
figure(8);
plot(1:
length(I2),I2,'
1:
length(Q2),Q2,'
\bfI2'
\bfQ2'
\fontsize{12}\bfI2与Q2路时域信号'
figure(9);
length(I3),I3,'
length(Q3),Q3,'
\bfI3'
\bfQ3'
\fontsize{12}\bfI3与Q3路时域信号'
figure(10);
length(I4),I4,'
length(Q4),Q4,'
\bfI4'
\bfQ4'
\fontsize{12}\bfI4与Q4路时域信号'
figure(11);
%画正交原图,看IQ通道是否正交
subplot(2,2,1),plot(I1,Q1),axisequal,grid;
\fontsize{12}\bf同相分量I1'
\fontsize{12}\bf正交分量Q1'
title('
\fontsize{12}\bfI1与Q1正交圆图'
subplot(2,2,2),plot(I2,Q2),axisequal,grid;
\fontsize{12}\bf同相分量I2'
\fontsize{12}\bf正交分量Q2'
\fontsize{12}\bfI2与Q2正交圆图'
subplot(2,2,3),plot(I3,Q3),axisequal,grid;
\fontsize{12}\bf同相分量I3'
\fontsize{12}\bf正交分量Q3'
\fontsize{12}\bfI3与Q3正交圆图'
subplot(2,2,4),plot(I4,Q4),axisequal,grid;
\fontsize{12}\bf同相分量I4'
\fontsize{12}\bf正交分量Q4'
\fontsize{12}\bfI4与Q4正交圆图'
%============================================================================%
%%%%鉴频%%%%
sig=2*I4+j*2*Q4;
an=abs(sig);
Pn0=atan2(Q4,I4);
Pn=unwrap(Pn0);
B5=[1,0,0,-1];
Fp=filter(B5,1,Pn);
Fn=filter(b2,1,Fp);
%鉴频低通滤波
ff=-diff(phi);
%理论值求解
ff=ff(1:
100:
length(ff));
figure(12);
subplot(2,1,1),plot(Fn,'
r'
\fontsize{12}\bf解调结果'
subplot(2,1,2),plot(ff,'
\fontsize{12}\bf理论解调结果'
第二题源程序:
%=================================================================%
%软件无线电第二题
%信道化接收
%================================================================%
clearall;
N=1120;
%低通滤波器长度
M=70;
%划分通道数
K=16;
%多相滤波器长度
fs=210e6;
%采样频率
Bw=85e6;
%信号带宽Bw=fh-fl
f0=262.5e6;
%信号中心频率
fl=f0-Bw/2;
%最低工作频率
fh=f0+Bw/2;
%最高工作频率
T=100e-6;
%LFM信号脉宽
Bs=3e6;
%每个信道带宽
N-1;
%------------------------输入信号------------------------
f1=242.5e6;
%chirp1中心频率
f2=267.5e6;
%chirp2中心频率
SNR=10;
%信噪比
t=-T/2:
1/fs:
T/2-1/fs;
w1=2*pi*f1;
BW1=600e3;
%chirp1带宽
w2=2*pi*f2;
BW2=300e3;
z1t=cos(w1*t+BW1/T*pi*t.*t);
%chirp1
z2t=cos(w2*t+BW2/T*pi*t.*t);
%chirp2
zt=z1t+z2t+(10.^(-10/20))*randn(1,length(t));
%信号+噪声
DFzt=fft(zt);
DFz1t=fft(z1t);
DFz2t=fft(z2t);
%信号的频谱
figure
(1)
subplot(3,2,1);
plot(0:
(fs*T-1),zt);
时域'
zt=chirp1+chirp2+noise'
subplot(3,2,2);
plot((-0.5:
1/(fs*T):
0.5-1/(fs*T))*fs,abs(DFzt));
频域'
subplot(3,2,3);
(fs*T-1),z1t);
chirp1'
subplot(3,2,4);
0.5-1/(fs*T))*fs,abs(DFz1t));
subplot(3,2,5);
(fs*T-1),z2t);
chirp2'
subplot(3,2,6);
0.5-1/(fs*T))*fs,abs(DFz2t));
%---------------------------设计低通滤波器------------------------------
f=[Bs,1.5*Bs]/fs;
m1=[1,0];
Rp=1;
Rs=40;
dat1=(10^(Rp/20)-1)/(10^(Rp/20)+1);
dat2=10^(-Rs/20);
rip=[dat1,dat2];
[MM,ff,mm,ww]=remezord(f,m1,rip);
hh=remez(N-1,ff,mm,ww);
figure
(2)
freqz(hh);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%----------------直接低通滤波---------------
fork=0:
M-1
z1h=conv(z1t,hh.*exp(-j*2*pi*k*n/M));
z2h=conv(z2t,hh.*exp(-j*2*pi*k*n/M));
zh=conv(zt,hh.*exp(-j*2*pi*k*n/M));
m=0:
size(t)-1;
W=exp(-j*2*pi/M*m*k);
z1w=z1h.*W;
z2w=z2h.*W;
zw=zh.*W;
z1k_temp=reshape([z1w,zeros(1,M-rem(length(z1w),M))],M,[]);
z2k_temp=reshape([z2w,zeros(1,M-rem(length(z2w),M))],M,[]);
zk_temp=reshape([zw,zeros(1,M-rem(length(z1w),M))],M,[]);
z1k(k+1,:
)=z1k_temp(k+1,:
z2k(k+1,:
)=z2k_temp(k+1,:
zk(k+1,:
)=zk_temp(k+1,:
ifk<
M/2
figure(3)
subplot(7,5,k+1);
plot([0:
length(zk_temp)-1]*pi/length(zk_temp),real(zk(k+1,:
)));
temp1=['
\fontsize{12}\bf通道'
num2str(k)];
axis([0pi-.6.6]);
title(temp1);
set(gca,'
figure(4)
length(zk_temp)-1]*pi/length(zk_temp),abs(fft(zk(k+1,:
))));
axis([0pi035]);
end
end
figure(5)
plot(real(zk(12,:
)))
通道11波形'
zk_fft=abs(fft(zk(12,:
)))/max(abs(fft(zk(12,:
figure(6)
plot(((0:
length(zk_fft)-1)*Bs/length(zk_fft)),fftshift(zk_fft))
频率'
归一化幅度'
通道11信号频谱'
figure(7)
plot(abs(xcorr(zk(12,:
),zk(12,:
))))
通道11脉压结果'
figure(8)
plot(real(zk(20,:
通道19波形'
zk_fft=abs(fft(zk(20,:
)))/max(abs(fft(zk(20,:
figure(9)
通道19信号频谱'
figure(10)
plot(abs(xcorr(zk(20,:
),zk(20,:
通道19脉压结果'
第三题源程序:
%========
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 完整 word 软件 无线电 习题