matlab编程界面显示和信号处理、滤波、分离线性调频信号分析源码程序 - matlab通信信号

下载频道> 资源分类> matlab源码> 通信信号> matlab编程界面显示和信号处理、滤波、分离线性调频信号分析源码程序

标题：matlab编程界面显示和信号处理、滤波、分离线性调频信号分析源码程序
分享到：

所属分类: 通信信号资源类型：程序源码文件大小: 4.47 KB 上传时间: 2019-12-01 21:39:50 下载次数: 6 资源积分：1分提供者: jiqiren 20191201094009208

内容：
% 线性调频信号分析：

% (1)Chirp复信号的建模及采样

% (2)Chirp复信号的时频特性研究

% (3)匹配滤波器实现Chirp复信号的脉冲压缩

% (4)窗函数对Chirp脉压的影响

% (5)匹配滤波器失配分析

%--------------------------------------------------

% 1、雷达仿真参数及初始化：

Tao=1.e-6; %s 雷达发射脉冲宽度

B=9.e+7; %Hz 雷达发射带宽

f0=1.e+10; %Hz 雷达工作频率,即信号载频

fs=2.e+8; %Hz A/D采样率

N=fs*Tao; % 采样点数

n=-N/2:N/2-1; % 采样点分配区间

t=n/fs; % 采样时间序列

%----------------------------------------------------------

figure(1);

%2、Chirp复信号的建模：

x=exp(j*pi*B*(t.^2)/Tao); % 线性调频信号建摸

subplot(2,3,1);

plot(t,x); % 时域信号图

title('一维Chirp信号');

xlabel('time[s]');

ylabel('scope');

%---------------------------------------------------------

% 3、信号FFT变换后分析频域特性：

X=fft(x); % 傅立叶变换

X1=fftshift(X); % 移动零点到频谱中心

magX=abs(X1/max(X1)); % 幅频特性

angX=angle(X); % 相频特性

k=(-length(X)/2:(length(X)-1)/2)/length(X); % 单位长度点数

f=fs*k; % 频率

subplot(2,3,2);

plot(f,magX); % 幅频图

axis([-B B 0 1.2]);

title('幅频谱');

xlabel('frequency[Hz]');

ylabel('scope');

subplot(2,3,3);

plot(f,unwrap(angX)); % 相频图

grid;

title('相频谱');

xlabel('frequency[Hz]');

ylabel('angle[drgree]');

%---------------------------------------------------------

% 4、匹配滤波器的建立及时频特性：

h=exp(-j*pi*B*(t.^2)/Tao); % 匹配滤波器建模

subplot(2,3,4);

plot(t,h); % 匹配滤波器时域图

title('匹配滤波器');

xlabel('time[s]');

ylabel('scope');

H=fft(h);

H1=fftshift(H);

magH=abs(H1/max(H1)); % 匹配滤波器幅频特性

angH=angle(H); % 匹配滤波器相频特性

subplot(2,3,5);

plot(f,magH); % 幅频图

axis([-B B 0 1.2]);

title('幅频谱');

xlabel('frequency[Hz]');

ylabel('scope');

subplot(2,3,6);

plot(f,unwrap(angH)); % 相频图

grid;

title('相频谱');

xlabel('frequency[Hz]');

ylabel('angle[drgree]');

%------------------------------------------------------------

figure(2)

% 5、Chirp信号通过匹配滤波器，实现脉冲压缩

Y=X1.*H1; % 通过滤波器的信号

magY=abs(Y/max(Y)); % 幅频特性

angY=angle(Y); % 相频特性

y=ifft(Y); % 变换到时域

y1=ifftshift(y);

y11=abs(y1)/max(abs(y1)); % 归一化

y111=20*log10(y11); % 用对数表示

subplot(2,3,1);

plot(t,y11);

grid;

axis([-Tao/10 Tao/10 0 1]);

title('信号通过滤波器');

xlabel('time[s]');

ylabel('scope');

subplot(2,3,2);

plot(f,magY);

axis([-B B 0 1.2]);

title('幅频谱');

xlabel('frequency[Hz]');

ylabel('scope');

subplot(2,3,3);

plot(f,angY);

title('相频谱');

xlabel('frequency[Hz]');

ylabel('angle[drgree]');

subplot(2,3,4);

plot(t,y111);

grid;

axis([-Tao/10 Tao/10 -30 0]);

title('信号通过滤波器对数图');

xlabel('time[s]');

ylabel('scope[dB]');

% 用对数形式表示，并插值

% 原理：频域补零相当于时域插值

