%****************************************************
%**********王程 仿真实验二 PAST**********************
%****************************************************
clear all ;
f = 10e3 ; %sample frequency
N = 4000 ; %sample number
lumda = 0.003 ; %预白化步长
beta = 0.99 ; %遗忘因子
%-------------准备工作--------------------------------
Wt = eye(5) ; %预白化矩阵初始化
W = eye(5); %分离矩阵初始化,对应PAST中的W
P = eye(5); %对应PAST算法中的P
% A = randn(5); %对应非奇异混合矩阵A
A = rand(5, 5) ;
s = zeros(5, N) ; %观测向量矩阵
z = zeros(5, N) ; %分离后的信号向量矩阵
for k = 1 : N %总采样点数
n = k / f ;
%--------------生成源信号向量-------------------------
s1 = sign(cos(2*pi*155*n)) ; %符号信号
s2 = sin(2*pi*800*n) ; %正弦波
s3 = sin(2*pi*90*n) ; %正弦波
s4 = sin(2*pi*300*n+6*cos(2*pi*60*n)) ; %调频信号
s5 = unifrnd(-1, 1) ; %均匀分布的噪声
s(:, k) = [s1; s2; s3; s4; s5] ; %源信号向量
x = A * s(:, k) ; %混合观测信号向量
%--------------预白化----------------------------------
y = Wt * x ;
Wt = Wt - lumda * [ y * y' - eye(5) ] * Wt ;
%--------------PAST------------------------------------
z(:, k) = tanh(W' * y) ; %g = tanh()
h = P * z(:, k) ;
m = h /(beta + z(:, k)' * h) ;
tri_P = (triu(P)-triu(m*h'))/beta;
P = tri_P + tril(tri_P',-1);
e = y - W * z(:, k) ;
W = W + e * m' ;
end
%------------显示分离后的盲信号-----------------------
% figure(1) ; %源信号向量
% hold on ;
%
% subplot(5, 1, 1) ;
% plot(1:N, s(1,:)) ;
% ylabel('s1') ;
%
% subplot(5, 1, 2) ;
% plot(1:N, s(2,:)) ;
% ylabel('s2') ;
%
% subplot(5, 1, 3) ;
% plot(1:N, s(3,:)) ;
% ylabel('s3') ;
%
% subplot(5, 1, 4) ;
% plot(1:N, s(4,:)) ;
% ylabel('s4') ;
%
% subplot(5, 1, 5) ;
% plot(1:N, s(5,:)) ;
% ylabel('s5') ;
%
% hold off ;
figure(2) ; %分离后得到的信号向量
hold on ;
subplot(5, 1, 1) ;
plot(1:N, z(1,:)) ;
ylabel('y1') ;
subplot(5, 1, 2) ;
plot(1:N, z(2,:)) ;
ylabel('y2') ;
subplot(5, 1, 3) ;
plot(1:N, z(3,:)) ;
ylabel('y3') ;
subplot(5, 1, 4) ;
plot(1:N, z(4,:)) ;
ylabel('y4') ;
subplot(5, 1, 5) ;
plot(1:N, z(5,:)) ;
ylabel('y5') ;
hold off ;
没有合适的资源?快使用搜索试试~ 我知道了~
资源推荐
资源详情
资源评论
![rar](https://img-home.csdnimg.cn/images/20210720083606.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![pdf](https://img-home.csdnimg.cn/images/20210720083512.png)
![pdf](https://img-home.csdnimg.cn/images/20210720083512.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![rar](https://img-home.csdnimg.cn/images/20210720083606.png)
![rar](https://img-home.csdnimg.cn/images/20210720083606.png)
![application/pdf](https://img-home.csdnimg.cn/images/20210720083512.png)
![application/pdf](https://img-home.csdnimg.cn/images/20210720083512.png)
![application/pdf](https://img-home.csdnimg.cn/images/20210720083512.png)
![application/pdf](https://img-home.csdnimg.cn/images/20210720083512.png)
![application/pdf](https://img-home.csdnimg.cn/images/20210720083512.png)
![doc](https://img-home.csdnimg.cn/images/20210720083327.png)
![pdf](https://img-home.csdnimg.cn/images/20210720083512.png)
![pdf](https://img-home.csdnimg.cn/images/20210720083512.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![pdf](https://img-home.csdnimg.cn/images/20210720083512.png)
收起资源包目录
![package](https://csdnimg.cn/release/downloadcmsfe/public/img/package.f3fc750b.png)
![file-type](https://csdnimg.cn/release/download/static_files/pc/images/minetype/TXT.png)
![file-type](https://csdnimg.cn/release/download/static_files/pc/images/minetype/UNKNOWN.png)
共 2 条
- 1
资源评论
![avatar-default](https://csdnimg.cn/release/downloadcmsfe/public/img/lazyLogo2.1882d7f4.png)
![avatar](https://profile-avatar.csdnimg.cn/dc78d2406d17417ca42db3bd43b9c72a_weixin_42652674.jpg!1)
御道御小黑
- 粉丝: 62
- 资源: 1万+
上传资源 快速赚钱
我的内容管理 展开
我的资源 快来上传第一个资源
我的收益
登录查看自己的收益我的积分 登录查看自己的积分
我的C币 登录后查看C币余额
我的收藏
我的下载
下载帮助
![voice](https://csdnimg.cn/release/downloadcmsfe/public/img/voice.245cc511.png)
![center-task](https://csdnimg.cn/release/downloadcmsfe/public/img/center-task.c2eda91a.png)
安全验证
文档复制为VIP权益,开通VIP直接复制
![dialog-icon](https://csdnimg.cn/release/downloadcmsfe/public/img/green-success.6a4acb44.png)