clc;clear;close all;
%% 采集信号
load m5
fs=10.24*10^3;
x=m5(:,1);x=x(500:3000);
N=length(x);%数据长度
x = awgn(x,7,'measured');%信噪比18dB
subplot(2,1,1),plot(x,'k.-','LineWidth',1,'MarkerSize',10);
% set(gca,'Fontname','Times New Roman','FontSize','12');
legend('5m处采样信号');xlim([0 1000]);xlabel('采样点数/n');ylabel('x(n)');%原始波形
%% 频谱分析
y=fft(x,N);%进行fft变换
mag=abs(y);%求幅值
f=(0:length(y)-1)'*fs/length(y)*16.5;%进行对应的频率转换
subplot(2,1,2),plot(f(1:N/2),mag(1:N/2),'ko-','LineWidth',1);%做频谱图%
xlim([0 30000]);xlabel('频率/Hz');ylabel('|X(k)|');legend('5m处信号频谱');
%% 低通测试
z=highp(x,300,500,0.2,20,fs);
z=lowp(z,1300,900,0.2,20,fs);
z=z+x*0.06;
z = awgn(z,18,'measured')/2.6;%信噪比18dB
figure;
subplot(211);plot(z,'k.-','LineWidth',1,'MarkerSize',10);xlim([0 1000])
s=fft(z,N);%进行fft变换
mag=(abs(s));%求幅值
f=(0:length(s)-1)'*fs/length(s)*16.5;%进行对应的频率转换
subplot(2,1,2),plot(f(1:N/2),mag(1:N/2),'ko-','LineWidth',1,'MarkerSize',5);%做频谱图%
xlim([0 30000]);xlabel('频率/Hz');ylabel('|X(k)|');legend('5m处信号频谱');