Z=X1.*H1;

t_add0=interp(t,16); % 时间插值

%加插值函数

[intp_sig,intp_logsig]=intp(Z); % 插值

z11=intp_sig;

z111=intp_logsig;

subplot(2,3,5);

plot(t_add0,z11);

grid;

axis([-Tao/10 Tao/10 0 1]);

title('信号通过滤波器/插值后');

xlabel('time[s]');

ylabel('scope');

subplot(2,3,6)

plot(t_add0,z111);

grid;

axis([-Tao/10 Tao/10 -30 0]);

title('信号通过滤波器对数图/插值后');

xlabel('time[s]');

ylabel('scope[dB]');

% 加评估函数

[O,width3dB,ratio,peak1,peak_1,peak2,peak_2,peak3,peak_3]=para(z11,t_add0); % 评估参数

width3dB

ratio

peak1

peak_1

peak2

peak_2

peak3

peak_3

%-------------------------------------------------------------

% 6、窗函数对Chirp脉压的影响（在频域中实现）

%---------------------------------------------------------------------

figure(3);

% 6.1 生成矩形窗函数

lengthY1=length(Y);

c1=zeros(1,lengthY1);

c1(71:130)=boxcar(60); % 矩形窗

subplot(2,2,1);

plot(f,c1);

axis([-B B 0 1.2]);

title('矩形窗');

xlabel('frequency[Hz]');

ylabel('scope');

% 矩形窗对Chirp脉压的影响

M1=abs(Y).*c1; % 对信号加矩形窗

subplot(2,2,2);

plot(f,M1/max(M1));

axis([-B B 0 1.2]);

title('矩形窗所截频域信号');

xlabel('frequency[Hz]');

ylabel('scope');

[intp_sig,intp_logsig]=intp(M1);

magD=intp_sig;

d11=intp_logsig;

subplot(2,2,3);

plot(t_add0,magD);

axis([-Tao/10 Tao/10 0 1]);

title('通过矩形窗信号');

xlabel('time[s]');

ylabel('scope');

subplot(2,2,4);

plot(t_add0,d11);

axis([-Tao/10 Tao/10 -30 0]);

grid;

title('通过矩形窗信号对数图');

xlabel('time[s]');

ylabel('scope[dB]');

% 加评估函数

[O,width3dB,ratio,peak1,peak_1,peak2,peak_2,peak3,peak_3]=para(magD,t_add0); % 评估参数

width3dB

ratio

peak1

peak_1

peak2

peak_2

peak3

peak_3

%----------------------------------------------------------------------------------

figure(4);

% 6.2 生成海明窗函数

c2=zeros(1,lengthY1);

c2(51:150)=hamming(100); % 海明窗

subplot(2,2,1);

plot(f,c2);

axis([-B B 0 1.2]);

title('海明窗');

xlabel('frequency[Hz]');

ylabel('scope');

% 海明窗对Chirp脉压的影响

M2=abs(Y).*c2; % 对信号加海明窗

subplot(2,2,2);

plot(f,M2/max(M2));

axis([-B B 0 1.2]);

title('海明窗所截频域信号');

xlabel('frequency[Hz]');

ylabel('scope');

[intp_sig,intp_logsig]=intp(M2);

magDm2=intp_sig;

dm211=intp_logsig;

subplot(2,2,3);

plot(t_add0,magDm2);

axis([-Tao/10 Tao/10 0 1]);

title('通过海明窗信号');

xlabel('time[s]');

ylabel('scope');

subplot(2,2,4);

plot(t_add0,dm211);

axis([-Tao/10 Tao/10 -40 0]);

grid;

title('通过海明窗信号对数图');

xlabel('time[s]');

ylabel('scope[dB]');

% 加评估函数

[O,width3dB,ratio,peak1,peak_1,peak2,peak_2,peak3,peak_3]=para(magDm2,t_add0); % 评估参数

width3dB

ratio

peak1

peak_1

peak2

peak_2

peak3

peak_3

%------------------------------------------------------------------------

figure(5);

% 6.3 生成三角窗函数

c3=zeros(1,lengthY1);

c3(61:140)=triang(80); % 三角窗

subplot(2,2,1);

plot(f,c3);

axis([-B B 0 1.2]);

title('三角窗');

xlabel('frequency[Hz]');

ylabel('scope');

% 三角窗对Chirp脉压的影响

M3=abs(Y).*c3; % 对信号加三角窗

subplot(2,2,2);

plot(f,M3/max(M3));

axis([-B B 0 1.2]);

title('三角窗所截频域信号');

