【干扰消除】最小二乘窄带干扰消除【含Matlab源码 2379期】

一、⛄获取代码方式

获取代码方式1：
完整代码已上传我的资源：【干扰消除】基于matlab最小二乘窄带干扰消除【含Matlab源码 2379期】
点击上面蓝色字体，直接付费下载，即可。

获取代码方式2：
付费专栏Matlab物理应用（初级版）

备注：
点击上面蓝色字体付费专栏Matlab物理应用（初级版），扫描上面二维码，付费29.9元订阅海神之光博客付费专栏Matlab物理应用（初级版），凭支付凭证，私信博主，可免费获得1份本博客上传CSDN资源代码（有效期为订阅日起，三天内有效）；
点击CSDN资源下载链接：1份本博客上传CSDN资源代码

⛄二、部分源代码

close all;clear;clc
%% 参数设定与信号生成
%-宽带信号的生成
tr=0.4;tf=2;arf=2.3;N=4096;
Fs=2;Ts=1/Fs;
point=(-2):Ts:(6);
g=1./(exp(-arfpoint/tr)+exp(arfpoint/tf));
s=g(2:end)-g(1:end-1);% s为有用宽带信号
plot(s);axis tight;
set(gca,‘XTick’,1:1.8:16);
set(gca,‘XTickLabel’,{‘-2’,‘-1’,‘0’,‘1’,‘2’,‘3’,‘4’,‘5’,‘6’});
xlabel(‘时间(ns)’);
ylabel(‘幅度’);
title(‘宽带信号’);
Ks=fft(s,1024);
Ks=fftshift(Ks);
AKs=abs(Ks);
PKs=phase(Ks);
figure;
subplot(211);
plot(AKs);axis tight;
set(gca,‘XTick’,1:102:1024);
set(gca,‘XTickLabel’,{‘-1’,‘-0.8’,‘-0.6’,‘-0.4’,‘-0.2’,‘0’,‘0.2’,‘0.4’,‘0.6’,‘0.8’,‘1’});
xlabel(‘频率(GHz)’);
ylabel(‘幅度’);
title(‘幅度’);
subplot(212);
plot(PKs);axis tight;
set(gca,‘XTick’,1:102:1024);
set(gca,‘XTickLabel’,{‘-1’,‘-0.8’,‘-0.6’,‘-0.4’,‘-0.2’,‘0’,‘0.2’,‘0.4’,‘0.6’,‘0.8’,‘1’});
xlabel(‘频率(GHz)’);
ylabel(‘相位(rad)’);
title(‘相位’);
%-窄带干扰信号的生成
F=0.1*[0.6,1,1.8,2.1,3,4.8,5.2,5.7,6.1,6.4,6.7,7,7.8,9.3]';
L=length(F);
om=2piF/Fs;
A=[0.5,1,1,0.5,0.1,0.3,0.5,1,1,1,0.5,0.3,1.5,0.5];
rand(‘seed’,1954);
phi=2pirand(L,1);
n=1:N;
y=Asin(omn+phiones(1,length(n)));
Ky=periodogram(y);
Ky=10log10(Ky);

⛄三、运行结果

⛄四、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1]刘萍.如何理解空气质量分指数（IAQI）计算公式并速算[J].黑龙江环境通报. 2014,38(02)

3 备注
简介此部分摘自互联网，仅供参考，若侵权，联系删除