gaofang.asv

雷达系统建模与仿真的matlab程序

function y=gaofang(s,fs,f0,f1)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%高频放大器
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[M,N]=size(s);
Vgain=3;                                         %高放增益
ws1=(f0-4*f1/3)*2*pi/fs;wp1=(f0-f1/3)*2*pi/fs;
ws2=(f0+7*f1/3)*2*pi/fs;wp2=(f0+4*f1/3)*2*pi/fs;
tr_width=min((wp1-ws1),(ws2-wp2));
M1=ceil(6.2*pi/tr_width)+1;
n=[0:M1-1];
wc1=(ws1+wp1)/2;wc2=(wp2+ws2)/2;
hd=ideal_lp(wc2,M1)-ideal_lp(wc1,M1);
w_ham=(hamming(M1))';
h=hd.*w_ham;
L=length(h);
x=Vgain*x;
z=zeros(M+L-1,N);
for i=1:N
    z(1:M+L-1,i)=conv(h,x(1:M,i)')';
end
y=zeros(M,N);
for i=1:N
    y(1:M,i)=z(1:M,i);
end
w=y(:)';
subplot(3,1,1);
i=0:length(w)-1;
plot(i,w);
subplot(3,1,2);
spectureplot(w);
subplot(3,1,3);
spectureplot(h);