xlabel('frequency[Hz]');

ylabel('scope');

[intp_sig,intp_logsig]=intp(M3);

magDm3=intp_sig;

dm311=intp_logsig;

subplot(2,2,3);

plot(t_add0,magDm3);

axis([-Tao/10 Tao/10 0 1]);

title('通过三角窗信号');

xlabel('time[s]');

ylabel('scope');

subplot(2,2,4);

plot(t_add0,dm311);

axis([-Tao/10 Tao/10 -45 0]);

grid;

title('通过三角窗信号对数图');

xlabel('time[s]');

ylabel('scope[dB]');

% 加评估函数

[O,width3dB,ratio,peak1,peak_1,peak2,peak_2,peak3,peak_3]=para(magDm3,t_add0); % 评估参数

width3dB

ratio

peak1

peak_1

peak2

peak_2

peak3

peak_3

%-------------------------------------------------------------------------

% 7、滤波起失配情况讨论:分四（或七）种情况

% 在s(t)=exp(j*2*pi(f0*t+u*t.^2+r(t)))中，令

% （1）f0不相等，u与r(t)相等

% （2）u不相等，其余相等

% （3）r(t)不相等，其余相等

% （4）f0,u,r(t)都不相等

% （5）加高频干扰

figure(6);

% 7.1

% f0?=0时Chirp复信号的建模：

s1=exp(j*pi*B*(t.^2)/Tao+j*2*pi*B/10*t); % 线性调频信号建摸

subplot(2,3,1);

plot(t,s1); % 时域信号图

title('f 0？=0时Chirp信号');

xlabel('time[s]');

ylabel('scope');

% 信号FFT变换后分析频域特性：

S1=fft(s1); % 傅立叶变换

angS1=angle(S1); % 相频特性

S11=fftshift(S1); % 移动零点到频谱中心

magS1=abs(S11/max(S11)); % 幅频特性

k=(-length(S1)/2:(length(S1)-1)/2)/length(S1); % 单位长度点数

f=fs*k; % 频率

subplot(2,3,2);

plot(f,magS1); % 幅频图

axis([-B B 0 1.2]);

title('幅频谱');

xlabel('frequency[Hz]');

ylabel('scope');

subplot(2,3,3);

plot(f,unwrap(angS1)); % 相频图

grid;

title('相频谱');

xlabel('frequency[Hz]');

ylabel('angle[drgree]');

% Chirp信号通过匹配滤波器，实现脉冲压缩

R1=S11.*H1; % 通过滤波器的信号

magR1=abs(R1/max(R1)); % 幅频特性

angR1=angle(R1); % 相频特性

r1=ifft(R1); % 变换到时域

r1=ifftshift(r1);

r11=abs(r1)/max(abs(r1)); % 归一化

r111=20*log10(r11); % 用对数表示

subplot(2,3,5);

plot(f,magR1);

axis([-B B 0 1.2]);

title('通过失配滤波器幅频谱');

xlabel('frequency[Hz]');

ylabel('scope');

% 用对数形式表示，并插值

% 原理：频域补零相当于时域插值

W1=S1.*H;

W1=fftshift(W1);

[intp_sig,intp_logsig]=intp(W1);

w11=intp_sig;

w111=intp_logsig;

subplot(2,3,4);

plot(t_add0,w11);

grid;

axis([-Tao/10-Tao/10 Tao/10 0 1]);

title('信号通过失配滤波器');

xlabel('time[s]');

ylabel('scope');

subplot(2,3,6)

plot(t_add0,w111);

grid;

axis([-Tao/10-Tao/10 Tao/10 -30 0]);

title('信号通过失配滤波器对数图');

xlabel('time[s]');

ylabel('scope[dB]');

% 加评估函数

[O,width3dB,ratio,peak1,peak_1,peak2,peak_2,peak3,peak_3]=para(w11,t_add0); % 评估参数

width3dB

ratio

peak1

peak_1

peak2

peak_2

peak3

peak_3

%-------------------------------------------------------------

figure(7);

% 7.2

% u1?=u时Chirp复信号的建模：

s2=exp(j*pi*(B+1.e+7)*(t.^2)/(Tao)); % 线性调频信号建摸

subplot(2,3,1);

plot(t,s2); % 时域信号图

title('u 1？=u时Chirp信号');

xlabel('time[s]');

ylabel('scope');

% 信号FFT变换后分析频域特性：

S2=fft(s2); % 傅立叶变换

angS2=angle(S2); % 相频特性

S22=fftshift(S2); % 移动零点到频谱中心

magS2=abs(S22/max(S22)); % 幅频特性

k=(-length(S2)/2:(length(S2)-1)/2)/length(S2); % 单位长度点数

f=fs*k; % 频率

subplot(2,3,2);

plot(f,magS2); % 幅频图

axis([-B B 0 1.2]);

title('幅频谱');

xlabel('frequency[Hz]');

ylabel('scope');

subplot(2,3,3);

plot(f,unwrap(angS2)); % 相频图

grid;

title('相频谱');

xlabel('frequency[Hz]');

ylabel('angle[drgree]');

% Chirp信号通过匹配滤波器，实现脉冲压缩

R2=S22.*H1; % 通过滤波器的信号

magR2=abs(R2/max(R2)); % 幅频特性

angR2=angle(R2); % 相频特性

r2=ifft(R2); % 变换到时域

r2=ifftshift(r2);

r21=abs(r2)/max(abs(r2)); % 归一化

r211=20*log10(r21); % 用对数表示

subplot(2,3,5);

plot(f,magR2);

axis([-B B 0 1.2]);

title('通过失配滤波器幅频谱');

xlabel('frequency[Hz]');

ylabel('scope');

% 用对数形式表示，并插值

% 原理：频域补零相当于时域插值

W2=S2.*H;

W2=fftshift(W2);

[intp_sig,intp_logsig]=intp(W2);

w21=intp_sig;

w211=intp_logsig;

subplot(2,3,4);

plot(t_add0,w21);

grid;

axis([-Tao/4 Tao/4 0 1]);

title('信号通过失配滤波器');

xlabel('time[s]');

ylabel('scope');

subplot(2,3,6)

plot(t_add0,w211);

grid;

axis([-Tao/4 Tao/4 -30 0]);

title('信号通过失配滤波器对数图');

xlabel('time[s]');

ylabel('scope[dB]');

% 加评估函数

[O,width3dB,ratio,peak1,peak_1,peak2,peak_2,peak3,peak_3]=para(w21,t_add0); % 评估参数

width3dB

ratio

peak1

peak_1

peak2

peak_2

peak3

peak_3

%-------------------------------------------------------------

figure(8);

% 7.3

% r(t)?=0时Chirp复信号的建模：

s3=exp(j*pi*B*(t.^2)/Tao+j*1.e+19*t.^3); % 线性调频信号建摸

subplot(2,3,1);

plot(t,s3); % 时域信号图

title('r(t)？=0时Chirp信号');

xlabel('time[s]');

ylabel('scope');

% 信号FFT变换后分析频域特性：

S3=fft(s3); % 傅立叶变换

angS3=angle(S3); % 相频特性

S33=fftshift(S3); % 移动零点到频谱中心

magS3=abs(S33/max(S33)); % 幅频特性

k=(-length(S3)/2:(length(S3)-1)/2)/length(S3); % 单位长度点数

f=fs*k; % 频率

subplot(2,3,2);

plot(f,magS3); % 幅频图

axis([-B B 0 1.2]);

title('幅频谱');

xlabel('frequency[Hz]');

ylabel('scope');

subplot(2,3,3);

plot(f,unwrap(angS3)); % 相频图

grid;

title('相频谱');

xlabel('frequency[Hz]');

ylabel('angle[drgree]');

% Chirp信号通过匹配滤波器，实现脉冲压缩

R3=S33.*H1; % 通过滤波器的信号

magR3=abs(R3/max(R3)); % 幅频特性

angR3=angle(R3); % 相频特性

r3=ifft(R3); % 变换到时域

r3=ifftshift(r3);

r31=abs(r3)/max(abs(r3)); % 归一化

r311=20*log10(r31); % 用对数表示

subplot(2,3,5);

plot(f,magR3);

axis([-B B 0 1.2]);

title('通过失配滤波器幅频谱');

xlabel('frequency[Hz]');

ylabel('scope');

% 用对数形式表示，并插值

% 原理：频域补零相当于时域插值

W3=S3.*H;

W3=fftshift(W3);

[intp_sig,intp_logsig]=intp(W3);

w31=intp_sig;

w311=intp_logsig;

subplot(2,3,4);

plot(t_add0,w31);

grid;

axis([-Tao/10 Tao/10 0 1]);

title('信号通过失配滤波器');

xlabel('time[s]');

ylabel('scope');

subplot(2,3,6)

plot(t_add0,w311);

grid;

axis([-Tao/10 Tao/10 -30 0]);

title('信号通过失配滤波器对数图');

xlabel('time[s]');

ylabel('scope[dB]');

% 加评估函数

[O,width3dB,ratio,peak1,peak_1,peak2,peak_2,peak3,peak_3]=para(w31,t_add0); % 评估参数

width3dB

ratio

peak1

peak_1

peak2

peak_2

peak3

peak_3

%-------------------------------------------------------------

figure(9);

% 7.3

% 高频干扰Chirp复信号的建模：

s4=exp(j*pi*B*(t.^2)/Tao+j*sin(2.e+10*t)); % 线性调频信号建摸

subplot(2,3,1);

plot(t,s4); % 时域信号图

title('高频干扰时Chirp信号');

xlabel('time[s]');

ylabel('scope');

% 信号FFT变换后分析频域特性：

S4=fft(s4); % 傅立叶变换

angS4=angle(S4); % 相频特性

S43=fftshift(S4); % 移动零点到频谱中心

magS4=abs(S43/max(S43)); % 幅频特性

k=(-length(S4)/2:(length(S4)-1)/2)/length(S4); % 单位长度点数

f=fs*k; % 频率

subplot(2,3,2);

plot(f,magS4); % 幅频图

axis([-B B 0 1.2]);

title('幅频谱');

xlabel('frequency[Hz]');

ylabel('scope');

subplot(2,3,3);

plot(f,unwrap(angS4)); % 相频图

grid;

title('相频谱');

xlabel('frequency[Hz]');

ylabel('angle[drgree]');

% Chirp信号通过匹配滤波器，实现脉冲压缩

R4=S43.*H1; % 通过滤波器的信号

magR4=abs(R4/max(R4)); % 幅频特性

angR4=angle(R4); % 相频特性

r4=ifft(R4); % 变换到时域

r4=ifftshift(r4);

r41=abs(r4)/max(abs(r4)); % 归一化

r411=20*log10(r41); % 用对数表示

subplot(2,3,5);

plot(f,magR4);

axis([-B B 0 1.2]);

title('通过失配滤波器幅频谱');

xlabel('frequency[Hz]');

ylabel('scope');

% 用对数形式表示，并插值

% 原理：频域补零相当于时域插值

W4=S4.*H;

W4=fftshift(W4);

[intp_sig,intp_logsig]=intp(W4);

w41=intp_sig;

w411=intp_logsig;

subplot(2,3,4);

plot(t_add0,w41);

grid;

%axis([-Tao/10 Tao/10 0 1]);

title('信号通过失配滤波器');

xlabel('time[s]');

ylabel('scope');

subplot(2,3,6)

plot(t_add0,w411);

grid;

%axis([-Tao/10 Tao/10 -30 0]);

title('信号通过失配滤波器对数图');

xlabel('time[s]');

ylabel('scope[dB]');

% 加评估函数

[O,width3dB,ratio,peak1,peak_1,peak2,peak_2,peak3,peak_3]=para(w41,t_add0); % 评估参数

width3dB

ratio

peak1

peak_1

peak2

peak_2

peak3

peak_3

%-------------------------------------------------------------

文件列表（点击上边下载按钮，如果是垃圾文件请在下面评价差评或者投诉）:

matlab编程界面显示和信号处理、滤波、分离线性调频信号分析源码程序/
matlab编程界面显示和信号处理、滤波、分离线性调频信号分析源码程序/code/
matlab编程界面显示和信号处理、滤波、分离线性调频信号分析源码程序/source/
matlab编程界面显示和信号处理、滤波、分离线性调频信号分析源码程序/结果/
matlab编程界面显示和信号处理、滤波、分离线性调频信号分析源码程序/code/main.m
matlab编程界面显示和信号处理、滤波、分离线性调频信号分析源码程序/程序说明.txt

关键词：信号处理 滤波 分离 线性 调频 信号 分析 源码 程序

matlab源码下载排行

回到顶部

友情链接

matlab源码网站源码编程语言网络编程手机通讯编程 KUNAG 官网 KUNAG 电商更多

: 联系方式| 版权声明| 招聘信息| 广告服务| 银行汇款| 法律顾问| 兼职技术| 付款方式| 关于我们|; 网站客服程序员兼职招聘 411241577@qq.com; 沪ICP备19040327号-3; 公安备案号：沪公网安备 31011802003874号; 库纳格流体控制系统（上海）有限公司版权所有; Copyright © 1999-2014, GUSUCODE.COM, All Rights Reserved

关键词： 信号处理 滤波 分离 线性 调频 信号 分析 源码 程序

相关推荐

matlab源码下载排行

友情链接

关键词：信号处理滤波分离线性调频信号分析源码程序